Adding UniversalAssembleBnd call to matrix multiply routines

library/MultiRegions/ContField2D.cpp

@@ -858,10 +858,10 @@ namespace Nektar
          *                      \f$N_{\mathrm{dof}}\f$.
          */
         void ContField2D::v_GeneralMatrixOp(
-                                       const GlobalMatrixKey             &gkey,
-                                       const Array<OneD,const NekDouble> &inarray,
-                                       Array<OneD,      NekDouble> &outarray,
-                                       CoeffState coeffstate)
+                const GlobalMatrixKey              &gkey,
+                const Array<OneD,const NekDouble>  &inarray,
+                      Array<OneD,      NekDouble>  &outarray,
+                      CoeffState                   coeffstate)
         {
             if(coeffstate == eGlobal)
             {
@@ -872,6 +872,7 @@ namespace Nektar
                 {
                     GlobalMatrixSharedPtr mat = GetGlobalMatrix(gkey);
                     mat->Multiply(inarray,outarray);
+                    m_locToGloMap->UniversalAssemble(outarray);
                 }
                 else
                 {

library/MultiRegions/ContField3D.cpp

@@ -542,6 +542,7 @@ namespace Nektar
               {
                   GlobalMatrixSharedPtr mat = GetGlobalMatrix(gkey);
                   mat->Multiply(inarray,outarray);
+                  m_locToGloMap->UniversalAssemble(outarray);
               }
               else
               {

library/MultiRegions/GlobalLinSysIterativeStaticCond.cpp

@@ -861,6 +861,9 @@ namespace Nektar
                 Array<OneD, NekDouble> out = pOutput+ nDir;
 
                 m_globalSchurCompl->Multiply(in,out);
+
+                m_locToGloMap->UniversalAssembleBnd(pOutput);
+
             }
             else
             {