from scipy import *
from core import Solver
import hdf5
test = zeros((5,5))
test[2,2] = 1
sol = Solver(test, (1,0), sol_method="trilinos", printing=3)
sol.converge()
sol.sync("Extern")
hdf5.write_solver_to_h5("mpitest-results.h5", sol)
sol.dbprint("Completed without error!")
# Check to see if a simulation has been done
if hdf5.has_dP_sim(save_path, dP) and options.force:
print "\tSimulation Detected! Results will be overwritten!"
elif hdf5.has_dP_sim(save_path, dP) and not options.force:
print "\tSimulation Detected! Skipping!"
continue
# If bigmode, just pass the filename instead of the solid
if options.bigmode:
S = save_path
else:
S = hdf5.get_S(save_path)
# Setup solver, printing full debug info, using the solver specified
sol = Solver(S, dP, printing=int(options.verb), sol_method=options.solver)
# setup is now implicit
# sol.setup()
#if we want to solve monolithically . . . or use the iterative solver
if options.mono:
sol.monolithic_solve()
else:
sol.converge(options.converge)
# Save dis shit.
print "Saving . . ."
hdf5.write_solver_to_h5(save_path, sol)
del sol
