ffc team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

Thanks for the report. It looks like there is a Numeric.vdot where
there should be a Numeric.dot. Changing from

    # Otherwise, use Numeric.vdot
    return Numeric.vdot(vec(v), vec(w))

to

    # Otherwise, use Numeric.dot
    return Numeric.dot(vec(v), vec(w))

in src/ffc/compilers/operators.py should do the trick.

If possible, could you try it out and see if it works?

/Anders

On Sun, Sep 25, 2005 at 02:30:50PM -0700, osman buyukisik wrote:
> On Sunday 25 September 2005 07:35, Anders Logg wrote:
> > Which version of Python and Python Numeric are you running?
> > I'm using Python 2.4.1 and Python Numeric 23.8-3 (from Debian)
> > and all forms compile cleanly on my system.
> >
> > This looks like a problem with using the correct dot function. The
> > Python Numeric dot function tries to take the conjugate of one of the
> > arguments. It would help if you commented out some of the exception
> > handling code in src/bin/ffc:
> >
> > #        except RuntimeError, exception:
> > #            print "*** " + str(exception)
> > #        except Exception, exception:
> > #            print "*** " + str(exception)
> >
> > which would print out a traceback to give some more information about
> > the error.
> 
> osman@linux:~/ffc-0.2.0/src/demo> ffc Heat.form
> This is FFC, the FEniCS Form Compiler, version 0.2.0.
> For further information, go to http://www/fenics.org/ffc/.
> 
> Parsing Heat.form
> Output written to Heat.py
> Traceback (most recent call last):
>   File "/usr/bin/ffc", line 94, in ?
>     main(sys.argv[1:])
>   File "/usr/bin/ffc", line 64, in main
>     execfile(outname, ns)
>   File "Heat.py", line 38, in ?
>     a = u1*v*dx + k*c*dot(grad(u1), grad(v))*dx
>   File "/usr/lib/python2.3/site-packages/ffc/compiler/operators.py", line 65, 
> in dot
>     return Numeric.vdot(vec(v), vec(w))
>   File "/usr/lib/python2.3/site-packages/Numeric/Numeric.py", line 353, in 
> vdot
>     return multiarray.matrixproduct(conjugate(ravel(a)),
> AttributeError: Product instance has no attribute 'conjugate'
> osman@linux:~/ffc-0.2.0/src/demo>       
> 
> _______________________________________________
> FFC-dev mailing list
> FFC-dev@xxxxxxxxxx
> http://www.fenics.org/cgi-bin/mailman/listinfo/ffc-dev
> 

-- 
Anders Logg
Research Assistant Professor
Toyota Technological Institute at Chicago
http://www.tti-c.org/logg/