1. Tests: with data Visit/2.13.0 failed opening data under /public/lchamand/Yisheng_chombo/out_143 small memory test with Visit/2.12.3: succeeded with pseuducolor plot rho, failed with px/rho big memory test: succeeded with pseducolor plot px/rho

Calculation using Expressions with scalar/vector mesh variables uses a lot computing resources..

Visit 2.12.3 on Bluehive bhc0198 with ~500GB memory crashes when plot "En/E_pot_rel_1"

AddPlot("Pseudocolor","En/E_pot_rel_1")
DrawPlots()
Query("Weighted Variable Sum")
erg_Epot_gas_par1_box = GetQueryOutputValue()

Definition structure of variable En/E_pot_rel_1

Temp/x1sc: 0.0000
Temp/y1sc: 0.0000
Temp/z1sc: 0.0000

=>

Pos/rVec: coord(Mesh)
Part/rVec_part1: {<Temp/x1sc>,<Temp/y1sc>,<Temp/z1sc>}

=>

Pos/rVec_rel_1: <Pos/rVec> -<Part/rVec_part1>

=>

Pos/r_rel_1: magnitude(<Pos/rVec_rel_1>)
Param/r_soft: 167600000000.000000
Param/G_Newton: 6.674e-8
Param/M_1: 729799999999999961087848627568640.000000

=>

Temp/u_1: <Pos/r_rel_1>/<Param/r_soft>
Pot/Phi0_rel_1: -1.0*<Param/G_Newton>*<Param/M_1>/<Pos/r_rel_1>
Temp/spline_A_1: -16./3*<Temp/u_1>^3 +48./5*<Temp/u_1>^5 -32./5*<Temp/u_1>^6 +14./5*<Temp/u_1>
Temp/spline_B_1:-1./15 -32./3*<Temp/u_1>^3 +16.*<Temp/u_1>^4 -48./5*<Temp/u_1>^5 +32./15*<Temp/u_1>^6 +48./15*<Temp/u_1>

=>

Temp/Phi_rel_1_A: if(lt(<Temp/u_1>,0.5), 1.0, 0.0)*<Pot/Phi0_rel_1>*<Temp/spline_A_1>
Temp/Phi_rel_1_B: if(and(ge(<Temp/u_1>,0.5),lt(<Temp/u_1>,1.)), 1.0, 0.0)*<Pot/Phi0_rel_1>*<Temp/spline_B_1>
Temp/Phi_rel_1_C: if(ge(<Temp/u_1>,1.), 1.0, 0.0)*<Pot/Phi0_rel_1>

=>

Pot/Phi_rel_1: <Temp/Phi_rel_1_A> +<Temp/Phi_rel_1_B> +<Temp/Phi_rel_1_C>

=>

En/E_pot_rel_1: rho*<Pot/Phi_rel_1>