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 PMeanPMetric.hpp 00028 * \brief 00029 * \author Jason Kraftcheck 00030 */ 00031 00032 #ifndef MSQ_PMEAN_PMETRIC_HPP 00033 #define MSQ_PMEAN_PMETRIC_HPP 00034 00035 #include "Mesquite.hpp" 00036 #include "Exponent.hpp" 00037 00038 #include "Vector3D.hpp" 00039 #include "Matrix3D.hpp" 00040 00041 #include <vector> 00042 00043 namespace MBMesquite 00044 { 00045 00046 class PatchData; 00047 class MsqError; 00048 class QualityMetric; 00049 00050 class PMeanPMetric 00051 { 00052 public: 00053 PMeanPMetric( double p ) : P( p ), P1( p - 1.0 ), P2( p - 2.0 ) {} 00054 00055 const Exponent& get_power() const 00056 { 00057 return P; 00058 } 00059 00060 protected: 00061 bool average( PatchData& pd, QualityMetric* qm, const std::vector< size_t >& qm_handles, double& value, MsqError& ); 00062 00063 bool average_with_indices( PatchData& pd, QualityMetric* qm, const std::vector< size_t >& qm_handles, double& value, 00064 std::vector< size_t >& indices, MsqError& err ); 00065 00066 bool average_with_gradient( PatchData& pd, QualityMetric* qm, const std::vector< size_t >& qm_handles, 00067 double& value, std::vector< size_t >& indices, std::vector< Vector3D >& gradient, 00068 MsqError& err ); 00069 00070 bool average_with_Hessian_diagonal( PatchData& pd, QualityMetric* metric, const std::vector< size_t >& qm_handles, 00071 double& value, std::vector< size_t >& indices, 00072 std::vector< Vector3D >& gradient, std::vector< SymMatrix3D >& Hessian_diagonal, 00073 MsqError& err ); 00074 00075 bool average_with_Hessian( PatchData& pd, QualityMetric* metric, const std::vector< size_t >& qm_handles, 00076 double& value, std::vector< size_t >& indices, std::vector< Vector3D >& gradient, 00077 std::vector< Matrix3D >& Hessian, MsqError& err ); 00078 00079 private: 00080 Exponent P; 00081 Exponent P1; 00082 Exponent P2; 00083 std::vector< size_t > mIndices, mOffsets; 00084 std::vector< Vector3D > mGrad; 00085 std::vector< Matrix3D > mHess; 00086 std::vector< SymMatrix3D > mDiag; 00087 std::vector< double > mValues; 00088 }; 00089 00090 } // namespace MBMesquite 00091 00092 #endif