Skip to content
Snippets Groups Projects
Commit 406b1eda authored by François Henrotte's avatar François Henrotte
Browse files

updates and correction

parent 87ce5349
No related branches found
No related tags found
No related merge requests found
Pipeline #5916 passed
Include "wires_common.pro"; Include "w2d_common.pro";
// global Gmsh options // global Gmsh options
Mesh.Optimize = 1; // optimize quality of tetrahedra Mesh.Optimize = 1; // optimize quality of tetrahedra
......
Include "wires_common.pro"; Include "w2d_common.pro";
Struct S::SOLVER [ Enum, mumps, gmres, fgmres,bcgs ]; Struct S::SOLVER [ Enum, mumps, gmres, fgmres,bcgs ];
DefineConstant[ DefineConstant[
...@@ -132,8 +132,8 @@ Function{ ...@@ -132,8 +132,8 @@ Function{
((R_2[]<rs) ? WI_2*Analytic_A[R_2[]] : 0 )); ((R_2[]<rs) ? WI_2*Analytic_A[R_2[]] : 0 ));
EndIf EndIf
If ( NumWires == 3 ) If ( NumWires == 3 )
Exact_a[] = WI_1*Analytic_A[R_1[]] + WI_2*Analytic_A[R_2[]] Exact_A[] = WI_1*Analytic_A[R_1[]] + WI_2*Analytic_A[R_2[]]
+ WI_3*Analytic_A[R_3[]] + (WI1_+WI_2+WI_3)*mu0/(2*Pi) * NumWires * Log[rout/rs] ; + WI_3*Analytic_A[R_3[]] + (WI_1+WI_2+WI_3)*mu0/(2*Pi) * NumWires * Log[rout/rs] ;
Correction_A[] = ((R_1[]<rs) ? WI_1*Analytic_A[R_1[]] : Correction_A[] = ((R_1[]<rs) ? WI_1*Analytic_A[R_1[]] :
((R_2[]<rs) ? WI_2*Analytic_A[R_2[]] : ((R_2[]<rs) ? WI_2*Analytic_A[R_2[]] :
((R_3[]<rs) ? WI_3*Analytic_A[R_3[]] : 0 ))); ((R_3[]<rs) ? WI_3*Analytic_A[R_3[]] : 0 )));
......
File moved
...@@ -578,11 +578,13 @@ PostOperation map UsingPost MagnetoDynamics { ...@@ -578,11 +578,13 @@ PostOperation map UsingPost MagnetoDynamics {
Print[ bs, OnElementsOf Vol_nu, File "bs.pos"]; Print[ bs, OnElementsOf Vol_nu, File "bs.pos"];
PrintGroup[ EdgesOf[ Vol_nu, ConnectedTo LWIRES ], PrintGroup[ EdgesOf[ Vol_nu, ConnectedTo LWIRES ],
In Vol_nu, File "group.pos"]; In Vol_nu, File "group.pos"];
PrintGroup[ _CO_Entity_71, In Vol_nu, File "Tree.pos"];
EndIf EndIf
PrintGroup[ EdgesOfTreeIn[ { Vol_Tree }, StartingOn { Sur_Tree, Lin_Tree } ],
In Vol_nu, File "Tree.pos"];
//PrintGroup[ _CO_Entity_44, In Vol_nu, File "Tree.pos"]; // PrintGroup[ EdgesOfTreeIn[ { Vol_Tree }, StartingOn { Sur_Tree, Lin_Tree } ],
// In Vol_nu, File "Tree.pos"];
} }
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment