46 MeshChecker::assert_geometric_element<QElementGeometricBase, ELEMENT>(2);
60 Node_pt.resize((1 + (Np - 1) * Nx) * (1 + (Np - 1) * Ny));
90 for (
unsigned l2 = 1;
l2 < Np;
l2++)
104 if ((Nx == 1) && (
l2 == (Np - 1)))
114 for (
unsigned l1 = 1;
l1 < Np;
l1++)
128 if ((Ny == 1) && (
l1 == (Np - 1)))
137 for (
unsigned l2 = 1;
l2 < Np;
l2++)
141 if (((Nx == 1) && (
l2 == (Np - 1))) || ((Ny == 1) && (
l1 == (Np - 1))))
158 if ((Nx == 1) &&
l2 == (Np - 1))
163 if ((Ny == 1) && (
l1 == (Np - 1)))
176 for (
unsigned j = 1;
j < (Nx - 1);
j++)
185 for (
unsigned l2 = 1;
l2 < Np;
l2++)
203 for (
unsigned l1 = 1;
l1 < Np;
l1++)
209 for (
unsigned l2 = 1;
l2 < Np;
l2++)
213 if ((Ny == 1) && (
l1 == (Np - 1)))
230 if ((Ny == 1) && (
l1 == (Np - 1)))
254 for (
unsigned l2 = 1;
l2 < (Np - 1);
l2++)
277 if (Xperiodic ==
true)
294 for (
unsigned l1 = 1;
l1 < Np;
l1++)
301 for (
unsigned l2 = 1;
l2 < (Np - 1);
l2++)
305 if ((Ny == 1) && (
l1 == (Np - 1)))
323 if ((Ny == 1) && (
l1 == (Np - 1)))
339 if (Xperiodic ==
true)
353 if ((Ny == 1) && (
l1 == (Np - 1)))
366 for (
unsigned i = 1;
i < (Ny - 1);
i++)
373 for (
unsigned l2 = 0;
l2 < Np;
l2++)
380 for (
unsigned l1 = 1;
l1 < Np;
l1++)
399 for (
unsigned l2 = 1;
l2 < Np;
l2++)
403 if ((Nx == 1) && (
l2 == (Np - 1)))
421 if ((Nx == 1) && (
l2 == (Np - 1)))
432 for (
unsigned j = 1;
j < (Nx - 1);
j++)
437 for (
unsigned l2 = 0;
l2 < Np;
l2++)
443 for (
unsigned l1 = 1;
l1 < Np;
l1++)
449 for (
unsigned l2 = 1;
l2 < Np;
l2++)
473 for (
unsigned l2 = 0;
l2 < Np;
l2++)
477 ->node_pt((Np - 1) * Np +
l2);
480 for (
unsigned l1 = 1;
l1 < Np;
l1++)
487 for (
unsigned l2 = 1;
l2 < (Np - 1);
l2++)
509 if (Xperiodic ==
true)
538 for (
unsigned l2 = 0;
l2 < Np;
l2++)
546 for (
unsigned l1 = 1;
l1 < (Np - 1);
l1++)
564 for (
unsigned l2 = 1;
l2 < Np;
l2++)
567 if ((Nx == 1) && (
l2 == Np - 1))
585 if ((Nx == 1) && (
l2 == Np - 1))
615 for (
unsigned l2 = 1;
l2 < Np;
l2++)
630 if ((Nx == 1) && (
l2 == Np - 1))
640 for (
unsigned j = 1;
j < (Nx - 1);
j++)
645 for (
unsigned l2 = 0;
l2 < Np;
l2++)
652 for (
unsigned l1 = 1;
l1 < (Np - 1);
l1++)
657 ->node_pt(Np *
l1 + (Np - 1));
660 for (
unsigned l2 = 1;
l2 < Np;
l2++)
680 ->node_pt(Np * (Np - 1) + (Np - 1));
682 for (
unsigned l2 = 1;
l2 < Np;
l2++)
708 for (
unsigned l2 = 0;
l2 < Np;
l2++)
712 ->node_pt((Np - 1) * Np +
l2);
716 for (
unsigned l1 = 1;
l1 < (Np - 1);
l1++)
721 ->node_pt(Np *
l1 + (Np - 1));
724 for (
unsigned l2 = 1;
l2 < (Np - 1);
l2++)
733 x_spacing_function(Nx - 1,
l2, Ny - 1,
l1);
735 y_spacing_function(Nx - 1,
l2, Ny - 1,
l1);
748 if (Xperiodic ==
true)
756 x_spacing_function(Nx - 1, Np - 1, Ny - 1,
l1);
758 y_spacing_function(Nx - 1, Np - 1, Ny - 1,
l1);
772 ->node_pt(Np * (Np - 1) + (Np - 1));
775 for (
unsigned l2 = 1;
l2 < (Np - 1);
l2++)
784 x_spacing_function(Nx - 1,
l2, Ny - 1, Np - 1);
786 y_spacing_function(Nx - 1,
l2, Ny - 1, Np - 1);
798 ->construct_boundary_node(
802 if (Xperiodic ==
true)
810 x_spacing_function(Nx - 1, Np - 1, Ny - 1, Np - 1);
812 y_spacing_function(Nx - 1, Np - 1, Ny - 1, Np - 1);
825 setup_boundary_element_info();