MOAB: Mesh Oriented datABase
(version 5.4.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 TSizeNB1.cpp 00028 * \brief 00029 * \author Jason Kraftcheck 00030 */ 00031 00032 #include "Mesquite.hpp" 00033 #include "TSizeNB1.hpp" 00034 #include "MsqMatrix.hpp" 00035 #include "TMPDerivs.hpp" 00036 #include "TMPCommon.hpp" 00037 00038 namespace MBMesquite 00039 { 00040 00041 std::string TSizeNB1::get_name() const 00042 { 00043 return "TSizeNB1"; 00044 } 00045 00046 TSizeNB1::~TSizeNB1() {} 00047 00048 template < unsigned DIM > 00049 static inline bool eval( const MsqMatrix< DIM, DIM >& T, double& result ) 00050 { 00051 double d1 = det( T ) - 1; 00052 result = d1 * d1; 00053 return true; 00054 } 00055 00056 template < unsigned DIM > 00057 static inline bool grad( const MsqMatrix< DIM, DIM >& T, double& result, MsqMatrix< DIM, DIM >& deriv ) 00058 { 00059 double d1 = det( T ) - 1; 00060 result = d1 * d1; 00061 deriv = 2 * d1 * transpose_adj( T ); 00062 return true; 00063 } 00064 00065 template < unsigned DIM > 00066 static inline bool hess( const MsqMatrix< DIM, DIM >& T, 00067 double& result, 00068 MsqMatrix< DIM, DIM >& deriv, 00069 MsqMatrix< DIM, DIM >* second ) 00070 { 00071 double d1 = det( T ) - 1; 00072 result = d1 * d1; 00073 const MsqMatrix< DIM, DIM > adjt = transpose_adj( T ); 00074 deriv = 2 * d1 * adjt; 00075 set_scaled_outer_product( second, 2, adjt ); 00076 pluseq_scaled_2nd_deriv_of_det( second, 2 * d1, T ); 00077 return true; 00078 } 00079 00080 TMP_T_TEMPL_IMPL_COMMON( TSizeNB1 ) 00081 00082 } // namespace MBMesquite