Zero out backsubstitution (#155)

jorgensd · jorgensd · commit f799a750f0f2 · 2025-03-03T21:49:07.000+01:00
diff --git a/cpp/MultiPointConstraint.h b/cpp/MultiPointConstraint.h
@@ -130,6 +130,9 @@ class MultiPointConstraint
   {
     for (auto slave : _slaves)
     {
+      // Zero out intial data in slave dof
+      vector[slave] = 0.0;
+      // Accumulate master contributions
       auto masters = _master_map->links(slave);
       auto coeffs = _coeff_map->links(slave);
       assert(masters.size() == coeffs.size());

Original file line number	Diff line number	Diff line change
`@@ -130,6 +130,9 @@ class MultiPointConstraint`
`130`	`130`	`{`
`131`	`131`	`for (auto slave : _slaves)`
`132`	`132`	`{`
	`133`	`+ // Zero out intial data in slave dof`
	`134`	`+ vector[slave] = 0.0;`
	`135`	`+ // Accumulate master contributions`
`133`	`136`	`auto masters = _master_map->links(slave);`
`134`	`137`	`auto coeffs = _coeff_map->links(slave);`
`135`	`138`	`assert(masters.size() == coeffs.size());`