Allow dirichlet functions to return NaN

Functionality like this was requested long ago; if I'd known it would be a one-line change I would have implemented it then.
colibius · Apr 22, 2014 · 313be2f · 313be2f
1 parent 1ea7b56
commit 313be2f
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/base/dof_map_constraints.C b/src/base/dof_map_constraints.C
@@ -831,7 +831,7 @@ private:
           for (unsigned int i = 0; i < n_dofs; i++)
             {
               DofConstraintRow empty_row;
-              if (dof_is_fixed[i])
+              if (dof_is_fixed[i] && !std::isnan(Ue(i)))
                 add_fn (dof_indices[i], empty_row, Ue(i));
             }
         }