<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="">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><blockquote type="cite" class=""><div class="">On Aug 13, 2015, at 9:34 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=""><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>
_______________________________________________<br class="">getdp mailing list<br class=""><a href="mailto:getdp@geuz.org" class="">getdp@geuz.org</a><br class="">http://www.geuz.org/mailman/listinfo/getdp<br class=""></div></blockquote></div><br class=""></div></body></html>