void L4Axisymm :: drawScalar(oofegGraphicContext &gc, TimeStep *tStep) { int i, indx, result = 0; WCRec p [ 4 ]; GraphicObj *tr; FloatArray v [ 4 ]; double s [ 4 ], defScale; if ( !gc.testElementGraphicActivity(this) ) { return; } EASValsSetLayer(OOFEG_VARPLOT_PATTERN_LAYER); if ( gc.giveIntVarMode() == ISM_recovered ) { for ( i = 1; i <= 4; i++ ) { result += this->giveInternalStateAtNode(v [ i - 1 ], gc.giveIntVarType(), gc.giveIntVarMode(), i, tStep); } if ( result != 4 ) { return; } indx = gc.giveIntVarIndx(); for ( i = 1; i <= 4; i++ ) { s [ i - 1 ] = v [ i - 1 ].at(indx); } if ( gc.getScalarAlgo() == SA_ISO_SURF ) { for ( i = 0; i < 4; i++ ) { if ( gc.getInternalVarsDefGeoFlag() ) { // use deformed geometry defScale = gc.getDefScale(); p [ i ].x = ( FPNum ) this->giveNode(i + 1)->giveUpdatedCoordinate(1, tStep, defScale); p [ i ].y = ( FPNum ) this->giveNode(i + 1)->giveUpdatedCoordinate(2, tStep, defScale); p [ i ].z = 0.; } else { p [ i ].x = ( FPNum ) this->giveNode(i + 1)->giveCoordinate(1); p [ i ].y = ( FPNum ) this->giveNode(i + 1)->giveCoordinate(2); p [ i ].z = 0.; } } //EASValsSetColor(gc.getYieldPlotColor(ratio)); gc.updateFringeTableMinMax(s, 4); tr = CreateQuadWD3D(p, s [ 0 ], s [ 1 ], s [ 2 ], s [ 3 ]); EGWithMaskChangeAttributes(LAYER_MASK, tr); EMAddGraphicsToModel(ESIModel(), tr); /* * } else if (gc.getScalarAlgo() == SA_ISO_LINE) { * * EASValsSetColor(context.getActiveCrackColor()); * EASValsSetLineWidth(OOFEG_ISO_LINE_WIDTH); * * for (i=0; i< 4; i++) { * if (gc.getInternalVarsDefGeoFlag()) { * // use deformed geometry * defScale = gc.getDefScale(); * p[i].x = (FPNum) this->giveNode(i+1)->giveUpdatedCoordinate(1,tStep,defScale); * p[i].y = (FPNum) this->giveNode(i+1)->giveUpdatedCoordinate(2,tStep,defScale); * p[i].z = 0.; * * } else { * p[i].x = (FPNum) this->giveNode(i+1)->giveCoordinate(1); * p[i].y = (FPNum) this->giveNode(i+1)->giveCoordinate(2); * p[i].z = 0.; * } * } * * // isoline implementation * oofeg_drawIsoLinesOnQuad (p, s); * */ } } else if ( gc.giveIntVarMode() == ISM_local ) { if ( numberOfGaussPoints != 4 ) { return; } IntArray ind(4); FloatArray *gpCoords; WCRec pp [ 9 ]; for ( i = 0; i < 4; i++ ) { if ( gc.getInternalVarsDefGeoFlag() ) { // use deformed geometry defScale = gc.getDefScale(); pp [ i ].x = ( FPNum ) this->giveNode(i + 1)->giveUpdatedCoordinate(1, tStep, defScale); pp [ i ].y = ( FPNum ) this->giveNode(i + 1)->giveUpdatedCoordinate(2, tStep, defScale); pp [ i ].z = 0.; } else { pp [ i ].x = ( FPNum ) this->giveNode(i + 1)->giveCoordinate(1); pp [ i ].y = ( FPNum ) this->giveNode(i + 1)->giveCoordinate(2); pp [ i ].z = 0.; } } for ( i = 0; i < 3; i++ ) { pp [ i + 4 ].x = 0.5 * ( pp [ i ].x + pp [ i + 1 ].x ); pp [ i + 4 ].y = 0.5 * ( pp [ i ].y + pp [ i + 1 ].y ); pp [ i + 4 ].z = 0.5 * ( pp [ i ].z + pp [ i + 1 ].z ); } pp [ 7 ].x = 0.5 * ( pp [ 3 ].x + pp [ 0 ].x ); pp [ 7 ].y = 0.5 * ( pp [ 3 ].y + pp [ 0 ].y ); pp [ 7 ].z = 0.5 * ( pp [ 3 ].z + pp [ 0 ].z ); pp [ 8 ].x = 0.25 * ( pp [ 0 ].x + pp [ 1 ].x + pp [ 2 ].x + pp [ 3 ].x ); pp [ 8 ].y = 0.25 * ( pp [ 0 ].y + pp [ 1 ].y + pp [ 2 ].y + pp [ 3 ].y ); pp [ 8 ].z = 0.25 * ( pp [ 0 ].z + pp [ 1 ].z + pp [ 2 ].z + pp [ 3 ].z ); for ( GaussPoint *gp: *this->giveDefaultIntegrationRulePtr() ) { gpCoords = gp->giveNaturalCoordinates(); if ( ( gpCoords->at(1) > 0. ) && ( gpCoords->at(2) > 0. ) ) { ind.at(1) = 0; ind.at(2) = 4; ind.at(3) = 8; ind.at(4) = 7; } else if ( ( gpCoords->at(1) < 0. ) && ( gpCoords->at(2) > 0. ) ) { ind.at(1) = 4; ind.at(2) = 1; ind.at(3) = 5; ind.at(4) = 8; } else if ( ( gpCoords->at(1) < 0. ) && ( gpCoords->at(2) < 0. ) ) { ind.at(1) = 5; ind.at(2) = 2; ind.at(3) = 6; ind.at(4) = 8; } else { ind.at(1) = 6; ind.at(2) = 3; ind.at(3) = 7; ind.at(4) = 8; } if ( giveIPValue(v [ 0 ], gp, gc.giveIntVarType(), tStep) == 0 ) { return; } indx = gc.giveIntVarIndx(); for ( i = 1; i <= 4; i++ ) { s [ i - 1 ] = v [ 0 ].at(indx); } for ( i = 0; i < 4; i++ ) { p [ i ].x = pp [ ind.at(i + 1) ].x; p [ i ].y = pp [ ind.at(i + 1) ].y; p [ i ].z = pp [ ind.at(i + 1) ].z; } gc.updateFringeTableMinMax(s, 4); tr = CreateQuadWD3D(p, s [ 0 ], s [ 1 ], s [ 2 ], s [ 3 ]); EGWithMaskChangeAttributes(LAYER_MASK, tr); EMAddGraphicsToModel(ESIModel(), tr); } } }
void QPlaneStress2d :: drawScalar(oofegGraphicContext &context) { int i, indx, n [ 4 ], result = 0; WCRec p [ 4 ], pp [ 9 ]; GraphicObj *tr; TimeStep *tStep = this->giveDomain()->giveEngngModel()->giveCurrentStep(); FloatArray v [ 8 ]; double s [ 9 ], ss [ 4 ], defScale; int ip; GaussPoint *gp; if ( !context.testElementGraphicActivity(this) ) { return; } EASValsSetLayer(OOFEG_VARPLOT_PATTERN_LAYER); if ( context.giveIntVarMode() == ISM_recovered ) { // ============ plot the recovered values (smoothed data) =============== for ( i = 1; i <= 8; i++ ) { result += this->giveInternalStateAtNode(v [ i - 1 ], context.giveIntVarType(), context.giveIntVarMode(), i, tStep); } if ( result != 8 ) { return; } indx = context.giveIntVarIndx(); for ( i = 1; i <= 8; i++ ) { s [ i - 1 ] = v [ i - 1 ].at(indx); } // auxiliary value at an added central node // computed as average of the values at all Gauss points s [ 8 ] = 0.; for ( ip = 1; ip <= integrationRulesArray [ 0 ]->giveNumberOfIntegrationPoints(); ip++ ) { gp = integrationRulesArray [ 0 ]->getIntegrationPoint(ip - 1); if ( giveIPValue(v [ 0 ], gp, context.giveIntVarType(), tStep) == 0 ) { return; } s [ 8 ] += v [ 0 ].at(indx); } s [ 8 ] /= integrationRulesArray [ 0 ]->giveNumberOfIntegrationPoints(); //s[8] = (s[4]+s[5]+s[6]+s[7])/4.; for ( i = 0; i < 8; i++ ) { if ( context.getInternalVarsDefGeoFlag() ) { // use deformed geometry defScale = context.getDefScale(); pp [ i ].x = ( FPNum ) this->giveNode(i + 1)->giveUpdatedCoordinate(1, tStep, defScale); pp [ i ].y = ( FPNum ) this->giveNode(i + 1)->giveUpdatedCoordinate(2, tStep, defScale); pp [ i ].z = 0.; } else { // use initial geometry pp [ i ].x = ( FPNum ) this->giveNode(i + 1)->giveCoordinate(1); pp [ i ].y = ( FPNum ) this->giveNode(i + 1)->giveCoordinate(2); pp [ i ].z = 0.; } } pp [ 8 ].x = ( pp [ 4 ].x + pp [ 5 ].x + pp [ 6 ].x + pp [ 7 ].x ) / 4.; pp [ 8 ].y = ( pp [ 4 ].y + pp [ 5 ].y + pp [ 6 ].y + pp [ 7 ].y ) / 4.; pp [ 8 ].z = 0.; for ( int t = 1; t <= 4; t++ ) { if ( t == 1 ) { n [ 0 ] = 0; n [ 1 ] = 4; n [ 2 ] = 8; n [ 3 ] = 7; } else if ( t == 2 ) { n [ 0 ] = 4; n [ 1 ] = 1; n [ 2 ] = 5; n [ 3 ] = 8; } else if ( t == 3 ) { n [ 0 ] = 5; n [ 1 ] = 2; n [ 2 ] = 6; n [ 3 ] = 8; } else { n [ 0 ] = 6; n [ 1 ] = 3; n [ 2 ] = 7; n [ 3 ] = 8; } ss [ 0 ] = s [ n [ 0 ] ]; ss [ 1 ] = s [ n [ 1 ] ]; ss [ 2 ] = s [ n [ 2 ] ]; ss [ 3 ] = s [ n [ 3 ] ]; for ( i = 0; i < 4; i++ ) { p [ i ].x = pp [ n [ i ] ].x; p [ i ].y = pp [ n [ i ] ].y; p [ i ].z = 0.; } if ( context.getScalarAlgo() == SA_ISO_SURF ) { /* * for ( i = 0; i < 4; i++ ) { * if ( context.getInternalVarsDefGeoFlag() ) { * // use deformed geometry * defScale = context.getDefScale(); * p [ i ].x = ( FPNum ) this->giveNode(n[i] + 1)->giveUpdatedCoordinate(1, tStep, defScale); * p [ i ].y = ( FPNum ) this->giveNode(n[i] + 1)->giveUpdatedCoordinate(2, tStep, defScale); * p [ i ].z = 0.; * } else { * // use initial geometry * p [ i ].x = ( FPNum ) this->giveNode(n[i] + 1)->giveCoordinate(1); * p [ i ].y = ( FPNum ) this->giveNode(n[i] + 1)->giveCoordinate(2); * p [ i ].z = 0.; * } * } */ //EASValsSetColor(gc.getYieldPlotColor(ratio)); context.updateFringeTableMinMax(ss, 4); tr = CreateQuadWD3D(p, ss [ 0 ], ss [ 1 ], ss [ 2 ], ss [ 3 ]); EGWithMaskChangeAttributes(LAYER_MASK, tr); EMAddGraphicsToModel(ESIModel(), tr); } else if ( ( context.getScalarAlgo() == SA_ZPROFILE ) || ( context.getScalarAlgo() == SA_COLORZPROFILE ) ) { //double landScale = context.getLandScale(); for ( i = 0; i < 4; i++ ) { /* * if ( context.getInternalVarsDefGeoFlag() ) { * // use deformed geometry * defScale = context.getDefScale(); * p [ i ].x = ( FPNum ) this->giveNode(i + 1)->giveUpdatedCoordinate(1, tStep, defScale); * p [ i ].y = ( FPNum ) this->giveNode(i + 1)->giveUpdatedCoordinate(2, tStep, defScale); * p [ i ].z = ss [ i ] * landScale; * } else { * // use initial geometry * p [ i ].x = ( FPNum ) this->giveNode(i + 1)->giveCoordinate(1); * p [ i ].y = ( FPNum ) this->giveNode(i + 1)->giveCoordinate(2); * p [ i ].z = ss [ i ] * landScale; * } */ // this fixes a bug in ELIXIR if ( fabs(ss [ i ]) < 1.0e-6 ) { ss [ i ] = 1.0e-6; } } if ( context.getScalarAlgo() == SA_ZPROFILE ) { EASValsSetColor( context.getDeformedElementColor() ); EASValsSetLineWidth(OOFEG_DEFORMED_GEOMETRY_WIDTH); tr = CreateQuad3D(p); EGWithMaskChangeAttributes(WIDTH_MASK | COLOR_MASK | LAYER_MASK, tr); } else { context.updateFringeTableMinMax(s, 4); tr = CreateQuadWD3D(p, ss [ 0 ], ss [ 1 ], ss [ 2 ], ss [ 3 ]); EGWithMaskChangeAttributes(LAYER_MASK, tr); } EMAddGraphicsToModel(ESIModel(), tr); } } } else if ( context.giveIntVarMode() == ISM_local ) { // ========== plot the local values (raw data) ===================== if ( numberOfGaussPoints != 4 ) { return; } IntArray ind(4); FloatArray *gpCoords; WCRec pp [ 9 ]; for ( i = 0; i < 8; i++ ) { if ( context.getInternalVarsDefGeoFlag() ) { // use deformed geometry defScale = context.getDefScale(); pp [ i ].x = ( FPNum ) this->giveNode(i + 1)->giveUpdatedCoordinate(1, tStep, defScale); pp [ i ].y = ( FPNum ) this->giveNode(i + 1)->giveUpdatedCoordinate(2, tStep, defScale); pp [ i ].z = 0.; } else { pp [ i ].x = ( FPNum ) this->giveNode(i + 1)->giveCoordinate(1); pp [ i ].y = ( FPNum ) this->giveNode(i + 1)->giveCoordinate(2); pp [ i ].z = 0.; } } pp [ 8 ].x = 0.25 * ( pp [ 0 ].x + pp [ 1 ].x + pp [ 2 ].x + pp [ 3 ].x ); pp [ 8 ].y = 0.25 * ( pp [ 0 ].y + pp [ 1 ].y + pp [ 2 ].y + pp [ 3 ].y ); pp [ 8 ].z = 0.; for ( ip = 1; ip <= integrationRulesArray [ 0 ]->giveNumberOfIntegrationPoints(); ip++ ) { gp = integrationRulesArray [ 0 ]->getIntegrationPoint(ip - 1); gpCoords = gp->giveCoordinates(); if ( ( gpCoords->at(1) > 0. ) && ( gpCoords->at(2) > 0. ) ) { ind.at(1) = 0; ind.at(2) = 4; ind.at(3) = 8; ind.at(4) = 7; } else if ( ( gpCoords->at(1) < 0. ) && ( gpCoords->at(2) > 0. ) ) { ind.at(1) = 4; ind.at(2) = 1; ind.at(3) = 5; ind.at(4) = 8; } else if ( ( gpCoords->at(1) < 0. ) && ( gpCoords->at(2) < 0. ) ) { ind.at(1) = 5; ind.at(2) = 2; ind.at(3) = 6; ind.at(4) = 8; } else { ind.at(1) = 6; ind.at(2) = 3; ind.at(3) = 7; ind.at(4) = 8; } if ( giveIPValue(v [ 0 ], gp, context.giveIntVarType(), tStep) == 0 ) { return; } indx = context.giveIntVarIndx(); for ( i = 1; i <= 4; i++ ) { s [ i - 1 ] = v [ 0 ].at(indx); } for ( i = 0; i < 4; i++ ) { p [ i ].x = pp [ ind.at(i + 1) ].x; p [ i ].y = pp [ ind.at(i + 1) ].y; p [ i ].z = pp [ ind.at(i + 1) ].z; } context.updateFringeTableMinMax(s, 4); tr = CreateQuadWD3D(p, s [ 0 ], s [ 1 ], s [ 2 ], s [ 3 ]); EGWithMaskChangeAttributes(LAYER_MASK, tr); EMAddGraphicsToModel(ESIModel(), tr); } } }
void QPlaneStrain :: drawScalar(oofegGraphicContext &context) { int i, indx; WCRec p [ 4 ]; GraphicObj *tr; TimeStep *tStep = this->giveDomain()->giveEngngModel()->giveCurrentStep(); FloatArray v [ 4 ]; double s [ 4 ], defScale; if ( !context.testElementGraphicActivity(this) ) { return; } EASValsSetLayer(OOFEG_VARPLOT_PATTERN_LAYER); if ( context.giveIntVarMode() == ISM_recovered ) { // ============ plot the recovered values (smoothed data) =============== /* * for ( i = 1; i <= 4; i++ ) { * result += this->giveInternalStateAtNode(v [ i - 1 ], context.giveIntVarType(), context.giveIntVarMode(), i, tStep); * } * * if ( result != 4 ) { * return; * } * * indx = context.giveIntVarIndx(); * * for ( i = 1; i <= 4; i++ ) { * s [ i - 1 ] = v [ i - 1 ].at(indx); * } * * if ( context.getScalarAlgo() == SA_ISO_SURF ) { * for ( i = 0; i < 4; i++ ) { * if ( context.getInternalVarsDefGeoFlag() ) { * // use deformed geometry * defScale = context.getDefScale(); * p [ i ].x = ( FPNum ) this->giveNode(i + 1)->giveUpdatedCoordinate(1, tStep, defScale); * p [ i ].y = ( FPNum ) this->giveNode(i + 1)->giveUpdatedCoordinate(2, tStep, defScale); * p [ i ].z = 0.; * } else { * p [ i ].x = ( FPNum ) this->giveNode(i + 1)->giveCoordinate(1); * p [ i ].y = ( FPNum ) this->giveNode(i + 1)->giveCoordinate(2); * p [ i ].z = 0.; * } * } * * //EASValsSetColor(gc.getYieldPlotColor(ratio)); * context.updateFringeTableMinMax(s, 4); * tr = CreateQuadWD3D(p, s [ 0 ], s [ 1 ], s [ 2 ], s [ 3 ]); * EGWithMaskChangeAttributes(LAYER_MASK, tr); * EMAddGraphicsToModel(ESIModel(), tr); * } else if ( ( context.getScalarAlgo() == SA_ZPROFILE ) || ( context.getScalarAlgo() == SA_COLORZPROFILE ) ) { * double landScale = context.getLandScale(); * * for ( i = 0; i < 4; i++ ) { * if ( context.getInternalVarsDefGeoFlag() ) { * // use deformed geometry * defScale = context.getDefScale(); * p [ i ].x = ( FPNum ) this->giveNode(i + 1)->giveUpdatedCoordinate(1, tStep, defScale); * p [ i ].y = ( FPNum ) this->giveNode(i + 1)->giveUpdatedCoordinate(2, tStep, defScale); * p [ i ].z = s [ i ] * landScale; * } else { * p [ i ].x = ( FPNum ) this->giveNode(i + 1)->giveCoordinate(1); * p [ i ].y = ( FPNum ) this->giveNode(i + 1)->giveCoordinate(2); * p [ i ].z = s [ i ] * landScale; * } * * // this fixes a bug in ELIXIR * if ( fabs(s [ i ]) < 1.0e-6 ) { * s [ i ] = 1.0e-6; * } * } * * if ( context.getScalarAlgo() == SA_ZPROFILE ) { * EASValsSetColor( context.getDeformedElementColor() ); * EASValsSetLineWidth(OOFEG_DEFORMED_GEOMETRY_WIDTH); * tr = CreateQuad3D(p); * EGWithMaskChangeAttributes(WIDTH_MASK | COLOR_MASK | LAYER_MASK, tr); * } else { * context.updateFringeTableMinMax(s, 4); * tr = CreateQuadWD3D(p, s [ 0 ], s [ 1 ], s [ 2 ], s [ 3 ]); * EGWithMaskChangeAttributes(LAYER_MASK, tr); * } * * EMAddGraphicsToModel(ESIModel(), tr); * } */ } else if ( context.giveIntVarMode() == ISM_local ) { // ========== plot the local values (raw data) ===================== if ( numberOfGaussPoints != 4 ) { return; } int ip; GaussPoint *gp; IntArray ind(4); FloatArray *gpCoords; WCRec pp [ 9 ]; for ( i = 0; i < 8; i++ ) { if ( context.getInternalVarsDefGeoFlag() ) { // use deformed geometry defScale = context.getDefScale(); pp [ i ].x = ( FPNum ) this->giveNode(i + 1)->giveUpdatedCoordinate(1, tStep, defScale); pp [ i ].y = ( FPNum ) this->giveNode(i + 1)->giveUpdatedCoordinate(2, tStep, defScale); pp [ i ].z = 0.; } else { pp [ i ].x = ( FPNum ) this->giveNode(i + 1)->giveCoordinate(1); pp [ i ].y = ( FPNum ) this->giveNode(i + 1)->giveCoordinate(2); pp [ i ].z = 0.; } } pp [ 8 ].x = 0.25 * ( pp [ 0 ].x + pp [ 1 ].x + pp [ 2 ].x + pp [ 3 ].x ); pp [ 8 ].y = 0.25 * ( pp [ 0 ].y + pp [ 1 ].y + pp [ 2 ].y + pp [ 3 ].y ); pp [ 8 ].z = 0.; for ( ip = 1; ip <= integrationRulesArray [ 0 ]->giveNumberOfIntegrationPoints(); ip++ ) { gp = integrationRulesArray [ 0 ]->getIntegrationPoint(ip - 1); gpCoords = gp->giveCoordinates(); if ( ( gpCoords->at(1) > 0. ) && ( gpCoords->at(2) > 0. ) ) { ind.at(1) = 0; ind.at(2) = 4; ind.at(3) = 8; ind.at(4) = 7; } else if ( ( gpCoords->at(1) < 0. ) && ( gpCoords->at(2) > 0. ) ) { ind.at(1) = 4; ind.at(2) = 1; ind.at(3) = 5; ind.at(4) = 8; } else if ( ( gpCoords->at(1) < 0. ) && ( gpCoords->at(2) < 0. ) ) { ind.at(1) = 5; ind.at(2) = 2; ind.at(3) = 6; ind.at(4) = 8; } else { ind.at(1) = 6; ind.at(2) = 3; ind.at(3) = 7; ind.at(4) = 8; } if ( giveIPValue(v [ 0 ], gp, context.giveIntVarType(), tStep) == 0 ) { return; } indx = context.giveIntVarIndx(); for ( i = 1; i <= 4; i++ ) { s [ i - 1 ] = v [ 0 ].at(indx); } for ( i = 0; i < 4; i++ ) { p [ i ].x = pp [ ind.at(i + 1) ].x; p [ i ].y = pp [ ind.at(i + 1) ].y; p [ i ].z = pp [ ind.at(i + 1) ].z; } context.updateFringeTableMinMax(s, 4); tr = CreateQuadWD3D(p, s [ 0 ], s [ 1 ], s [ 2 ], s [ 3 ]); EGWithMaskChangeAttributes(LAYER_MASK, tr); EMAddGraphicsToModel(ESIModel(), tr); } } }
void Quad1_ht :: drawScalar(oofegGraphicContext &context) { int i, indx, result = 0; WCRec p [ 4 ]; GraphicObj *tr; TimeStep *tStep = this->giveDomain()->giveEngngModel()->giveCurrentStep(); double s [ 4 ]; FloatArray v [ 4 ]; InternalStateType itype = context.giveIntVarType(); if ( !context.testElementGraphicActivity(this) ) { return; } EASValsSetLayer(OOFEG_VARPLOT_PATTERN_LAYER); if ( itype == IST_HydrationDegree ) { for ( i = 1; i <= 4; i++ ) { result += this->giveInternalStateAtNode(v [ i - 1 ], context.giveIntVarType(), context.giveIntVarMode(), i, tStep); } if ( result != 4 ) { return; } indx = context.giveIntVarIndx(); for ( i = 1; i <= 4; i++ ) { s [ i - 1 ] = v [ i - 1 ].at(indx); } if ( context.getScalarAlgo() == SA_ISO_SURF ) { for ( i = 0; i < 4; i++ ) { p [ i ].x = ( FPNum ) this->giveNode(i + 1)->giveCoordinate(1); p [ i ].y = ( FPNum ) this->giveNode(i + 1)->giveCoordinate(2); p [ i ].z = 0.; } context.updateFringeTableMinMax(s, 4); tr = CreateQuadWD3D(p, s [ 0 ], s [ 1 ], s [ 2 ], s [ 3 ]); EGWithMaskChangeAttributes(LAYER_MASK, tr); EMAddGraphicsToModel(ESIModel(), tr); } } else if ( ( ( ( emode == HeatTransferEM ) || ( emode == HeatMass1TransferEM ) ) && ( itype == IST_Temperature ) ) || ( ( emode == HeatMass1TransferEM ) && ( itype == IST_MassConcentration_1 ) ) ) { IntArray dofMask(1); if ( itype == IST_Temperature ) { dofMask.at(1) = T_f; } else { dofMask.at(1) = C_1; } FloatArray r; for ( i = 0; i < 4; i++ ) { this->giveNode(i + 1)->giveUnknownVector(r, dofMask, VM_Total, tStep); s [ i ] = r.at(1); p [ i ].x = ( FPNum ) this->giveNode(i + 1)->giveCoordinate(1); p [ i ].y = ( FPNum ) this->giveNode(i + 1)->giveCoordinate(2); p [ i ].z = 0.; } context.updateFringeTableMinMax(s, 4); tr = CreateQuadWD3D(p, s [ 0 ], s [ 1 ], s [ 2 ], s [ 3 ]); EGWithMaskChangeAttributes(LAYER_MASK, tr); EMAddGraphicsToModel(ESIModel(), tr); } }