Added a simple check of symmetries obeyed by the interactions.

fpm.toml

@@ -70,6 +70,11 @@ name = "test_periodictable"
 source-dir="test"
 main = "test_periodictable.f90"
 
+[[test]]
+name = "check_interactions_symmetries"
+source-dir="test"
+main = "check_interactions_symmetries.f90"
+
 [[test]]
 name = "bte_regression"
 source-dir="test"

src/interactions.f90

@@ -47,7 +47,8 @@ module interactions
        calculate_echimp_interaction_ibzk, calculate_el_rta_rates, &
        calculate_bound_scatt_rates, calculate_thinfilm_scatt_rates, &
        calculate_4ph_rta_rates, calculate_coarse_grained_3ph_vertex, &
-       calculate_W_fromcgV2, calculate_W3ph_OTF, calculate_Y_OTF
+       calculate_W_fromcgV2, calculate_W3ph_OTF, calculate_Y_OTF, &
+       Vm2_3ph
 
   !external chdir, system
 

test/3C-SiC/fpm_run_symcheck_caf.sh

@@ -0,0 +1,23 @@
+#!/bin/bash
+
+# Download data from public repo
+#curl -X GET "https://nomad-lab.eu/prod/v1/api/v1/uploads/5b-c4JYRSG6qitMFjFUUzQ/raw/?offset=0&length=-1&compress=true" -o 3C-SiC.zip
+
+#Extract data
+#unzip 3C-SiC.zip
+#rm 3C-SiC.zip
+
+workdir="./3C-SiC_test_output"
+datadir="./3C-SiC/elphbolt_input_data"
+inputdir="./test/3C-SiC"
+
+#mkdir -p $workdir
+#cp $datadir/* $workdir
+#cp $inputdir/input.nml $workdir
+#rm -r 3C-SiC
+cd $workdir
+
+#gcc+opencoarrays
+cafrun -n 1 ../build/caf_*/test/check_interactions_symmetries | tee 3C_SiC_test.output
+
+cd ..