Sacado_LFad_LogicalSparse.hpp

Go to the documentation of this file.

// $Id: Sacado_LFad_LogicalSparse.hpp,v 1.2 2008/06/16 21:57:46 etphipp Exp $ 
// $Source: /space/CVS/Trilinos/packages/sacado/src/Sacado_LFad_LogicalSparse.hpp,v $ 
// @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:

      typedef ValT value_type;

      typedef LogT logical_type;


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


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


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


      LogicalSparseImp(const int sz, const int i, const value_type & x) : 
  s_(sz, x) { 
  s_.dx_[i]=logical_type(1); 
      }

      LogicalSparseImp(const LogicalSparseImp& x) : 
  s_(x.s_) {}

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

      ~LogicalSparseImp() {}


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


      void resize(int sz) { s_.resize(sz); }

      void zero() { s_.zero(); }



      const value_type& val() const { return s_.val_;}

      value_type& val() { return s_.val_;}



      int size() const { return s_.size();}

      bool hasFastAccess() const { return s_.size()!=0;}

      bool isPassive() const { return s_.size()!=0;}
      
      void setIsConstant(bool is_const) { 
  if (is_const && s_.size()!=0)
    s_.resize(0);
      }

      const logical_type* dx() const { return &(s_.dx_[0]);}

      logical_type dx(int i) const { 
  return s_.size() ? s_.dx_[i] : logical_type(0); }
    
      logical_type& fastAccessDx(int i) { return s_.dx_[i];}

      const logical_type& fastAccessDx(int i) const { return s_.dx_[i];}
    


      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);


    protected:

      Fad::DynamicStorage<value_type,logical_type> s_;

    }; // 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:



      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& val) {
  LogicalSparseImp< ValT,LogT >::operator=(val);
  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

---