fixed burgers test case

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;