Time_Def.hpp

Go to the documentation of this file.

//*****************************************************************//
//    Albany 2.0:  Copyright 2012 Sandia Corporation               //
//    This Software is released under the BSD license detailed     //
//    in the file "license.txt" in the top-level Albany directory  //
//*****************************************************************//

#include <fstream>
#include "Teuchos_TestForException.hpp"
#include "Phalanx_DataLayout.hpp"
#include "Sacado_ParameterRegistration.hpp"
#include "Albany_Utils.hpp"

namespace LCM {

template<typename EvalT, typename Traits>
Time<EvalT, Traits>::
Time(Teuchos::ParameterList& p) :
  time      (p.get<std::string>("Time Name"),
       p.get<Teuchos::RCP<PHX::DataLayout> >("Workset Scalar Data Layout")),
  deltaTime (p.get<std::string>("Delta Time Name"),
       p.get<Teuchos::RCP<PHX::DataLayout> >("Workset Scalar Data Layout")),
  timeValue(0.0)
{
  if (p.isType<bool>("Disable Transient"))
    enableTransient = !p.get<bool>("Disable Transient");
  else enableTransient = true;

  // Add Time as a Sacado-ized parameter
  Teuchos::RCP<ParamLib> paramLib = 
    p.get< Teuchos::RCP<ParamLib> >("Parameter Library", Teuchos::null);
  new Sacado::ParameterRegistration<EvalT, SPL_Traits>("Time", this, paramLib);

  this->addEvaluatedField(time);
  this->addEvaluatedField(deltaTime);

  timeName = p.get<std::string>("Time Name")+"_old";;
  this->setName("Time"+PHX::TypeString<EvalT>::value);
}

// **********************************************************************
template<typename EvalT, typename Traits>
void Time<EvalT, Traits>::
postRegistrationSetup(typename Traits::SetupData d,
                      PHX::FieldManager<Traits>& fm)
{
  this->utils.setFieldData(time,fm);
  this->utils.setFieldData(deltaTime,fm);
}

// **********************************************************************
template<typename EvalT, typename Traits>
void Time<EvalT, Traits>::
evaluateFields(typename Traits::EvalData workset)
{
  time(0) = workset.current_time;

  Albany::MDArray timeOld = (*workset.stateArrayPtr)[timeName];
  deltaTime(0) = time(0) - timeOld(0);
  //std::cout << "time = " << time(0) << ", timeOld = " << timeOld(0) << "\n";
}

// **********************************************************************
template<typename EvalT,typename Traits>
typename Time<EvalT,Traits>::ScalarT& 
Time<EvalT,Traits>::getValue(const std::string &n)
{
  return timeValue;
}

// **********************************************************************
// **********************************************************************
}