// saves full information stored in this status
// temporary variables are NOT stored
contextIOResultType
MisesMatStatus :: saveContext(DataStream &stream, ContextMode mode, void *obj)
{
contextIOResultType iores;
// save parent class status
if ( ( iores = StructuralMaterialStatus :: saveContext(stream, mode, obj) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
// write raw data
// write plastic strain (vector)
if ( ( iores = plasticStrain.storeYourself(stream) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
// write cumulative plastic strain (scalar)
if ( !stream.write(kappa) ) {
THROW_CIOERR(CIO_IOERR);
}
// write damage (scalar)
if ( !stream.write(damage) ) {
THROW_CIOERR(CIO_IOERR);
}
return CIO_OK;
}
contextIOResultType
FiberedCrossSection :: restoreIPContext(DataStream &stream, ContextMode mode, GaussPoint *masterGp)
//
// restores full material context (saves state variables, that completely describe
// current state)
//
// restores also slaves of master gp
//
{
contextIOResultType iores;
if ( ( iores = CrossSection :: restoreIPContext(stream, mode, masterGp) ) != CIO_OK ) {
THROW_CIOERR(iores); // saved masterGp
}
for ( int i = 1; i <= numberOfFibers; i++ ) {
// creates also slaves if they don't exists
GaussPoint *slaveGP = this->giveSlaveGaussPoint(masterGp, i - 1);
StructuralMaterial *mat = dynamic_cast< StructuralMaterial * >( domain->giveMaterial( fiberMaterials.at(i) ) );
if ( ( iores = mat->restoreIPContext(stream, mode, slaveGP) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
}
return CIO_OK;
}
contextIOResultType
RCSDMaterialStatus :: saveContext(DataStream *stream, ContextMode mode, void *obj)
//
// saves full information stored in this Status
// no temp variables stored
//
{
contextIOResultType iores;
// save parent class status
if ( ( iores = RCM2MaterialStatus :: saveContext(stream, mode, obj) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
// write a raw data
if ( !stream->write(& maxEquivStrain, 1) ) {
THROW_CIOERR(CIO_IOERR);
}
if ( !stream->write(& damageCoeff, 1) ) {
THROW_CIOERR(CIO_IOERR);
}
if ( !stream->write(& mode, 1) ) {
THROW_CIOERR(CIO_IOERR);
}
if ( ( iores = Ds0.storeYourself(stream, mode) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
return CIO_OK;
}
contextIOResultType XfemManager :: restoreContext(DataStream *stream, ContextMode mode, void *obj)
{
contextIOResultType iores;
if ( mode & CM_Definition ) {
if ( !stream->read(& this->numberOfEnrichmentItems, 1) ) {
THROW_CIOERR(CIO_IOERR);
}
this->enrichmentItemList.resize(this->numberOfEnrichmentItems);
}
for ( int i = 1; i <= this->numberOfEnrichmentItems; i++ ) {
EnrichmentItem *obj;
if ( mode & CM_Definition ) {
std :: string name;
if ( !stream->read(name) ) {
THROW_CIOERR(CIO_IOERR);
}
std :: unique_ptr< EnrichmentItem > ei( classFactory.createEnrichmentItem(name.c_str(), i, this, this->domain) );
obj = ei.get();
enrichmentItemList.insert(enrichmentItemList.begin() + i-1, std :: move(ei));
} else {
obj = this->giveEnrichmentItem(i);
}
if ( ( iores = obj->restoreContext(stream, mode) ) != CIO_OK ) {
THROW_CIOERR(CIO_IOERR);
}
}
return CIO_OK;
}
contextIOResultType Set :: restoreContext(DataStream &stream, ContextMode mode, void *obj)
{
contextIOResultType iores;
if ( ( iores = FEMComponent :: restoreContext(stream, mode, obj) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
if ( mode & CM_Definition ) {
if ( ( iores = elements.restoreYourself(stream) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
if ( ( iores = elementBoundaries.restoreYourself(stream) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
if ( ( iores = elementEdges.restoreYourself(stream) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
if ( ( iores = nodes.restoreYourself(stream) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
}
this->totalNodes.clear();
return CIO_OK;
}
contextIOResultType SimpleSlaveDof :: saveContext(DataStream &stream, ContextMode mode, void *obj)
//
// saves full node context (saves state variables, that completely describe
// current state)
//
{
contextIOResultType iores;
if ( ( iores = Dof :: saveContext(stream, mode, obj) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
if ( mode & CM_Definition ) {
if ( mode & CM_DefinitionGlobal ) {
int _masterGlobNum = dofManager->giveDomain()->giveDofManager(masterDofMngr)->giveGlobalNumber();
if ( !stream.write(_masterGlobNum) ) {
THROW_CIOERR(CIO_IOERR);
}
} else {
if ( !stream.write(masterDofMngr) ) {
THROW_CIOERR(CIO_IOERR);
}
}
}
return CIO_OK;
}
contextIOResultType
RheoChainMaterialStatus :: restoreContext(DataStream *stream, ContextMode mode, void *obj)
//
// restore the state variables from a stream
//
{
contextIOResultType iores;
if ( ( iores = StructuralMaterialStatus :: restoreContext(stream, mode, obj) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
// read raw data
for ( int i = 0; i < nUnits; i++ ) {
if ( ( iores = hiddenVars [ i ].restoreYourself(stream, mode) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
}
if ( ( iores = shrinkageStrain.restoreYourself(stream, mode) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
return CIO_OK;
}
contextIOResultType
TransportMaterialStatus :: restoreContext(DataStream &stream, ContextMode mode, void *obj)
//
// restores full material context (saves state variables, that completely describe
// current state)
//
{
contextIOResultType iores;
if ( ( iores = MaterialStatus :: restoreContext(stream, mode, obj) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
if ( ( iores = gradient.restoreYourself(stream) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
if ( ( iores = field.restoreYourself(stream) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
if ( ( iores = flux.restoreYourself(stream) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
return CIO_OK;
}
contextIOResultType StaticStructural :: restoreContext(DataStream *stream, ContextMode mode, void *obj)
{
contextIOResultType iores;
int closeFlag = 0;
int istep, iversion;
FILE *file = NULL;
this->resolveCorrespondingStepNumber(istep, iversion, obj);
if ( stream == NULL ) {
if ( !this->giveContextFile(& file, istep, iversion, contextMode_read) ) {
THROW_CIOERR(CIO_IOERR); // override
}
stream = new FileDataStream(file);
closeFlag = 1;
}
if ( ( iores = StructuralEngngModel :: restoreContext(stream, mode, obj) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
if ( ( iores = this->field->restoreContext(*stream, mode) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
if ( closeFlag ) {
fclose(file);
delete stream;
stream = NULL;
}
return CIO_OK;
}
contextIOResultType MasterDof :: saveContext(DataStream &stream, ContextMode mode, void *obj)
//
// saves full node context (saves state variables, that completely describe
// current state)
//
{
contextIOResultType iores;
if ( ( iores = Dof :: saveContext(stream, mode, obj) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
if ( mode & CM_Definition ) {
if ( !stream.write(bc) ) {
THROW_CIOERR(CIO_IOERR);
}
if ( !stream.write(ic) ) {
THROW_CIOERR(CIO_IOERR);
}
}
// store equation number of receiver
if ( !stream.write(equationNumber) ) {
THROW_CIOERR(CIO_IOERR);
}
if ( ( mode & CM_UnknownDictState ) || ( dofManager->giveDomain()->giveEngngModel()->requiresUnknownsDictionaryUpdate() ) ) {
if ( ( iores = unknowns.saveContext(stream, mode, obj) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
}
return CIO_OK;
}
contextIOResultType
PrimaryField :: restoreContext(DataStream &stream, ContextMode mode)
{
contextIOResultType iores(CIO_IOERR);
if ( !stream.read(actualStepNumber) ) {
THROW_CIOERR(CIO_IOERR);
}
if ( !stream.read(actualStepIndx) ) {
THROW_CIOERR(CIO_IOERR);
}
for ( int i = 0; i <= nHistVectors; i++ ) {
if ( ( iores = solutionVectors[i].restoreYourself(stream) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
}
for ( int i = 0; i <= nHistVectors; i++ ) {
solStepList[i] = TimeStep(emodel);
if ( ( iores = solStepList[i].restoreContext(stream, mode) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
}
return CIO_OK;
}
contextIOResultType
PatchIntegrationRule :: saveContext(DataStream *stream, ContextMode mode, void *obj)
{
//
// saves full context (saves state variables, that completely describe
// current state)
//
// save parent data
contextIOResultType iores;
if ( ( iores = IntegrationRule :: saveContext(stream, mode, obj) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
// save patch data
if ( this->patch ) {
// store patch type
int _type = this->patch->givePatchType();
if ( !stream->write(& _type, 1) ) {
THROW_CIOERR(CIO_IOERR);
}
patch->saveContext(stream, mode, obj);
} else {
OOFEM_ERROR("saveContex : can't store NULL patch");
}
return CIO_OK;
}
contextIOResultType
PatchIntegrationRule :: restoreContext(DataStream *stream, ContextMode mode, void *obj)
{
//
// restores full element context (saves state variables, that completely describe
// current state)
//
contextIOResultType iores;
if ( stream == NULL ) {
OOFEM_ERROR("restoreContex : can't write into NULL stream");
}
if ( ( iores = IntegrationRule :: restoreContext(stream, mode, obj) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
// restore patch data
if ( this->patch ) {
delete this->patch;
}
int _ptype;
if ( !stream->read(& _ptype, 1) ) {
THROW_CIOERR(CIO_IOERR);
}
// create new patch
this->patch = CreateUsrDefPatch( ( Patch :: PatchType ) _ptype, this->giveElement() );
this->patch->restoreContext(stream, mode, obj);
return CIO_OK;
}
contextIOResultType
MisesMatStatus :: restoreContext(DataStream &stream, ContextMode mode, void *obj)
//
// restores full information stored in stream to this Status
//
{
contextIOResultType iores;
// read parent class status
if ( ( iores = StructuralMaterialStatus :: restoreContext(stream, mode, obj) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
// read plastic strain (vector)
if ( ( iores = plasticStrain.restoreYourself(stream) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
// read cumulative plastic strain (scalar)
if ( !stream.read(kappa) ) {
THROW_CIOERR(CIO_IOERR);
}
// read damage (scalar)
if ( !stream.read(damage) ) {
THROW_CIOERR(CIO_IOERR);
}
return CIO_OK; // return succes
}
contextIOResultType IncrementalLinearStatic :: saveContext(DataStream *stream, ContextMode mode, void *obj)
{
int closeFlag = 0;
contextIOResultType iores;
FILE *file = NULL;
if ( stream == NULL ) {
if ( !this->giveContextFile(& file, this->giveCurrentStep()->giveNumber(),
this->giveCurrentStep()->giveVersion(), contextMode_write) ) {
THROW_CIOERR(CIO_IOERR);
}
stream = new FileDataStream(file);
closeFlag = 1;
}
if ( ( iores = StructuralEngngModel :: saveContext(stream, mode) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
if ( closeFlag ) {
fclose(file);
delete stream;
stream = NULL;
} // ensure consistent records
return CIO_OK;
}
contextIOResultType
BinghamFluidMaterial2Status :: saveContext(DataStream *stream, ContextMode mode, void *obj)
//
// saves full ms context (saves state variables, that completely describe
// current state)
// saving the data in TDictionary is left to material (yield crit. level).
{
contextIOResultType iores;
if ( stream == NULL ) {
_error("saveContex : can't write into NULL stream");
}
if ( ( iores = FluidDynamicMaterialStatus :: saveContext(stream, mode, obj) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
if ( !stream->write(& devStrainMagnitude, 1) ) {
THROW_CIOERR(CIO_IOERR);
}
if ( !stream->write(& devStressMagnitude, 1) ) {
THROW_CIOERR(CIO_IOERR);
}
if ( ( iores = deviatoricStrainVector.storeYourself(stream, mode) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
return CIO_OK;
}
contextIOResultType IncrementalLinearStatic :: restoreContext(DataStream *stream, ContextMode mode, void *obj)
{
int closeFlag = 0, istep, iversion;
contextIOResultType iores;
FILE *file = NULL;
this->resolveCorrespondingStepNumber(istep, iversion, obj);
if ( stream == NULL ) {
if ( !this->giveContextFile(& file, istep, iversion, contextMode_read) ) {
THROW_CIOERR(CIO_IOERR);
}
stream = new FileDataStream(file);
closeFlag = 1;
}
if ( ( iores = StructuralEngngModel :: restoreContext(stream, mode, obj) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
if ( closeFlag ) {
fclose(file);
delete stream;
stream = NULL;
} // ensure consistent records
return CIO_OK;
}
contextIOResultType
BinghamFluidMaterial2Status :: restoreContext(DataStream *stream, ContextMode mode, void *obj)
//
// restores full material context (saves state variables, that completely describe
// current state)
//
{
// FloatArray *s;
contextIOResultType iores;
if ( stream == NULL ) {
_error("saveContex : can't write into NULL stream");
}
if ( ( iores = FluidDynamicMaterialStatus :: restoreContext(stream, mode, obj) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
if ( !stream->read(& devStrainMagnitude, 1) ) {
THROW_CIOERR(CIO_IOERR);
}
if ( !stream->read(& devStressMagnitude, 1) ) {
THROW_CIOERR(CIO_IOERR);
}
if ( ( iores = deviatoricStrainVector.restoreYourself(stream, mode) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
return CIO_OK;
}
contextIOResultType
RheoChainMaterialStatus :: saveContext(DataStream *stream, ContextMode mode, void *obj)
//
// saves full information stored in this Status
//
{
contextIOResultType iores;
if ( stream == NULL ) {
OOFEM_ERROR("can't write into NULL stream");
}
if ( ( iores = StructuralMaterialStatus :: saveContext(stream, mode, obj) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
// write raw data
for ( int i = 0; i < nUnits; i++ ) {
if ( ( iores = hiddenVars [ i ].storeYourself(stream, mode) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
}
if ( ( iores = shrinkageStrain.storeYourself(stream, mode) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
return CIO_OK;
}
contextIOResultType
TransientTransportProblem :: saveContext(DataStream *stream, ContextMode mode, void *obj)
{
contextIOResultType iores;
int closeFlag = 0;
FILE *file = NULL;
if ( stream == NULL ) {
if ( !this->giveContextFile(& file, this->giveCurrentStep()->giveNumber(),
this->giveCurrentStep()->giveVersion(), contextMode_write) ) {
THROW_CIOERR(CIO_IOERR); // override
}
stream = new FileDataStream(file);
closeFlag = 1;
}
if ( ( iores = EngngModel :: saveContext(stream, mode) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
if ( ( iores = field->saveContext(*stream, mode) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
if ( closeFlag ) {
fclose(file);
delete stream;
stream = NULL;
}
// ensure consistent records
return CIO_OK;
}
contextIOResultType XfemManager :: saveContext(DataStream *stream, ContextMode mode, void *obj)
{
contextIOResultType iores;
if ( mode & CM_Definition ) {
if ( !stream->write(& this->numberOfEnrichmentItems, 1) ) {
THROW_CIOERR(CIO_IOERR);
}
}
for ( int i = 1; i <= this->numberOfEnrichmentItems; i++ ) {
EnrichmentItem *obj = this->giveEnrichmentItem(i);
if ( ( mode & CM_Definition ) ) {
if ( !stream->write( obj->giveInputRecordName() ) ) {
THROW_CIOERR(CIO_IOERR);
}
}
if ( ( iores = obj->saveContext(stream, mode) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
}
return CIO_OK;
}
contextIOResultType
MicroplaneMaterial :: saveIPContext(DataStream *stream, ContextMode mode, GaussPoint* gp)
{
contextIOResultType iores;
int mPlaneIndex;
Microplane *mPlane;
IntegrationPointStatus *status;
if ( stream == NULL ) {
_error("saveContex : can't write into NULL stream");
}
if ( gp == NULL ) {
THROW_CIOERR(CIO_BADOBJ);
}
// save master
if ( ( iores = StructuralMaterial :: saveIPContext(stream, mode, gp) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
// save microplanes
for ( mPlaneIndex = 0; mPlaneIndex < numberOfMicroplanes; mPlaneIndex++ ) {
mPlane = this->giveMicroplane(mPlaneIndex, gp);
status = this->giveMicroplaneStatus(mPlane);
if ( status ) {
if ( ( iores = status->saveContext(stream, mode, gp) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
}
}
return CIO_OK;
}
contextIOResultType
Material :: restoreIPContext(DataStream *stream, ContextMode mode, GaussPoint *gp)
//
// restores full material status (saves state variables, that completely describe
// current state) stored in gp->matstatusDict with key = (int)this->giveClassID()
// restoring of corresponding context if it is defined for current material in
// gp status dictionary should be performed here by overloading this function.
// (such code should invoke also corresponding function for yield conditions,
// submaterials and so on)
//
//
{
contextIOResultType iores;
if ( gp == NULL ) {
THROW_CIOERR(CIO_BADOBJ);
}
// read raw data - context
MaterialStatus *status = this->giveStatus(gp);
if ( status ) {
if ( ( iores = status->restoreContext(stream, mode, gp) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
}
return CIO_OK;
}
contextIOResultType
PlasticMaterialStatus :: saveContext(DataStream &stream, ContextMode mode, void *obj)
//
// saves full information stored in this Status
// no temp variables stored
//
{
contextIOResultType iores;
// save parent class status
if ( ( iores = StructuralMaterialStatus :: saveContext(stream, mode, obj) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
// write a raw data
if ( ( iores = plasticStrainVector.storeYourself(stream) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
if ( ( iores = strainSpaceHardeningVarsVector.storeYourself(stream) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
if ( !stream.write(state_flag) ) {
THROW_CIOERR(CIO_IOERR);
}
if ( !stream.write(gamma) ) {
THROW_CIOERR(CIO_IOERR);
}
return CIO_OK;
}
contextIOResultType
MicroplaneMaterial :: restoreIPContext(DataStream &stream, ContextMode mode, GaussPoint *gp)
{
contextIOResultType iores;
int mPlaneIndex;
Microplane *mPlane;
IntegrationPointStatus *status;
// corresponding gp is passed in obj
if ( gp == NULL ) {
THROW_CIOERR(CIO_BADOBJ);
}
// save master
if ( ( iores = StructuralMaterial :: restoreIPContext(stream, mode, gp) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
// save microplanes
for ( mPlaneIndex = 0; mPlaneIndex < numberOfMicroplanes; mPlaneIndex++ ) {
mPlane = this->giveMicroplane(mPlaneIndex, gp);
status = this->giveMicroplaneStatus(mPlane);
if ( status ) {
if ( ( iores = status->restoreContext(stream, mode, gp) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
}
}
return CIO_OK;
}
contextIOResultType
PlasticMaterialStatus :: restoreContext(DataStream &stream, ContextMode mode, void *obj)
//
// restores full information stored in stream to this Status
//
{
contextIOResultType iores;
// read parent class status
if ( ( iores = StructuralMaterialStatus :: restoreContext(stream, mode, obj) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
if ( ( iores = plasticStrainVector.restoreYourself(stream) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
if ( ( iores = strainSpaceHardeningVarsVector.restoreYourself(stream) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
if ( !stream.read(state_flag) ) {
THROW_CIOERR(CIO_IOERR);
}
if ( !stream.read(gamma) ) {
THROW_CIOERR(CIO_IOERR);
}
return CIO_OK; // return success
}
contextIOResultType FreeWarping :: saveContext(DataStream *stream, ContextMode mode, void *obj)
//
// saves state variable - displacement vector
//
{
contextIOResultType iores;
int closeFlag = 0;
FILE *file = NULL;
if ( stream == NULL ) {
if ( !this->giveContextFile(& file, this->giveCurrentStep()->giveNumber(),
this->giveCurrentStep()->giveVersion(), contextMode_write) ) {
THROW_CIOERR(CIO_IOERR); // override
}
stream = new FileDataStream(file);
closeFlag = 1;
}
if ( ( iores = StructuralEngngModel :: saveContext(stream, mode) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
if ( ( iores = displacementVector.storeYourself(*stream) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
if ( closeFlag ) {
fclose(file);
delete stream;
stream = NULL;
}
return CIO_OK;
}
contextIOResultType
StructuralMaterialStatus :: restoreContext(DataStream *stream, ContextMode mode, void *obj)
//
// restores full material context (saves state variables, that completely describe
// current state)
//
{
contextIOResultType iores;
if ( stream == NULL ) {
OOFEM_ERROR("can't write into NULL stream");
}
if ( ( iores = MaterialStatus :: restoreContext(stream, mode, obj) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
if ( ( iores = strainVector.restoreYourself(stream, mode) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
if ( ( iores = stressVector.restoreYourself(stream, mode) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
return CIO_OK;
}
contextIOResultType
RCSDMaterialStatus :: restoreContext(DataStream *stream, ContextMode mode, void *obj)
//
// restores full information stored in stream to this Status
//
{
contextIOResultType iores;
// read parent class status
if ( ( iores = RCM2MaterialStatus :: restoreContext(stream, mode, obj) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
// read raw data
if ( !stream->read(& maxEquivStrain, 1) ) {
THROW_CIOERR(CIO_IOERR);
}
if ( !stream->read(& damageCoeff, 1) ) {
THROW_CIOERR(CIO_IOERR);
}
if ( !stream->read(& mode, 1) ) {
THROW_CIOERR(CIO_IOERR);
}
if ( ( iores = Ds0.restoreYourself(stream, mode) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
return CIO_OK; // return succes
}
contextIOResultType LinearStability :: restoreContext(DataStream *stream, ContextMode mode, void *obj)
//
// restore state variable - displacement vector
//
{
int activeVector, version;
int istep = 1, iversion = 1;
int closeFlag = 0;
contextIOResultType iores;
FILE *file = NULL;
this->resolveCorrespondingStepNumber(activeVector, version, obj);
if ( eigVal.isEmpty() ) { // not restored before
if ( stream == NULL ) {
if ( !this->giveContextFile(& file, istep, iversion, contextMode_read) ) {
THROW_CIOERR(CIO_IOERR); // override
}
stream = new FileDataStream(file);
closeFlag = 1;
}
if ( ( iores = StructuralEngngModel :: restoreContext(stream, mode, ( void * ) & istep) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
if ( ( iores = displacementVector.restoreYourself(stream, mode) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
if ( ( iores = eigVal.restoreYourself(stream, mode) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
if ( ( iores = eigVec.restoreYourself(stream, mode) ) != CIO_OK ) {
THROW_CIOERR(iores);
}
if ( closeFlag ) {
fclose(file);
delete stream;
stream = NULL;
} // ensure consistent records
}
// if (istep > numberOfRequiredEigenValues) istep = numberOfRequiredEigenValues ;
// printf( "Restoring - corresponding index is %d, EigenValue is %lf\n",
// istep,eigVal.at(istep));
// setActiveVector (istep);
if ( activeVector > numberOfRequiredEigenValues ) {
activeVector = numberOfRequiredEigenValues;
}
OOFEM_LOG_INFO( "Restoring - corresponding index is %d, EigenValue is %f\n",
activeVector, eigVal.at(activeVector) );
this->giveCurrentStep()->setTime( ( double ) activeVector );
return CIO_OK;
}
