MOAB: Mesh Oriented datABase  (version 5.4.0)
ElementPMeanP.cpp
Go to the documentation of this file.
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,
00069                                             size_t handle,
00070                                             double& value,
00071                                             std::vector< size_t >& indices,
00072                                             std::vector< Vector3D >& gradient,
00073                                             MsqError& err )
00074 {
00075     ElemSampleQM* qm = get_quality_metric();
00076     mHandles.clear();
00077     qm->get_element_evaluations( pd, handle, mHandles, err );
00078     MSQ_ERRFALSE( err );
00079     bool result = average_with_gradient( pd, qm, mHandles, value, indices, gradient, err );
00080     return !MSQ_CHKERR( err ) && result;
00081 }
00082 
00083 bool ElementPMeanP::evaluate_with_Hessian( PatchData& pd,
00084                                            size_t handle,
00085                                            double& value,
00086                                            std::vector< size_t >& indices,
00087                                            std::vector< Vector3D >& gradient,
00088                                            std::vector< Matrix3D >& Hessian,
00089                                            MsqError& err )
00090 {
00091     ElemSampleQM* qm = get_quality_metric();
00092     mHandles.clear();
00093     qm->get_element_evaluations( pd, handle, mHandles, err );
00094     MSQ_ERRFALSE( err );
00095     bool result = average_with_Hessian( pd, qm, mHandles, value, indices, gradient, Hessian, err );
00096     return !MSQ_CHKERR( err ) && result;
00097 }
00098 
00099 bool ElementPMeanP::evaluate_with_Hessian_diagonal( PatchData& pd,
00100                                                     size_t handle,
00101                                                     double& value,
00102                                                     std::vector< size_t >& indices,
00103                                                     std::vector< Vector3D >& gradient,
00104                                                     std::vector< SymMatrix3D >& diagonal,
00105                                                     MsqError& err )
00106 {
00107     ElemSampleQM* qm = get_quality_metric();
00108     mHandles.clear();
00109     qm->get_element_evaluations( pd, handle, mHandles, err );
00110     MSQ_ERRFALSE( err );
00111     bool result = average_with_Hessian_diagonal( pd, qm, mHandles, value, indices, gradient, diagonal, err );
00112     return !MSQ_CHKERR( err ) && result;
00113 }
00114 
00115 }  // namespace MBMesquite
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Defines