fixed ispc & removed AVX2 kernels

still need to implement the trefftz cd kernel rip avx2
samayala22 · Feb 14, 2024 · 655d7e7 · 655d7e7
1 parent 0713929
commit 655d7e7
Show file tree

Hide file tree

Showing 13 changed files with 184 additions and 352 deletions.
diff --git a/mesh/infinite_rectangular_2x8.x b/mesh/infinite_rectangular_2x8.x
@@ -0,0 +1,23 @@
+       1
+       3        9        1
+     0.0000000000000000     0.50000000000000000      1.0000000000000000      0.0000000000000000
+    0.50000000000000000      1.0000000000000000      0.0000000000000000     0.50000000000000000
+     1.0000000000000000      0.0000000000000000     0.50000000000000000      1.0000000000000000
+     0.0000000000000000     0.50000000000000000      1.0000000000000000      0.0000000000000000
+    0.50000000000000000      1.0000000000000000      0.0000000000000000     0.50000000000000000
+     1.0000000000000000      0.0000000000000000     0.50000000000000000      1.0000000000000000
+     0.0000000000000000     0.50000000000000000      1.0000000000000000
+     0.0000000000000000      0.0000000000000000      0.0000000000000000      3750.0000000000000
+     3750.0000000000005      3750.0000000000000      7500.0000000000000      7500.0000000000009
+     7500.0000000000000      11250.000000000000      11250.000000000000      11250.000000000000
+     15000.000000000000      15000.000000000002      15000.000000000000      18750.000000000000
+     18750.000000000000      18750.000000000000      17500.000000000000      17500.000000000000
+     17500.000000000000      13750.000000000000      13750.000000000000      13750.000000000000
+     10000.000000000000      10000.000000000000      10000.000000000000
+     0.0000000000000000      0.0000000000000000      0.0000000000000000      0.0000000000000000
+     0.0000000000000000      0.0000000000000000      0.0000000000000000      0.0000000000000000
+     0.0000000000000000      0.0000000000000000      0.0000000000000000      0.0000000000000000
+     0.0000000000000000      0.0000000000000000      0.0000000000000000      0.0000000000000000
+     0.0000000000000000      0.0000000000000000      0.0000000000000000      0.0000000000000000
+     0.0000000000000000      0.0000000000000000      0.0000000000000000      0.0000000000000000
+     0.0000000000000000      0.0000000000000000      0.0000000000000000
diff --git a/tests/test_non_linear.cpp b/tests/test_non_linear.cpp
@@ -64,8 +64,8 @@ int main(int argc, char** argv) {
     const std::vector<std::string> backends = {"cpu"};
     std::vector<std::pair<std::string, std::unique_ptr<LiftCurveFunctor>>> lift_curves;
     lift_curves.emplace_back(std::make_pair("spallart1", std::make_unique<SpallartLiftCurve>(1.2f, 0.28f, 0.02f, 2.f*PI_f, 2.f*PI_f)));
-    lift_curves.emplace_back(std::make_pair("spallart2", std::make_unique<SpallartLiftCurve>(0.72f, 0.28f, 0.04f, 2.f*PI_f, 1.5f*PI_f)));
-    lift_curves.emplace_back(std::make_pair("polar", std::make_unique<ThinAirfoilPolarLiftCurve>()));
+    // lift_curves.emplace_back(std::make_pair("spallart2", std::make_unique<SpallartLiftCurve>(0.72f, 0.28f, 0.04f, 2.f*PI_f, 1.5f*PI_f)));
+    // lift_curves.emplace_back(std::make_pair("polar", std::make_unique<ThinAirfoilPolarLiftCurve>()));
 
     std::vector<f32> test_alphas = {0, 5, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20};
     std::transform(test_alphas.begin(), test_alphas.end(), test_alphas.begin(), to_radians);

diff --git a/vlm/backends/cpu/include/vlm_backend_cpu.hpp b/vlm/backends/cpu/include/vlm_backend_cpu.hpp
@@ -5,16 +5,17 @@
 
 namespace vlm {
 
-class BackendGeneric : public Backend {
+class BackendCPU : public Backend {
     public:
         std::vector<f32> lhs;
         std::vector<f32> rhs;
         std::vector<i32> ipiv;
         std::vector<f32> gamma;
         std::vector<f32> delta_gamma;
+        std::vector<f32> trefftz_buffer;
 
-        BackendGeneric(Mesh& mesh);
-        ~BackendGeneric();
+        BackendCPU(Mesh& mesh);
+        ~BackendCPU();
         void reset() override;
         void compute_lhs(const FlowData& flow) override;
         void compute_rhs(const FlowData& flow) override;

diff --git a/vlm/backends/cpu/include/vlm_backend_cpu_kernels.hpp b/vlm/backends/cpu/include/vlm_backend_cpu_kernels.hpp
@@ -9,4 +9,5 @@ void kernel_influence(
     f32 normalx[], f32 normaly[], f32 normalz[],
     u32 ia, u32 lidx, f32 sigma_p4
     );
-}
+}
+