v0.9.1 |

gauss_points_on_quad.cpp

Testing gauss points coordinates and weights for a quad face

/** \file gauss_points_on_quad.cpp

\example gauss_points_on_quad.cpp

\brief Testing gauss points coordinates and weights for a quad face

*/

/* This file is part of MoFEM.

* MoFEM is free software: you can redistribute it and/or modify it under

* the terms of the GNU Lesser General Public License as published by the

* Free Software Foundation, either version 3 of the License, or (at your

* option) any later version.

*

* MoFEM is distributed in the hope that it will be useful, but WITHOUT

* ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or

* FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public

* License for more details.

*

* You should have received a copy of the GNU Lesser General Public

* License along with MoFEM. If not, see <http://www.gnu.org/licenses/>. */

#include <MoFEM.hpp>

using namespace MoFEM;

try {

moab::Core mb_instance;

moab::Interface &moab = mb_instance;

ParallelComm *pcomm = ParallelComm::get_pcomm(&moab, MYPCOMM_INDEX);

if (pcomm == NULL)

pcomm = new ParallelComm(&moab, PETSC_COMM_WORLD);

MoFEM::Core core(moab);

MatrixDouble pts_quad;

int rule_ksi = 6;

int rule_eta = 8;

rule_eta);

int nb_gauss_pts = pts_quad.size2();

double sum_coords = 0, sum_gauss_pts = 0;

sum_gauss_pts += pts_quad(2, i);

}

if (fabs(20.0 - sum_coords) > eps) {

}

if (fabs(1.0 - sum_gauss_pts) > eps) {

SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY, "wrong result");

}

}

}

Generated on Mon Feb 3 2020 15:26:17 for MoFEM by Doxygen 1.8.15 and hosted at