ffc team mailing list archive

[Alexander Jarosch <alexanj@xxxxx>]

Yups, did it and get this:

Python 2.4.2 (#1, May  1 2006, 12:35:13)
[GCC 3.4.5 (Gentoo 3.4.5, ssp-3.4.5-1.0, pie-8.7.9)] on linux2
Type "help", "copyright", "credits" or "license" for more information.
>>> from FIAT.Lagrange import*
>>> from FIAT.shapes import *
>>> element = VectorLagrange(TRIANGLE, 2)
Traceback (most recent call last):
 File "<stdin>", line 1, in ?
 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
>>>

so it seems that it is not your problem, well shall I post it on the 
FIAT list or does Rob get it through this one?

Alex

Anders Logg wrote:

> 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
> >    
>
>  


--
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/