float *
avtR2Fvariance::FinalizePass(int pass)
{
float *rv = NULL;
if (pass == 0)
{
double *rt2 = new double[nBins];
int *cnt2 = new int[nBins];
SumIntArrayAcrossAllProcessors(count, cnt2, nBins);
SumDoubleArrayAcrossAllProcessors(running_total_ave, rt2, nBins);
for (int i = 0 ; i < nBins ; i++)
{
if (cnt2[i] > 0)
running_total_ave[i] = rt2[i] / cnt2[i];
else
running_total_ave[i] = undefinedVal;
}
delete [] rt2;
delete [] cnt2;
}
else
{
rv = new float[nBins];
double *rt2 = new double[nBins];
int *cnt2 = new int[nBins];
SumIntArrayAcrossAllProcessors(count, cnt2, nBins);
SumDoubleArrayAcrossAllProcessors(running_total_variance, rt2, nBins);
for (int i = 0 ; i < nBins ; i++)
{
if (cnt2[i] > 0)
rv[i] = rt2[i] / cnt2[i];
else
rv[i] = undefinedVal;
}
delete [] rt2;
delete [] cnt2;
}
this->pass++;
return rv;
}
void
avtConnComponentsWeightedVariableQuery::PostExecute(void)
{
// sum vals across all processors
double *sum_res_dbl = new double[nComps];
SumDoubleArrayAcrossAllProcessors(&sumPerComp[0], sum_res_dbl, nComps);
memcpy(&sumPerComp[0],sum_res_dbl,nComps * sizeof(double));
delete [] sum_res_dbl;
// create output message
if(PAR_Rank() == 0)
{
std::string msg = "";
char buff[2048];
if(nComps == 1)
{SNPRINTF(buff,2048,"Found %d connected component\n",nComps);}
else
{SNPRINTF(buff,2048,"Found %d connected components\n",nComps);}
msg += buff;
std::string format = "Component %d Weighted Sum = ("
+ queryAtts.GetFloatFormat() +")\n";
for(int i=0;i<nComps;i++)
{
SNPRINTF(buff,1024,
format.c_str(),
i,
sumPerComp[i]);
msg += buff;
}
// set output message
SetResultMessage(msg);
// set output values
SetResultValues(sumPerComp);
}
}
float *
avtR2Frms::FinalizePass(int pass)
{
float *rv = new float[nBins];
double *rt2 = new double[nBins];
int *cnt2 = new int[nBins];
SumIntArrayAcrossAllProcessors(count, cnt2, nBins);
SumDoubleArrayAcrossAllProcessors(running_total, rt2, nBins);
for (int i = 0 ; i < nBins ; i++)
{
if (cnt2[i] > 0)
rv[i] = sqrt(rt2[i] / cnt2[i]);
else
rv[i] = undefinedVal;
}
delete [] rt2;
delete [] cnt2;
return rv;
}
void
avtIntegrationRF::OutputRawValues(const char *filename)
{
int nvals = windowSize[0]*windowSize[1];
double *out_vals = new double[nvals];
SumDoubleArrayAcrossAllProcessors(vals, out_vals, nvals);
if (PAR_Rank() == 0)
{
ofstream ofile(filename);
if (ofile.fail())
{
avtCallback::IssueWarning("The integration option for the ray"
" casted volume rendering could not output a file, because"
" it could not open the file.");
}
ofile << windowSize[0] << " " << windowSize[1] << endl;
for (int i = 0 ; i < nvals ; i++)
ofile << out_vals[i] << endl;
}
if (!issuedWarning)
{
char msg[1024];
SNPRINTF(msg, 1024, "The integration was outputted to the file \"%s\". "
"If you are running client/server, the file is located where "
"the remote server is running. This message will only be issued"
" once per VisIt session. Further renderings will overwrite"
" previous versions of this file.", filename);
avtCallback::IssueWarning(msg);
issuedWarning = true;
}
delete [] out_vals;
}
void
avtConnComponentsCentroidQuery::PostExecute(void)
{
// get # of cells per component (from all processors)
int *sum_res_int = new int[nComps];
SumIntArrayAcrossAllProcessors(&nCellsPerComp[0], sum_res_int, nComps);
memcpy(&nCellsPerComp[0],sum_res_int,nComps * sizeof(int));
delete [] sum_res_int;
// get centroid values (from all processors)
double *sum_res_dbl = new double[nComps];
SumDoubleArrayAcrossAllProcessors(&xCentroidPerComp[0],
sum_res_dbl,
nComps);
memcpy(&xCentroidPerComp[0],sum_res_dbl,nComps * sizeof(double));
SumDoubleArrayAcrossAllProcessors(&yCentroidPerComp[0],
sum_res_dbl,
nComps);
memcpy(&yCentroidPerComp[0],sum_res_dbl,nComps * sizeof(double));
SumDoubleArrayAcrossAllProcessors(&zCentroidPerComp[0],
sum_res_dbl,
nComps);
memcpy(&zCentroidPerComp[0],sum_res_dbl,nComps * sizeof(double));
delete [] sum_res_dbl;
// create output message
if(PAR_Rank() == 0)
{
std::string msg = "";
char buff[2048];
if(nComps == 1)
{SNPRINTF(buff,2048,"Found %d connected component\n",nComps);}
else
{SNPRINTF(buff,2048,"Found %d connected components\n",nComps);}
msg += buff;
// pack values into a a single vector for query output
std::vector<double> result_vec(nComps *3);
for(int i=0;i<nComps;i++)
{
// get number of cells for current component
double n_comp_cells = (double)nCellsPerComp[i];
// calculate centriod values for the current component
xCentroidPerComp[i] /= n_comp_cells;
yCentroidPerComp[i] /= n_comp_cells;
zCentroidPerComp[i] /= n_comp_cells;
// pack into result vector
result_vec[i*3 + 0] = xCentroidPerComp[i];
result_vec[i*3 + 1] = yCentroidPerComp[i];
result_vec[i*3 + 2] = zCentroidPerComp[i];
}
std::string format = "Component %d [%d cells] Centroid = ("
+ queryAtts.GetFloatFormat() +","
+ queryAtts.GetFloatFormat() +","
+ queryAtts.GetFloatFormat() +")\n";
// prepare the output message
for(int i=0;i<nComps;i++)
{
SNPRINTF(buff,1024,
format.c_str(),
i,
nCellsPerComp[i],
xCentroidPerComp[i],
yCentroidPerComp[i],
zCentroidPerComp[i]);
msg += buff;
}
// set result message
SetResultMessage(msg);
// set result values
SetResultValues(result_vec);
}
}
void
avtMassDistributionQuery::PostExecute(void)
{
int i;
avtWeightedVariableSummationQuery summer;
avtDataObject_p dob = GetInput();
summer.SetInput(dob);
QueryAttributes qa;
summer.PerformQuery(&qa);
double totalMass = qa.GetResultsValue()[0];
bool didVolume = false;
bool didRVolume = false;
bool didSA = false;
if (totalMass == 0.)
{
if (dob->GetInfo().GetAttributes().GetTopologicalDimension() == 3)
{
avtTotalVolumeQuery tvq;
avtDataObject_p dob = GetInput();
tvq.SetInput(dob);
QueryAttributes qa;
tvq.PerformQuery(&qa);
totalMass = qa.GetResultsValue()[0];
didVolume = true;
}
else if (dob->GetInfo().GetAttributes().GetMeshCoordType() == AVT_RZ
|| dob->GetInfo().GetAttributes().GetMeshCoordType() == AVT_ZR)
{
avtTotalRevolvedVolumeQuery rvq;
avtDataObject_p dob = GetInput();
rvq.SetInput(dob);
QueryAttributes qa;
rvq.PerformQuery(&qa);
totalMass = qa.GetResultsValue()[0];
didRVolume = true;
}
else
{
avtTotalSurfaceAreaQuery saq;
avtDataObject_p dob = GetInput();
saq.SetInput(dob);
QueryAttributes qa;
saq.PerformQuery(&qa);
totalMass = qa.GetResultsValue()[0];
didSA = true;
}
}
int times = 0;
char name[1024];
sprintf(name, "md%d.ult", times++);
if (PAR_Rank() == 0)
{
bool lookingForUnused = true;
while (lookingForUnused)
{
ifstream ifile(name);
if (ifile.fail())
lookingForUnused = false;
else
sprintf(name, "md%d.ult", times++);
}
}
char msg[1024];
const char *mass_string = (didVolume ? "volume" :
(didRVolume ? "revolved volume" :
(didSA ? "area" : "mass")));
std::string format = "The %s distribution has been outputted as an "
"Ultra file (%s), which can then be imported into VisIt. The"
" total %s considered was " + queryAtts.GetFloatFormat() +"\n";
SNPRINTF(msg,1024,format.c_str(),
mass_string, name, mass_string, totalMass);
SetResultMessage(msg);
SetResultValue(0.);
double *m2 = new double[numBins];
SumDoubleArrayAcrossAllProcessors(mass, m2, numBins);
delete [] mass;
mass = m2;
double totalMassFromLines = 0.;
for (i = 0 ; i < numBins ; i++)
totalMassFromLines += mass[i];
if (PAR_Rank() == 0)
{
if (totalMassFromLines == 0.)
{
sprintf(msg, "The mass distribution could not be calculated "
"becuase none of the lines intersected the data set."
" If you have used a fairly large number of lines, then "
"this may be indicative of an error state.");
SetResultMessage(msg);
return;
}
ofstream ofile(name);
if (ofile.fail())
{
sprintf(msg, "Unable to write out file containing distribution.");
SetResultMessage(msg);
}
if (!ofile.fail())
ofile << "# Mass distribution" << endl;
MapNode result_node;
doubleVector curve;
double binWidth = (maxLength-minLength) / numBins;
for (int i = 0 ; i < numBins ; i++)
{
double x1 = minLength + (i)*binWidth;
double x2 = minLength + (i+1)*binWidth;
double y = (totalMass*mass[i]) / (totalMassFromLines*binWidth);
curve.push_back(x1);
curve.push_back(y);
curve.push_back(x2);
curve.push_back(y);
if (!ofile.fail())
{
ofile << x1 << " " << y << endl;
ofile << x2 << " " << y << endl;
}
}
result_node["mass_distribution"] = curve;
SetXmlResult(result_node.ToXML());
SetResultValues(curve);
}
}
