dolfin team mailing list archive

["Garth N. Wells" <gnw20@xxxxxxxxx>]

Robert Kirby wrote:
> I didn't see it as a sub of cpp.  thanks for pointing that out.  no_prec 
> is a better name
> given what "None" means in python.

I find the enums is C++ clumsy, so is there a disadvantage in just using 
strings for setting solvers, preconditioners, etc?

Garth

> On Mon, Apr 13, 2009 at 3:00 PM, Johan Hake <hake@xxxxxxxxx 
> <mailto:hake@xxxxxxxxx>> wrote:
>
>     On Monday 13 April 2009 21:32:22 Robert Kirby wrote:
>      > Aha.  I forgot about the inheritance issue.This should die on
>     forming an
>      > ILU since I don't have matrix values.
>      >
>      > doflin::none is not exposed to PyDOLFIN, so I can't set the
>     preconditioner
>      > to do nothing rather than attempt to compute an ILU.
>
>     It is exposed to dolfin.cpp. You can for example do:
>
>      dolfin.uBLASKrylovSolver( dolfin.gmres, dolfin.cpp.none )
>
>     We can allways expose it to dolfin too but however I am not sure
>     "none" is a
>     good name. Maybe "no_prec" is better?
>
>     Johan
>
>      > Thanks,
>      > Rob
>      >
>      > On Mon, Apr 13, 2009 at 2:30 PM, Johan Hake <hake@xxxxxxxxx
>     <mailto:hake@xxxxxxxxx>> wrote:
>      > > On Wednesday 08 April 2009 23:15:42 Robert Kirby wrote:
>      > > > Hi all,I've gotten a prototype working where I construct
>      > > > a uBLASSparseMatrix and get it to compute the same result as my
>      > > > matrix-free uBLASKrylovMatrix.  Both of these are constructed
>      > > > in PyDOLFIN.  However, I can feed the uBLASSparseMatrix
>      > > > I construct into a dolfin.uBLASKrylovSolver solve method,
>      > > > but I get a type error when I try to put my
>      > > > uBLASKrylovMatrix with mult and dim implemented in Python
>      > > > into the solver.
>      > > >
>      > > > In fact, I get a
>      > > >
>      > > > TypeError: in method 'uBLASKrylovSolver_solve', argument 2 of
>     type
>      > > > 'dolfin::uBLASKrylovMatrix const &'
>      > >
>      > > You need to call:
>      > >
>      > >  dolfin.uBLASKrylovMatrix.__init__(self)
>      > >
>      > > to initialize the super class.
>      > >
>      > > When I did this I got a bit further. Now an assertion is triggered:
>      > >
>      > > *** Assertion (_x.size() == _M.size1())
>      > >    [at dolfin/la/uBLASILUPreconditioner.cpp:41 in solve()]
>      > >
>      > > Johan
>      > >
>      > > > Not very revealing.  I've attached a horribly ugly source
>     code for
>      > > > anyone interested in doing a post-mortem on this.
>      > > >
>      > > > Thanks,
>      > > > Rob
>
>
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> DOLFIN-dev mailing list
> DOLFIN-dev@xxxxxxxxxx
> http://www.fenics.org/mailman/listinfo/dolfin-dev