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 TSizeB1.cpp 00028 * \brief 00029 * \author Jason Kraftcheck 00030 */ 00031 00032 #include "Mesquite.hpp" 00033 #include "TSizeB1.hpp" 00034 #include "MsqMatrix.hpp" 00035 #include "MsqError.hpp" 00036 #include "TMPDerivs.hpp" 00037 #include "TMPCommon.hpp" 00038 00039 namespace MBMesquite 00040 { 00041 00042 std::string TSizeB1::get_name() const 00043 { 00044 return "TSizeB1"; 00045 } 00046 00047 TSizeB1::~TSizeB1() {} 00048 00049 bool TSizeB1::evaluate( const MsqMatrix< 2, 2 >& T, double& result, MsqError& err ) 00050 { 00051 double d = det( T ); 00052 if( TMetric::invalid_determinant( d ) ) 00053 { 00054 MSQ_SETERR( err )( barrier_violated_msg, MsqError::BARRIER_VIOLATED ); 00055 return false; 00056 } 00057 result = d + 1.0 / d - 2.0; 00058 return true; 00059 } 00060 00061 bool TSizeB1::evaluate_with_grad( const MsqMatrix< 2, 2 >& T, double& result, MsqMatrix< 2, 2 >& deriv_wrt_T, 00062 MsqError& err ) 00063 { 00064 double d = det( T ); 00065 if( TMetric::invalid_determinant( d ) ) 00066 { 00067 MSQ_SETERR( err )( barrier_violated_msg, MsqError::BARRIER_VIOLATED ); 00068 return false; 00069 } 00070 result = d + 1.0 / d - 2.0; 00071 deriv_wrt_T = ( 1 - 1 / ( d * d ) ) * transpose_adj( T ); 00072 return true; 00073 } 00074 00075 bool TSizeB1::evaluate_with_hess( const MsqMatrix< 2, 2 >& T, double& result, MsqMatrix< 2, 2 >& deriv_wrt_T, 00076 MsqMatrix< 2, 2 > second_wrt_T[3], MsqError& err ) 00077 { 00078 double d = det( T ); 00079 if( TMetric::invalid_determinant( d ) ) 00080 { 00081 MSQ_SETERR( err )( barrier_violated_msg, MsqError::BARRIER_VIOLATED ); 00082 return false; 00083 } 00084 result = d + 1.0 / d - 2.0; 00085 MsqMatrix< 2, 2 > adjt = transpose_adj( T ); 00086 const double f = 1 - 1 / ( d * d ); 00087 deriv_wrt_T = f * adjt; 00088 00089 set_scaled_outer_product( second_wrt_T, 2 / ( d * d * d ), adjt ); 00090 pluseq_scaled_2nd_deriv_of_det( second_wrt_T, f, T ); 00091 return true; 00092 } 00093 00094 bool TSizeB1::evaluate( const MsqMatrix< 3, 3 >& T, double& result, MsqError& err ) 00095 { 00096 double d = det( T ); 00097 if( TMetric::invalid_determinant( d ) ) 00098 { 00099 MSQ_SETERR( err )( barrier_violated_msg, MsqError::BARRIER_VIOLATED ); 00100 return false; 00101 } 00102 result = d + 1.0 / d - 2.0; 00103 return true; 00104 } 00105 00106 bool TSizeB1::evaluate_with_grad( const MsqMatrix< 3, 3 >& T, double& result, MsqMatrix< 3, 3 >& deriv_wrt_T, 00107 MsqError& err ) 00108 { 00109 double d = det( T ); 00110 if( TMetric::invalid_determinant( d ) ) 00111 { 00112 MSQ_SETERR( err )( barrier_violated_msg, MsqError::BARRIER_VIOLATED ); 00113 return false; 00114 } 00115 result = d + 1.0 / d - 2.0; 00116 deriv_wrt_T = ( 1 - 1 / ( d * d ) ) * transpose_adj( T ); 00117 return true; 00118 } 00119 00120 bool TSizeB1::evaluate_with_hess( const MsqMatrix< 3, 3 >& T, double& result, MsqMatrix< 3, 3 >& deriv_wrt_T, 00121 MsqMatrix< 3, 3 > second_wrt_T[6], MsqError& err ) 00122 { 00123 double d = det( T ); 00124 if( TMetric::invalid_determinant( d ) ) 00125 { 00126 MSQ_SETERR( err )( barrier_violated_msg, MsqError::BARRIER_VIOLATED ); 00127 return false; 00128 } 00129 result = d + 1.0 / d - 2.0; 00130 MsqMatrix< 3, 3 > adjt = transpose_adj( T ); 00131 const double f = 1 - 1 / ( d * d ); 00132 deriv_wrt_T = f * adjt; 00133 00134 set_scaled_outer_product( second_wrt_T, 2 / ( d * d * d ), adjt ); 00135 pluseq_scaled_2nd_deriv_of_det( second_wrt_T, f, T ); 00136 return true; 00137 } 00138 00139 } // namespace MBMesquite