Rheolef  7.1
an efficient C++ finite element environment
vortex_position.cc
Go to the documentation of this file.
1 #include "rheolef.h"
26 using namespace rheolef;
27 int main (int argc, char** argv) {
28  environment rheolef (argc, argv);
29  check_macro (communicator().size() == 1, "please, use sequentially");
30  field psi_h;
31  din >> psi_h;
32  size_t idof_min = 0;
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);
37  idof_min = idof;
38  }
39  const disarray<point>& xdof = psi_h.get_space().get_xdofs();
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;
43 }
rheolef::space_numbering::ndof
size_type ndof(const basis_basic< T > &b, const geo_size &gs, size_type map_dim)
Definition: space_numbering.cc:28
check_macro
check_macro(expr1.have_homogeneous_space(Xh1), "dual(expr1,expr2); expr1 should have homogeneous space. HINT: use dual(interpolate(Xh, expr1),expr2)")
main
int main(int argc, char **argv)
Definition: vortex_position.cc:27
field
see the field page for the full documentation
rheolef.h
rheolef - reference manual
rheolef::din
idiststream din
see the diststream page for the full documentation
Definition: diststream.h:427
rheolef::environment
see the environment page for the full documentation
Definition: environment.h:115
rheolef
This file is part of Rheolef.
Definition: compiler_eigen.h:37
Float
see the Float page for the full documentation
point
see the point page for the full documentation
rheolef::disarray
see the disarray page for the full documentation
Definition: disarray.h:459
rheolef::dout
odiststream dout(cout)
see the diststream page for the full documentation
Definition: diststream.h:430