rheolef  6.6
reference_element2_tst.cc
Go to the documentation of this file.
1 #include "rheolef/reference_element.h"
2 //#include "rheolef/geo_element.h"
3 using namespace rheolef;
4 using namespace std;
5 
6 void
7 put_prism (ostream& out, size_t my_order)
8 {
10  typedef point_basic<size_type> ilat;
11  for (size_type k = 0; k <= my_order; k++) {
12  for (size_type j = 0; j <= my_order; j++) {
13  for (size_type i = 0; i+j <= my_order; i++) {
14  size_type loc_inod = reference_element_P::ilat2loc_inod (my_order, ilat(i, j, k));
15  out << " loc_inod("<<i<<","<<j<<","<<k<<") = " << loc_inod << endl;
16  }
17  }
18  }
19 }
20 int main() {
22  cout << "P1(P):"<<endl; put_prism (cout, 1);
23  cout << "P2(P):"<<endl; put_prism (cout, 2);
24  cout << "P3(P):"<<endl; put_prism (cout, 3);
25  return 0;
26 }
27 
static void put_prism(ostream &vtk, const geo_element &K, const numbering< T, sequential > &my_numb, const geo_basic< T, sequential > &omega, const array< point_basic< Float >, sequential > &my_node)
static size_type ilat2loc_inod(size_type order, const point_basic< size_type > &ilat)
point - vertex of a mesh
Definition: point.h:22
STL namespace.
std::vector< int >::size_type size_type
irheostream, orheostream - large data streams
Definition: compiler.h:7
int main()
reference_element_H::size_type size_type