From 902c0032ecda2391819dda0d234612338b75a582 Mon Sep 17 00:00:00 2001 From: Geoffroy Lesur Date: Fri, 21 Feb 2025 10:35:00 +0100 Subject: [PATCH 1/5] add the job to the test suite --- .github/workflows/idefix-ci-jobs.yml | 2 ++ 1 file changed, 2 insertions(+) diff --git a/.github/workflows/idefix-ci-jobs.yml b/.github/workflows/idefix-ci-jobs.yml index dfb89e0f..b95385e4 100644 --- a/.github/workflows/idefix-ci-jobs.yml +++ b/.github/workflows/idefix-ci-jobs.yml @@ -177,6 +177,8 @@ jobs: run: scripts/ci/run-tests $IDEFIX_DIR/test/MHD/sphBragTDiffusion -all $TESTME_OPTIONS - name: Spherical anisotropic viscosity run: scripts/ci/run-tests $IDEFIX_DIR/test/MHD/sphBragViscosity -all $TESTME_OPTIONS + - name: Collisionless thermal conduction + run: scripts/ci/run-tests $IDEFIX_DIR/test/MHD/clessTDiffusion -all $TESTME_OPTIONS Examples: needs: [Fargo, Dust, Planet, ShearingBox, SelfGravity] From 24a62e49eba5c49af1393490dacef8918497ef16 Mon Sep 17 00:00:00 2001 From: Geoffroy Lesur Date: Fri, 21 Feb 2025 10:37:17 +0100 Subject: [PATCH 2/5] update reference commit # --- reference | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/reference b/reference index 5ecb9588..c3289420 160000 --- a/reference +++ b/reference @@ -1 +1 @@ -Subproject commit 5ecb95884f845a23fa0baf45492770c30e22b4fc +Subproject commit c328942083b618a85b84eb16ce9fd35e67c00597 From c34b8190030ffb12ae6598597fe2d364e0a0bf57 Mon Sep 17 00:00:00 2001 From: Geoffroy Lesur Date: Fri, 21 Feb 2025 23:21:28 +0100 Subject: [PATCH 3/5] clean up test --- test/MHD/clessTDiffusion/setup.cpp | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/test/MHD/clessTDiffusion/setup.cpp b/test/MHD/clessTDiffusion/setup.cpp index 651a81f5..dceafd12 100644 --- a/test/MHD/clessTDiffusion/setup.cpp +++ b/test/MHD/clessTDiffusion/setup.cpp @@ -45,7 +45,7 @@ void UserDefBoundary(Hydro *hydro, int dir, BoundarySide side, real t) { IdefixArray4D Vs = hydro->Vs; IdefixArray1D x1 = data->x[IDIR]; - real rc,vc,vwind0; + real rc,vwind0; real cs=cs_vescGlob*sqrt(2.); real va_vesc = va_vescGlob; real mu = va_vesc * sqrt(2.); @@ -96,8 +96,8 @@ void Setup::InitFlow(DataBlock &data) { // Create a host copy DataBlockHost d(data); - real r,th,rl; - real PonRho, vwind0, rc, vc; + real r,rl; + real PonRho, vwind0, rc; real cs=cs_vescGlob*sqrt(2.); @@ -110,7 +110,7 @@ void Setup::InitFlow(DataBlock &data) { for(int j = 0; j < d.np_tot[JDIR] ; j++) { for(int i = 0; i < d.np_tot[IDIR] ; i++) { r=d.x[IDIR](i); - th=d.x[JDIR](j); + real vwind; vwind = ParkerWind(r/rc) * cs; From 6f3bedd7fc7b425ad6cc86a0a537c8a7027b952c Mon Sep 17 00:00:00 2001 From: Geoffroy Lesur Date: Sat, 22 Feb 2025 11:33:27 +0100 Subject: [PATCH 4/5] too much cleaning --- test/MHD/clessTDiffusion/setup.cpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/test/MHD/clessTDiffusion/setup.cpp b/test/MHD/clessTDiffusion/setup.cpp index dceafd12..4e369b2d 100644 --- a/test/MHD/clessTDiffusion/setup.cpp +++ b/test/MHD/clessTDiffusion/setup.cpp @@ -45,7 +45,7 @@ void UserDefBoundary(Hydro *hydro, int dir, BoundarySide side, real t) { IdefixArray4D Vs = hydro->Vs; IdefixArray1D x1 = data->x[IDIR]; - real rc,vwind0; + real rc,vc,vwind0; real cs=cs_vescGlob*sqrt(2.); real va_vesc = va_vescGlob; real mu = va_vesc * sqrt(2.); From dbdd28d60a7f9f7ce605ba381b55a74c0f7e87ed Mon Sep 17 00:00:00 2001 From: Geoffroy Lesur Date: Sat, 22 Feb 2025 16:13:33 +0100 Subject: [PATCH 5/5] final cleanup --- test/MHD/clessTDiffusion/setup.cpp | 3 +-- 1 file changed, 1 insertion(+), 2 deletions(-) diff --git a/test/MHD/clessTDiffusion/setup.cpp b/test/MHD/clessTDiffusion/setup.cpp index 4e369b2d..c540c0de 100644 --- a/test/MHD/clessTDiffusion/setup.cpp +++ b/test/MHD/clessTDiffusion/setup.cpp @@ -45,7 +45,7 @@ void UserDefBoundary(Hydro *hydro, int dir, BoundarySide side, real t) { IdefixArray4D Vs = hydro->Vs; IdefixArray1D x1 = data->x[IDIR]; - real rc,vc,vwind0; + real rc,vwind0; real cs=cs_vescGlob*sqrt(2.); real va_vesc = va_vescGlob; real mu = va_vesc * sqrt(2.); @@ -53,7 +53,6 @@ void UserDefBoundary(Hydro *hydro, int dir, BoundarySide side, real t) { PonRho = cs*cs; rc = 0.25 / (cs_vescGlob*cs_vescGlob); - vc = cs; vwind0 = ParkerWind(1./rc) * cs; hydro->boundary->BoundaryFor("UserDefBoundary", dir, side,