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

[Michael . Walker]

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