Commits
1735 1735 | } |
1736 1736 | } |
1737 1737 | if (status==GSL_SUCCESS || status==GSL_EMAXITER) { |
1738 1738 | // Current policy is to assume that exceeding max |
1739 1739 | // number of iterations is not a deal-breaker, leave it |
1740 1740 | // to SNR calculation to decide if the results are |
1741 1741 | // useful. |
1742 1742 | casa_param(3*icor + 0, iant) = gsl_vector_get(res, iparam+0); |
1743 1743 | casa_param(3*icor + 1, iant) = gsl_vector_get(res, iparam+1); |
1744 1744 | casa_param(3*icor + 2, iant) = gsl_vector_get(res, iparam+2); |
1745 - | // Sssshhh! We leave the first SNR as the FFT one. |
1745 + | // Shhh! We leave the first SNR as the FFT one. |
1746 1746 | for (size_t i=1; i!=3; i++) { |
1747 1747 | casa_snr(3*icor + i, iant) = gsl_vector_get(snr_vector, iparam+0); |
1748 1748 | } |
1749 1749 | } else { // gsl solver failed; flag data |
1750 1750 | logSink << "Least-squares solver failed to converge; flagging" << endl; |
1751 1751 | casa_flags(3*icor + 0, iant) = false; |
1752 1752 | casa_flags(3*icor + 1, iant) = false; |
1753 1753 | casa_flags(3*icor + 2, iant) = false; |
1754 1754 | } |
1755 1755 | } |