Hi Marc,
I propose the following (or a similar) modification to the code in qtbwave.cpp, in order to avoid bad values in output.
Ciao!
MODIFIED CODE:
-----------------------
for(int iMode = 0;iMode<*nModes;iMode++) {
const QtbValue * mode = dispersion.mode( iMode );
for (int iSample = 0; iSample<*nSamples; iSample++) {
if ( mode[iSample].isValid() )
*(slowness++) = mode[iSample].value();
else
*(slowness++) = -1;
}
}
ORIGINAL CODE
-----------------------
for(int iMode = 0;iMode<*nModes;iMode++) {
const QtbValue * mode = dispersion.mode( iMode );
for (int iSample = 0; iSample<*nSamples; iSample++) {
*(slowness++) = mode[iSample].value();
}
}
