Fix the pipelines and remove unecessary test cases

pyTurbulence/solver.py

@@ -22,7 +22,7 @@ def solve(user_params)->dict:
         - "k0" (int): Peak wavenumber (default: 4).
         - "domain_size" (list, tuple, np.ndarray): Size of the simulation domain, should be a 3D vector (default: [2π, 2π, 2π]).
         - "domain_resolution" (list, tuple, np.ndarray): Resolution of the domain, should be a 3D vector (default: [64, 64, 64]).
-        - "Spectrum" (str): Type of spectrum for turbulence generation. Options: "PassotPouquet", "Gaussian" (default: "PassotPouquet").
+        - "Spectrum" (str): Type of spectrum for turbulence generation. Options: "PassotPouquet" (default: "PassotPouquet").
         - "gamma" (float): Heat capacity ratio, typically between 1 and 2 (default: 1.4).
         - "case" (int): Type of turbulent fields to generate. Options: 0, 1 (default: 1). 0: incompressible - 1: compressible.
         - "output_dir" (str): Directory to store output results (default: "results").
@@ -210,7 +210,7 @@ def _validate_params(params) -> dict:
 
     # Define allowed values for specific parameters
     allowed_values = {
-        "Spectrum": {"PassotPouquet", "Gaussian"},  # Allowed spectrum types
+        "Spectrum": {"PassotPouquet"},  # Allowed spectrum types
         "case": {0, 1},  # Allowed integer cases
     }
 

pyTurbulence/spectrum.py

@@ -26,8 +26,6 @@ def energy_spectrum(spectrum: str, k: np.ndarray, urms: float, k0: float) -> np.
     E_k = np.empty_like(k)
     if spectrum == "PassotPouquet":
         E_k = (urms**2 * 16. * np.sqrt(2. / np.pi) * k**4 / k0**5) * np.exp(-2. * (k / k0)**2)
-    if spectrum == "Gaussian":
-        E_k = 0.00013*k**4*np.exp(-2.*k**2/k0**2)
     return E_k
 
 def compute_tke_spectrum(u : np.ndarray,v : np.ndarray,w : np.ndarray,

test/test_Spectrum.py

@@ -14,15 +14,4 @@ def test_energy_spectrum_passot_pouquet():
     assert np.isclose(k0_spectrum, k0,rtol=1e-2), "Test failed: PassotPouquet spectrum does not match expected k0"
     assert np.isclose(TKE_spectrum, TKE, rtol=1e-2), "Test failed: PassotPouquet spectrum does not conserve TKE"
 
-def test_energy_spectrum_constant():
-    k = np.linspace(0,1024,1024)
-    TKE = 0.5
-    k0 = 4.0
-    spectrum = energy_spectrum("Constant", k, TKE, k0)
-    TKE_spectrum = np.trapz(spectrum, k)
-    idxMax = np.argmax(spectrum)
-    k0_spectrum = k[idxMax]
-    assert np.isclose(k0_spectrum, k0,1e-2), "Test failed: Constant spectrum does not match expected k0"
-
 test_energy_spectrum_passot_pouquet()
-test_energy_spectrum_constant()
\ No newline at end of file