/usr/include/trilinos/Sacado_ELRFad_Expression.hpp is in libtrilinos-dev 10.4.0.dfsg-1ubuntu2.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 | // $Id$
// $Source$
// @HEADER
// ***********************************************************************
//
// Sacado Package
// Copyright (2006) Sandia Corporation
//
// Under the terms of Contract DE-AC04-94AL85000 with Sandia Corporation,
// the U.S. Government retains certain rights in this software.
//
// This library is free software; you can redistribute it and/or modify
// it under the terms of the GNU Lesser General Public License as
// published by the Free Software Foundation; either version 2.1 of the
// License, or (at your option) any later version.
//
// This library is distributed in the hope that it will be useful, but
// WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
// Lesser General Public License for more details.
//
// You should have received a copy of the GNU Lesser General Public
// License along with this library; if not, write to the Free Software
// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307
// USA
// Questions? Contact David M. Gay (dmgay@sandia.gov) or Eric T. Phipps
// (etphipp@sandia.gov).
//
// ***********************************************************************
//
// The forward-mode AD classes in Sacado are a derivative work of the
// expression template classes in the Fad package by Nicolas Di Cesare.
// The following banner is included in the original Fad source code:
//
// ************ DO NOT REMOVE THIS BANNER ****************
//
// Nicolas Di Cesare <Nicolas.Dicesare@ann.jussieu.fr>
// http://www.ann.jussieu.fr/~dicesare
//
// CEMRACS 98 : C++ courses,
// templates : new C++ techniques
// for scientific computing
//
//********************************************************
//
// A short implementation ( not all operators and
// functions are overloaded ) of 1st order Automatic
// Differentiation in forward mode (FAD) using
// EXPRESSION TEMPLATES.
//
//********************************************************
// @HEADER
#ifndef SACADO_ELRFAD_EXPRESSION_HPP
#define SACADO_ELRFAD_EXPRESSION_HPP
#include "Sacado_Traits.hpp"
namespace Sacado {
namespace ELRFad {
//! Base template specification for %ExprPromote
/*!
* The %ExprPromote classes provide a mechanism for computing the
* promoted expression-type of a binary operation.
*/
template <typename A, typename B> struct ExprPromote {};
//! Specialization of %ExprPromote for a single type
template <typename A> struct ExprPromote<A,A> {
typedef A type;
};
//! Wrapper for a generic expression template
/*!
* This template class serves as a wrapper for all Fad expression
* template classes.
*/
template <typename ExprT>
class Expr {};
//! Constant expression template
/*!
* This template class represents a constant expression.
*/
template <typename ConstT>
class ConstExpr {
public:
//! Typename of argument values
typedef ConstT value_type;
//! Typename of base-expressions
typedef ConstT base_expr_type;
//! Number of arguments
static const int num_args = 0;
//! Constructor
ConstExpr(const ConstT& constant) : constant_(constant) {}
//! Return size of the derivative array of the operation
int size() const { return 0; }
//! Return value of operation
value_type val() const { return constant_; }
//! Return partials w.r.t. arguments
void computePartials(const value_type& bar,
value_type partials[]) const {}
//! Rturn tangent component \c i of arguments
void getTangents(int i, value_type dots[]) const {}
//! Return tangent component \c i of argument \c Arg
template <int Arg>
value_type getTangent(int i) const { return 0.0; }
//! Return whether argument is active
template <int Arg>
bool isActive() const { return false; }
protected:
//! The constant
const ConstT& constant_;
}; // class ConstExpr
template <typename T> struct ExprConstRef {
typedef const T& type;
};
template <typename T> struct ExprConstRef< ConstExpr<T> > {
typedef const ConstExpr<T> type;
};
} // namespace ELRFad
} // namespace Sacado
#endif // SACADO_ELRFAD_EXPRESSION_HPP
|