Mentions légales du service

Skip to content
Snippets Groups Projects
Commit 036273ac authored by Martin Genet's avatar Martin Genet
Browse files

assemble_ener now assemble energy by energy

parent 197f1141
No related branches found
No related tags found
No related merge requests found
......@@ -181,14 +181,22 @@ class ImageRegistrationProblem(Problem):
def assemble_ener(self):
ener_form = 0.
ener = 0.
for energy in self.energies:
ener_form += dolfin.Constant(energy.w) * energy.ener_form
ener = dolfin.assemble(
ener_form)
ener_ = dolfin.assemble(
energy.ener_form)
self.printer.print_var("ener_"+ener.name,ener_)
ener += energy.w * ener_
#self.printer.print_var("ener",ener)
# ener_form = 0.
# for energy in self.energies:
# ener_form += dolfin.Constant(energy.w) * energy.ener_form
#
# ener = dolfin.assemble(
# ener_form)
# #self.printer.print_var("ener",ener)
return ener
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment