26 #ifndef OOMPH_ANNULAR_MESH_TEMPLATE_CC
27 #define OOMPH_ANNULAR_MESH_TEMPLATE_CC
36 template<
class ELEMENT>
40 const double& azimuthal_fraction,
50 const unsigned n_node = this->nnode();
51 for (
unsigned n = 0; n < n_node; n++)
54 Node* nod_pt = this->node_pt(n);
59 xi[0] = (1.0 - 1.0e-10) * (-azimuthal_fraction * nod_pt->
x(0)) * 2.0 *
68 double w = nod_pt->
x(1);
74 double norm = sqrt(base[0] * base[0] + base[1] * base[1]);
75 N[0] = base[0] / norm;
76 N[1] = base[1] / norm;
79 nod_pt->
x(0) = base[0] + w * (h + a - norm) *
N[0];
80 nod_pt->
x(1) = base[1] + w * (h + a - norm) *
N[1];
88 xi_bound[0] = atan2(nod_pt->
x(1), nod_pt->
x(0));
95 xi_bound[0] = sqrt(pow(nod_pt->
x(0), 2) + pow(nod_pt->
x(1), 2));
102 xi_bound[0] = atan2(nod_pt->
x(1), nod_pt->
x(0));
109 xi_bound[0] = sqrt(pow(nod_pt->
x(0), 2) + pow(nod_pt->
x(1), 2));
114 this->Boundary_coordinate_exists[0] =
true;
115 this->Boundary_coordinate_exists[1] =
true;
116 this->Boundary_coordinate_exists[2] =
true;
117 this->Boundary_coordinate_exists[3] =
true;
////////////////////////////////////////////////////////////////////
void position(const Vector< double > &zeta, Vector< double > &r) const
Position Vector at Lagrangian coordinate zeta.
Nodes are derived from Data, but, in addition, have a definite (Eulerian) position in a space of a gi...
double & x(const unsigned &i)
Return the i-th nodal coordinate.
virtual bool is_on_boundary() const
Test whether the Node lies on a boundary. The "bulk" Node cannot lie on a boundary,...
virtual void set_coordinates_on_boundary(const unsigned &b, const unsigned &k, const Vector< double > &boundary_zeta)
Set the vector of the k-th generalised boundary coordinates on mesh boundary b. Broken virtual interf...
void wrap_into_annular_shape(const double &a, const double &h, const double &azimuthal_fraction, const double &phi)
Wrap mesh into annular shape.
const double Pi
50 digits from maple
//////////////////////////////////////////////////////////////////// ////////////////////////////////...