Skip to content

Commit

Permalink
Update tests
Browse files Browse the repository at this point in the history
  • Loading branch information
@pghysels
pghysels committed Jun 4, 2024
1 parent 8e49ca1 commit a4e5ea6
Show file tree
Hide file tree
Showing 10 changed files with 582 additions and 131 deletions.
7 changes: 5 additions & 2 deletions examples/dense/CMakeLists.txt
View file
Original file line number Diff line number Diff line change
Expand Up @@ -17,8 +17,11 @@ target_link_libraries(testCovariance strumpack)
target_link_libraries(testQC strumpack)
target_link_libraries(testFront strumpack)

target_include_directories(testCovariance PRIVATE /home/pieterg/LBL/summer2022/random/mfem/mfem-4.4)
target_link_libraries(testCovariance /home/pieterg/LBL/summer2022/random/mfem/mfem-4.4/libmfem.a)
# target_include_directories(testCovariance PRIVATE /home/pieterg/LBL/summer2022/random/mfem/mfem-4.4)
# target_link_libraries(testCovariance /home/pieterg/LBL/summer2022/random/mfem/mfem-4.4/libmfem.a)
target_include_directories(testCovariance PRIVATE /global/homes/p/pghysels/perlmutter/mfem-4.7)
target_link_libraries(testCovariance /global/homes/p/pghysels/perlmutter/mfem-4.7/libmfem.a)


add_dependencies(examples
KernelRegression
Expand Down
53 changes: 30 additions & 23 deletions examples/dense/test2DSIE.cpp
View file
Original file line number Diff line number Diff line change
Expand Up @@ -102,33 +102,40 @@ int main(int argc, char* argv[]) {
TaskTimer tassmbly("");
tassmbly.start();
DenseMatrix<std::complex<double>> Lop(N, N);

std::string fname("Lop_" + std::to_string(N) + ".bin");
std::ifstream f(fname.c_str());
if (f.good())
Lop = DenseMatrix<std::complex<double>>::read(fname);
else {
#pragma omp parallel for
for (int i=0; i<N; i++) {
double p[] = {xyz(0,i), xyz(1,i)};
double k = w * n_num(p[0], p[1]);
for (int j=0; j<N; j++) {
Lop(i, j) = 0.;
if (i == j)
Lop(i, j) = 1. / (2. * M_PI) *
(dl[j] - dl[j] * std::log(dl[j] / 2.));
for (int aa=0; aa<nquad; aa++) {
auto nq = (aa - 0.5) / nquad;
double q[] = {pn0(0,j) + nq * (pn1(0,j)-pn0(0,j)),
pn0(1,j) + nq * (pn1(1,j)-pn0(1,j))};
double rvec[] = {p[0]-q[0], p[1]-q[1]};
auto r = norm(rvec, 2);
auto G = std::complex<double>(0, 1./4.) * BesselH0(k * r);
if (std::abs(i-j) > 0)
Lop(i, j) += dl[j] / nquad * G;
else {
auto G0 = -1. / (2. * M_PI) * std::log(r);
Lop(i, j) += dl[j] / nquad * (G - G0);
}
for (int i=0; i<N; i++) {
double p[] = {xyz(0,i), xyz(1,i)};
double k = w * n_num(p[0], p[1]);
for (int j=0; j<N; j++) {
Lop(i, j) = 0.;
if (i == j)
Lop(i, j) = 1. / (2. * M_PI) *
(dl[j] - dl[j] * std::log(dl[j] / 2.));
for (int aa=0; aa<nquad; aa++) {
auto nq = (aa - 0.5) / nquad;
double q[] = {pn0(0,j) + nq * (pn1(0,j)-pn0(0,j)),
pn0(1,j) + nq * (pn1(1,j)-pn0(1,j))};
double rvec[] = {p[0]-q[0], p[1]-q[1]};
auto r = norm(rvec, 2);
auto G = std::complex<double>(0, 1./4.) * BesselH0(k * r);
if (std::abs(i-j) > 0)
Lop(i, j) += dl[j] / nquad * G;
else {
auto G0 = -1. / (2. * M_PI) * std::log(r);
Lop(i, j) += dl[j] / nquad * (G - G0);
}
}
}
}
std::cout << "# SIE assembly time: " << tassmbly.elapsed() << std::endl;
Lop.write(fname);
}
std::cout << "# SIE assembly time: " << tassmbly.elapsed() << std::endl;


// strumpack::HSS::HSSOptions<std::complex<double>> hss_opts;
// hss_opts.set_from_command_line(argc, argv);
Expand Down
121 changes: 121 additions & 0 deletions experiments/pp_cov.sh
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,121 @@
#!/bin/bash

out=out_cov
comp=stable

dim=3

for k in 9 19 29; do
for tol in 1e-2 1e-4 1e-6; do
rm -rf tmp*
for nnz in 1 2 4 8; do
grep "errors:" ${out}/out_dim${dim}_k${k}*tol${tol}*${comp}_SJLT_nnz${nnz} | awk '{print $2}' > tmp_err_med_nnz${nnz}
grep "errors:" ${out}/out_dim${dim}_k${k}*tol${tol}*${comp}_SJLT_nnz${nnz} | awk '{print $3}' > tmp_err_min_nnz${nnz}
grep "errors:" ${out}/out_dim${dim}_k${k}*tol${tol}*${comp}_SJLT_nnz${nnz} | awk '{print $4}' > tmp_err_max_nnz${nnz}
done
grep "errors:" ${out}/out_dim${dim}_k${k}*tol${tol}*${comp}_Gaussian | awk '{print $2}' >> tmp_err_med_gaussian
grep "errors:" ${out}/out_dim${dim}_k${k}*tol${tol}*${comp}_Gaussian | awk '{print $3}' >> tmp_err_min_gaussian
grep "errors:" ${out}/out_dim${dim}_k${k}*tol${tol}*${comp}_Gaussian | awk '{print $4}' >> tmp_err_max_gaussian
grep "errors:" ${out}/out_dim${dim}_k${k}*tol${tol}*${comp}_SRHT | awk '{print $2}' >> tmp_err_med_srht
grep "errors:" ${out}/out_dim${dim}_k${k}*tol${tol}*${comp}_SRHT | awk '{print $3}' >> tmp_err_min_srht
grep "errors:" ${out}/out_dim${dim}_k${k}*tol${tol}*${comp}_SRHT | awk '{print $4}' >> tmp_err_max_srht
case $tol in
1e-2)
color='red'
mark='';;
1e-4)
color='green'
mark='square';;
1e-6)
color='blue'
mark='triangle';;
esac
printf '\t\\addplot[color=%s, mark=%s*, only marks, mark size=2pt, error bars/.cd, y dir=both,y explicit]\n' $color $mark
printf '\tcoordinates { %% %s, %s\n' $tol $k
printf '\t\t(G, %s) -= (0, %s) += (0, %s)\n' `cat tmp_err_med_gaussian` `cat tmp_err_min_gaussian` `cat tmp_err_max_gaussian`
printf '\t\t(S1, %s) -= (0, %s) += (0, %s)\n' `cat tmp_err_med_nnz1` `cat tmp_err_min_nnz1` `cat tmp_err_max_nnz1`
printf '\t\t(S2, %s) -= (0, %s) += (0, %s)\n' `cat tmp_err_med_nnz2` `cat tmp_err_min_nnz2` `cat tmp_err_max_nnz2`
printf '\t\t(S4, %s) -= (0, %s) += (0, %s)\n' `cat tmp_err_med_nnz4` `cat tmp_err_min_nnz4` `cat tmp_err_max_nnz4`
printf '\t\t(S8, %s) -= (0, %s) += (0, %s)\n' `cat tmp_err_med_nnz8` `cat tmp_err_min_nnz8` `cat tmp_err_max_nnz8`
printf '\t\t(H, %s) -= (0, %s) += (0, %s)\n' `cat tmp_err_med_srht` `cat tmp_err_min_srht` `cat tmp_err_max_srht`
printf '\t};\n'
echo ""
done
done


for k in 9 19 29; do
for tol in 1e-2 1e-4 1e-6; do
rm -rf tmp*
for nnz in 1 2 4 8; do
grep "ranks:" ${out}/out_dim${dim}_k${k}*tol${tol}*${comp}_SJLT_nnz${nnz} | awk '{print $2}' > tmp_ranks_med_nnz${nnz}
grep "ranks:" ${out}/out_dim${dim}_k${k}*tol${tol}*${comp}_SJLT_nnz${nnz} | awk '{print $3}' > tmp_ranks_min_nnz${nnz}
grep "ranks:" ${out}/out_dim${dim}_k${k}*tol${tol}*${comp}_SJLT_nnz${nnz} | awk '{print $4}' > tmp_ranks_max_nnz${nnz}
done
grep "ranks:" ${out}/out_dim${dim}_k${k}*tol${tol}*${comp}_Gaussian | awk '{print $2}' >> tmp_ranks_med_gaussian
grep "ranks:" ${out}/out_dim${dim}_k${k}*tol${tol}*${comp}_Gaussian | awk '{print $3}' >> tmp_ranks_min_gaussian
grep "ranks:" ${out}/out_dim${dim}_k${k}*tol${tol}*${comp}_Gaussian | awk '{print $4}' >> tmp_ranks_max_gaussian
grep "ranks:" ${out}/out_dim${dim}_k${k}*tol${tol}*${comp}_SRHT | awk '{print $2}' >> tmp_ranks_med_srht
grep "ranks:" ${out}/out_dim${dim}_k${k}*tol${tol}*${comp}_SRHT | awk '{print $3}' >> tmp_ranks_min_srht
grep "ranks:" ${out}/out_dim${dim}_k${k}*tol${tol}*${comp}_SRHT | awk '{print $4}' >> tmp_ranks_max_srht
case $tol in
1e-2)
color='red'
mark='';;
1e-4)
color='green'
mark='square';;
1e-6)
color='blue'
mark='triangle';;
esac
# case $k in
# 10000) mark='';;
# 20000) mark='square';;
# 40000) mark='triangle';;
# esac
# case $tol in
# 1e-2) color='red';;
# 1e-4) color='green';;
# 1e-6) color='blue';;
# esac
printf '\t\\addplot[color=%s, mark=%s*, only marks, mark size=2pt, error bars/.cd, y dir=both,y explicit]\n' $color $mark
printf '\tcoordinates { %% %s, %s\n' $tol $k
printf '\t\t(G, %s) -= (0, %s) += (0, %s)\n' `cat tmp_ranks_med_gaussian` `cat tmp_ranks_min_gaussian` `cat tmp_ranks_max_gaussian`
printf '\t\t(S1, %s) -= (0, %s) += (0, %s)\n' `cat tmp_ranks_med_nnz1` `cat tmp_ranks_min_nnz1` `cat tmp_ranks_max_nnz1`
printf '\t\t(S2, %s) -= (0, %s) += (0, %s)\n' `cat tmp_ranks_med_nnz2` `cat tmp_ranks_min_nnz2` `cat tmp_ranks_max_nnz2`
printf '\t\t(S4, %s) -= (0, %s) += (0, %s)\n' `cat tmp_ranks_med_nnz4` `cat tmp_ranks_min_nnz4` `cat tmp_ranks_max_nnz4`
printf '\t\t(S8, %s) -= (0, %s) += (0, %s)\n' `cat tmp_ranks_med_nnz8` `cat tmp_ranks_min_nnz8` `cat tmp_ranks_max_nnz8`
printf '\t\t(H, %s) -= (0, %s) += (0, %s)\n' `cat tmp_ranks_med_srht` `cat tmp_ranks_min_srht` `cat tmp_ranks_max_srht`
printf '\t};\n'
echo ""

done
done

for tol in 1e-2 1e-4 1e-6; do
echo $tol
for k in 9 19 29; do
rm -rf tmp*
for nnz in 1 2 4 8; do
grep "A\*S time\|AT\*S time" ${out}/out_dim${dim}_k${k}*tol${tol}*${comp}_SJLT_nnz${nnz} | awk '{sum+=$5} END {print sum/3.0/1000.0}' > tmp_sample_nnz${nnz}
grep "times:" ${out}/out_dim${dim}_k${k}*tol${tol}*${comp}_SJLT_nnz${nnz} | awk '{print $2/1000.0}' > tmp_constr_nnz${nnz}
done

grep "A\*S time\|AT\*S time" ${out}/out_dim${dim}_k${k}*tol${tol}*${comp}_Gaussian | awk '{sum+=$5} END {print sum/3.0/1000.0}' >> tmp_sample_gaussian
grep "times:" ${out}/out_dim${dim}_k${k}*tol${tol}*${comp}_Gaussian | awk '{print $2/1000.0}' >> tmp_constr_gaussian

grep "A\*S time\|AT\*S time" ${out}/out_dim${dim}_k${k}*tol${tol}*${comp}_SRHT | awk '{sum+=$5} END {print sum/3.0/1000.0}' >> tmp_sample_srht
grep "times:" ${out}/out_dim${dim}_k${k}*tol${tol}*${comp}_SRHT | awk '{print $2/1000.0}' >> tmp_constr_srht

grep "Number of unknowns:" ${out}/out_dim${dim}_k${k}*tol${tol}*${comp}_Gaussian | awk '{print $4/1000.0}' >> tmp_N_gaussian
grep "\% of dense" ${out}/out_dim${dim}_k${k}*tol${tol}*${comp}_Gaussian | awk '{print $6}' | tail -n 1 | sed 's/%//'>> tmp_compr
# printf '& & %s & %5.0f & %5.0f & %5.0f & %5.0f & %5.0f & %5.0f & %5.0f & %5.0f & %5.0f & %5.0f & %5.0f & %5.0f & %5.1f \\\\ \n' $k `cat tmp_sample_gaussian` `cat tmp_sample_nnz1` `cat tmp_sample_nnz2` `cat tmp_sample_nnz4` `cat tmp_sample_nnz8` `cat tmp_sample_srht` `cat tmp_constr_gaussian` `cat tmp_constr_nnz1` `cat tmp_constr_nnz2` `cat tmp_constr_nnz4` `cat tmp_constr_nnz8` `cat tmp_constr_srht` `cat tmp_compr`

# printf '& & %s & %5.3f & %5.3f & %5.3f & %5.3f & %5.3f & %5.3f & %5.3f & %5.3f & %5.3f & %5.3f & %5.3f & %5.3f & %5.1f \\\\ \n' $k `cat tmp_sample_gaussian` `cat tmp_sample_nnz1` `cat tmp_sample_nnz2` `cat tmp_sample_nnz4` `cat tmp_sample_nnz8` `cat tmp_sample_srht` `cat tmp_constr_gaussian` `cat tmp_constr_nnz1` `cat tmp_constr_nnz2` `cat tmp_constr_nnz4` `cat tmp_constr_nnz8` `cat tmp_constr_srht` `cat tmp_compr`
printf '& & %s & %3.3g & %3.3g & %3.3g & %3.3g & %3.3g & %3.3g & %3.3g & %3.3g & %3.3g & %3.3g & %3.3g & %3.3g & %5.1f \\\\ \n' $k `cat tmp_sample_gaussian` `cat tmp_sample_nnz1` `cat tmp_sample_nnz2` `cat tmp_sample_nnz4` `cat tmp_sample_nnz8` `cat tmp_sample_srht` `cat tmp_constr_gaussian` `cat tmp_constr_nnz1` `cat tmp_constr_nnz2` `cat tmp_constr_nnz4` `cat tmp_constr_nnz8` `cat tmp_constr_srht` `cat tmp_compr`
done
echo ""
done

rm tmp*
116 changes: 116 additions & 0 deletions experiments/pp_front.sh
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,116 @@
#!/bin/bash

out=out_front
comp=stable

for k in 100 150 200; do
for tol in 1e-2 1e-4 1e-6; do
rm -rf tmp*
for nnz in 1 2 4 8; do
grep "errors:" ${out}/out_dim_k${k}*tol${tol}*${comp}_SJLT_nnz${nnz} | awk '{print $2}' > tmp_err_med_nnz${nnz}
grep "errors:" ${out}/out_dim_k${k}*tol${tol}*${comp}_SJLT_nnz${nnz} | awk '{print $3}' > tmp_err_min_nnz${nnz}
grep "errors:" ${out}/out_dim_k${k}*tol${tol}*${comp}_SJLT_nnz${nnz} | awk '{print $4}' > tmp_err_max_nnz${nnz}
done
grep "errors:" ${out}/out_dim_k${k}*tol${tol}*${comp}_Gaussian | awk '{print $2}' >> tmp_err_med_gaussian
grep "errors:" ${out}/out_dim_k${k}*tol${tol}*${comp}_Gaussian | awk '{print $3}' >> tmp_err_min_gaussian
grep "errors:" ${out}/out_dim_k${k}*tol${tol}*${comp}_Gaussian | awk '{print $4}' >> tmp_err_max_gaussian
grep "errors:" ${out}/out_dim_k${k}*tol${tol}*${comp}_SRHT | awk '{print $2}' >> tmp_err_med_srht
grep "errors:" ${out}/out_dim_k${k}*tol${tol}*${comp}_SRHT | awk '{print $3}' >> tmp_err_min_srht
grep "errors:" ${out}/out_dim_k${k}*tol${tol}*${comp}_SRHT | awk '{print $4}' >> tmp_err_max_srht
case $tol in
1e-2)
color='red'
mark='';;
1e-4)
color='green'
mark='square';;
1e-6)
color='blue'
mark='triangle';;
esac
printf '\t\\addplot[color=%s, mark=%s*, only marks, mark size=2pt, error bars/.cd, y dir=both,y explicit]\n' $color $mark
printf '\tcoordinates { %% %s, %s\n' $tol $k
printf '\t\t(G, %s) -= (0, %s) += (0, %s)\n' `cat tmp_err_med_gaussian` `cat tmp_err_min_gaussian` `cat tmp_err_max_gaussian`
printf '\t\t(S1, %s) -= (0, %s) += (0, %s)\n' `cat tmp_err_med_nnz1` `cat tmp_err_min_nnz1` `cat tmp_err_max_nnz1`
printf '\t\t(S2, %s) -= (0, %s) += (0, %s)\n' `cat tmp_err_med_nnz2` `cat tmp_err_min_nnz2` `cat tmp_err_max_nnz2`
printf '\t\t(S4, %s) -= (0, %s) += (0, %s)\n' `cat tmp_err_med_nnz4` `cat tmp_err_min_nnz4` `cat tmp_err_max_nnz4`
printf '\t\t(S8, %s) -= (0, %s) += (0, %s)\n' `cat tmp_err_med_nnz8` `cat tmp_err_min_nnz8` `cat tmp_err_max_nnz8`
printf '\t\t(H, %s) -= (0, %s) += (0, %s)\n' `cat tmp_err_med_srht` `cat tmp_err_min_srht` `cat tmp_err_max_srht`
printf '\t};\n'
echo ""
done
done


for k in 100 150 200; do
for tol in 1e-2 1e-4 1e-6; do
rm -rf tmp*
for nnz in 1 2 4 8; do
grep "ranks:" ${out}/out_dim_k${k}*tol${tol}*${comp}_SJLT_nnz${nnz} | awk '{print $2}' > tmp_ranks_med_nnz${nnz}
grep "ranks:" ${out}/out_dim_k${k}*tol${tol}*${comp}_SJLT_nnz${nnz} | awk '{print $3}' > tmp_ranks_min_nnz${nnz}
grep "ranks:" ${out}/out_dim_k${k}*tol${tol}*${comp}_SJLT_nnz${nnz} | awk '{print $4}' > tmp_ranks_max_nnz${nnz}
done
grep "ranks:" ${out}/out_dim_k${k}*tol${tol}*${comp}_Gaussian | awk '{print $2}' >> tmp_ranks_med_gaussian
grep "ranks:" ${out}/out_dim_k${k}*tol${tol}*${comp}_Gaussian | awk '{print $3}' >> tmp_ranks_min_gaussian
grep "ranks:" ${out}/out_dim_k${k}*tol${tol}*${comp}_Gaussian | awk '{print $4}' >> tmp_ranks_max_gaussian
grep "ranks:" ${out}/out_dim_k${k}*tol${tol}*${comp}_SRHT | awk '{print $2}' >> tmp_ranks_med_srht
grep "ranks:" ${out}/out_dim_k${k}*tol${tol}*${comp}_SRHT | awk '{print $3}' >> tmp_ranks_min_srht
grep "ranks:" ${out}/out_dim_k${k}*tol${tol}*${comp}_SRHT | awk '{print $4}' >> tmp_ranks_max_srht
case $tol in
1e-2)
color='red'
mark='';;
1e-4)
color='green'
mark='square';;
1e-6)
color='blue'
mark='triangle';;
esac
# case $k in
# 10000) mark='';;
# 20000) mark='square';;
# 40000) mark='triangle';;
# esac
# case $tol in
# 1e-2) color='red';;
# 1e-4) color='green';;
# 1e-6) color='blue';;
# esac
printf '\t\\addplot[color=%s, mark=%s*, only marks, mark size=2pt, error bars/.cd, y dir=both,y explicit]\n' $color $mark
printf '\tcoordinates { %% %s, %s\n' $tol $k
printf '\t\t(G, %s) -= (0, %s) += (0, %s)\n' `cat tmp_ranks_med_gaussian` `cat tmp_ranks_min_gaussian` `cat tmp_ranks_max_gaussian`
printf '\t\t(S1, %s) -= (0, %s) += (0, %s)\n' `cat tmp_ranks_med_nnz1` `cat tmp_ranks_min_nnz1` `cat tmp_ranks_max_nnz1`
printf '\t\t(S2, %s) -= (0, %s) += (0, %s)\n' `cat tmp_ranks_med_nnz2` `cat tmp_ranks_min_nnz2` `cat tmp_ranks_max_nnz2`
printf '\t\t(S4, %s) -= (0, %s) += (0, %s)\n' `cat tmp_ranks_med_nnz4` `cat tmp_ranks_min_nnz4` `cat tmp_ranks_max_nnz4`
printf '\t\t(S8, %s) -= (0, %s) += (0, %s)\n' `cat tmp_ranks_med_nnz8` `cat tmp_ranks_min_nnz8` `cat tmp_ranks_max_nnz8`
printf '\t\t(H, %s) -= (0, %s) += (0, %s)\n' `cat tmp_ranks_med_srht` `cat tmp_ranks_min_srht` `cat tmp_ranks_max_srht`
printf '\t};\n'
echo ""

done
done

for tol in 1e-2 1e-4 1e-6; do
echo $tol
for k in 100 150 200; do
rm -rf tmp*
for nnz in 1 2 4 8; do
grep "A\*S time\|AT\*S time" ${out}/out_dim${dim}_k${k}*tol${tol}*${comp}_SJLT_nnz${nnz} | awk '{sum+=$5} END {print sum/3.0/1000.0}' > tmp_sample_nnz${nnz}
grep "times:" ${out}/out_dim${dim}_k${k}*tol${tol}*${comp}_SJLT_nnz${nnz} | awk '{print $2/1000.0}' > tmp_constr_nnz${nnz}
done

grep "A\*S time\|AT\*S time" ${out}/out_dim${dim}_k${k}*tol${tol}*${comp}_Gaussian | awk '{sum+=$5} END {print sum/3.0/1000.0}' >> tmp_sample_gaussian
grep "times:" ${out}/out_dim${dim}_k${k}*tol${tol}*${comp}_Gaussian | awk '{print $2/1000.0}' >> tmp_constr_gaussian

grep "A\*S time\|AT\*S time" ${out}/out_dim${dim}_k${k}*tol${tol}*${comp}_SRHT | awk '{sum+=$5} END {print sum/3.0/1000.0}' >> tmp_sample_srht
grep "times:" ${out}/out_dim${dim}_k${k}*tol${tol}*${comp}_SRHT | awk '{print $2/1000.0}' >> tmp_constr_srht

grep "Number of unknowns:" ${out}/out_dim${dim}_k${k}*tol${tol}*${comp}_Gaussian | awk '{print $4/1000.0}' >> tmp_N_gaussian
grep "\% of dense" ${out}/out_dim${dim}_k${k}*tol${tol}*${comp}_Gaussian | awk '{print $6}' | tail -n 1 | sed 's/%//'>> tmp_compr
printf '& & %s & %3.4g & %3.4g & %3.4g & %3.4g & %3.4g & %3.4g & %3.4g & %3.4g & %3.4g & %3.4g & %3.4g & %3.4g & %5.1f \\\\ \n' $k `cat tmp_sample_gaussian` `cat tmp_sample_nnz1` `cat tmp_sample_nnz2` `cat tmp_sample_nnz4` `cat tmp_sample_nnz8` `cat tmp_sample_srht` `cat tmp_constr_gaussian` `cat tmp_constr_nnz1` `cat tmp_constr_nnz2` `cat tmp_constr_nnz4` `cat tmp_constr_nnz8` `cat tmp_constr_srht` `cat tmp_compr`
done
echo ""
done

rm tmp*
Loading

0 comments on commit a4e5ea6

Please sign in to comment.