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 AWSizeB1.cpp 00028 * \brief 00029 * \author Jason Kraftcheck 00030 */ 00031 00032 #include "Mesquite.hpp" 00033 #include "AWSizeB1.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 AWSizeB1::get_name() const 00043 { 00044 return "AWSizeB1"; 00045 } 00046 00047 AWSizeB1::~AWSizeB1() {} 00048 00049 bool AWSizeB1::evaluate( const MsqMatrix< 2, 2 >& A, const MsqMatrix< 2, 2 >& W, double& result, MsqError& err ) 00050 { 00051 const double alpha = det( A ); 00052 const double omega = det( W ); 00053 const double prod = alpha * omega; 00054 if( AWMetric::invalid_determinant( prod ) ) 00055 { 00056 MSQ_SETERR( err )( barrier_violated_msg_aw, MsqError::BARRIER_VIOLATED ); 00057 return false; 00058 } 00059 00060 result = ( alpha - omega ); 00061 result *= result; 00062 result /= prod; 00063 return true; 00064 } 00065 00066 bool AWSizeB1::evaluate_with_grad( const MsqMatrix< 2, 2 >& A, const MsqMatrix< 2, 2 >& W, double& result, 00067 MsqMatrix< 2, 2 >& deriv_wrt_A, MsqError& err ) 00068 { 00069 const double alpha = det( A ); 00070 const double omega = det( W ); 00071 const double prod = alpha * omega; 00072 if( AWMetric::invalid_determinant( prod ) ) 00073 { 00074 MSQ_SETERR( err )( barrier_violated_msg_aw, MsqError::BARRIER_VIOLATED ); 00075 return false; 00076 } 00077 00078 result = ( alpha - omega ); 00079 result *= result; 00080 result /= prod; 00081 deriv_wrt_A = transpose_adj( A ); 00082 deriv_wrt_A *= ( alpha * alpha - omega * omega ) / ( alpha * prod ); 00083 return true; 00084 } 00085 00086 bool AWSizeB1::evaluate( const MsqMatrix< 3, 3 >& A, const MsqMatrix< 3, 3 >& W, double& result, MsqError& err ) 00087 { 00088 const double alpha = det( A ); 00089 const double omega = det( W ); 00090 const double prod = alpha * omega; 00091 if( AWMetric::invalid_determinant( prod ) ) 00092 { 00093 MSQ_SETERR( err )( barrier_violated_msg_aw, MsqError::BARRIER_VIOLATED ); 00094 return false; 00095 }; 00096 00097 result = ( alpha - omega ); 00098 result *= result; 00099 result /= prod; 00100 return true; 00101 } 00102 00103 bool AWSizeB1::evaluate_with_grad( const MsqMatrix< 3, 3 >& A, const MsqMatrix< 3, 3 >& W, double& result, 00104 MsqMatrix< 3, 3 >& deriv_wrt_A, MsqError& err ) 00105 { 00106 const double alpha = det( A ); 00107 const double omega = det( W ); 00108 const double prod = alpha * omega; 00109 if( AWMetric::invalid_determinant( prod ) ) 00110 { 00111 MSQ_SETERR( err )( barrier_violated_msg_aw, MsqError::BARRIER_VIOLATED ); 00112 return false; 00113 } 00114 00115 result = ( alpha - omega ); 00116 result *= result; 00117 result /= prod; 00118 deriv_wrt_A = transpose_adj( A ); 00119 deriv_wrt_A *= ( alpha * alpha - omega * omega ) / ( alpha * prod ); 00120 return true; 00121 } 00122 00123 } // namespace MBMesquite