Resolution 
{
	{
		Name analysis; 
		System 
		{
			{ Name resol_getdpThermal; NameOfFormulation getdpThermal_formulation; } 
			{ Name resol_getdpEMGeneral2D; NameOfFormulation getdpEMGeneral2D_formulation; Type ComplexValue; Frequency Freq;} 
		} 
		Operation 
		{ 
			DeleteFile["../getdpresults/times.dat"]; 
			InitSolution[resol_getdpThermal] ; 
			SaveSolution[resol_getdpThermal] ; 
			SetTime[ 0.000000e+00]; 
			Print[{$Time}, File "../getdpresults/times.dat"]; 
			InitSolution[resol_getdpEMGeneral2D] ; 
        	TimeLoopAdaptive[ 0.000000e+00, 5.000000e+00, 1.000000e+00, 1.000000e-08, 1.000000e+05, "Euler", List[Breakpoints], PostOperation { { adaptivePostOp_getdpThermal, 5.000000e-04, 1.000000e+00, MeanL1Norm  }  { adaptivePostOp_getdpEMGeneral2D, 5.000000e-04, 1.000000e+00, MeanL1Norm  } } ] 
 		   	{ 
				Test[ ($TimeStep<=1)] 
				{ 
					UpdateConstraint[resol_getdpEMGeneral2D, Region[All], Assign ]; 
					Print ["Solving formulation [getdpEMGeneral2D]"]; 
					Generate[resol_getdpEMGeneral2D] ;
					Solve[resol_getdpEMGeneral2D];
					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_getdpEMGeneral2D, Region[All], Assign ]; 
						Print ["Solving formulation [getdpEMGeneral2D]"]; 
						Generate[resol_getdpEMGeneral2D] ;
						Solve[resol_getdpEMGeneral2D];
						Evaluate[$value4 = 0];
						Evaluate[$value2 = 0];
					} 
					{ 
						Generate[resol_getdpEMGeneral2D] ;
						Print["Skipping getdpEMGeneral2D..."];
					} 
				} 

				Print ["Solving formulation [getdpThermal]"]; 
				IterativeLoop [ 50, 0.001000, 0.400000]
				{ 
					UpdateConstraint[resol_getdpThermal, Region[All], Assign ]; 
					Print ["Solving formulation [getdpThermal]"]; 
					GenerateJac[resol_getdpThermal];
					SolveJac[resol_getdpThermal];
				} 

 		   	} 
 		   	{ 
				SaveSolution[resol_getdpThermal]; 
				Print[{$Time}, File "../getdpresults/times.dat"]; 
				SaveSolution[resol_getdpEMGeneral2D]; 
				PostOperation[postOp_getdpThermal]; 
				PostOperation[postOp_getdpEMGeneral2D]; 

        	} 
		} 
	} 
}