diff --git a/pub/lib/import.cpp b/pub/lib/import.cpp index 921a052ff347e88ef28ef7f3cf533eb1ace4ebb0..aebe61b6114f801d3d9b6a933695650d2eea0a97 100644 --- a/pub/lib/import.cpp +++ b/pub/lib/import.cpp @@ -141,7 +141,7 @@ void NImport::read_tab(string qualif) } if (!horizontal) { s = PSpec(new CSpec); - s->z.resize(nz); + s->z.resize(nzdat); for (int iz = 0; iz < nzdat; ++iz) s->z[iz] = RObjDbl(new CObjDbl(zdat[iz])); } diff --git a/pub/trivia/integrate.cpp b/pub/trivia/integrate.cpp index 3993cfb16bd72710eb7a79751441d0e8e4311a7f..f5ce5ed81df9155052e3439910e47c6955167756 100644 --- a/pub/trivia/integrate.cpp +++ b/pub/trivia/integrate.cpp @@ -52,6 +52,8 @@ double triv::integrate( case 6: rule = GSL_INTEG_GAUSS15; break; + default: + throw "Impossible case"; } // QAG adaptive integration gsl_integration_qag(&F, low, hig, epsabs, epsrel, nwork, rule, work, &val, &err);