From 594ef3d9f16a9f94badd0126a5f83556006ba617 Mon Sep 17 00:00:00 2001
From: abtin98 <abtin.ameri@yahoo.com>
Date: Wed, 11 Sep 2019 19:51:22 -0400
Subject: [PATCH] fixed burgers test case

---
 src/numerical_flux/split_form_numerical_flux.cpp | 8 ++++----
 1 file changed, 4 insertions(+), 4 deletions(-)

diff --git a/src/numerical_flux/split_form_numerical_flux.cpp b/src/numerical_flux/split_form_numerical_flux.cpp
index c190d2d15..dfa30347b 100644
--- a/src/numerical_flux/split_form_numerical_flux.cpp
+++ b/src/numerical_flux/split_form_numerical_flux.cpp
@@ -18,17 +18,17 @@ std::array<real, nstate> SplitFormNumFlux<dim,nstate,real>::evaluate_flux(
 	    conv_phys_split_flux = pde_physics->convective_numerical_split_flux (soln_int,soln_ext);
 	    //std::cout << "done evaluating the conv num split flux" <<std::endl;
 
-	   // const real conv_max_eig_int = pde_physics->max_convective_eigenvalue(soln_int);
+	    const real conv_max_eig_int = pde_physics->max_convective_eigenvalue(soln_int);
 	   // std::cout << "1st eig" << std::endl;
-	   // const real conv_max_eig_ext = pde_physics->max_convective_eigenvalue(soln_ext);
+	    const real conv_max_eig_ext = pde_physics->max_convective_eigenvalue(soln_ext);
 	    //std::cout << "2nd eig" << std::endl;
-	  //  const real conv_max_eig = std::max(conv_max_eig_int, conv_max_eig_ext);
+	    const real conv_max_eig = std::max(conv_max_eig_int, conv_max_eig_ext);
 
 	   // std::cout << "obtained the max eig" <<std::endl;
 	    // Scalar dissipation
 	    std::array<real, nstate> numerical_flux_dot_n;
 	    for (int s=0; s<nstate; s++) {
-	        numerical_flux_dot_n[s] = conv_phys_split_flux[s]*normal_int /*- 0.5 * conv_max_eig * (soln_ext[s]-soln_int[s])*/;
+	        numerical_flux_dot_n[s] = conv_phys_split_flux[s]*normal_int - 0.5 * conv_max_eig * (soln_ext[s]-soln_int[s]);
 	    }
 	   // std::cout << "about to return split num flux" <<std::endl;
 	    return numerical_flux_dot_n;