<div dir="ltr"><div><div>Dear all<br><br></div>This addition is not harmful in my point
 of view. I just want to comment a bit on that. Basically in
 the case of structural solver, the new SolveSolutionStep function will comprise of:<br><br></div>Clear the linear system<br><br>loop until convergence:<br>  scheme.InitializeNonLinIteration<br>  builder_and_solver.Build<br>  linear_solver.Solve<br><div>  scheme.FinalizeNonLinIteration<br><br></div><div>We wrap it under the name ExecuteIteration anyway (see uzawa_contact_strategy.py). Therefore it&#39;s just a name change.<br><br><br></div><div>Cheers,<br></div><div>Bui<br></div><div>P/S:
 Is there a bounding volume tree implementation for efficient contact 
search already in the Kratos kernel/ else where in the application? I 
just want to know if I can use a better implementation, otherwise I will
 contribute my one to the repository.<br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Oct 7, 2015 at 11:57 AM, Riccardo Rossi <span dir="ltr">&lt;<a href="mailto:rrossi@cimne.upc.edu" target="_blank">rrossi@cimne.upc.edu</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div><div><div><div><div>Dear Nelson,<br><br></div> the point is the following:<br><br></div>imagine that you do a FSI simulation and your structure is elastoplastic: you first have a  wrong guess (too high)  of the forces acting on the structure and you call &quot;solve&quot;. The structure becomes plastic and after reaching convergence it calls FinalizeSolutionStep where the constitutive law stores the internal variables for the next step. <br><br></div></div>then you correct the forces (now much smaller) and you call solve again. The structure should not become plastic anymore HOWEVER it already stored the internal variables at the end of the previous step, so the results would wrongly believe that the structural behaviour is plastic, while really should have stayed elastic<br><br></div>still... note that if you call &quot;Solve&quot; it will behave only as before. The difference is that you can choose when your Internal variables have to be stored (at the end of the complete fsi loop, not at the end of each solve!!)<br><br></div>cheers<br></div><div>Riccardo<br></div><div><div><div><div><br><br></div></div></div></div></div><div class="gmail_extra"><div><div class="h5"><br><div class="gmail_quote">On Wed, Oct 7, 2015 at 9:40 AM,  <span dir="ltr">&lt;<a href="mailto:maireni@cimne.upc.edu" target="_blank">maireni@cimne.upc.edu</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Riccardo,<br>
<br>
I do not understrand yet why it is not necessary to call the initialize,<br>
predict and finalize functions.<br>
This methods are escencially important in structural and solid<br>
application where some internal variables<br>
or other need to be updated or reseted (mainly  in constitutive law).<br>
If bassically  does not want to call this method i prefere in my opinion<br>
to include a new member in base<br>
class of strategy where those calls can be avoided.<br>
<br>
It really cost calls this function in each time step?<br>
<br>
Regards<br>
<br>
Nelson<br>
<div><div><br>
<br>
<br>
El 2015-10-07 15:05, Riccardo Rossi escribió:<br>
&gt; Dear All,<br>
&gt;<br>
&gt;   i believe there is a weakness in the current design of the solving<br>
&gt; strategies, and i would like to propose a small improvement to the<br>
&gt; current interface:<br>
&gt;<br>
&gt; the point is that as of today the strategy has essentially the<br>
&gt; following methods:<br>
&gt;<br>
&gt; Initialize (to be called once)<br>
&gt; InitializeSolutionStep<br>
&gt; Predict<br>
&gt; FinalizeSolutionStep<br>
&gt;<br>
&gt; and<br>
&gt;<br>
&gt; Solve --&gt; which combines everything in one single call.<br>
&gt;<br>
&gt; the problem is that is someone calls solve the InitializeSolutionStep<br>
&gt; and FinalizeSolutionStep are inevitably called, which is a problem if<br>
&gt; for example one<br>
&gt; would like to do repeated solves within a single time step.<br>
&gt;<br>
&gt; the typical scenario is that of a FSI simulation, but the same problem<br>
&gt; is also found for example in the RVE solver.<br>
&gt;<br>
&gt; My proposal would be to add a new function<br>
&gt;<br>
&gt; &quot;SolveSolutionStep&quot;<br>
&gt;<br>
&gt; which ONLY does the solution but WITHOUT CALLING<br>
&gt; InitializeSolutionStep, Predict, FinalizeSolutionStep<br>
&gt;<br>
&gt; The point of doing this is that the &quot;solve&quot; call would be:<br>
&gt;<br>
&gt; Solve()<br>
&gt;     InitializeSolutionStep<br>
&gt;     Predict<br>
&gt;     SolveSolutionStep<br>
&gt;     FinalizeSolutionStep<br>
&gt;<br>
&gt; but one would be able to call the very same thing from the python (or<br>
&gt; from a higher level strategy)<br>
&gt;<br>
&gt; so that for example a FSI loop would become<br>
&gt;<br>
&gt;     structural_solver.InitializeSolutionStep()<br>
&gt;     fluid_solver.InitializeSolutionStep()<br>
&gt;<br>
&gt;     structural_solver.Predict()<br>
&gt;<br>
&gt;     while(not converged)<br>
&gt;         structural_solver.SolveSolutionStep()<br>
&gt;<br>
&gt;         apply_velocities_to_fluid()<br>
&gt;<br>
&gt;         fluid_solver.SolveSolutionStep()<br>
&gt;<br>
&gt;         apply_forces_to_structure()<br>
&gt;<br>
&gt;         check_convergence()<br>
&gt;<br>
&gt;     structural_solver.FinalizeSolutionStep()<br>
&gt;     fluid_solver.FinalizeSolutionStep()<br>
&gt;    <br>
&gt;<br>
&gt; Of course if someone uses the current Solve() nothing changes and<br>
&gt; current behaviour is unmodified...<br>
&gt;<br>
&gt; any voice against such change? if not i&#39;ll go ahead...<br>
&gt;<br>
&gt; ciao<br>
&gt; Riccardo<br>
&gt;<br>
&gt;        <br>
&gt;       <br>
&gt;<br>
&gt; --<br>
&gt;<br>
&gt; Riccardo Rossi<br>
&gt;<br>
&gt;  PhD, Civil Engineer<br>
&gt;<br>
</div></div>&gt; member of the Kratos Team: <a href="http://www.cimne.com/kratos" rel="noreferrer" target="_blank">www.cimne.com/kratos</a> [1]<br>
<span>&gt;<br>
&gt; lecturer at Universitat Politècnica de Catalunya, BarcelonaTech (UPC)<br>
&gt;<br>
&gt; Research fellow at International Center for Numerical Methods in<br>
&gt; Engineering (CIMNE)<br>
&gt;<br>
&gt; C/ Gran Capità, s/n, Campus Nord UPC, Ed. C1, Despatx C9<br>
&gt;<br>
</span>&gt;  08034 – Barcelona – Spain – <a href="http://www.cimne.com" rel="noreferrer" target="_blank">www.cimne.com</a> [2]  - <br>
&gt;<br>
&gt;  T.<a href="tel:%28%2B34%29%2093%20401%2056%2096" value="+34934015696" target="_blank">(+34) 93 401 56 96</a> skype: ROUGERED4<br>
&gt;<br>
&gt;  <br>
&gt;<br>
&gt;  [3]<br>
&gt;<br>
&gt;  [4] [5] [6] [7] [8] [9]<br>
<span>&gt;<br>
&gt; Les dades personals contingudes en aquest missatge són tractades amb<br>
&gt; la finalitat de mantenir el contacte professional entre CIMNE i voste.<br>
&gt; Podra exercir els drets d&#39;accés, rectificació, cancel·lació i<br>
&gt; oposició, dirigint-se a <a href="mailto:cimne@cimne.upc.edu" target="_blank">cimne@cimne.upc.edu</a>. La utilització de la<br>
&gt; seva adreça de correu electronic per part de CIMNE queda subjecte a<br>
&gt; les disposicions de la Llei 34/2002, de Serveis de la Societat de la<br>
&gt; Informació i el Comerç Electronic.<br>
&gt;<br>
</span>&gt;  Imprimiu aquest missatge, només si és estrictament necessari. [3]<br>
&gt;<br>
&gt; Links:<br>
&gt; ------<br>
&gt; [1] <a href="http://www.cimne.com/kratos" rel="noreferrer" target="_blank">http://www.cimne.com/kratos</a><br>
&gt; [2] <a href="http://www.cimne.com" rel="noreferrer" target="_blank">http://www.cimne.com</a><br>
&gt; [3] <a href="http://www.cimne.com/" rel="noreferrer" target="_blank">http://www.cimne.com/</a><br>
&gt; [4] <a href="https://www.facebook.com/cimne" rel="noreferrer" target="_blank">https://www.facebook.com/cimne</a><br>
&gt; [5] <a href="http://blog.cimne.com/" rel="noreferrer" target="_blank">http://blog.cimne.com/</a><br>
&gt; [6] <a href="http://vimeo.com/cimne" rel="noreferrer" target="_blank">http://vimeo.com/cimne</a><br>
&gt; [7] <a href="http://www.youtube.com/user/CIMNEvideos" rel="noreferrer" target="_blank">http://www.youtube.com/user/CIMNEvideos</a><br>
&gt; [8] <a href="http://www.linkedin.com/company/cimne" rel="noreferrer" target="_blank">http://www.linkedin.com/company/cimne</a><br>
&gt; [9] <a href="https://twitter.com/cimne" rel="noreferrer" target="_blank">https://twitter.com/cimne</a><br>
&gt;<br>
&gt; _______________________________________________<br>
&gt; Kratos mailing list<br>
&gt; <a href="mailto:Kratos@listas.cimne.upc.edu" target="_blank">Kratos@listas.cimne.upc.edu</a><br>
&gt; <a href="http://listas.cimne.upc.edu/cgi-bin/mailman/listinfo/kratos" rel="noreferrer" target="_blank">http://listas.cimne.upc.edu/cgi-bin/mailman/listinfo/kratos</a><br>
_______________________________________________<br>
Kratos mailing list<br>
<a href="mailto:Kratos@listas.cimne.upc.edu" target="_blank">Kratos@listas.cimne.upc.edu</a><br>
<a href="http://listas.cimne.upc.edu/cgi-bin/mailman/listinfo/kratos" rel="noreferrer" target="_blank">http://listas.cimne.upc.edu/cgi-bin/mailman/listinfo/kratos</a><br>
</blockquote></div><br><br clear="all"><br>-- <br></div></div><div><div dir="ltr"><p></p><div><div dir="ltr"><div><div class="h5"><p style="padding:0px;margin:0px;font-family:Arial,Helvetica,sans-serif;color:rgb(0,73,150);line-height:18px"><b>Riccardo Rossi<br></b></p><p style="padding:0px;margin:0px;font-family:Arial,Helvetica,sans-serif;color:rgb(120,120,120);font-size:12px;line-height:18px">

PhD, Civil Engineer</p><p style="padding:0px;margin:0px;font-family:Arial,Helvetica,sans-serif;color:rgb(120,120,120);font-size:12px;line-height:18px"><br></p><p style="padding:0px;margin:0px;font-family:Arial,Helvetica,sans-serif;color:rgb(120,120,120);font-size:12px;line-height:18px">member of the Kratos Team: <a href="http://www.cimne.com/kratos" target="_blank">www.cimne.com/kratos</a><br></p><p style="padding:0px;margin:0px"><span style="color:rgb(120,120,120);font-family:Arial,Helvetica,sans-serif;font-size:12px;line-height:18px">lecturer at Universitat Politècnica de Catalunya, BarcelonaTech (UPC)</span><br>

</p><p style="padding:0px;margin:0px;font-family:Arial,Helvetica,sans-serif;color:rgb(120,120,120);font-size:12px;line-height:18px">Research fellow at International Center for Numerical Methods in Engineering (CIMNE)</p>
<p style="padding:0px;margin:0px;font-family:Arial,Helvetica,sans-serif;color:rgb(120,120,120);font-size:12px;line-height:18px">

<br></p><p style="padding:0px;margin:0px;font-family:Arial,Helvetica,sans-serif;color:rgb(120,120,120);font-size:12px;line-height:18px">C/ Gran Capità, s/n, Campus Nord UPC, Ed. C1, Despatx C9</p></div></div><p style="padding:0px;margin:0px;font-family:Arial,Helvetica,sans-serif;color:rgb(120,120,120);font-size:12px;line-height:18px">

08034 – Barcelona – Spain – <a href="http://www.cimne.com" target="_blank">www.cimne.com</a>  - </p><span class=""><p style="padding:0px;margin:0px;font-family:Arial,Helvetica,sans-serif;color:rgb(120,120,120);font-size:12px;line-height:18px">

<span>T</span>.<a value="+34934010794">(+34) 93 401 56 96</a> skype: <b>rougered4</b><br></p><p style="padding:0px;margin:0px;font-family:Arial,Helvetica,sans-serif;color:rgb(120,120,120);font-size:12px;line-height:18px"> </p><p style="padding:0px;margin:0px;font-family:Arial,Helvetica,sans-serif;color:rgb(120,120,120);font-size:12px;line-height:11px">

<a href="http://www.cimne.com/" target="_blank"><img src="https://web.cimne.upc.edu/groups/publicacions/signatura/logo_cimne_negre.png" height="17" width="80" border="0"></a></p><p style="font-size:medium;font-family:&#39;Times New Roman&#39;">

<a href="https://www.facebook.com/cimne" target="_blank"><img src="https://web.cimne.upc.edu/groups/publicacions/signatura/facebook-logo.png" height="24" width="24" border="0"></a><a href="http://blog.cimne.com/" target="_blank"><img src="https://web.cimne.upc.edu/groups/publicacions/signatura/wordpress.png" height="24" width="24" border="0"></a><a href="http://vimeo.com/cimne" target="_blank"><img src="https://web.cimne.upc.edu/groups/publicacions/signatura/vimeo.png" height="24" width="24" border="0"></a><a href="http://www.youtube.com/user/CIMNEvideos" target="_blank"><img src="https://web.cimne.upc.edu/groups/publicacions/signatura/you-tube1.png" height="24" width="24" border="0"></a><a href="http://www.linkedin.com/company/cimne" target="_blank"><img src="https://web.cimne.upc.edu/groups/publicacions/signatura/linkedin-logo.png" height="24" width="24" border="0"></a><a href="https://twitter.com/cimne" target="_blank"><img src="https://web.cimne.upc.edu/groups/publicacions/signatura/twitter.png" height="24" width="24" border="0"></a></p>

</span><p style="padding:0px;margin:0px;font-family:Arial,Helvetica,sans-serif;color:rgb(155,155,155);font-size:10px;line-height:14px">Les
 dades personals contingudes en aquest missatge són tractades amb la 
finalitat de mantenir el contacte professional entre CIMNE i voste. Podra exercir els drets d&#39;accés, rectificació, cancel·lació i oposició, dirigint-se a <a href="mailto:cimne@cimne.upc.edu" target="_blank">cimne@cimne.upc.edu</a>. La utilització de la seva adreça de correu electronic
 per part de CIMNE queda subjecte a les disposicions de la Llei 34/2002,
 de Serveis de la Societat de la Informació i el Comerç Electronic.</p>

<p style="padding:0px;margin:0px;font-family:Arial,Helvetica,sans-serif;color:rgb(4,155,2);font-size:10px;line-height:18px"><img src="https://web.cimne.upc.edu/groups/publicacions/signatura/image002.png" height="17" width="20"> Imprimiu aquest missatge, només si és estrictament necessari.</p>

</div></div><span lang="ES"><a href="http://www.cimne.com/" target="_blank"><span style="font-size:10.0pt;font-family:&quot;Arial&quot;,&quot;sans-serif&quot;"></span></a></span></div></div>
</div>
<br>_______________________________________________<br>
Kratos mailing list<br>
<a href="mailto:Kratos@listas.cimne.upc.edu">Kratos@listas.cimne.upc.edu</a><br>
<a href="http://listas.cimne.upc.edu/cgi-bin/mailman/listinfo/kratos" rel="noreferrer" target="_blank">http://listas.cimne.upc.edu/cgi-bin/mailman/listinfo/kratos</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature">With Best Regards !<br>Giang Bui<div>To learn and to excel</div></div>
</div>