ffc team mailing list archive

[Alexander Jarosch <alexanj@xxxxx>]

Hello everybody,

I just pulled the recent ffc and tried to compile a form, but I get the 
output:

# ffc -d 3 Stokes.form
This is FFC, the FEniCS Form Compiler, version 0.3.2.
For further information, go to http://www/fenics.org/ffc/.
Parsing Stokes.form
Output written to Stokes.py
Traceback (most recent call last):
 File "/usr/bin/ffc", line 125, in ?
   main(sys.argv[1:])
 File "/usr/bin/ffc", line 82, in main
   execfile(outname, ns)
 File "Stokes.py", line 21, in ?
   P2 = FiniteElement("Vector Lagrange", "triangle", 2)
 File "/usr/lib/python2.4/site-packages/ffc/compiler/finiteelement.py", 
line 
73                                                                                                                                

, in __init__
   self.element = VectorLagrange(self.fiat_shape, degree, num_components)
 File "/usr/lib/python2.4/site-packages/FIAT/Lagrange.py", line 108, in 
__init_                                                                                                                                

_
   U = polynomial.OrthogonalPolynomialArraySet( shape , n , nc )
 File "/usr/lib/python2.4/site-packages/FIAT/polynomial.py", line 394, 
in 
Ortho                                                                                                                                

gonalPolynomialArraySet
   b = PolynomialBase( element_shape , degree )
 File "/usr/lib/python2.4/site-packages/FIAT/polynomial.py", line 29, 
in 
__init                                                                                                                                

__
   vinv = LinearAlgebra.inverse( v )
 File "/usr/lib64/python2.4/site-packages/Numeric/LinearAlgebra.py", 
line 
108,                                                                                                                                 

in inverse
   return solve_linear_equations(a, Numeric.identity(a.shape[0]))
 File "/usr/lib64/python2.4/site-packages/Numeric/Numeric.py", line 
604, in 
ide                                                                                                                                

ntity
   return resize(array([1]+n*[0],typecode=typecode), (n,n))
 File "/usr/lib64/python2.4/site-packages/Numeric/Numeric.py", line 
398, in 
res                                                                                                                                

ize
   return reshape(a, new_shape)
ValueError: total size of new array must be unchanged

looks like a numeric package problem, but I have to say, I'm stuck here. 
Does anybody have an idea what to do? I use python 2.4 and the numeric 
package 23.7

thanks in advance

Alex

--
Alexander H. Jarosch

Jarðvísindastofnun Háskólans
Institute of Earth Sciences, University of Iceland
Náttúrufræðahús, Askja
Building of Natural Sciences, Askja
Sturlugata 7
IS - 101 Reykjavík
Iceland

Tel.: +354 525 4906
http://raunvis.hi.is/~jarosch/