Match boundary conditions of 2D FD test

Author: stevendargaville

tests/adv_diff_cg_supg.c

@@ -74,10 +74,10 @@ static inline PetscErrorCode ComputeSUPGStabilization(PetscInt dim, PetscReal h,
   PetscFunctionReturn(PETSC_SUCCESS);
 }
 
-// Dirichlet BC: u = 1.0 (for inflow boundaries)
+// Dirichlet BC: u = 0.0 (for inflow boundaries)
 static PetscErrorCode dirichlet_bc(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, void *ctx)
 {
-  u[0] = 1.0;
+  u[0] = 0.0;
   return PETSC_SUCCESS;
 }
 
@@ -404,7 +404,7 @@ int main(int argc, char **argv)
   PetscCall(SetupSUPG(dm));
 
   PetscCall(DMCreateGlobalVector(dm, &u));
-  PetscCall(VecSet(u, 0.0));
+  PetscCall(VecSet(u, 1.0));
   PetscCall(PetscObjectSetName((PetscObject)u, "adv_diff"));
   PetscCall(DMPlexSetSNESLocalFEM(dm, PETSC_FALSE, &options));
   // Only compute the jacobian once in case we're doing a second solve
@@ -423,7 +423,7 @@ int main(int argc, char **argv)
   // gpu implementation and we don't want a copy occuring back to the cpu
   if (second_solve)
   {
-   PetscCall(VecSet(u, 0.0));
+   PetscCall(VecSet(u, 1.0));
    PetscCall(SNESSolve(snes, NULL, u));
   }