Skip to content
Snippets Groups Projects
Commit ecba5846 authored by Matthew Webster's avatar Matthew Webster
Browse files

fixed unknown operation

parent 4c7f2c96
No related branches found
No related tags found
No related merge requests found
......@@ -194,19 +194,22 @@ int inputParser(int argc, char *argv[], short output_dt, bool forceOutputType=fa
for (int j=0;j<nsize;j++) column_volume.value(j,0,0)=temp_volume(x+j*(xoff!=1),y+j*(yoff!=1),z+j*(zoff!=1),t+j*(toff!=1));
//This goes along the appropriate axis (non unitary offset variable) and fills a "column" volume with data
if (string(argv[i]+2) == "max") inputVolume.value(x,y,z,t)=column_volume.max();
if (string(argv[i]+2) == "min") inputVolume.value(x,y,z,t)=column_volume.min();
if (string(argv[i]+2) == "mean") inputVolume.value(x,y,z,t)=(T)column_volume.mean();
if (string(argv[i]+2) == "std") inputVolume.value(x,y,z,t)=(T)column_volume.stddev();
if (string(argv[i]+2) == "maxn") inputVolume.value(x,y,z,t)=column_volume.maxcoordx();
if (string(argv[i]+2) == "median") inputVolume.value(x,y,z,t)=column_volume.percentile(0.5);
if (string(argv[i]+2) == "perc") inputVolume.value(x,y,z,t)=column_volume.percentile(atof(argv[i+1])/100.0);
if (string(argv[i]+2) == "ar1")
{
else if (string(argv[i]+2) == "min") inputVolume.value(x,y,z,t)=column_volume.min();
else if (string(argv[i]+2) == "mean") inputVolume.value(x,y,z,t)=(T)column_volume.mean();
else if (string(argv[i]+2) == "std") inputVolume.value(x,y,z,t)=(T)column_volume.stddev();
else if (string(argv[i]+2) == "maxn") inputVolume.value(x,y,z,t)=column_volume.maxcoordx();
else if (string(argv[i]+2) == "median") inputVolume.value(x,y,z,t)=column_volume.percentile(0.5);
else if (string(argv[i]+2) == "perc") inputVolume.value(x,y,z,t)=column_volume.percentile(atof(argv[i+1])/100.0);
else if (string(argv[i]+2) == "ar1") {
column_volume-=(T)column_volume.mean();
double sumsq=column_volume.sumsquares();
inputVolume(x,y,z,t)=0;
if(sumsq!=0) for (int k=1;k<nsize;k++) inputVolume(x,y,z,t)+=(T)(column_volume(k,0,0)*column_volume(k-1,0,0)/sumsq);
}
else {
cerr << "Error unknown Dimensionality operation: " << string(argv[i]) << endl;
return(1);
}
}
if (string(argv[i]+2) == "perc") i++;
}
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment