#include <varstep.h>

Inheritance diagram for VariableStepFileDrive:

Collaboration diagram for VariableStepFileDrive:

Public Member Functions
	VariableStepFileDrive (unsigned int uL, const DriveHandler pDH, const char const sFileName, integer id, bool bl, bool pz, Drive::Bailout bo)

virtual	~VariableStepFileDrive (void)

virtual std::ostream &	Restart (std::ostream &out) const

virtual void	ServePending (const doublereal &t)

Public Member Functions inherited from FileDrive
	FileDrive (unsigned int uL, const DriveHandler *pDH, const std::string &s, integer nd, const std::vector< doublereal > &v0)

virtual	~FileDrive (void)

virtual Drive::Type	GetDriveType (void) const

virtual integer	iGetNumDrives (void) const

virtual doublereal	dGet (const doublereal &t, int i=1) const

Public Member Functions inherited from Drive
	Drive (unsigned int uL, const DriveHandler *pDH)

virtual	~Drive (void)

Public Member Functions inherited from WithLabel
	WithLabel (unsigned int uL=0, const std::string &sN="")

virtual	~WithLabel (void)

void	PutLabel (unsigned int uL)

void	PutName (const std::string &sN)

unsigned int	GetLabel (void) const

const std::string &	GetName (void) const

Protected Attributes
integer	iNumSteps

integer	iCurrStep

bool	bLinear

bool	bPadZeroes

Bailout	boWhen

doublereal *	pd

doublereal **	pvd

Protected Attributes inherited from FileDrive
std::string	sFileName

integer	iNumDrives

doublereal *	pdVal

Protected Attributes inherited from Drive
const DriveHandler *	pDrvHdl

Protected Attributes inherited from WithLabel
unsigned int	uLabel

std::string	sName

Additional Inherited Members
Public Types inherited from Drive
enum	Type { UNKNOWN = -1, FILEDRIVE = 0, LASTDRIVETYPE }

enum	Bailout { BO_NONE = 0x0, BO_UPPER = 0x1, BO_LOWER = 0x2, BO_ANY = (BO_UPPER \| BO_LOWER) }

Static Protected Attributes inherited from Drive
static doublereal	dReturnValue = 0.

Detailed Description

Definition at line 41 of file varstep.h.

Constructor & Destructor Documentation

VariableStepFileDrive::VariableStepFileDrive	(	unsigned int	uL,
		const DriveHandler *	pDH,
		const char *const	sFileName,
		integer	id,
		bool	bl,
		bool	pz,
		Drive::Bailout	bo
	)

Definition at line 48 of file varstep.cc.

References ASSERT, buf, c, DriveHandler::dGetTime(), iCurrStep, FileDrive::iNumDrives, iNumSteps, MBDYN_EXCEPT_ARGS, pd, pvd, SAFENEWARR, ServePending(), and STRLENOF.

 : FileDrive(uL, pDH, sFileName, ind, v0),
 iNumSteps(-1), iCurrStep(-1),
 bLinear(bl), bPadZeroes(pz), boWhen(bo), pd(0), pvd(0)
 {
         ASSERT(iNumDrives > 0);
         ASSERT(sFileName != NULL);
 
         std::ifstream in(sFileName);
         if (!in) {
                 silent_cerr("can't open file \""
                         << sFileName << "\"" << std::endl);
                 throw ErrGeneric(MBDYN_EXCEPT_ARGS);
         }
 
         /*
          * Mangia gli eventuali commenti iniziali
          */
         char c = '\0';
         while (in.get(c), c == '#') {
                 char buf[1024];
 
                 do {
                         in.getline(buf, sizeof(buf));
                 } while (strlen(buf) == STRLENOF(buf) && buf[STRLENOF(buf)] != '\n');
         }
 
         if (c != '#') {
                 in.putback(c);
         }
 
         if (iNumSteps == -1) {
                 std::streampos pos = in.tellg();
 
                 integer ins;
                 for (ins = 0; !in.eof(); ins++) {
                         char buf[1024];
 
                         do {
                                 in.getline(buf, sizeof(buf));
                         } while (strlen(buf) == STRLENOF(buf) && buf[STRLENOF(buf)] != '\n');
                 }
                 iNumSteps = --ins;
 
                 in.clear();
                 in.seekg(pos);
 
                 silent_cout("VariableStepFileDrive(" << uL << "): "
                         "counted " << ins << " steps" << std::endl);
         }
 
         SAFENEWARR(pd, doublereal, (1 + iNumDrives)*iNumSteps);
         SAFENEWARR(pvd, doublereal*, 1 + iNumDrives);
 
         for (integer i = iNumDrives + 1; i-- > 0; ) {
                 pvd[i] = pd + i*iNumSteps;
         }
 
         for (integer j = 0; j < iNumSteps; j++) {
                 // 0 -> iNumDrives to account for time
                 for (integer i = 0; i <= iNumDrives; i++) {
                         in >> pvd[i][j];
                         if (in.eof()) {
                                 silent_cerr("unexpected end of file '"
                                         << sFileName << '\'' << std::endl);
                                 throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                         }
 
                         int c;
                         for (c = in.get(); isspace(c); c = in.get()) {
                                 if (c == '\n') {
                                         if (i == 0) {
                                                 silent_cerr("unexpected end of line #" << j + 1 << " after time=" << pvd[0][j] << ", column #" << i + 1 << " of file '"
                                                         << sFileName << '\'' << std::endl);
                                                 throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                                         }
                                         if (i != iNumDrives) {
                                                 silent_cerr("unexpected end of line #" << j + 1 << ", time=" << pvd[0][j] << " after channel #" << i << ", column #" << i + 1 << " of file '"
                                                         << sFileName << '\'' << std::endl);
                                                 throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                                         }
                                         break;
                                 }
                         }
 
                         if (i == iNumDrives && c != '\n') {
                                 silent_cerr("missing end-of-line at line #" << j + 1 << ", time=" << pvd[0][j] << " of file '"
                                         << sFileName << '\'' << std::endl);
                                 throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                         }
 
                         in.putback(c);
                 }
 
                 if (j > 0) {
                         if (pvd[0][j] <= pvd[0][j - 1]) {
                                 silent_cerr("time[" << j << "]=" << pvd[0][j]
                                         << " <= time[" << j - 1 << "]=" << pvd[0][j - 1]
                                         << " in file '" << sFileName << "'" << std::endl);
                                 throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                         }
                 }
         }
 
         // All data is available, so initialize the buffer accordingly
         // use bisection to initialize iCurrStep
         doublereal dTime = pDH->dGetTime();
         iCurrStep = bisec(pvd[0], dTime, 0, iNumSteps - 1);
         if (iCurrStep < 0) {
                 iCurrStep++;
         }
 
         ServePending(dTime);
 }

Here is the call graph for this function:

VariableStepFileDrive::~VariableStepFileDrive ( void )

virtual

Definition at line 166 of file varstep.cc.

References pd, pvd, and SAFEDELETEARR.

 {
         SAFEDELETEARR(pd);
         SAFEDELETEARR(pvd);
 }

Member Function Documentation

std::ostream & VariableStepFileDrive::Restart ( std::ostream & out ) const

virtual

Implements FileDrive.

Definition at line 175 of file varstep.cc.

 {
         return out << "0. /* VariableStepFileDrive: not implemented yet! */"
                 << std::endl;
 }

void VariableStepFileDrive::ServePending ( const doublereal & t )

virtual

Implements Drive.

Definition at line 182 of file varstep.cc.

References ASSERT, bLinear, Drive::BO_LOWER, Drive::BO_UPPER, boWhen, bPadZeroes, iCurrStep, FileDrive::iNumDrives, iNumSteps, MBDYN_EXCEPT_ARGS, FileDrive::pdVal, and pvd.

Referenced by VariableStepFileDrive().

 {
         if (t <= pvd[0][0]) {
                 if (t < pvd[0][0] && (boWhen & Drive::BO_LOWER)) {
                         throw Solver::EndOfSimulation(EXIT_SUCCESS,
                                 MBDYN_EXCEPT_ARGS,
                                 "A variable step file drive lower bound is halting the simulation");
                 }
 
                 if (bPadZeroes) {
                         for (int i = 1; i <= iNumDrives; i++) {
                                 pdVal[i] = 0.;
                         }
 
                 } else {
                         for (int i = 1; i <= iNumDrives; i++) {
                                 pdVal[i] = pvd[i][0];
                         }
                 }
 
         } else if (t >= pvd[0][iNumSteps - 1]) {
                 if (t > pvd[0][iNumSteps - 1] && (boWhen & Drive::BO_UPPER)) {
                         throw Solver::EndOfSimulation(EXIT_SUCCESS,
                                 MBDYN_EXCEPT_ARGS,
                                 "A variable step file drive upper bound is halting the simulation");
                 }
 
                 if (bPadZeroes) {
                         for (int i = 1; i <= iNumDrives; i++) {
                                 pdVal[i] = 0.;
                         }
 
                 } else {
                         for (int i = 1; i <= iNumDrives; i++) {
                                 pdVal[i] = pvd[i][iNumSteps - 1];
                         }
                 }
 
         } else {
                 // look for step exactly before
                 // note: use linear search, under the assumption
                 // that we always start from a relatively close guess
                 while (pvd[0][iCurrStep] > t) {
                         iCurrStep--;
                 }
 
                 // no need to check for under/uverflow
                 ASSERT(iCurrStep >= 0);
 
                 while (pvd[0][iCurrStep + 1] <= t) {
                         iCurrStep++;
                 }
                 ASSERT(iCurrStep < iNumSteps);
 
                 integer j1 = iCurrStep;
                 if (bLinear) {
                         integer j2 = j1 + 1;
                         doublereal dt1 = pvd[0][j1];
                         doublereal dt2 = pvd[0][j2];
                         doublereal dw1 = (dt2 - t)/(dt2 - dt1);
                         doublereal dw2 = (t - dt1)/(dt2 - dt1);
 
                         for (int i = 1; i <= iNumDrives; i++) {
                                 pdVal[i] = pvd[i][j2]*dw2 + pvd[i][j1]*dw1;
                         }
 
                 } else {
                         for (int i = 1; i <= iNumDrives; i++) {
                                 pdVal[i] = pvd[i][j1];
                         }
                 }
         }
 }

Member Data Documentation

bool VariableStepFileDrive::bLinear

protected

Definition at line 45 of file varstep.h.

Referenced by ServePending().

Bailout VariableStepFileDrive::boWhen

protected

Definition at line 47 of file varstep.h.

Referenced by ServePending().

bool VariableStepFileDrive::bPadZeroes

protected

Definition at line 46 of file varstep.h.

Referenced by ServePending().

integer VariableStepFileDrive::iCurrStep

protected

Definition at line 44 of file varstep.h.

Referenced by ServePending(), and VariableStepFileDrive().

integer VariableStepFileDrive::iNumSteps

protected

Definition at line 43 of file varstep.h.

Referenced by ServePending(), and VariableStepFileDrive().

doublereal* VariableStepFileDrive::pd

protected

Definition at line 49 of file varstep.h.

Referenced by VariableStepFileDrive(), and ~VariableStepFileDrive().

doublereal** VariableStepFileDrive::pvd

protected

Definition at line 50 of file varstep.h.

Referenced by ServePending(), VariableStepFileDrive(), and ~VariableStepFileDrive().

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

Public Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation