#include <fixedstep.h>

Inheritance diagram for FixedStepFileDrive:

Collaboration diagram for FixedStepFileDrive:

Public Member Functions
	FixedStepFileDrive (unsigned int uL, const DriveHandler pDH, const char const sFileName, integer is, integer id, doublereal t0, doublereal dt, bool bl, bool pz, Drive::Bailout bo)

virtual	~FixedStepFileDrive (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
doublereal	dT0

doublereal	dDT

integer	iNumSteps

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 fixedstep.h.

Constructor & Destructor Documentation

FixedStepFileDrive::FixedStepFileDrive	(	unsigned int	uL,
		const DriveHandler *	pDH,
		const char *const	sFileName,
		integer	is,
		integer	id,
		doublereal	t0,
		doublereal	dt,
		bool	bl,
		bool	pz,
		Drive::Bailout	bo
	)

Definition at line 48 of file fixedstep.cc.

References ASSERT, c, dDT, dFromFile, DriveHandler::dGetTime(), dT0, FIXED_BUFSIZE, FileDrive::iNumDrives, iNumSteps, MBDYN_EXCEPT_ARGS, pd, pvd, SAFENEWARR, ServePending(), and STRLENOF.

 : FileDrive(uL, pDH, sFileName, ind, v0),
 dT0(t0), dDT(dt), iNumSteps(ins),
 bLinear(bl), bPadZeroes(pz), boWhen(bo), pd(0), pvd(0)
 {
         ASSERT(iNumDrives > 0);
         ASSERT(sFileName != NULL);
         ASSERT(dDT > 0.);
 
         std::ifstream in(sFileName);
         if (!in) {
                 silent_cerr("FixedStepFileDrive(" << uL << "): "
                         "can't open file \"" << sFileName << "\""
                         << std::endl);
                 throw ErrGeneric(MBDYN_EXCEPT_ARGS);
         }
 
         /*
          * Mangia gli eventuali commenti iniziali
          */
         unsigned uLineNo = 0;
         char cTmp = '\0';
         while (in.get(cTmp), cTmp == '#') {
                 // increase this to enlarge the size of the buffer
 #define FIXED_BUFSIZE (8192)
                 char tmpbuf[FIXED_BUFSIZE];
 
                 do {
                         in.getline(tmpbuf, sizeof(tmpbuf));
                         ++uLineNo;
                         if (in.fail()) {
                                 silent_cerr("FixedStepFileDrive(" << uL << "): "
                                         "can't get line " << uLineNo << " "
                                         "from file \"" << sFileName << "\"" << std::endl);
                                 throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                         }
                         int idx = 0;
                         while (isspace(tmpbuf[idx])) {
                                 idx++;
                         }
 
                         if (dT0 == dFromFile && strncasecmp(&tmpbuf[idx], "initial time:", STRLENOF("initial time:")) == 0) {
                                 double d;
                                 if (sscanf(&tmpbuf[idx + STRLENOF("initial time:")], "%le", &d) != 1) {
                                         silent_cerr("FixedStepFileDrive(" << uL << "): "
                                                 "can't parse \"initial time\" line \"" << &tmpbuf[idx] << "\" "
                                                 "from file \"" << sFileName << "\"" << std::endl);
                                         throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                                 }
                                 dT0 = d;
 
                         } else if (dDT == dFromFile && strncasecmp(&tmpbuf[idx], "time step:", STRLENOF("time step:")) == 0) {
                                 double d;
                                 if (sscanf(&tmpbuf[idx + STRLENOF("time step:")], "%le", &d) != 1) {
                                         silent_cerr("FixedStepFileDrive(" << uL << "): "
                                                 "can't parse \"time step\" line \"" << &tmpbuf[idx] << "\" "
                                                 "from file \"" << sFileName << "\"" << std::endl);
                                         throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                                 }
                                 if (d <= 0.) {
                                         silent_cerr("FixedStepFileDrive(" << uL << "): "
                                                 "invalid time step value " << d << " "
                                                 "in file \"" << sFileName << "\"" << std::endl);
                                         throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                                 }
 
                                 dDT = d;
                         }
                 } while (strlen(tmpbuf) == STRLENOF(tmpbuf) && tmpbuf[STRLENOF(tmpbuf)] != '\n');
         }
 
         if (dT0 == dFromFile) {
                 silent_cerr("FixedStepFileDrive(" << uL << "): "
                         "expecting \"initial time\" line in file \"" << sFileName << "\""
                         << std::endl);
                 throw ErrGeneric(MBDYN_EXCEPT_ARGS);
         }
 
         if (dDT == dFromFile) {
                 silent_cerr("FixedStepFileDrive(" << uL << "): "
                         "expecting \"time step\" line in file \"" << sFileName << "\""
                         << std::endl);
                 throw ErrGeneric(MBDYN_EXCEPT_ARGS);
         }
 
         if (cTmp != '#') {
                 in.putback(cTmp);
         }
 
         if (ins == -1) {
                 std::streampos pos = in.tellg();
 
                 for (ins = 0; !in.eof(); ins++) {
                         char tmpbuf[FIXED_BUFSIZE];
 
                         do {
                                 in.getline(tmpbuf, sizeof(tmpbuf));
                                 ++uLineNo;
                                 if (in.eof()) {
                                         break;
                                 }
 
                                 if (in.fail()) {
                                         silent_cerr("FixedStepFileDrive(" << uL << "): "
                                                 "can't get line " << uLineNo << " "
                                                 "from file \"" << sFileName << "\"" << std::endl);
                                         throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                                 }
                         } while (strlen(tmpbuf) == STRLENOF(tmpbuf) && tmpbuf[STRLENOF(tmpbuf)] != '\n');
                 }
                 iNumSteps = --ins;
 
                 in.clear();
                 in.seekg(pos);
 
                 silent_cout("FixedStepFileDrive(" << uL << "): "
                         "counted " << ins << " steps" << std::endl);
         }
 
         SAFENEWARR(pd, doublereal, iNumDrives*iNumSteps);
         SAFENEWARR(pvd, doublereal*, iNumDrives + 1);
 
         /* Attenzione: il primo puntatore e' vuoto
          * (ne e' stato allocato uno in piu'),
          * cosi' i drives possono essere numerati da 1 a n */
         for (integer i = iNumDrives; i-- > 0; ) {
                 pvd[i + 1] = pd + i*iNumSteps;
         }
 
         for (integer j = 0; j < iNumSteps; j++) {
                 for (integer i = 1; 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 != iNumDrives) {
                                                 silent_cerr("unexpected end of line #" << j + 1 << " (" << uLineNo << ") after channel #" << i << ", column #" << i << " of file '"
                                                         << sFileName << '\'' << std::endl);
                                                 throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                                         }
 
                                         ++uLineNo;
                                         break;
                                 }
                         }
 
                         if (i == iNumDrives && c != '\n') {
                                 silent_cerr("missing end-of-line at line #" << j + 1 << " (" << uLineNo << ") of file '"
                                         << sFileName << '\'' << std::endl);
                                 throw ErrGeneric(MBDYN_EXCEPT_ARGS);
                         }
 
                         ++uLineNo;
                         in.putback(c);
                 }
         }
 
         /* All data is available, so initialize the buffer accordingly */
         ServePending(pDH->dGetTime());
 }

Here is the call graph for this function:

FixedStepFileDrive::~FixedStepFileDrive ( void )

virtual

Definition at line 220 of file fixedstep.cc.

References pd, pvd, and SAFEDELETEARR.

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

Member Function Documentation

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

virtual

Implements FileDrive.

Definition at line 229 of file fixedstep.cc.

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

void FixedStepFileDrive::ServePending ( const doublereal & t )

virtual

Implements Drive.

Definition at line 236 of file fixedstep.cc.

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

Referenced by FixedStepFileDrive().

 {
         doublereal tt = t - dT0;
 
         if (tt < 0) {
                 if (boWhen & Drive::BO_LOWER) {
                         throw Solver::EndOfSimulation(EXIT_SUCCESS,
                                 MBDYN_EXCEPT_ARGS,
                                 "A fixed 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 (tt > dDT*(iNumSteps - 1)) {
                 if (boWhen & Drive::BO_UPPER) {
                         throw Solver::EndOfSimulation(EXIT_SUCCESS,
                                 MBDYN_EXCEPT_ARGS,
                                 "A fixed 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 {
                 integer j1 = integer(floor(tt/dDT));
                 if (bLinear) {
                         if (j1 == iNumSteps - 1) {
                                 for (int i = 1; i <= iNumDrives; i++) {
                                         pdVal[i] = pvd[i][j1];
                                 }
 
                         } else {
                                 integer j2 = j1 + 1;
                                 doublereal dt1 = dT0 + j1*dDT;
                                 doublereal dt2 = dt1 + dDT;
                                 doublereal dw1 = (dt2 - t)/dDT;
                                 doublereal dw2 = (t - dt1)/dDT;
 
                                 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 FixedStepFileDrive::bLinear

protected

Definition at line 46 of file fixedstep.h.

Referenced by ServePending().

Bailout FixedStepFileDrive::boWhen

protected

Definition at line 48 of file fixedstep.h.

Referenced by ServePending().

bool FixedStepFileDrive::bPadZeroes

protected

Definition at line 47 of file fixedstep.h.

Referenced by ServePending().

doublereal FixedStepFileDrive::dDT

protected

Definition at line 44 of file fixedstep.h.

Referenced by FixedStepFileDrive(), and ServePending().

doublereal FixedStepFileDrive::dT0

protected

Definition at line 43 of file fixedstep.h.

Referenced by FixedStepFileDrive(), and ServePending().

integer FixedStepFileDrive::iNumSteps

protected

Definition at line 45 of file fixedstep.h.

Referenced by FixedStepFileDrive(), and ServePending().

doublereal* FixedStepFileDrive::pd

protected

Definition at line 50 of file fixedstep.h.

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

doublereal** FixedStepFileDrive::pvd

protected

Definition at line 51 of file fixedstep.h.

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

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