> Can anyone help me with this error, my approach or a better way to impose > neumann bcs.
The VectorTools::compute_no_normal_flux function can't deal with the Raviart-Thomas element you are using. That said, since the only degrees of freedom on the faces of cells are the normal fluxes, you could use VectorTools::interpolate_boundary_values with a ZeroFunction object, or alternatively DoFTools::make_zero_boundary_constraints This will have the same effect. Best W. -- ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
