<html><head><meta http-equiv="Content-Type" content="text/html charset=us-ascii"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">I actually do not know. I am not yet fully proficient with Getup!<div class=""><br class=""></div><div class=""><br class=""></div><div class="">Best,</div><div class=""><br class=""></div><div class="">Frederic<br class=""><div class=""><br class=""></div><div class=""><br class=""><div><blockquote type="cite" class=""><div class="">On Aug 13, 2015, at 9:47 AM, Peter Kis <<a href="mailto:md2z34@gmail.com" class="">md2z34@gmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div dir="ltr" class=""><div class=""><div class=""><div class="">Dear Frederic, <br class=""><br class=""></div>my final goal is to calculate the current (I = int_coil(js - sigma[]*Dt[{Aphi})) at every time instance and update the current density in order to achieve a certain current value. <br class=""><br class="">Regarding your suggestion, is it possible to define such constraint which depends on the solution (Aphi)?<br class=""><br class=""></div>Best Regards<br class=""><br class=""></div>Peter <br class=""></div><div class="gmail_extra"><br class=""><div class="gmail_quote">On Thu, Aug 13, 2015 at 4:40 PM, Frederic Trillaud Pighi <span dir="ltr" class=""><<a href="mailto:ftrillaudp@pumas.iingen.unam.mx" target="_blank" class="">ftrillaudp@pumas.iingen.unam.mx</a>></span> wrote:<br class=""><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word" class="">Dear Peter,<div class=""><br class=""></div><div class="">You may adjust your TimeFunction to your need in the constraint. It is interpolated during the resolution.</div><div class=""><br class=""></div><div class="">Best,</div><div class=""><br class=""></div><div class="">Frederic </div><div class=""><br class=""><div class=""><blockquote type="cite" class=""><div class=""><div class="h5"><div class="">On Aug 13, 2015, at 9:34 AM, Peter Kis <<a href="mailto:md2z34@gmail.com" target="_blank" class="">md2z34@gmail.com</a>> wrote:</div><br class=""></div></div><div class=""><div class=""><div class="h5"><div dir="ltr" class=""><div class=""><div class=""><div class=""><div class=""><div class=""><div class="">Dear List Members,<br class=""><br class=""></div>I'm just wondering if it's possible to scale the current density from the time iteration loop. <br class=""><br class="">Here I introduced a parameter <b class="">current_factor</b>, <br class=""><br class=""><span style="font-family:monospace,monospace" class="">Function{ <br class=""> ......<br class=""> <b class="">current_factor</b> = 1.;
<br class="">}
</span><br class=""><br class=""></div>which appears in the formulation as a scaling factor for the current density:<br class=""><br class=""><span style="font-family:monospace,monospace" class="">Formulation{ // OK
<br class=""> {Name VectorPotentialFormulation; Type FemEquation;
<br class=""> Quantity{
<br class=""> {Name Aphi; Type Local; NameOfSpace Afs;}
<br class=""> {Name js; Type Local; NameOfSpace jfs;}
<br class=""> }
<br class=""> Equation{
<br class=""> Galerkin{ [1/mu0*Dof{Curl Aphi},{Curl Aphi}];
<br class=""> In Omega; Jacobian JVol; Integration I1;}
<br class=""> <br class=""> Galerkin{ DtDof[sigma[]*Dof{Aphi},{Aphi}];
<br class=""> In Omega; Jacobian JVol; Integration I1;}
<br class=""> <br class=""> Galerkin{ [-<b class="">current_factor</b>*Dof{js},{Aphi}];
<br class=""> In Omega_coil; Jacobian JVol; Integration I1;}
<br class=""> }
<br class=""> }
<br class="">}
</span><br class=""><br class=""></div>and I'd like to update this parameter during the time iteration:<br class=""><br class=""><span style="font-family:monospace,monospace" class="">Resolution{
<br class=""> {Name Solution;
<br class=""> System{
<br class=""> {Name Syst; NameOfFormulation VectorPotentialFormulation;}
<br class=""> }
<br class=""> Operation{
<br class=""> InitSolution[Syst];SaveSolution[Syst];
<br class=""> TimeLoopTheta[t0,tmax,dt,theta] {<br class=""> <b class="">current_factor</b> = 1/1.7;
// This has no effect!!<br class=""> Generate[Syst]; Solve[Syst]; SaveSolution[Syst];
<br class=""> }
<br class=""> }
<br class=""> }
<br class="">}</span><br class=""><br class=""></div>but this has no impact on the solution. Apparently the value of current_vector is still considered one as it's defined in Functions. <br class=""><br class="">Is there any way to update this parameter from the time loop? <br class=""><br class=""></div>Thanks for your help in advance.<br class=""><br class=""><br class=""></div>Peter<br class=""><div class=""><div class=""><br class=""></div></div></div></div></div>
_______________________________________________<br class="">getdp mailing list<br class=""><a href="mailto:getdp@geuz.org" target="_blank" class="">getdp@geuz.org</a><br class=""><a href="http://www.geuz.org/mailman/listinfo/getdp" target="_blank" class="">http://www.geuz.org/mailman/listinfo/getdp</a><br class=""></div></blockquote></div><br class=""></div></div></blockquote></div><br class=""></div>
</div></blockquote></div><br class=""></div></div></body></html>