Interpolate_Central Class Reference

#include <Interpolate.h>

Inheritance diagram for Interpolate_Central:

Collaboration diagram for Interpolate_Central:

Public Member Functions
	Interpolate_Central (Agent *ag, Agent *bkag, const std::string attr, const std::string alpsuf="_alp")
void	run (double t, double dt, double alpha)
Public Member Functions inherited from InterpolateBase
	InterpolateBase (Agent *ag, Agent *bkag, int cond=0)
virtual void	init (double t)
void	backup ()
Public Member Functions inherited from Action
	Action (void *p=0, char *name=NULL)
	Action (int n, const char *at[], int *i=NULL, void *p=0, char *name=NULL)
	Action (int n, const std::string at[], int *i=NULL, void *p=0, char *name=NULL)
virtual	~Action ()
virtual void	declare (Scheduler &)
virtual void	finalize ()
virtual char *	name ()
void	set_name (const char *name)
virtual void	print (FILE *f)
virtual void	print_toposort (FILE *f)
virtual void	schedule ()

Additional Inherited Members
Protected Types inherited from Action
enum	{ IN =1, OUT =2, INOUT =3 }
Protected Member Functions inherited from InterpolateBase
void	extrapolate_Linear (double dt, double dt_old, double time_old, int a_old, double time_new, int a_new, double time_out, int a_out, int a_gra=-100)
Protected Member Functions inherited from Action
int	get_attribute_handle (int i)
int	get_attribute_handle_const (int i)
int	get_attribute_handle (const std::string str)
void	set_attr (int n, const std::string at[], int *id=NULL)
void	set_attr (int n, const char *at[], int *id=NULL)
void	set_io (int n, const int *io)
void	set_io (const char *io)
int	get_io (int i)
Protected Attributes inherited from InterpolateBase
Agent *	agent
Agent *	bkagent
int	attr_hdls [4]
int	bkup_hdls [3]
int	conditional
int	on
Protected Attributes inherited from Action
char *	action_name
char **	attr
int *	idx
int	count
void *	usr_ptr
std::vector< int >	inout

Detailed Description

Definition at line 76 of file Interpolate.h.

Constructor & Destructor Documentation

Interpolate_Central	(	Agent *	ag,
		Agent *	bkag,
		const std::string	attr,
		const std::string	alpsuf = "_alp"
	)

Definition at line 345 of file Interpolate.C.

References Action::action_name, Action::attr, Action::IN, Action::OUT, Action::set_attr(), and Action::set_io().

                                                                                                            : InterpolateBase(ag, bkag)
{
  action_name = (char *)"Interpolate_Central";
    // INIT_INTERP_HANDLES() in "man_basic.f90"
  std::string atts[4];
  atts[0] = attr;
  atts[1] = attr + "_old";
  atts[2] = attr +  alpsuf;
  atts[3] = attr + "_grad";

  set_attr(4, atts);

  // Define I/O of the arguments.
  int io[] = {IN, IN, OUT, IN};
  set_io( 4, io);

/*
    // register attributes
  std::string surf_win = agent->get_surface_window();
  std::string::size_type pos = attr.find(".");
  std::string wname, aname;
  if ( pos == std::string::npos) {
    COM_assertion_msg(0, "ERROR: Here!");
  }
  else {
    wname = attr.substr( 0, pos);
    aname = attr.substr( pos, attr.size()-pos);
  }

  if (attr_action[0] == 'n')
    agent->register_new_attribute( wname, aname, 'e', COM_DOUBLE, 1, "kg/(m^2 s)");

  if (attr_action[1] == 'n')
    agent->register_new_attribute( wname, aname+"_old", 'e', COM_DOUBLE, 1, "kg/(m^2 s)");
  else if (attr_action[1] == 'c')
    agent->register_clone_attribute( 0, wname, aname+"_old", surf_win, aname);   // attr_old

  if (attr_action[2] == 'c')
    agent->register_clone_attribute( 0, wname, aname+alpsuf, surf_win, aname);  // attr_grad

  if (attr_action[3] == 'n')
    agent->register_new_attribute( wname, aname+"_grad", 'e', COM_DOUBLE, 1, "kg/(m^2 s)");
  else if (attr_action[3] == 'c')
    agent->register_clone_attribute( 0, wname, aname+"_grad", surf_win, aname);  // attr_grad
*/
}

Here is the call graph for this function:

Member Function Documentation

void run ( double  t, double  dt, double  alpha  )

virtual

Implements InterpolateBase.

Definition at line 392 of file Interpolate.C.

References InterpolateBase::agent, Action::attr, InterpolateBase::attr_hdls, COM_assertion_msg, COM_call_function(), RocBlas::copy, InterpolateBase::extrapolate_Linear(), Agent::get_coupling(), Agent::get_old_dt(), Coupling::get_rocmancontrol_param(), MAN_DEBUG, Coupling::new_start(), and RocmanControl_parameters::order.

{
  if (attr_hdls[2] <= 0) return;

  MAN_DEBUG(3, ("Interpolate_Central::run (%s) called with t:%e dt:%e alpha:%e.\n", attr[0], t, dt, alpha));

  COM_assertion_msg(alpha>=0.0, "ERROR: Extrapolate_Central called with invalid alpha!");

  COM_assertion_msg(alpha>=-1.e-6 && alpha <= 1+1.e-6, "ERROR: Abort!");

  int order = agent->get_coupling()->get_rocmancontrol_param()->order;
  if (order == 0 || agent->get_coupling()->new_start(t)) {
    COM_call_function( RocBlas::copy, &attr_hdls[0], &attr_hdls[2]);
    return;
  }

  double base = -0.5;

  // linear interpolation.
  extrapolate_Linear(dt, agent->get_old_dt(), base, attr_hdls[1], base+1.0, attr_hdls[0], alpha, attr_hdls[2], attr_hdls[3]);
}

Here is the call graph for this function:

---

The documentation for this class was generated from the following files:

• Interpolate.h
• Interpolate.C