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 (2010) kraftche@cae.wisc.edu 00024 00025 ***************************************************************** */ 00026 00027 /** \file AWUntangleBeta.hpp 00028 * \brief 00029 * \author Jason Kraftcheck 00030 */ 00031 00032 #ifndef MSQ_AW_UNTANGLE_BETA_HPP 00033 #define MSQ_AW_UNTANGLE_BETA_HPP 00034 00035 #include "Mesquite.hpp" 00036 #include "AWMetric.hpp" 00037 00038 namespace MBMesquite 00039 { 00040 00041 /** \f$ \left\{ |\alpha - \gamma \, \omega| - ( \alpha - \gamma \, \omega ) \right\}^p \f$ */ 00042 class AWUntangleBeta : public AWMetric 00043 { 00044 private: 00045 double mGamma; 00046 00047 public: 00048 MESQUITE_EXPORT 00049 AWUntangleBeta( double gamma = 0.5 ) : mGamma( gamma ) {} 00050 00051 MESQUITE_EXPORT virtual ~AWUntangleBeta(); 00052 00053 MESQUITE_EXPORT virtual std::string get_name() const; 00054 00055 MESQUITE_EXPORT virtual bool evaluate( const MsqMatrix< 2, 2 >& A, const MsqMatrix< 2, 2 >& W, double& result, 00056 MsqError& err ); 00057 00058 MESQUITE_EXPORT virtual bool evaluate( const MsqMatrix< 3, 3 >& A, const MsqMatrix< 3, 3 >& W, double& result, 00059 MsqError& err ); 00060 00061 MESQUITE_EXPORT virtual bool evaluate_with_grad( const MsqMatrix< 2, 2 >& A, const MsqMatrix< 2, 2 >& W, 00062 double& result, MsqMatrix< 2, 2 >& deriv_wrt_A, MsqError& err ); 00063 00064 MESQUITE_EXPORT virtual bool evaluate_with_grad( const MsqMatrix< 3, 3 >& A, const MsqMatrix< 3, 3 >& W, 00065 double& result, MsqMatrix< 3, 3 >& deriv_wrt_A, MsqError& err ); 00066 /* 00067 MESQUITE_EXPORT virtual 00068 bool evaluate_with_hess( const MsqMatrix<2,2>& A, 00069 const MsqMatrix<2,2>& W, 00070 double& result, 00071 MsqMatrix<2,2>& deriv_wrt_A, 00072 MsqMatrix<2,2> second_wrt_A[3], 00073 MsqError& err ); 00074 00075 MESQUITE_EXPORT virtual 00076 bool evaluate_with_hess( const MsqMatrix<3,3>& A, 00077 const MsqMatrix<3,3>& W, 00078 double& result, 00079 MsqMatrix<3,3>& deriv_wrt_A, 00080 MsqMatrix<3,3> second_wrt_A[6], 00081 MsqError& err ); 00082 */ 00083 00084 private: 00085 template < unsigned D > 00086 inline bool eval( const MsqMatrix< D, D >& A, const MsqMatrix< D, D >& W, double& result ); 00087 template < unsigned D > 00088 inline bool grad( const MsqMatrix< D, D >& A, const MsqMatrix< D, D >& W, double& result, 00089 MsqMatrix< D, D >& first ); 00090 // template <unsigned D> inline 00091 // bool hess( const MsqMatrix<D,D>& T, double& result, MsqMatrix<D,D>& first, MsqMatrix<D,D>* 00092 // second ); 00093 }; 00094 00095 } // namespace MBMesquite 00096 00097 #endif