<html>
<body>
Hello Christophe,<br><br>
Thank you for your help!<br><br>
<blockquote type=cite class=cite cite>A correct definition of the
discrete vector potential involves the edges<br>
of the _complement_ of the spanning tree. The easiest way to do
this<br>
is to define the vector potential everywhere, and use the constraint
to<br>
fix the dofs to zero on the edges of the tree (so that only the edges
of<br>
the complement lead to actual unknowns).</blockquote><br>
If I understand well, this means:<br><br>
<font face="Courier New, Courier">FunctionSpace {<br>
{ Name Hcurl_a_Gauge; Type Form1;<br>
BasisFunction {<br>
// a = a s <br>
// e e
<br>
{ Name se; NameOfCoef ae; Function
BF_Edge;<br>
Support DomainC_Mag ; Entity
EdgesOf[All]; }<br>
}<br>
Constraint {<br>
//gauge condition ae = 0 on tree edges and on
some boundaries<br>
{ NameOfCoef ae;<br>
EntityType EdgesOfTreeIn;
EntitySubType StartingOn;<br>
NameOfConstraint
GaugeCondition_a_Mag_3D; }<br>
}<br>
}<br><br>
Constraint {<br>
{ Name GaugeCondition_a_Mag_3D; Type Assign;<br>
Case {<br>
{ Region DomainC_Mag ; SubRegion Surf;
Value 0.; }<br>
}<br>
}<br>
</font>Is this correct? If yes, probably something else is wrong
because<br>
unfortunately with this the problem does not converge:<br><br>
GetDP : 998 1.0825175e+06 9.4973969e-01<br>
GetDP : 999 1.0559604e+06 9.2643995e-01<br>
GetDP : Warning : Iterative solver has iterated too many
times<br>
GetDP : 1000 1.0548821e+06 9.2549394e-01<br>
GetDP : Solver : 1000 Iterations / Residual:
1.21005e+06<br><br>
What can I do?<br><br>
Florin<br><br>
<br>
<x-sigsep><p></x-sigsep>
Florin Ciuprina, Ph.D.<br>
Associate Professor<br>
----------------------------------------------------------------------------<br>
POLITEHNICA University of Bucharest<br>
Electrical Engineering Faculty<br>
Electrical Materials Laboratory<br>
Spl.Independentei 313, 77206, Bucharest, Romania<br>
Tel: +40.21.402 92 91<br>
Fax: +40.21.410 43 55<br>
e-mail: florin@elmat.pub.ro<br>
----------------------------------------------------------------------------</body>
</html>