WIP

chaos-polymtl · Sep 24, 2024 · 1adb1c5 · 1adb1c5
1 parent 915cefd
commit 1adb1c5
Showing 1 changed file with 5 additions and 4 deletions.
diff --git a/source/solvers/tracer.cc b/source/solvers/tracer.cc
@@ -234,7 +234,7 @@ Tracer<dim>::assemble_system_matrix_dg()
         const unsigned int                                   &face_no,
         TracerScratchData<dim>                               &scratch_data,
         StabilizedMethodsCopyData                            &copy_data) {
-      double beta = 1;
+      double beta = 10;
       beta *= 1 / compute_cell_diameter<dim>(cell->measure(), 1);
 
       // Identify which boundary condition corresponds to the boundary id. If
@@ -330,7 +330,7 @@ Tracer<dim>::assemble_system_matrix_dg()
         const unsigned int                                   &nsf,
         TracerScratchData<dim>                               &scratch_data,
         StabilizedMethodsCopyData                            &copy_data) {
-      double beta = 0;
+      double beta = 10;
       beta *= 1 / compute_cell_diameter<dim>(cell->measure(), 1);
       FEInterfaceValues<dim> &fe_iv = scratch_data.fe_interface_values_tracer;
       fe_iv.reinit(cell, f, sf, ncell, nf, nsf);
@@ -602,7 +602,7 @@ Tracer<dim>::assemble_system_rhs_dg()
         const unsigned int                                   &face_no,
         TracerScratchData<dim>                               &scratch_data,
         StabilizedMethodsCopyData                            &copy_data) {
-      double beta = 1;
+      double beta = 10;
 
       beta *= 1 / compute_cell_diameter<dim>(cell->measure(), 1);
 
@@ -730,7 +730,8 @@ Tracer<dim>::assemble_system_rhs_dg()
         TracerScratchData<dim>                               &scratch_data,
         StabilizedMethodsCopyData                            &copy_data) {
       // TODO refactor and put inside a parameter formally
-      const double            beta  = 0.;
+      double beta = 10;
+      beta *= 1 / compute_cell_diameter<dim>(cell->measure(), 1);
       FEInterfaceValues<dim> &fe_iv = scratch_data.fe_interface_values_tracer;
       fe_iv.reinit(cell, f, sf, ncell, nf, nsf);