Branch data     Line data    Source code

       1                 :            : /* *****************************************************************
       2                 :            :     MESQUITE -- The Mesh Quality Improvement Toolkit
       3                 :            : 
       4                 :            :     Copyright 2006 Sandia National Laboratories.  Developed at the
       5                 :            :     University of Wisconsin--Madison under SNL contract number
       6                 :            :     624796.  The U.S. Government and the University of Wisconsin
       7                 :            :     retain certain rights to this software.
       8                 :            : 
       9                 :            :     This library is free software; you can redistribute it and/or
      10                 :            :     modify it under the terms of the GNU Lesser General Public
      11                 :            :     License as published by the Free Software Foundation; either
      12                 :            :     version 2.1 of the License, or (at your option) any later version.
      13                 :            : 
      14                 :            :     This library is distributed in the hope that it will be useful,
      15                 :            :     but WITHOUT ANY WARRANTY; without even the implied warranty of
      16                 :            :     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
      17                 :            :     Lesser General Public License for more details.
      18                 :            : 
      19                 :            :     You should have received a copy of the GNU Lesser General Public License
      20                 :            :     (lgpl.txt) along with this library; if not, write to the Free Software
      21                 :            :     Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
      22                 :            : 
      23                 :            :     (2006) [email protected]
      24                 :            : 
      25                 :            :   ***************************************************************** */
      26                 :            : 
      27                 :            : /** \file TSizeNB1.cpp
      28                 :            :  *  \brief
      29                 :            :  *  \author Jason Kraftcheck
      30                 :            :  */
      31                 :            : 
      32                 :            : #include "Mesquite.hpp"
      33                 :            : #include "TSizeNB1.hpp"
      34                 :            : #include "MsqMatrix.hpp"
      35                 :            : #include "TMPDerivs.hpp"
      36                 :            : #include "TMPCommon.hpp"
      37                 :            : 
      38                 :            : namespace MBMesquite
      39                 :            : {
      40                 :            : 
      41                 :          7 : std::string TSizeNB1::get_name() const
      42                 :            : {
      43         [ +  - ]:          7 :     return "TSizeNB1";
      44                 :            : }
      45                 :            : 
      46         [ -  + ]:         40 : TSizeNB1::~TSizeNB1() {}
      47                 :            : 
      48                 :            : template < unsigned DIM >
      49                 :    2341956 : static inline bool eval( const MsqMatrix< DIM, DIM >& T, double& result )
      50                 :            : {
      51                 :    2341956 :     double d1 = det( T ) - 1;
      52                 :    2341956 :     result    = d1 * d1;
      53                 :    2341956 :     return true;
      54                 :            : }
      55                 :            : 
      56                 :            : template < unsigned DIM >
      57                 :     938836 : static inline bool grad( const MsqMatrix< DIM, DIM >& T, double& result, MsqMatrix< DIM, DIM >& deriv )
      58                 :            : {
      59                 :     938836 :     double d1 = det( T ) - 1;
      60                 :     938836 :     result    = d1 * d1;
      61   [ #  #  +  - ]:     938836 :     deriv     = 2 * d1 * transpose_adj( T );
      62                 :     938836 :     return true;
      63                 :            : }
      64                 :            : 
      65                 :            : template < unsigned DIM >
      66                 :          0 : static inline bool hess( const MsqMatrix< DIM, DIM >& T, double& result, MsqMatrix< DIM, DIM >& deriv,
      67                 :            :                          MsqMatrix< DIM, DIM >* second )
      68                 :            : {
      69 [ #  # ][ #  # ]:          0 :     double d1                        = det( T ) - 1;
      70                 :          0 :     result                           = d1 * d1;
      71 [ #  # ][ #  # ]:          0 :     const MsqMatrix< DIM, DIM > adjt = transpose_adj( T );
      72 [ #  # ][ #  # ]:          0 :     deriv                            = 2 * d1 * adjt;
      73 [ #  # ][ #  # ]:          0 :     set_scaled_outer_product( second, 2, adjt );
      74 [ #  # ][ #  # ]:          0 :     pluseq_scaled_2nd_deriv_of_det( second, 2 * d1, T );
      75                 :          0 :     return true;
      76                 :            : }
      77                 :            : 
      78                 :    6561584 : TMP_T_TEMPL_IMPL_COMMON( TSizeNB1 )
      79                 :            : 
      80                 :            : }  // namespace MBMesquite