fiat team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

Can you try the following (in Python) and see if you get the same
error. In that case, we can nail it down to a FIAT problem, which
would be good since then it would be Rob's problem and not mine. :-)

from FIAT.Lagrange import *
from FIAT.shapes import *

element = VectorLagrange(TRIANGLE, 2)

/Anders


On Tue, May 02, 2006 at 09:48:40AM +0000, Alexander Jarosch wrote:
> 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/
> 
> 
> _______________________________________________
> FFC-dev mailing list
> FFC-dev@xxxxxxxxxx
> http://www.fenics.org/cgi-bin/mailman/listinfo/ffc-dev