Re: [Getfem-users] solvers won't solve

[Ronan Perrussel]

Dear Michael,

if you use the "stable" version of the web site, I am quite sure it will 
not work with numpy (this version seems to be anterior to all the 
modifications made by Luis).
You should use the latest svn version (some information here: 
https://gna.org/svn/?group=getfem).

Best regards,
Ronan

address@hidden a écrit :
> Dear Luis,
> I have removed getfem from my system and tried to install the latest
> version. I am having the same problem I had back in March, namely that the
> gmake goes looking for numarray and calls an error because I have numpy,
> which it doesn't find. The first time I installed it Yves Renard emailed
> me an archive which was what I built the installation from. So that is the
> reason that my installation was not the stable release.
>
> So, what do I need to do to make the latest stable release install and
> work properly with numpy, rather than numarray?
>
> Regards
>
> Michael
>
>   
> > Michael,
> >
> > the 2931 release is not the stable release, and its python
> > interface contains a bug (fixed in current release).
> > Could you use the latest version and try again?
> >
> > regards,
> > Luis
> >
> >
> > address@hidden escribió:
> >     
> > > Thanks Luis,
> > > According to the getfem.py file it is  __version__ = "$Revision: 2931
> > > $".
> > > I'm pretty sure its version 3.1, the latest stable release when I
> > > downloaded it in early March of this year.
> > > Below is the relevant part of the code. I have left out the construction
> > > of the mesh m, which is working properly.
> > > thanks
> > > Michael
> > >
> > > fem = MeshFem(m,1)
> > > fem.set_fem(Fem('FEM_PK(2,1)'))
> > >
> > > ##################
> > > #
> > > print "integration method"
> > > #
> > > ##################
> > > mim = MeshIm(m)
> > > mim.set_integ(Integ('IM_TRIANGLE(1)'))
> > >
> > > ##################
> > > #
> > > print "bricks"
> > > #
> > > ##################
> > > #B = MdBrick('bilaplacian',mim,fem)
> > > B = MdBrick('generic_elliptic',mim,fem)
> > >
> > > m.set_region(41,m.faces_from_cvid(face_trans[50]))
> > > m.set_region(42,m.faces_from_cvid(face_trans[51]))
> > >
> > > C = MdBrick('dirichlet',B,41,fem,'penalized')
> > > C.set_param('R',array([0.]) )
> > >
> > > D = MdBrick('dirichlet',C,42,fem,'penalized')
> > > D.set_param('R',array([0.]) )
> > >
> > >
> > > mds = MdState('real')
> > > D.solve(mds)
> > > state = mds.state()
> > >
> > >
> > >       
> > > > Hello Michael,
> > > >
> > > > Could you show me the code of your python-script and your version of
> > > > getfem?
> > > >
> > > > regards,
> > > > Luis.
> > > >
> > > > address@hidden escribió:
> > > >
> > > >         
> > > > > Hello again,
> > > > > I've been poking around the workings of the python interface and found
> > > > > that the solvers that getfem couldn't find aren't supported by getfem,
> > > > > which is fine but those it can find
> > > > > a) just won't solve a simple diffusion equation
> > > > > b) but solve it without difficulty on a colleague's computer. The main
> > > > > differences between his system and mine are 1) numarray on his and
> > > > > numpy
> > > > > on mine (yes I did make the necessary changes to the import statement)
> > > > > and
> > > > > 2) ubuntu on his and fedora on mine.
> > > > >
> > > > > I have yet to verify if the solvers work outside of getfem, but in the
> > > > > meantime, does anyone have any ideas about what may cause a solver to
> > > > > fail?
> > > > >
> > > > > superlu - fails to solve and program crashes with error 'SuperLU
> > > > > failed
> > > > > to
> > > > > solve info=9
> > > > > cg/ildlt - doesn't crash but doesn't solve either
> > > > > gmres/ilu[t[p]] - program hangs
> > > > >
> > > > >
> > > > > Michael Walker
> > > > >
> > > > > Plant Modelling Group
> > > > > CIRAD, Montpellier
> > > > > 04 67 61 57 27
> > > > >
> > > > > _______________________________________________
> > > > > Getfem-users mailing list
> > > > > address@hidden
> > > > > https://mail.gna.org/listinfo/getfem-users
> > > > >
> > > > >
> > > > >           
> > > > _______________________________________________
> > > > Getfem-users mailing list
> > > > address@hidden
> > > > https://mail.gna.org/listinfo/getfem-users
> > > >
> > > >
> > > >         
> > > Michael Walker
> > >
> > > Plant Modelling Group
> > > CIRAD, Montpellier
> > > 04 67 61 57 27
> > >
> > >       
> >     
>
>
> Michael Walker
>
> Plant Modelling Group
> CIRAD, Montpellier
> 04 67 61 57 27
>
> _______________________________________________
> Getfem-users mailing list
> address@hidden
> https://mail.gna.org/listinfo/getfem-users
>

ronan_perrussel.vcf
Description: Vcard