adding test using dummy_vis to verify fourfit4 can process at least 256 channels

Author: jpbarrett

source/cpp_src/Applications/dummy_vis.cc

@@ -326,6 +326,8 @@ int main(int argc, char** argv)
     std::string baseline_shortname = "AB";
     std::string ref_station = "A";
     std::string rem_station = "B";
+    std::string ref_station_code = "Aa";
+    std::string rem_station_code = "Bb";
     std::string ref_station_name = "DUMMY_A";
     std::string rem_station_name = "DUMMY_B";
     std::string root_code = "abcdef";
@@ -563,17 +565,17 @@ int main(int argc, char** argv)
     // =========================================================
 
 
-    std::cout<<"baseline_name"<<baseline_name<<std::endl;
-    std::cout<<"ref station"<<ref_station<<std::endl;
-    std::cout<<"rem station"<<rem_station<<std::endl;
+    std::cout<<"baseline_name: "<<baseline_name<<std::endl;
+    std::cout<<"ref station: "<<ref_station<<std::endl;
+    std::cout<<"rem station: "<<rem_station<<std::endl;
 
     vis.Insert(std::string("name"), std::string("visibilities"));
     vis.Insert(std::string("baseline"), baseline_name);
     vis.Insert(std::string("baseline_shortname"), baseline_shortname);
-    vis.Insert(std::string("reference_station"), ref_station);
-    vis.Insert(std::string("remote_station"), rem_station);
-    vis.Insert(std::string("reference_station_name"), ref_station);
-    vis.Insert(std::string("remote_station_name"), rem_station);
+    vis.Insert(std::string("reference_station"), ref_station_code);
+    vis.Insert(std::string("remote_station"), rem_station_code);
+    vis.Insert(std::string("reference_station_name"), ref_station_name);
+    vis.Insert(std::string("remote_station_name"), rem_station_name);
     vis.Insert(std::string("reference_station_mk4id"), ref_station);
     vis.Insert(std::string("remote_station_mk4id"), rem_station);
     vis.Insert(std::string("correlation_date"), start_time);
@@ -622,10 +624,10 @@ int main(int argc, char** argv)
     wt.Insert(std::string("name"), std::string("weights"));
     wt.Insert(std::string("baseline"), baseline_name);
     wt.Insert(std::string("baseline_shortname"), baseline_shortname);
-    wt.Insert(std::string("reference_station"), ref_station);
-    wt.Insert(std::string("remote_station"), rem_station);
-    wt.Insert(std::string("reference_station_name"), ref_station);
-    wt.Insert(std::string("remote_station_name"), rem_station);
+    wt.Insert(std::string("reference_station"), ref_station_code);
+    wt.Insert(std::string("remote_station"), rem_station_code);
+    wt.Insert(std::string("reference_station_name"), ref_station_name);
+    wt.Insert(std::string("remote_station_name"), rem_station_name);
     wt.Insert(std::string("reference_station_mk4id"), ref_station);
     wt.Insert(std::string("remote_station_mk4id"), rem_station);
     wt.Insert(std::string("correlation_date"), std::string("2025-01-01T00:00:00"));
@@ -651,11 +653,11 @@ int main(int argc, char** argv)
     tags.Insert(std::string("root_code"), root_code);
     tags.Insert(std::string("baseline"), baseline_name);
     tags.Insert(std::string("baseline_shortname"), baseline_shortname);
-    tags.Insert(std::string("reference_station"), ref_station);
-    tags.Insert(std::string("remote_station"), rem_station);
+    tags.Insert(std::string("reference_station"), ref_station_code);
+    tags.Insert(std::string("remote_station"), rem_station_code);
     tags.Insert(std::string("reference_station_name"), ref_station_name);
-    tags.Insert(std::string("remote_station_name"), rem_station);
-    tags.Insert(std::string("reference_station_mk4id"), ref_station_name);
+    tags.Insert(std::string("remote_station_name"), rem_station_name);
+    tags.Insert(std::string("reference_station_mk4id"), ref_station);
     tags.Insert(std::string("remote_station_mk4id"), rem_station);
     tags.Insert(std::string("correlation_date"), std::string("2025-01-01T00:00:00"));
     tags.Insert(std::string("origin"), std::string("simulated"));
@@ -744,6 +746,8 @@ int main(int argc, char** argv)
     std::string vex_string = GetVexString();
     mho_json vex_obj = mho_json::parse(vex_string);
 
+
+
     //write out the 'ovex'/'vex' json object as a json file
     std::string root_file = "./dummy.";
     root_file += root_code + "." + "root.json";

source/cpp_src/Applications/test/CMakeLists.txt

@@ -35,6 +35,7 @@ if(BASH_PROGRAM)
             chk_adhoc_phase_poly
             chk_adhoc_phase_sinewave
             chk_adhoc_flagging
+            chk_highnchannels
         )
     endif(PY_EXE)
 

source/cpp_src/Applications/test/chk_highnchannels.sh

@@ -0,0 +1,55 @@
+#!/bin/bash
+verb=false
+[ -n "$testverb" ] && verb=true
+
+[ -d "$srcdir" ] || { echo srcdir not set; exit 1; }
+${HOPS_SETUP-'false'} || . $srcdir/chk_env.sh
+
+if [ ! -d "$srcdir/dummy" ]; then
+  mkdir -p "$srcdir/dummy"
+fi
+
+export DATADIR=`cd $srcdir/dummy; pwd`
+RET_VAL=0
+EXP_DIR=$DATADIR
+cd $EXP_DIR
+SNR_VAL=200
+NCHAN=256
+CHANW=64
+NAPS=30
+NSPECPTS=128
+
+#run dummy_vis to generate fake visibilty data to test if we can process n_channels > 128 (here we test 256)
+dummy_vis --snr ${SNR_VAL} -c ${NCHAN} -w ${CHANW} -n ${NAPS} -s ${NSPECPTS}
+
+echo "Running: fourfit4 -m 4 -b AB -P XX ./"
+output_file=$(fourfit4 -m 4 -b AB -P XX ./ 2>&1 | awk '{print $NF}')
+echo "fourfit4 output file: $output_file"
+
+#convert the fringe file to json
+hops2json ${output_file}
+
+#use jq (json query) to extract the plot_data element and pipe to file
+echo "jq '.[].tags.parameters.fringe.snr | select( . != null )' "${output_file}.json" > ./fdump.json"
+jq '.[].tags.parameters.fringe.snr | select( . != null )' "${output_file}.json" > ./fdump.json
+
+#use jq (json query) to extract the plot_data element and pipe to file
+echo "jq '.[].tags.parameters.config.nchannels | select( . != null )' "${output_file}.json" > ./nchan.json"
+jq '.[].tags.parameters.config.nchannels | select( . != null )' "${output_file}.json" > ./nchan.json
+
+NCHAN_USED=$(cat ./nchan.json)
+echo "fourfit4 used: $NCHAN_USED channels."
+
+#determine if the SNR is near the specified value of SNR_VAL
+SNR=$(cat ./fdump.json)
+
+result=$(echo "sqrt(($SNR - $SNR_VAL)^2) < 1.0" | bc -l)
+echo "The calculated SNR of the simulated fringe = $SNR, expected: $SNR_VAL"
+echo "Result ok? $result"  # prints 1 if true, 0 if false
+
+#invert the logic for the return code (0 is normal, 1 is fail)
+if [ "$result" -eq 1 ]; then
+    exit 0
+else
+    exit 1
+fi