def solve_all(self, problems, *, n_jobs=1, **kwargs):
"""Solve several problems.
Optional keyword arguments are passed to `BEMSolver.solve`.
Parameters
----------
problems: list of LinearPotentialFlowProblem
several problems to be solved
n_jobs: int, optional (default: 1)
the number of jobs to run in parallel using the optional dependency `joblib`
By defaults: do not use joblib and solve sequentially.
Returns
-------
list of LinearPotentialFlowResult
the solved problems
"""
if n_jobs == 1: # force sequential resolution
return [self.solve(pb, **kwargs) for pb in sorted(problems)]
else:
joblib = silently_import_optional_dependency("joblib")
if joblib is None:
raise ImportError(
f"Setting the `n_jobs` argument to {n_jobs} requires the missing optional dependency 'joblib'."
)
groups_of_problems = LinearPotentialFlowProblem._group_for_parallel_resolution(
problems)
groups_of_results = joblib.Parallel(n_jobs=n_jobs)(
joblib.delayed(self.solve_all)(grp, n_jobs=1, **kwargs)
for grp in groups_of_problems)
results = [res for grp in groups_of_results
for res in grp] # flatten the nested list
return results
def test_results():
assert isinstance(LinearPotentialFlowProblem().make_results_container(), LinearPotentialFlowResult)
pb = DiffractionProblem(g=10, rho=1023, free_surface=np.infty)
res = DiffractionResult(pb)
assert res.g == pb.g == 10
assert "DiffractionResult" in str(res)
pb = RadiationProblem(g=10, rho=1023, free_surface=np.infty)
res = RadiationResult(pb)
assert res.g == pb.g == 10
assert "RadiationResult" in str(res)
def test_LinearPotentialFlowProblem():
# Without a body
pb = LinearPotentialFlowProblem(omega=1.0)
assert pb.omega == 1.0
assert pb.period == 2*np.pi
assert pb.wavenumber == 1.0/9.81
assert pb.wavelength == 9.81*2*np.pi
assert LinearPotentialFlowProblem(free_surface=np.infty, sea_bottom=-np.infty).depth == np.infty
assert LinearPotentialFlowProblem(free_surface=0.0, sea_bottom=-np.infty).depth == np.infty
pb = LinearPotentialFlowProblem(free_surface=0.0, sea_bottom=-1.0, omega=1.0)
assert pb.depth == 1.0
assert np.isclose(pb.omega**2, pb.g*pb.wavenumber*np.tanh(pb.wavenumber*pb.depth))
assert pb.dimensionless_wavenumber == pb.wavenumber*1.0
with pytest.raises(NotImplementedError):
LinearPotentialFlowProblem(free_surface=2.0)
with pytest.raises(NotImplementedError):
LinearPotentialFlowProblem(free_surface=np.infty, sea_bottom=0.0)
with pytest.raises(ValueError):
LinearPotentialFlowProblem(free_surface=0.0, sea_bottom=1.0)
with pytest.raises(TypeError):
LinearPotentialFlowProblem(wave_direction=1.0)
with pytest.raises(TypeError):
LinearPotentialFlowProblem(radiating_dof="Heave")
# With a body
sphere = Sphere(center=(0, 0, -2.0))
sphere.add_translation_dof(direction=(0, 0, 1), name="Heave")
pb = LinearPotentialFlowProblem(body=sphere, omega=1.0)
pb.boundary_condition = sphere.mesh.faces_normals @ (1, 1, 1)
assert list(pb.influenced_dofs.keys()) == ['Heave']
pb2 = LinearPotentialFlowProblem(body=sphere, omega=2.0)
pb2.boundary_condition = sphere.mesh.faces_normals @ (1, 1, 1)
assert pb < pb2
# Test transformation to result class
res = pb.make_results_container()
assert isinstance(res, LinearPotentialFlowResult)
assert res.problem is pb
assert res.omega == pb.omega
assert res.period == pb.period
assert res.body is pb.body