Codeformat

Author: KatrinKoesler

src/model/ValveTanh.cpp

@@ -45,24 +45,21 @@ void ValveTanh::update_solution(
   double Rmax = parameters[global_param_ids[ParamId::RMAX]];
   double steep = parameters[global_param_ids[ParamId::STEEPNESS]];
 
-  //Helper functions
+  // Helper functions
   double fun_tanh = tanh(steep * (p_out - p_in));
   double fun_cosh = 0.5 * steep / pow(cosh(steep * (p_in - p_out)), 2);
-  
+
   // Nonlinear terms
-  system.C(global_eqn_ids[0]) =
-      -0.5 * q_in * (Rmax - Rmin) * fun_tanh;
+  system.C(global_eqn_ids[0]) = -0.5 * q_in * (Rmax - Rmin) * fun_tanh;
   system.C(global_eqn_ids[2]) = -0.5 * (1 + fun_tanh);
 
   // Derivatives of non-linear terms
   system.dC_dy.coeffRef(global_eqn_ids[0], global_var_ids[0]) =
-      0.5 * q_in * (Rmax - Rmin) * steep * (1.0 - pow(fun_tanh,2));
+      0.5 * q_in * (Rmax - Rmin) * steep * (1.0 - pow(fun_tanh, 2));
   system.dC_dy.coeffRef(global_eqn_ids[0], global_var_ids[1]) =
       -0.5 * (Rmax - Rmin) * fun_tanh;
   system.dC_dy.coeffRef(global_eqn_ids[0], global_var_ids[2]) =
-      -0.5 * q_in * (Rmax - Rmin) * steep * (1.0 - pow(fun_tanh,2));
-  system.dC_dy.coeffRef(global_eqn_ids[2], global_var_ids[0]) =
-      fun_cosh;
-  system.dC_dy.coeffRef(global_eqn_ids[2], global_var_ids[2]) =
-      -fun_cosh;
+      -0.5 * q_in * (Rmax - Rmin) * steep * (1.0 - pow(fun_tanh, 2));
+  system.dC_dy.coeffRef(global_eqn_ids[2], global_var_ids[0]) = fun_cosh;
+  system.dC_dy.coeffRef(global_eqn_ids[2], global_var_ids[2]) = -fun_cosh;
 }

src/model/ValveTanh.h

@@ -39,7 +39,7 @@
  * \f[
  * Q_{in}-Q_{out}=0
  * \f]
- *  
+ *
  * \f[
  * \text{valve\_status} = \frac{1}{2}\left(1+tanh\{k(P_{out}-P{in})\}\right)
  * \f]