#include "clitkImageUncertainty_ggo.h"
#include "clitkImageCommon.h"
#include "clitkCommon.h"
+#include "clitkPortability.h"
// itk include
#include "itkImageRegionConstIterator.h"
#include "itkImageRegionIterator.h"
-#include <cmath> // for isfinite
-
//====================================================================
int main(int argc, char * argv[]) {
double squared = pii.Get();
double mean = pi.Get();
double uncert = sqrt((NumberOfEvents*squared - mean*mean) / ((NumberOfEvents-1)*(mean*mean)));
- if (!std::isnormal(uncert)) uncert = 1.;
+ if (!IsNormal(uncert)) uncert = 1.;
po.Set(uncert);
++pi;
++pii;