04 2022 档案
摘要:fortran对于浮点数异常默认是忽略的,如果要开启,有两种方法,一种可以通过调用c函数去开启这部分异常处理,并且需要将函数打包成库,链接到fortran中,另一个是通过编译选项设置。 第一种方法的代码: trapfpe.c /* * ref: * https://gcc.gnu.org/onlin
阅读全文
摘要:2022.5.31更新: 发现这是WSL 1特有的问题,遇到除0的情况,不会报错,只会直接终止程序,并且用文中提到的编译参数: gfortran -g main2.f90 -ffpe-trap=invalid,zero,overflow 也不会给出报错信息,导致无法排查错误。换成WSL 2,编译的时
阅读全文