Sacado_LFad_LogicalSparse.hpp

Go to the documentation of this file.

// $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).
// 
// ***********************************************************************
// @HEADER

#ifndef SACADO_LFAD_LOGICALSPARSE_HPP
#define SACADO_LFAD_LOGICALSPARSE_HPP

#include "Sacado_LFad_LogicalSparseTraits.hpp"
#include "Sacado_LFad_ExpressionTraits.hpp"
#include "Sacado_Fad_DynamicStorage.hpp"

namespace Sacado {

  namespace LFad {


    template <typename ExprT> class Expr {};

    template <typename ValT, typename LogT>
    class LogicalSparseImp : 
      public Fad::DynamicStorage<ValT,LogT> {

      typedef Fad::DynamicStorage<ValT,LogT> Storage;

    public:

      typedef ValT value_type;

      typedef LogT logical_type;


      LogicalSparseImp() : Storage(value_type(0)) {}


      LogicalSparseImp(const value_type & x) : Storage(x) {}


      LogicalSparseImp(const int sz, const value_type & x) : Storage(sz, x) {}


      LogicalSparseImp(const int sz, const int i, const value_type & x) : 
  Storage(sz, x) { 
  this->fastAccessDx(i)=logical_type(1); 
      }

      LogicalSparseImp(const LogicalSparseImp& x) : 
  Storage(x) {}

      template <typename S> LogicalSparseImp(const Expr<S>& x);

      ~LogicalSparseImp() {}


      void diff(const int ith, const int n);

      template <typename S>
      bool isEqualTo(const Expr<S>& x) const {
  typedef IsEqual<value_type> IE;
  if (x.size() != this->size()) return false;
  bool eq = IE::eval(x.val(), this->val());
  for (int i=0; i<this->size(); i++)
    eq = eq && IE::eval(x.dx(i), this->dx(i));
  return eq;
      }



      bool hasFastAccess() const { return this->size()!=0;}

      bool isPassive() const { return this->size()!=0;}
      
      void setIsConstant(bool is_const) { 
  if (is_const && this->size()!=0)
    this->resize(0);
      }
    


      LogicalSparseImp& operator=(const value_type& val);

      LogicalSparseImp& 
      operator=(const LogicalSparseImp& x);

      template <typename S> 
      LogicalSparseImp& operator=(const Expr<S>& x); 



      LogicalSparseImp& operator += (const value_type& x);

      LogicalSparseImp& operator -= (const value_type& x);

      LogicalSparseImp& operator *= (const value_type& x);

      LogicalSparseImp& operator /= (const value_type& x);

      template <typename S> 
      LogicalSparseImp& operator += (const Expr<S>& x);

      template <typename S> 
      LogicalSparseImp& operator -= (const Expr<S>& x);
  
      template <typename S> 
      LogicalSparseImp& operator *= (const Expr<S>& x);

      template <typename S> 
      LogicalSparseImp& operator /= (const Expr<S>& x);


    }; // class LogicalSparseImp

    template <typename ValT, typename LogT> 
    class Expr< LogicalSparseImp<ValT,LogT> > : 
      public LogicalSparseImp<ValT,LogT> {

    public:

      Expr() : 
  LogicalSparseImp<ValT,LogT>() {}


      Expr(const ValT & x) : 
  LogicalSparseImp<ValT,LogT>(x) {}


      Expr(const int sz, const ValT & x) : 
  LogicalSparseImp<ValT,LogT>(sz,x) {}


      Expr(const int sz, const int i, const ValT & x) : 
  LogicalSparseImp<ValT,LogT>(sz,i,x) {}

      Expr(const Expr& x) : 
  LogicalSparseImp<ValT,LogT>(x) {}

      template <typename S> Expr(const Expr<S>& x) :
  LogicalSparseImp<ValT,LogT>(x) {}

      ~Expr() {}

    }; // class Expr<LogicalSparseImp>

    template <typename ValT, typename LogT >
    class LogicalSparse : public Expr< LogicalSparseImp<ValT,LogT > > {

    public:

      template <typename T, typename U = LogT> 
      struct apply {
  typedef LogicalSparse<T,U> type;
      };



      LogicalSparse() : 
  Expr< LogicalSparseImp< ValT,LogT > >() {}


      LogicalSparse(const ValT& x) : 
  Expr< LogicalSparseImp< ValT,LogT > >(x) {}


      LogicalSparse(const int sz, const ValT& x) : 
  Expr< LogicalSparseImp< ValT,LogT > >(sz,x) {}


      LogicalSparse(const int sz, const int i, const ValT & x) : 
  Expr< LogicalSparseImp< ValT,LogT > >(sz,i,x) {}

      LogicalSparse(const LogicalSparse& x) : 
  Expr< LogicalSparseImp< ValT,LogT > >(x) {}

      template <typename S> LogicalSparse(const Expr<S>& x) : 
  Expr< LogicalSparseImp< ValT,LogT > >(x) {}


      ~LogicalSparse() {}

      LogicalSparse& operator=(const ValT& v) {
  LogicalSparseImp< ValT,LogT >::operator=(v);
  return *this;
      }

      LogicalSparse& operator=(const LogicalSparse& x) {
  LogicalSparseImp< ValT,LogT >::operator=(static_cast<const LogicalSparseImp< ValT,LogT >&>(x));
  return *this;
      }

      template <typename S> LogicalSparse& operator=(const Expr<S>& x) 
      {
  LogicalSparseImp< ValT,LogT >::operator=(x);
  return *this;
      }
  
    }; // class LogicalSparse<ValT,LogT>

  } // namespace LFad

} // namespace Sacado

#include "Sacado_LFad_LogicalSparseImp.hpp"
#include "Sacado_LFad_LogicalSparseOps.hpp"

#endif // SACADO_LFAD_LOGICALSPARSE_HPP