<html><head></head><body><div style="color:#000; background-color:#fff; font-family:verdana, helvetica, sans-serif;font-size:13px"><div id="yui_3_16_0_1_1455113157552_21799">Thanks for your reply.</div><div id="yui_3_16_0_1_1455113157552_21801"><br></div><div id="yui_3_16_0_1_1455113157552_21803" dir="ltr">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.</div><div dir="ltr"><br></div><div id="yui_3_16_0_1_1455113157552_21806" class="qtdSeparateBR"><div>Bests,</div><div>FM<br></div><br></div><div style="display: block;" id="yui_3_16_0_1_1455113157552_21811" class="yahoo_quoted"> <div id="yui_3_16_0_1_1455113157552_21810" style="font-family: verdana, helvetica, sans-serif; font-size: 13px;"> <div id="yui_3_16_0_1_1455113157552_21809" style="font-family: HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 16px;"> <div id="yui_3_16_0_1_1455113157552_21808" dir="ltr"> <font id="yui_3_16_0_1_1455113157552_21807" face="Arial" size="2"> <hr id="yui_3_16_0_1_1455113157552_21935" size="1"> <b id="yui_3_16_0_1_1455113157552_21938"><span id="yui_3_16_0_1_1455113157552_21937" style="font-weight:bold;">From:</span></b> Vincent Huber <vincent.huber@cemosis.fr><br> <b><span style="font-weight: bold;">To:</span></b> Farshid Mossaiby <mossaiby@yahoo.com> <br><b><span style="font-weight: bold;">Cc:</span></b> PETSc Users List <petsc-users@mcs.anl.gov><br> <b><span style="font-weight: bold;">Sent:</span></b> Thursday, February 11, 2016 1:33 PM<br> <b><span style="font-weight: bold;">Subject:</span></b> Re: [petsc-users] Parallel FEM code using PETSc<br> </font> </div> <div id="yui_3_16_0_1_1455113157552_22068" class="y_msg_container"><br><div id="yiv6496152878"><div id="yui_3_16_0_1_1455113157552_22067"><div id="yui_3_16_0_1_1455113157552_22066" dir="ltr"><div id="yui_3_16_0_1_1455113157552_22065" class="yiv6496152878markdown-here-wrapper" style="font-family:'Courier New';"><div style="margin:1.2em 0px!important;">Hello,</div>
<div style="margin:1.2em 0px!important;">I’m not convinced the following will be find <code style="font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;font-size:inherit;line-height:inherit;font-family:inherit;margin:0px 0.15em;padding:0px 0.15em;white-space:pre-wrap;border:1px solid rgb(234,234,234);border-radius:3px;display:inline;background-color:rgb(248,248,248);">simple</code> enough but <a rel="nofollow" shape="rect" target="_blank" href="http://www.feelpp.org/">Feel++</a> match the other expectations.</div>
<div style="margin:1.2em 0px!important;">In the example on the main page of the web site, you only have to change the line</div>
<pre style="font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;font-size:inherit;line-height:inherit;font-family:inherit;font-size:1em;line-height:1.2em;margin:1.2em 0px;"><code style="font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;font-size:inherit;line-height:inherit;font-family:inherit;margin:0px 0.15em;padding:0px 0.15em;white-space:pre-wrap;border:1px solid rgb(234,234,234);border-radius:3px;display:inline;background-color:rgb(248,248,248);white-space:pre;overflow:auto;border-radius:3px;border:1px solid rgb(204,204,204);padding:0.5em 0.7em;display:block!important;display:block;padding:0.5em;color:rgb(51,51,51);background:rgb(248,248,248);"> a+=on(_range=boundaryfaces(mesh), _rhs=l, _element=u, _expr=cst(<span class="yiv6496152878hljs-number" style="color:rgb(0,128,128);">0.</span>) );
</code></pre>
<div style="margin:1.2em 0px!important;">to </div>
<pre style="font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;font-size:inherit;line-height:inherit;font-family:inherit;font-size:1em;line-height:1.2em;margin:1.2em 0px;"><code style="font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;font-size:inherit;line-height:inherit;font-family:inherit;margin:0px 0.15em;padding:0px 0.15em;white-space:pre-wrap;border:1px solid rgb(234,234,234);border-radius:3px;display:inline;background-color:rgb(248,248,248);white-space:pre;overflow:auto;border-radius:3px;border:1px solid rgb(204,204,204);padding:0.5em 0.7em;display:block!important;display:block;padding:0.5em;color:rgb(51,51,51);background:rgb(248,248,248);"> a+=on(_range=boundaryfaces(mesh), _rhs=l, _element=u, _expr=expr(soption(<span class="yiv6496152878hljs-string" style="color:rgb(221,17,68);">"functions.g"</span>)) );
</code></pre>
<div id="yui_3_16_0_1_1455113157552_22074" style="margin:1.2em 0px!important;">and run the code with</div>
<pre id="yui_3_16_0_1_1455113157552_22077" style="font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;font-size:inherit;line-height:inherit;font-family:inherit;font-size:1em;line-height:1.2em;margin:1.2em 0px;"><code id="yui_3_16_0_1_1455113157552_22076" class="yiv6496152878hljs yiv6496152878language-sh" style="font-style:inherit;font-variant:inherit;font-weight:inherit;font-stretch:inherit;font-size:inherit;line-height:inherit;font-family:inherit;margin:0px 0.15em;padding:0px 0.15em;white-space:pre-wrap;border:1px solid rgb(234,234,234);border-radius:3px;display:inline;background-color:rgb(248,248,248);white-space:pre;overflow:auto;border-radius:3px;border:1px solid rgb(204,204,204);padding:0.5em 0.7em;display:block!important;display:block;padding:0.5em;color:rgb(51,51,51);background:rgb(248,248,248);">mpirun -np xxx ./feelpp_theApp --functions.g=<span id="yui_3_16_0_1_1455113157552_22075" class="yiv6496152878hljs-string" style="color:rgb(221,17,68);">"2*x*y=x=y"</span>
</code></pre>
<div style="margin:1.2em 0px!important;">To change the bc you want to impose.</div>
<div style="margin:1.2em 0px!important;">I (highly) recommend Ubuntu 14.04 to test the library.<br clear="none">The documentation is <a rel="nofollow" shape="rect" target="_blank" href="http://book.feelpp.org/">here</a></div>
<div id="yui_3_16_0_1_1455113157552_22078" style="margin:1.2em 0px!important;">You may also have a look to <a rel="nofollow" shape="rect" target="_blank" href="http://www.freefem.org/">FreeFem++</a></div>
<div id="yui_3_16_0_1_1455113157552_22079" style="margin:1.2em 0px!important;">All the best,</div>
<div style="margin:1.2em 0px!important;">VH</div>
<div title-off="" title="MDH:PGRpdiBjbGFzcz0iZ21haWxfZGVmYXVsdCIgc3R5bGU9ImZvbnQtZmFtaWx5OiAnY291cmllciBu ZXcnLCBtb25vc3BhY2U7Ij5IZWxsbyw8L2Rpdj48ZGl2IGNsYXNzPSJnbWFpbF9kZWZhdWx0IiBz dHlsZT0iZm9udC1mYW1pbHk6ICdjb3VyaWVyIG5ldycsIG1vbm9zcGFjZTsiPjxicj48L2Rpdj48 ZGl2IGNsYXNzPSJnbWFpbF9kZWZhdWx0IiBzdHlsZT0iZm9udC1mYW1pbHk6ICdjb3VyaWVyIG5l dycsIG1vbm9zcGFjZTsiPkknbSBub3QgY29udmluY2VkIHRoZSBmb2xsb3dpbmcgd2lsbCBiZSBm aW5kIGBzaW1wbGVgIGVub3VnaCBidXQgW0ZlZWwrK10od3d3LmZlZWxwcC5vcmcpIG1hdGNoIHRo ZSBvdGhlciBleHBlY3RhdGlvbnMuPGJyPjwvZGl2PjxkaXYgY2xhc3M9ImdtYWlsX2RlZmF1bHQi IHN0eWxlPSJmb250LWZhbWlseTogJ2NvdXJpZXIgbmV3JywgbW9ub3NwYWNlOyI+PGJyPjwvZGl2 PjxkaXYgY2xhc3M9ImdtYWlsX2RlZmF1bHQiIHN0eWxlPSJmb250LWZhbWlseTogJ2NvdXJpZXIg bmV3JywgbW9ub3NwYWNlOyI+SW4gdGhlIGV4YW1wbGUgb24gdGhlIG1haW4gcGFnZSBvZiB0aGUg d2ViIHNpdGUsIHlvdSBvbmx5IGhhdmUgdG8gY2hhbmdlIHRoZSBsaW5lPC9kaXY+PGRpdiBjbGFz cz0iZ21haWxfZGVmYXVsdCIgc3R5bGU9ImZvbnQtZmFtaWx5OiAnY291cmllciBuZXcnLCBtb25v c3BhY2U7Ij5gYGBjKys8L2Rpdj48ZGl2IGNsYXNzPSJnbWFpbF9kZWZhdWx0IiBzdHlsZT0iIj48 ZGl2IGNsYXNzPSJnbWFpbF9kZWZhdWx0IiBzdHlsZT0iZm9udC1mYW1pbHk6ICdjb3VyaWVyIG5l dycsIG1vbm9zcGFjZTsiPiZuYnNwOyAmbmJzcDsgYSs9b24oX3JhbmdlPWJvdW5kYXJ5ZmFjZXMo bWVzaCksIF9yaHM9bCwgX2VsZW1lbnQ9dSwgX2V4cHI9Y3N0KDAuKSApOzwvZGl2PjxkaXYgc3R5 bGU9ImZvbnQtZmFtaWx5OiAnY291cmllciBuZXcnLCBtb25vc3BhY2U7Ij5gYGA8L2Rpdj48ZGl2 IHN0eWxlPSJmb250LWZhbWlseTogJ2NvdXJpZXIgbmV3JywgbW9ub3NwYWNlOyI+dG8mbmJzcDs8 L2Rpdj48ZGl2IHN0eWxlPSJmb250LWZhbWlseTogJ2NvdXJpZXIgbmV3JywgbW9ub3NwYWNlOyI+ YGBgYysrPC9kaXY+PGRpdiBzdHlsZT0iIj48ZGl2IHN0eWxlPSIiPjxmb250IGZhY2U9ImNvdXJp ZXIgbmV3LCBtb25vc3BhY2UiPiZuYnNwOyAmbmJzcDsgYSs9b24oX3JhbmdlPWJvdW5kYXJ5ZmFj ZXMobWVzaCksIF9yaHM9bCwgX2VsZW1lbnQ9dSwgX2V4cHI9ZXhwcihzb3B0aW9uKCJmdW5jdGlv bnMuZyIpKSApOzwvZm9udD48L2Rpdj48L2Rpdj48ZGl2IHN0eWxlPSJmb250LWZhbWlseTogJ2Nv dXJpZXIgbmV3JywgbW9ub3NwYWNlOyI+YGBgPC9kaXY+PGRpdiBzdHlsZT0iZm9udC1mYW1pbHk6 ICdjb3VyaWVyIG5ldycsIG1vbm9zcGFjZTsiPmFuZCBydW4gdGhlIGNvZGUgd2l0aDwvZGl2Pjxk aXYgc3R5bGU9ImZvbnQtZmFtaWx5OiAnY291cmllciBuZXcnLCBtb25vc3BhY2U7Ij48YnI+PC9k aXY+PGRpdiBzdHlsZT0iZm9udC1mYW1pbHk6ICdjb3VyaWVyIG5ldycsIG1vbm9zcGFjZTsiPmBg YHNoPC9kaXY+PGRpdiBzdHlsZT0iZm9udC1mYW1pbHk6ICdjb3VyaWVyIG5ldycsIG1vbm9zcGFj ZTsiPm1waXJ1biAtbnAgeHh4IC4vZmVlbHBwX3RoZUFwcCAtLWZ1bmN0aW9ucy5nPSIyKngqeT14 PXkiPC9kaXY+PGRpdiBzdHlsZT0iZm9udC1mYW1pbHk6ICdjb3VyaWVyIG5ldycsIG1vbm9zcGFj ZTsiPmBgYDwvZGl2PjxkaXYgc3R5bGU9ImZvbnQtZmFtaWx5OiAnY291cmllciBuZXcnLCBtb25v c3BhY2U7Ij48YnI+PC9kaXY+PGRpdiBzdHlsZT0iZm9udC1mYW1pbHk6ICdjb3VyaWVyIG5ldycs IG1vbm9zcGFjZTsiPlRvIGNoYW5nZSB0aGUgYmMgeW91IHdhbnQgdG8gaW1wb3NlLjwvZGl2Pjxk aXYgc3R5bGU9ImZvbnQtZmFtaWx5OiAnY291cmllciBuZXcnLCBtb25vc3BhY2U7Ij48YnI+PC9k aXY+PGRpdiBzdHlsZT0iZm9udC1mYW1pbHk6ICdjb3VyaWVyIG5ldycsIG1vbm9zcGFjZTsiPkkg KGhpZ2hseSkgcmVjb21tZW5kIFVidW50dSAxNC4wNCB0byB0ZXN0IHRoZSBsaWJyYXJ5LjwvZGl2 PjxkaXYgc3R5bGU9ImZvbnQtZmFtaWx5OiAnY291cmllciBuZXcnLCBtb25vc3BhY2U7Ij5UaGUg ZG9jdW1lbnRhdGlvbiBpcyBbaGVyZV0oaHR0cDovL2Jvb2suZmVlbHBwLm9yZy8pPC9kaXY+PGRp diBzdHlsZT0iZm9udC1mYW1pbHk6ICdjb3VyaWVyIG5ldycsIG1vbm9zcGFjZTsiPjxicj48L2Rp dj48ZGl2IHN0eWxlPSJmb250LWZhbWlseTogJ2NvdXJpZXIgbmV3JywgbW9ub3NwYWNlOyI+WW91 IG1heSBhbHNvIGhhdmUgYSBsb29rIHRvIFtGcmVlRmVtKytdKGh0dHA6Ly93d3cuZnJlZWZlbS5v cmcvKTwvZGl2PjxkaXYgc3R5bGU9ImZvbnQtZmFtaWx5OiAnY291cmllciBuZXcnLCBtb25vc3Bh Y2U7Ij48YnI+PC9kaXY+PGRpdiBzdHlsZT0iZm9udC1mYW1pbHk6ICdjb3VyaWVyIG5ldycsIG1v bm9zcGFjZTsiPkFsbCB0aGUgYmVzdCw8L2Rpdj48ZGl2IHN0eWxlPSJmb250LWZhbWlseTogJ2Nv dXJpZXIgbmV3JywgbW9ub3NwYWNlOyI+PGJyPjwvZGl2PjxkaXYgc3R5bGU9ImZvbnQtZmFtaWx5 OiAnY291cmllciBuZXcnLCBtb25vc3BhY2U7Ij5WSDwvZGl2PjwvZGl2Pg==" style="height:0;width:0;max-height:0;max-width:0;overflow:hidden;font-size:0em;padding:0;margin:0;"></div></div><div class="yiv6496152878gmail_extra"><br clear="none"><div class="yiv6496152878gmail_quote">2016-02-11 10:18 GMT+01:00 Farshid Mossaiby <span dir="ltr"><<a rel="nofollow" shape="rect" ymailto="mailto:mossaiby@yahoo.com" target="_blank" href="mailto:mossaiby@yahoo.com">mossaiby@yahoo.com</a>></span>:<br clear="none"><div class="yiv6496152878yqt1386111815" id="yiv6496152878yqt41736"><blockquote class="yiv6496152878gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div><div style="color:#000;background-color:#fff;font-family:verdana, helvetica, sans-serif;font-size:13px;"><div>Hi all,</div><div><br clear="none"></div><div dir="ltr">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.</div><div dir="ltr"><br clear="none"></div><div dir="ltr">TIA<br clear="none"></div></div></div></blockquote></div></div><br clear="none"><br clear="all"><div><br clear="none"></div>-- <br clear="none"><div class="yiv6496152878gmail_signature"><div dir="ltr"><font face="courier new, monospace">Docteur Ingénieur de recherche</font><div><font face="courier new, monospace"><a rel="nofollow" shape="rect" target="_blank" href="http://www.cemosis.fr/">CeMoSiS</a> - <a rel="nofollow" shape="rect" ymailto="mailto:vincent.huber@cemosis.fr" target="_blank" href="mailto:vincent.huber@cemosis.fr">vincent.huber@cemosis.fr</a></font></div><div><font face="courier new, monospace">Tel: +33 (0)3 68 8<b>5 02 06</b></font></div><div><font face="courier new, monospace">IRMA - 7, rue René Descartes</font></div><div><font face="courier new, monospace">67 000 Strasbourg</font></div><div><br clear="none"></div><div><img data-id="aef3fac3-ef4f-0219-3adc-a43590192fc9" src="http://www.google.com/a/cpanel/cemosis.fr/images/logo.gif"></div></div></div>
</div></div></div></div><br><br></div> </div> </div> </div></div></body></html>