Commit 980c5e18 authored by Laurent Belcour's avatar Laurent Belcour

Updating moment analysis

parent 51c5fa4d
......@@ -151,6 +151,17 @@ int main(int argc, char** argv)
}
// Get the minimum element of the dataset
vec _d_min = d->get(0);
for(int i=1; i<d->size(); ++i)
{
const vec y = d->get(i);
for(int j=0; j<nY; ++j)
{
_d_min[nX + j] = std::min(_d_min[nX + j], y[nX + j]);
}
}
// Compute the volume in which we integrate and then compute the
// dt to apply for each element.
double dt = 1.0;
......@@ -219,6 +230,8 @@ int main(int argc, char** argv)
vec y = d->value(x);
for(int k=0; k<nY; ++k)
{
y[k] -= _d_min[nX + k];
double val = y[k] * dt;
m_0[k] += val ;
......
......@@ -22,7 +22,7 @@
<parameter name="param" value="RETRO_TL_TVL_PROJ_DPHI" />
<parameter name="min" value="[0.2, -2, -2]" />
<parameter name="max" value="[0.3, 2, 2]" />
<parameter name="ymin" value="[5.0]" />
<parameter name="ymin" value="[2.5]" />
</action>
<!-- Compute the moments on the 15 degrees data slice -->
......@@ -52,7 +52,7 @@
<parameter name="param" value="RETRO_TL_TVL_PROJ_DPHI" />
<parameter name="min" value="[0.3, -2, -2]" />
<parameter name="max" value="[0.8, 2, 2]" />
<parameter name="ymin" value="[5.0]" />
<parameter name="ymin" value="[2.5]" />
</action>
<!-- Compute the moments on the 30 degrees data slice -->
......@@ -82,7 +82,7 @@
<parameter name="param" value="RETRO_TL_TVL_PROJ_DPHI" />
<parameter name="min" value="[0.9, -2, -2]" />
<parameter name="max" value="[1.1, 2, 2]" />
<parameter name="ymin" value="[5.0]" />
<parameter name="ymin" value="[2.5]" />
</action>
<!-- Compute the moments on the 60 degrees data slice -->
......@@ -113,7 +113,7 @@
<parameter name="param" value="SCHLICK_TL_TK_PROJ_DPHI" />
<parameter name="min" value="[0.2, -2, -2]" />
<parameter name="max" value="[0.3, 2, 2]" />
<parameter name="ymin" value="[5.0]" />
<parameter name="ymin" value="[2.5]" />
</action>
<!-- Compute the moments on the 15 degrees data slice -->
......@@ -143,7 +143,7 @@
<parameter name="param" value="SCHLICK_TL_TK_PROJ_DPHI" />
<parameter name="min" value="[0.5, -2, -2]" />
<parameter name="max" value="[0.6, 2, 2]" />
<parameter name="ymin" value="[5.0]" />
<parameter name="ymin" value="[2.5]" />
</action>
<!-- Compute the moments on the 15 degrees data slice -->
......@@ -173,7 +173,7 @@
<parameter name="param" value="SCHLICK_TL_TK_PROJ_DPHI" />
<parameter name="min" value="[0.9, -2, -2]" />
<parameter name="max" value="[1.1, 2, 2]" />
<parameter name="ymin" value="[5.0]" />
<parameter name="ymin" value="[2.5]" />
</action>
<!-- Compute the moments on the 15 degrees data slice -->
......
......@@ -22,7 +22,7 @@
<parameter name="param" value="RETRO_TL_TVL_PROJ_DPHI" />
<parameter name="min" value="[0.2, -2, -2]" />
<parameter name="max" value="[0.3, 2, 2]" />
<parameter name="ymin" value="[1.0]" />
<parameter name="ymin" value="[3.0]" />
</action>
<!-- Compute the moments on the 15 degrees data slice -->
......@@ -52,7 +52,7 @@
<parameter name="param" value="RETRO_TL_TVL_PROJ_DPHI" />
<parameter name="min" value="[0.3, -2, -2]" />
<parameter name="max" value="[0.8, 2, 2]" />
<parameter name="ymin" value="[1.0]" />
<parameter name="ymin" value="[5.0]" />
</action>
<!-- Compute the moments on the 30 degrees data slice -->
......@@ -82,7 +82,7 @@
<parameter name="param" value="RETRO_TL_TVL_PROJ_DPHI" />
<parameter name="min" value="[0.9, -2, -2]" />
<parameter name="max" value="[1.1, 2, 2]" />
<parameter name="ymin" value="[1.0]" />
<parameter name="ymin" value="[0.5]" />
</action>
<!-- Compute the moments on the 60 degrees data slice -->
......@@ -113,7 +113,7 @@
<parameter name="param" value="SCHLICK_TL_TK_PROJ_DPHI" />
<parameter name="min" value="[0.2, -2, -2]" />
<parameter name="max" value="[0.3, 2, 2]" />
<parameter name="ymin" value="[1.0]" />
<parameter name="ymin" value="[5.0]" />
</action>
<!-- Compute the moments on the 15 degrees data slice -->
......@@ -143,7 +143,7 @@
<parameter name="param" value="SCHLICK_TL_TK_PROJ_DPHI" />
<parameter name="min" value="[0.5, -2, -2]" />
<parameter name="max" value="[0.6, 2, 2]" />
<parameter name="ymin" value="[1.0]" />
<parameter name="ymin" value="[3.0]" />
</action>
<!-- Compute the moments on the 15 degrees data slice -->
......@@ -173,7 +173,7 @@
<parameter name="param" value="SCHLICK_TL_TK_PROJ_DPHI" />
<parameter name="min" value="[0.9, -2, -2]" />
<parameter name="max" value="[1.1, 2, 2]" />
<parameter name="ymin" value="[1.0]" />
<parameter name="ymin" value="[0.5]" />
</action>
<!-- Compute the moments on the 15 degrees data slice -->
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment