def __init__(self, problem, solver_params_fluid, solver_params_flow0d):
self.pb = problem
self.solver_params_fluid = solver_params_fluid
self.solver_params_flow0d = solver_params_flow0d
# initialize nonlinear solver class
self.solnln = solver_nonlin.solver_nonlinear_constraint_monolithic(
self.pb, self.pb.pbs.V_v, self.pb.pbs.V_p,
self.solver_params_fluid, self.solver_params_flow0d)
def __init__(self, problem, solver_params_solid, solver_params_flow0d):
self.pb = problem
self.solver_params_solid = solver_params_solid
self.solver_params_flow0d = solver_params_flow0d
# initialize nonlinear solver class
self.solnln = solver_nonlin.solver_nonlinear_constraint_monolithic(
self.pb, self.pb.pbs.V_u, self.pb.pbs.V_p,
self.solver_params_solid, self.solver_params_flow0d)
if self.pb.pbs.prestress_initial:
# add coupling work to prestress weak form
self.pb.pbs.weakform_prestress_u -= self.pb.work_coupling_prestr
# initialize solid mechanics solver
self.solverprestr = SolidmechanicsSolver(self.pb.pbs,
self.solver_params_solid)