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 AWShapeOrientNB1.cpp 00028 * \brief 00029 * \author Jason Kraftcheck 00030 */ 00031 00032 #include "Mesquite.hpp" 00033 #include "AWShapeOrientNB1.hpp" 00034 #include "MsqMatrix.hpp" 00035 #include "TMPDerivs.hpp" 00036 #include "TMPCommon.hpp" 00037 00038 namespace MBMesquite 00039 { 00040 00041 std::string AWShapeOrientNB1::get_name() const 00042 { 00043 return "AWShapeOrientNB1"; 00044 } 00045 00046 AWShapeOrientNB1::~AWShapeOrientNB1() {} 00047 00048 template < unsigned DIM > 00049 static inline bool eval( const MsqMatrix< DIM, DIM >& A, const MsqMatrix< DIM, DIM >& W, double& result ) 00050 { 00051 result = std::sqrt( sqr_Frobenius( A ) * sqr_Frobenius( W ) ); 00052 result -= A % W; 00053 result *= result; 00054 return true; 00055 } 00056 00057 template < unsigned DIM > 00058 static inline bool grad( const MsqMatrix< DIM, DIM >& A, const MsqMatrix< DIM, DIM >& W, double& result, 00059 MsqMatrix< DIM, DIM >& deriv ) 00060 { 00061 const double nsW = sqr_Frobenius( W ); 00062 const double nsA = sqr_Frobenius( A ); 00063 const double nW = std::sqrt( nsW ); 00064 const double nA = std::sqrt( nsA ); 00065 const double dot = A % W; 00066 result = nA * nW - dot; 00067 result *= result; 00068 00069 deriv = nA * W; 00070 double tmp; 00071 if( divide( dot, nA, tmp ) ) deriv += tmp * A; 00072 deriv *= -nW; 00073 deriv += nsW * A; 00074 deriv += dot * W; 00075 deriv *= 2.0; 00076 return true; 00077 } 00078 00079 /* 00080 template <unsigned DIM> static inline 00081 bool hess( const MsqMatrix<DIM,DIM>& A, 00082 const MsqMatrix<DIM,DIM>& W, 00083 double& result, 00084 MsqMatrix<DIM,DIM>& deriv, 00085 MsqMatrix<DIM,DIM>* second ) 00086 { 00087 } 00088 */ 00089 00090 TMP_AW_TEMPL_IMPL_COMMON_NO2ND( AWShapeOrientNB1 ) 00091 00092 } // namespace MBMesquite