Commit e72ca507 authored by Spencer Sherwin's avatar Spencer Sherwin
Browse files

Updated v_Integral method so that default weights in the collapsed coordinate...

Updated v_Integral method so that default weights in the collapsed coordinate direction are the given weights scaled by the appropriate Jacobian. This implies that the provided weights are of for a standard Legendre inner product unless otherwise specfiied by another switch case.
parent d7510ee8
......@@ -90,22 +90,18 @@ namespace Nektar
switch(m_base[1]->GetPointsType())
{
case LibUtilities::eGaussLobattoLegendre: // Legendre inner product
{
for(i = 0; i < nquad1; ++i)
{
w1_tmp[i] = 0.5*(1-z1[i])*w1[i];
}
break;
}
case LibUtilities::eGaussRadauMAlpha1Beta0: // (0,1) Jacobi Inner product
case LibUtilities::eGaussRadauMAlpha1Beta0: // (0,1) Jacobi Inner product
{
Vmath::Smul(nquad1, 0.5, w1, 1, w1_tmp,1);
break;
}
default:
default:
{
ASSERTL0(false, "populate swith for this point type");
// include jacobian factor on whatever coordinates are defined.
for(i = 0; i < nquad1; ++i)
{
w1_tmp[i] = 0.5*(1-z1[i])*w1[i];
}
break;
}
}
......
Supports Markdown
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