MOAB: Mesh Oriented datABase
(version 5.2.1)
|
00001 /* ***************************************************************** 00002 MESQUITE -- The Mesh Quality Improvement Toolkit 00003 00004 Copyright 2006 Sandia National Laboratories. Developed at the 00005 University of Wisconsin--Madison under SNL contract number 00006 624796. The U.S. Government and the University of Wisconsin 00007 retain certain rights to this software. 00008 00009 This library is free software; you can redistribute it and/or 00010 modify it under the terms of the GNU Lesser General Public 00011 License as published by the Free Software Foundation; either 00012 version 2.1 of the License, or (at your option) any later version. 00013 00014 This library is distributed in the hope that it will be useful, 00015 but WITHOUT ANY WARRANTY; without even the implied warranty of 00016 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU 00017 Lesser General Public License for more details. 00018 00019 You should have received a copy of the GNU Lesser General Public License 00020 (lgpl.txt) along with this library; if not, write to the Free Software 00021 Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA 00022 00023 (2006) kraftche@cae.wisc.edu 00024 00025 ***************************************************************** */ 00026 00027 /** \file ElementPMeanP.cpp 00028 * \brief 00029 * \author Jason Kraftcheck 00030 */ 00031 00032 #include "Mesquite.hpp" 00033 #include "ElementPMeanP.hpp" 00034 #include "ElemSampleQM.hpp" 00035 #include "MsqError.hpp" 00036 #include "PatchData.hpp" 00037 00038 namespace MBMesquite 00039 { 00040 00041 ElementPMeanP::ElementPMeanP( double p, ElemSampleQM* metric ) : PMeanPMetric( p ), mMetric( metric ) {} 00042 00043 ElementPMeanP::~ElementPMeanP() {} 00044 00045 std::string ElementPMeanP::get_name() const 00046 { 00047 std::string result( "ElementPMeanP(" ); 00048 result += mMetric->get_name(); 00049 result += ")"; 00050 return result; 00051 } 00052 00053 int ElementPMeanP::get_negate_flag() const 00054 { 00055 return get_quality_metric()->get_negate_flag(); 00056 } 00057 00058 bool ElementPMeanP::evaluate( PatchData& pd, size_t handle, double& value, MsqError& err ) 00059 { 00060 ElemSampleQM* qm = get_quality_metric(); 00061 mHandles.clear(); 00062 qm->get_element_evaluations( pd, handle, mHandles, err ); 00063 MSQ_ERRFALSE( err ); 00064 bool result = average( pd, qm, mHandles, value, err ); 00065 return !MSQ_CHKERR( err ) && result; 00066 } 00067 00068 bool ElementPMeanP::evaluate_with_gradient( PatchData& pd, size_t handle, double& value, std::vector< size_t >& indices, 00069 std::vector< Vector3D >& gradient, MsqError& err ) 00070 { 00071 ElemSampleQM* qm = get_quality_metric(); 00072 mHandles.clear(); 00073 qm->get_element_evaluations( pd, handle, mHandles, err ); 00074 MSQ_ERRFALSE( err ); 00075 bool result = average_with_gradient( pd, qm, mHandles, value, indices, gradient, err ); 00076 return !MSQ_CHKERR( err ) && result; 00077 } 00078 00079 bool ElementPMeanP::evaluate_with_Hessian( PatchData& pd, size_t handle, double& value, std::vector< size_t >& indices, 00080 std::vector< Vector3D >& gradient, std::vector< Matrix3D >& Hessian, 00081 MsqError& err ) 00082 { 00083 ElemSampleQM* qm = get_quality_metric(); 00084 mHandles.clear(); 00085 qm->get_element_evaluations( pd, handle, mHandles, err ); 00086 MSQ_ERRFALSE( err ); 00087 bool result = average_with_Hessian( pd, qm, mHandles, value, indices, gradient, Hessian, err ); 00088 return !MSQ_CHKERR( err ) && result; 00089 } 00090 00091 bool ElementPMeanP::evaluate_with_Hessian_diagonal( PatchData& pd, size_t handle, double& value, 00092 std::vector< size_t >& indices, std::vector< Vector3D >& gradient, 00093 std::vector< SymMatrix3D >& diagonal, MsqError& err ) 00094 { 00095 ElemSampleQM* qm = get_quality_metric(); 00096 mHandles.clear(); 00097 qm->get_element_evaluations( pd, handle, mHandles, err ); 00098 MSQ_ERRFALSE( err ); 00099 bool result = average_with_Hessian_diagonal( pd, qm, mHandles, value, indices, gradient, diagonal, err ); 00100 return !MSQ_CHKERR( err ) && result; 00101 } 00102 00103 } // namespace MBMesquite