Resolution { { Name analysis; System { { Name resol_getdpThermal; NameOfFormulation getdpThermal_formulation; } { Name resol_getdpEMGeneral3D; NameOfFormulation getdpEMGeneral3D_formulation; Type ComplexValue; Frequency Freq;} } Operation { DeleteFile["../getdpresults/times.dat"]; InitSolution[resol_getdpThermal] ; InitSolution[resol_getdpEMGeneral3D] ; Test[ ($TimeStep<=1)] { UpdateConstraint[resol_getdpEMGeneral3D, Region[All], Assign ]; Print ["Solving formulation [getdpEMGeneral3D]"]; Generate[resol_getdpEMGeneral3D] ; Solve[resol_getdpEMGeneral3D]; Evaluate[$value4 = 0]; Evaluate[$value2 = 0]; } { Evaluate[$value4 = $value4 + #4]; Evaluate[$value2 = $value2 + #2]; Print[{#4*100}, Format 'Maximal change of Magnetic Permeability is %f ...']; Print[{$value4*100}, Format "Accumulated change of Magnetic Permeability is %f ..."]; Print[{#2*100}, Format 'Maximal change of Electrical Conductivity is %f ...']; Print[{$value2*100}, Format "Accumulated change of Electrical Conductivity is %f ..."]; Test[($value4 > 1.000000e-01) || ($value2 > 1.000000e-01) ] { UpdateConstraint[resol_getdpEMGeneral3D, Region[All], Assign ]; Print ["Solving formulation [getdpEMGeneral3D]"]; Generate[resol_getdpEMGeneral3D] ; Solve[resol_getdpEMGeneral3D]; Evaluate[$value4 = 0]; Evaluate[$value2 = 0]; } { Generate[resol_getdpEMGeneral3D] ; Print["Skipping getdpEMGeneral3D..."]; } } Print ["Solving formulation [getdpThermal]"]; IterativeLoop [ 50, 0.001000, 0.200000] { UpdateConstraint[resol_getdpThermal, Region[All], Assign ]; Print ["Solving formulation [getdpThermal]"]; GenerateJac[resol_getdpThermal]; SolveJac[resol_getdpThermal]; } SaveSolution[resol_getdpThermal]; SaveSolution[resol_getdpEMGeneral3D]; Print[{$Time}, File "../getdpresults/times.dat"]; PostOperation[postOp_getdpThermal]; PostOperation[postOp_getdpEMGeneral3D]; } } }