OOF.File.Load.Data(filename='../33676.skel')
OOF.Mesh.New(name='mesh', skeleton='MeshSubject_02_2dec.pgm:skeleton', element_types=['T3_3', 'Q4_4'])
OOF.Subproblem.Field.Define(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', field=Displacement)
OOF.Subproblem.Field.Activate(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', field=Displacement)
OOF.Subproblem.Equation.Activate(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', equation=Force_Balance)
OOF.Subproblem.Equation.Activate(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', equation=Plane_Stress)
OOF.Mesh.Boundary_Conditions.New(name='bc', mesh='MeshSubject_02_2dec.pgm:skeleton:mesh', condition=DirichletBC(field=Displacement,field_component='y',equation=Force_Balance,eqn_component='y',profile=ConstantProfile(value=3),boundary='top'))
OOF.Mesh.Boundary_Conditions.New(name='bc<2>', mesh='MeshSubject_02_2dec.pgm:skeleton:mesh', condition=DirichletBC(field=Displacement,field_component='y',equation=Force_Balance,eqn_component='y',profile=ConstantProfile(value=0),boundary='bottom'))
OOF.Solver.Solve(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', solver=LinearDriver(method=CGSolver(max_iterations=1000,tolerance=1e-13,preconditioner=ILUPreconditioner())))
OOF.Solver.Solve(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', solver=LinearDriver(method=CGSolver(max_iterations=1000,tolerance=1e-13,preconditioner=ILUPreconditioner())))
OOF.Solver.Solve(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', solver=LinearDriver(method=CGSolver(max_iterations=3000,tolerance=1e-13,preconditioner=ILUPreconditioner())))
OOF.Solver.Solve(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', solver=LinearDriver(method=GMRESSolver(max_iterations=1000,krylov_dimension=200,tolerance=1e-13,preconditioner=ILUPreconditioner())))
OOF.Subproblem.Equation.Activate(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', equation=Heat_Eqn)
OOF.Subproblem.Equation.Activate(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', equation=Plane_Heat_Flux)
OOF.Subproblem.Equation.Deactivate(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', equation=Force_Balance)
OOF.Subproblem.Equation.Deactivate(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', equation=Plane_Stress)
OOF.Subproblem.Field.Deactivate(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', field=Displacement)
OOF.Subproblem.Field.Define(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', field=Temperature)
OOF.Subproblem.Field.Activate(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', field=Temperature)
OOF.Mesh.Boundary_Conditions.New(name='bc<3>', mesh='MeshSubject_02_2dec.pgm:skeleton:mesh', condition=DirichletBC(field=Temperature,field_component='',equation=Heat_Eqn,eqn_component='',profile=ConstantProfile(value=0.0),boundary='bottom'))
OOF.Mesh.Boundary_Conditions.New(name='bc<4>', mesh='MeshSubject_02_2dec.pgm:skeleton:mesh', condition=DirichletBC(field=Temperature,field_component='',equation=Heat_Eqn,eqn_component='',profile=ConstantProfile(value=1),boundary='top'))
OOF.Mesh.Boundary_Conditions.Disable(mesh='MeshSubject_02_2dec.pgm:skeleton:mesh', name='bc')
OOF.Mesh.Boundary_Conditions.Disable(mesh='MeshSubject_02_2dec.pgm:skeleton:mesh', name='bc<2>')
OOF.Solver.Solve(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', solver=LinearDriver(method=CGSolver(max_iterations=1000,tolerance=1e-13,preconditioner=ILUPreconditioner())))
OOF.Subproblem.Equation.Deactivate(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', equation=Heat_Eqn)
OOF.Subproblem.Equation.Deactivate(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', equation=Plane_Heat_Flux)
OOF.Subproblem.Equation.Activate(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', equation=Force_Balance)
OOF.Subproblem.Equation.Activate(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', equation=Plane_Stress)
OOF.Subproblem.Field.Deactivate(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', field=Temperature)
OOF.Subproblem.Field.Activate(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', field=Displacement)
OOF.Mesh.Boundary_Conditions.Enable(mesh='MeshSubject_02_2dec.pgm:skeleton:mesh', name='bc<2>')
OOF.Mesh.Boundary_Conditions.Enable(mesh='MeshSubject_02_2dec.pgm:skeleton:mesh', name='bc')
OOF.Mesh.Boundary_Conditions.Disable(mesh='MeshSubject_02_2dec.pgm:skeleton:mesh', name='bc<3>')
OOF.Mesh.Boundary_Conditions.Disable(mesh='MeshSubject_02_2dec.pgm:skeleton:mesh', name='bc<4>')
OOF.Solver.Solve(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', solver=LinearDriver(method=CGSolver(max_iterations=1000,tolerance=1e-13,preconditioner=ILUPreconditioner())))
OOF.File.Save.Mesh(filename='/Users/ccarter/People/Neri/MeshSubject_02/vertical_loads/33676_not_converged.msh', mode='w', format='ascii', mesh='MeshSubject_02_2dec.pgm:skeleton:mesh')
OOF.Subproblem.Field.Deactivate(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', field=Displacement)
OOF.Subproblem.Field.Activate(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', field=Temperature)
OOF.Subproblem.Equation.Activate(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', equation=Heat_Eqn)
OOF.Subproblem.Equation.Activate(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', equation=Plane_Heat_Flux)
OOF.Subproblem.Equation.Deactivate(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', equation=Force_Balance)
OOF.Subproblem.Equation.Deactivate(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', equation=Plane_Stress)
OOF.Mesh.Boundary_Conditions.Enable(mesh='MeshSubject_02_2dec.pgm:skeleton:mesh', name='bc<4>')
OOF.Mesh.Boundary_Conditions.Enable(mesh='MeshSubject_02_2dec.pgm:skeleton:mesh', name='bc<3>')
OOF.Mesh.Boundary_Conditions.Disable(mesh='MeshSubject_02_2dec.pgm:skeleton:mesh', name='bc')
OOF.Mesh.Boundary_Conditions.Disable(mesh='MeshSubject_02_2dec.pgm:skeleton:mesh', name='bc<2>')
OOF.Solver.Solve(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', solver=LinearDriver(method=CGSolver(max_iterations=1000,tolerance=1e-13,preconditioner=ILUPreconditioner())))
OOF.Subproblem.Equation.Activate(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', equation=Force_Balance)
OOF.Subproblem.Equation.Activate(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', equation=Plane_Stress)
OOF.Subproblem.Field.Activate(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', field=Displacement)
OOF.Mesh.Boundary_Conditions.Enable(mesh='MeshSubject_02_2dec.pgm:skeleton:mesh', name='bc<2>')
OOF.Mesh.Boundary_Conditions.Enable(mesh='MeshSubject_02_2dec.pgm:skeleton:mesh', name='bc')
OOF.Solver.Solve(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', solver=LinearDriver(method=GMRESSolver(max_iterations=1000,krylov_dimension=200,tolerance=1e-13,preconditioner=ILUPreconditioner())))
OOF.File.Save.Mesh(filename='/Users/ccarter/People/Neri/MeshSubject_02/vertical_loads/33676_not_converged.msh', mode='w', format='ascii', mesh='MeshSubject_02_2dec.pgm:skeleton:mesh')
OOF.Solver.Solve(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', solver=LinearDriver(method=CGSolver(max_iterations=1000,tolerance=1e-13,preconditioner=ILUPreconditioner())))
OOF.Solver.Solve(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', solver=LinearDriver(method=GMRESSolver(max_iterations=1000,krylov_dimension=200,tolerance=1e-13,preconditioner=ILUPreconditioner())))
OOF.Solver.Solve(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', solver=LinearDriver(method=GMRESSolver(max_iterations=3000,krylov_dimension=200,tolerance=1e-13,preconditioner=ILUPreconditioner())))
OOF.Solver.Solve(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', solver=LinearDriver(method=CGSolver(max_iterations=1000,tolerance=1e-13,preconditioner=ILUPreconditioner())))
OOF.Solver.Solve(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', solver=LinearDriver(method=GMRESSolver(max_iterations=3000,krylov_dimension=200,tolerance=1e-13,preconditioner=ILUPreconditioner())))
OOF.Solver.Solve(subproblem='MeshSubject_02_2dec.pgm:skeleton:mesh:default', solver=LinearDriver(method=GMRESSolver(max_iterations=3000,krylov_dimension=200,tolerance=1e-13,preconditioner=ILUPreconditioner())))
OOF.File.Save.Mesh(filename='/Users/ccarter/People/Neri/MeshSubject_02/vertical_loads/33676_converged.msh', mode='w', format='ascii', mesh='MeshSubject_02_2dec.pgm:skeleton:mesh')
OOF.Mesh.Analyze.Direct_Output(mesh='MeshSubject_02_2dec.pgm:skeleton:mesh', data=getOutput('Flux:Invariant',invariant=Magnitude(),flux=Heat_Flux), domain=EntireMesh(), sampling=GridSampleSet(x_points=10,y_points=10,show_x=True,show_y=True), destination='/Users/ccarter/People/Neri/MeshSubject_02/vertical_loads/jmag.cvs', comment_character='#', separator='comma')
OOF.Mesh.Analyze.Direct_Output(mesh='MeshSubject_02_2dec.pgm:skeleton:mesh', data=getOutput('Flux:Invariant',invariant=Magnitude(),flux=Heat_Flux), domain=EntireMesh(), sampling=GridSampleSet(x_points=200,y_points=200,show_x=True,show_y=True), destination='/Users/ccarter/People/Neri/MeshSubject_02/vertical_results/jmag.cvs', comment_character='#', separator='comma')
OOF.Mesh.Analyze.Direct_Output(mesh='MeshSubject_02_2dec.pgm:skeleton:mesh', data=getOutput('Flux:Invariant',invariant=MatrixTrace(),flux=Stress), domain=EntireMesh(), sampling=GridSampleSet(x_points=200,y_points=200,show_x=True,show_y=True), destination='/Users/ccarter/People/Neri/MeshSubject_02/vertical_results/stress.cvs', comment_character='#', separator='comma')
OOF.Mesh.Analyze.Direct_Output(mesh='MeshSubject_02_2dec.pgm:skeleton:mesh', data=getOutput('Energy',etype='Total'), domain=EntireMesh(), sampling=GridSampleSet(x_points=200,y_points=200,show_x=True,show_y=True), destination='/Users/ccarter/People/Neri/MeshSubject_02/vertical_results/energy.cvs', comment_character='#', separator='comma')
OOF.Mesh.Analyze.Direct_Output(mesh='MeshSubject_02_2dec.pgm:skeleton:mesh', data=getOutput('Field:Component',component='z',field=Displacement_z), domain=EntireMesh(), sampling=GridSampleSet(x_points=200,y_points=200,show_x=True,show_y=True), destination='/Users/ccarter/People/Neri/MeshSubject_02/vertical_results/z.cvs', comment_character='#', separator='comma')
OOF.Mesh.Analyze.Direct_Output(mesh='MeshSubject_02_2dec.pgm:skeleton:mesh', data=getOutput('Strain:Component',component='yy',type=GeometricStrain()), domain=EntireMesh(), sampling=GridSampleSet(x_points=200,y_points=200,show_x=True,show_y=True), destination='/Users/ccarter/People/Neri/MeshSubject_02/vertical_results/z.cvs', comment_character='#', separator='comma')
OOF.Mesh.Analyze.Direct_Output(mesh='MeshSubject_02_2dec.pgm:skeleton:mesh', data=getOutput('Field:Component',component='z',field=Displacement_z), domain=EntireMesh(), sampling=GridSampleSet(x_points=200,y_points=200,show_x=True,show_y=True), destination='/Users/ccarter/People/Neri/MeshSubject_02/vertical_results/z.cvs', comment_character='#', separator='comma')
OOF.Mesh.Analyze.Direct_Output(mesh='MeshSubject_02_2dec.pgm:skeleton:mesh', data=getOutput('Strain:Component',component='yy',type=GeometricStrain()), domain=EntireMesh(), sampling=GridSampleSet(x_points=200,y_points=200,show_x=True,show_y=True), destination='/Users/ccarter/People/Neri/MeshSubject_02/vertical_results/strain.cvs', comment_character='#', separator='comma')