On Fri, Oct 09, 2009 at 09:31:27AM +0200, Johan Hake wrote: > On Friday 09 October 2009 09:18:30 Johan Hake wrote: > > On Friday 09 October 2009 09:06:19 azumi.hayak...@kuroda-electric.com wrote: > > > Hello list, > > > > > > I would like to use get_eigenpair function from within > > > python programmes. The 1st thing I tried is modifying > > > demo/la/eigensolver/python/demo.py as follows: > > > ---- original > > > (lr, lc) = esolver.get_eigenvalue(i) > > > > > > ---- modified > > > r = PETScVector(A.size(0)) > > > c = PETScVector(A.size(0)) > > > esolver.get_eigenpair(lr, lc, r, c, i) > > > > There is a bug in the Python wrapper (thanks for noticing!). > > > > Try using > > > > lr, lc = esolver.getEigenpair(r, c, i) > > > > if you are on stable. If you pull the latest code you will soon be able to > > use: > > > > lr, lc = esolver.get_eigenpair(r, c, i) > > While at it shouldn't the PETScVectors also be return arguments? Leaving us > with: > > r, c, lr, lc = esolver.get_eigenpair(i)
Yes, or rather lr, lc, r, c = esolver.get_eigenpair(i) as in C++. We should also add get_eigenvector to just return r and c. -- Anders > Johan > > > Johan > > > > > This gave me an error below: > > > Traceback (most recent call last): > > > File "./demo.py", line 46, in <module> > > > esolver.get_eigenpair(lr, lc, r, c, i) > > > TypeError: in method 'SLEPcEigenSolver_get_eigenpair', argument 2 of > > > type 'double &' > > > > > > Could some one give some advice on how to use > > > get_eigenpair function? Incidentally, cpp demo > > > ran without an error, and seemed to work fine. > > > > > > Thanks in advance, > > > -- Azumi Hayakawa > > > _______________________________________________ > > > FEniCS-users mailing list > > > fenics-us...@fenics.org > > > http://fenics.org/mailman/listinfo/fenics-users > > > > _______________________________________________ > > FEniCS-users mailing list > > fenics-us...@fenics.org > > http://fenics.org/mailman/listinfo/fenics-users > > > _______________________________________________ > DOLFIN-dev mailing list > DOLFIN-dev@fenics.org > http://www.fenics.org/mailman/listinfo/dolfin-dev
signature.asc
Description: Digital signature
_______________________________________________ DOLFIN-dev mailing list DOLFIN-dev@fenics.org http://www.fenics.org/mailman/listinfo/dolfin-dev
