Commit 20a7bd05 authored by Baptiste Durand's avatar Baptiste Durand

Correction of FullScaleModel

Handling of simulations without imposed force.
parent 0ff60840
......@@ -265,7 +265,8 @@ class FullScaleModel(object):
L_terms += contrib_list
L = sum(L_terms)
else:
L = 0.0
zero = fe.Constant(np.zeros(shape=self.v.ufl_shape))
L = fe.dot(zero, self.v) * self.measures[self.part.dim]
logger.info("Assembling system...")
K, res = fe.assemble_system(self.a, L, self.Dirichlet_bc)
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment