[petsc-users] Parallel FEM code using PETSc

Farshid Mossaiby mossaiby at yahoo.com
Thu Feb 11 05:22:37 CST 2016

Thanks for your reply.
Indeed I want to expand the code to non-linear cases and use it as a study tool for some other researches. This requires plain C++ code to be available so that I can change different parts of it easily. Using a complex library unfortunately does not fit my needs for this research.

      From: Vincent Huber <vincent.huber at cemosis.fr>
 To: Farshid Mossaiby <mossaiby at yahoo.com> 
Cc: PETSc Users List <petsc-users at mcs.anl.gov>
 Sent: Thursday, February 11, 2016 1:33 PM
 Subject: Re: [petsc-users] Parallel FEM code using PETSc
Hello,I’m not convinced the following will be find simple enough but Feel++ match the other expectations.In the example on the main page of the web site, you only have to change the line    a+=on(_range=boundaryfaces(mesh), _rhs=l, _element=u, _expr=cst(0.) );
to     a+=on(_range=boundaryfaces(mesh), _rhs=l, _element=u, _expr=expr(soption("functions.g")) );
and run the code withmpirun -np xxx ./feelpp_theApp --functions.g="2*x*y=x=y"
To change the bc you want to impose.I (highly) recommend Ubuntu 14.04 to test the library.
The documentation is hereYou may also have a look to FreeFem++All the best,VH​
2016-02-11 10:18 GMT+01:00 Farshid Mossaiby <mossaiby at yahoo.com>:

Hi all,
Can somebody point me to a simple, linear, parallel FEM code based on PETSc, implementing efficient application of non-homogeneous Dirichlet BC? Domain decomposition code will is also welcome.

Docteur Ingénieur de rechercheCeMoSiS - vincent.huber at cemosis.frTel: +33 (0)3 68 85 02 06IRMA - 7, rue René Descartes67 000 Strasbourg

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20160211/c93bfc76/attachment-0001.html>

More information about the petsc-users mailing list