nan_check.hpp

Go to the documentation of this file.

#include <Kokkos_Core.hpp>
#include <string>
#include "space_settings.hpp"

// Loop over 4D Host View and print the index of any detected NaNs
inline void nan_check_view(const Kokkos::View<double****,HostType>& view, std::string name){
#ifdef COL_F_NAN_CHECK
    printf("\nChecking for NaNs in %s",name.c_str());
    for(int i=0;i<view.extent(0);i++){
        for(int j=0;j<view.extent(1);j++){
            for(int k=0;k<view.extent(2);k++){
                for(int l=0;l<view.extent(3);l++){
                    if(isnan(view(i,j,k,l))){
                        printf("\nNaN found in %s at (%d,%d,%d,%d)",view.label().c_str(),i,j,k,l);
                    }
                }
            }
        }
    }
    fflush(stdout);
#endif
}