Alex, > reinitialization_constraints.add_entry(i,0, 0.0 ); > > reinitialization_constraints.set_inhomogeneity(i,old_reinitialization_so > lution(i));
you can remove the first line -- it leads to the constraint x_i = 0 * x_0 + old_reinitialization_solution(i) The first part is of course entirely irrelevant. > For the moment this approach does not seem to be doing what I meant to > do, although I need to look into that more thoroughly. The usual question is of course what exactly is happening. There is little we can guess that may go wrong unless we know what exactly is happening. > Do I have to use the relevant VectorTools::interpolate_boundary_values > to enforce the inhomogeneity in the constraint? No. W. ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
