27int main (
int argc,
char** argv) {
33 Float psi_min = std::numeric_limits<Float>::max();
34 for (
size_t idof = 0, ndof = psi_h.ndof(); idof < ndof; idof++) {
35 if (psi_h.dof(idof) >= psi_min)
continue;
36 psi_min = psi_h.dof(idof);
40 point xmin = xdof [idof_min];
41 dout <<
"xc\t\tyc\t\tpsi" << std::endl
42 << xmin[0] <<
"\t" << xmin[1] <<
"\t" << psi_min << std::endl;
see the Float page for the full documentation
see the communicator page for the full documentation
see the field page for the full documentation
see the point page for the full documentation
see the disarray page for the full documentation
see the environment page for the full documentation
check_macro(expr1.have_homogeneous_space(Xh1), "dual(expr1,expr2); expr1 should have homogeneous space. HINT: use dual(interpolate(Xh, expr1),expr2)")
This file is part of Rheolef.
rheolef - reference manual