MOAB: Mesh Oriented datABase
(version 5.4.1)
|
00001 /* ***************************************************************** 00002 MESQUITE -- The Mesh Quality Improvement Toolkit 00003 00004 Copyright 2009 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 (2009) kraftche@cae.wisc.edu 00024 00025 ***************************************************************** */ 00026 00027 /** \file TSum.cpp 00028 * \brief 00029 * \author Jason Kraftcheck 00030 */ 00031 00032 #include "Mesquite.hpp" 00033 #include "TMetricBarrier.hpp" 00034 #include "TSum.hpp" 00035 #include "MsqMatrix.hpp" 00036 #include "MsqError.hpp" 00037 #include "TMPCommon.hpp" 00038 00039 namespace MBMesquite 00040 { 00041 00042 std::string TSum::get_name() const 00043 { 00044 return mu1->get_name() + '+' + mu2->get_name(); 00045 } 00046 00047 TSum::~TSum() {} 00048 00049 template < unsigned DIM > 00050 inline bool TSum::eval( const MsqMatrix< DIM, DIM >& T, double& result, MsqError& err ) 00051 { 00052 double val2; 00053 bool rval = mu1->evaluate( T, result, err ); 00054 MSQ_ERRZERO( err ); 00055 bool rval2 = mu2->evaluate( T, val2, err ); 00056 MSQ_ERRZERO( err ); 00057 result += val2; 00058 return rval && rval2; 00059 } 00060 00061 template < unsigned DIM > 00062 inline bool TSum::grad( const MsqMatrix< DIM, DIM >& T, double& result, MsqMatrix< DIM, DIM >& deriv, MsqError& err ) 00063 { 00064 double val2; 00065 MsqMatrix< DIM, DIM > grad2; 00066 bool rval = mu1->evaluate_with_grad( T, result, deriv, err ); 00067 MSQ_ERRZERO( err ); 00068 bool rval2 = mu2->evaluate_with_grad( T, val2, grad2, err ); 00069 MSQ_ERRZERO( err ); 00070 result += val2; 00071 deriv += grad2; 00072 return rval && rval2; 00073 } 00074 00075 template < unsigned DIM > 00076 inline bool TSum::hess( const MsqMatrix< DIM, DIM >& T, 00077 double& result, 00078 MsqMatrix< DIM, DIM >& deriv_wrt_T, 00079 MsqMatrix< DIM, DIM >* second_wrt_T, 00080 MsqError& err ) 00081 { 00082 const int HL = ( DIM * ( DIM + 1 ) ) / 2; 00083 double val2; 00084 MsqMatrix< DIM, DIM > grad2, hess2[HL]; 00085 00086 bool rval = mu1->evaluate_with_hess( T, result, deriv_wrt_T, second_wrt_T, err ); 00087 MSQ_ERRZERO( err ); 00088 bool rval2 = mu2->evaluate_with_hess( T, val2, grad2, hess2, err ); 00089 MSQ_ERRZERO( err ); 00090 result += val2; 00091 deriv_wrt_T += grad2; 00092 for( int i = 0; i < HL; ++i ) 00093 second_wrt_T[i] += hess2[i]; 00094 return rval && rval2; 00095 } 00096 00097 TMP_T_TEMPL_IMPL_COMMON_ERR( TSum ) 00098 00099 } // namespace MBMesquite