void
POIExportModule :: exportIntVarAs(InternalStateType valID, FILE *stream, TimeStep *tStep)
{
int i, region;
IntArray toMap(1);
Domain *d = emodel->giveDomain(1);
FloatArray poiCoords(3);
FloatArray val;
toMap.at(1) = ( int ) valID;
// loop over POIs
for ( auto &poi: POIList ) {
poiCoords.at(1) = poi.x;
poiCoords.at(2) = poi.y;
poiCoords.at(3) = poi.z;
region = poi.region;
this->giveMapper()->__init(d, toMap, poiCoords, * d->giveSet(region), tStep);
if ( !this->giveMapper()->__mapVariable(val, poiCoords, valID, tStep) ) {
OOFEM_WARNING("Failed to map variable");
val.clear();
}
fprintf(stream, "%10d ", poi.id);
for ( i = 1; i <= val.giveSize(); i++ ) {
fprintf( stream, " %15e", val.at(i) );
}
fprintf(stream, "\n");
}
}
void
POIExportModule :: exportIntVarAs(InternalStateType valID, FILE *stream, TimeStep *tStep)
{
int i, region;
IntArray toMap(1);
Domain *d = emodel->giveDomain(1);
FloatArray poiCoords(3);
FloatArray val;
toMap.at(1) = ( int ) valID;
// loop over POIs
std::list< POI_dataType > :: iterator PoiIter;
for ( PoiIter = POIList.begin(); PoiIter != POIList.end(); ++PoiIter ) {
poiCoords.at(1) = ( * PoiIter ).x;
poiCoords.at(2) = ( * PoiIter ).y;
poiCoords.at(3) = ( * PoiIter ).z;
region = ( * PoiIter ).region;
this->giveMapper()->__init(d, toMap, poiCoords, region, tStep);
if ( !this->giveMapper()->__mapVariable(val, poiCoords, valID, tStep) ) {
OOFEM_WARNING("POIExportModule :: exportIntVarAs - Failed to map variable");
val.resize(0);
}
fprintf(stream, "%10d ", ( * PoiIter ).id);
for ( i = 1; i <= val.giveSize(); i++ ) {
fprintf( stream, " %15e", val.at(i) );
}
fprintf(stream, "\n");
}
}
void
POIExportModule :: exportIntVars(FILE *stream, TimeStep *tStep)
{
int i, n = internalVarsToExport.giveSize();
InternalStateType type;
FloatArray poiCoords(3);
if ( n == 0 ) {
return;
}
// loop over POIs
POI_dataType &poi = *POIList.begin();
poiCoords.at(1) = poi.x;
poiCoords.at(2) = poi.y;
poiCoords.at(3) = poi.z;
//int region = poi.region;
for ( i = 1; i <= n; i++ ) {
type = ( InternalStateType ) internalVarsToExport.at(i);
fprintf(stream, "\n\nPOI_INTVAR_DATA %d\n", type);
this->exportIntVarAs(type, stream, tStep);
}
this->giveMapper()->finish(tStep);
}
void
POIExportModule :: exportIntVars(FILE *stream, TimeStep *tStep)
{
int i, region, n = internalVarsToExport.giveSize();
Domain *d = emodel->giveDomain(1);
InternalStateType type;
FloatArray poiCoords(3);
if ( n == 0 ) {
return;
}
// loop over POIs
std::list< POI_dataType > :: iterator PoiIter = POIList.begin();
poiCoords.at(1) = ( * PoiIter ).x;
poiCoords.at(2) = ( * PoiIter ).y;
poiCoords.at(3) = ( * PoiIter ).z;
region = ( * PoiIter ).region;
this->giveMapper()->__init(d, internalVarsToExport, poiCoords, region, tStep);
for ( i = 1; i <= n; i++ ) {
type = ( InternalStateType ) internalVarsToExport.at(i);
fprintf(stream, "\n\nPOI_INTVAR_DATA %d\n", type);
this->exportIntVarAs(type, stream, tStep);
}
this->giveMapper()->finish(tStep);
}
