Commit e1c397fa authored by Chris Cantwell's avatar Chris Cantwell

Merge branch 'fix/DirectFullCoeffInitZeroFix' into 'master'

Added zero initialization for Coeff in DirectFull, added test for same

See merge request !898
parents 20f598c6 584782b2
......@@ -94,6 +94,7 @@ v4.4.2
- Fix deadlock in Hdf5 with homogeneous expansions (!858)
- Fix a crash when Interpolator is called on an empty field (!869)
- Fix petsc compile without MPI (!873)
- Fix uninitialised coefficients in DirectFull solver (!898)
**NekMesh**
- Fix missing periodic boundary meshing and boundary layer mesh adjustment
......
......@@ -128,6 +128,7 @@ namespace Nektar
}
else
{
Vmath::Zero(nGlobDofs-nDirDofs, &pOutput[nDirDofs], 1);
// Calculate Dirichlet forcing and subtract it from the rhs
m_expList.lock()->GeneralMatrixOp(
m_linSysKey, pOutput, tmp, eGlobal);
......
......@@ -182,6 +182,7 @@ IF( NEKTAR_SOLVER_ADR )
ADD_NEKTAR_TEST(ImDiffusion_m12)
ADD_NEKTAR_TEST(ImDiffusion_m6)
ADD_NEKTAR_TEST(ImDiffusion_VarCoeff)
ADD_NEKTAR_TEST(ImDiffusion_VarCoeff_DirectFull)
ADD_NEKTAR_TEST(ImDiffusion_Quad_Periodic_m7)
ADD_NEKTAR_TEST(ImDiffusion_Hex_Periodic_m5)
ADD_NEKTAR_TEST(SVV_Quad)
......
<?xml version="1.0" encoding="utf-8"?>
<test>
<description> 2D unsteady CG implicit diffusion, variable coeffs. </description>
<executable>ADRSolver</executable>
<parameters>-I GlobalSysSoln=DirectFull ImDiffusion_VarCoeff.xml</parameters>
<files>
<file description="Session File"> ImDiffusion_VarCoeff.xml </file>
</files>
<metrics>
<metric type="L2" id="1">
<value variable="u" tolerance="1e-12"> 1.68029e-08 </value>
</metric>
<metric type="Linf" id="2">
<value variable="u" tolerance="1e-12"> 3.25511e-08 </value>
</metric>
</metrics>
</test>
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