如何使用FMZM高精度fortran打印变量
我想使用FM软件包的功能来研究数论中的一个问题。 我安装了包,编译了库,并运行了提供的两个测试套件,所有这些都没有任何问题 我写了如何使用FMZM高精度fortran打印变量,fortran,arbitrary-precision,fmzm,Fortran,Arbitrary Precision,Fmzm,我想使用FM软件包的功能来研究数论中的一个问题。 我安装了包,编译了库,并运行了提供的两个测试套件,所有这些都没有任何问题 我写了test.f: use fmzm TYPE (IM), SAVE :: n n = 0 WRITE(*,*) n end 并使用 gfortran -c -O3 -Wall test.f gfortran fmsave.o fm.o fmzm90.o test.o -o test 没有返回任何错误或警告。但发现了结果却让人失望:
test.f
:
use fmzm
TYPE (IM), SAVE :: n
n = 0
WRITE(*,*) n
end
并使用
gfortran -c -O3 -Wall test.f
gfortran fmsave.o fm.o fmzm90.o test.o -o test
没有返回任何错误或警告。但发现了结果却让人失望:
./test
200000
数字200000似乎是变量n的内存位置。我做了一点实验,如果我把n改为complextype(ZM)
它会输出200000,199999
。类似地,如果我声明并初始化两个变量而不是一个
如果我更改类型(IM),将::n
保存为整数n
,并按照上面的方式编译,我将得到预期的0
作为输出
如果我将代码替换为
n = 0
do
n=n+1
if (n < 10) WRITE(*,*) n
end do
n=0
做
n=n+1
如果(n<10)写入(*,*)n
结束
然后输出
200000
重复9次,然后停止。因此,WRITE
函数只查找位置,而不查找值<代码>打印也一样。正如Vladimir指出的,除非你写出一个标准类型的值,否则你似乎必须首先使用一种策略转换成字符串,例如
str = im_format('i100', n)
str = adjustl(str)
str = trim(str)
print*, str
中间两行可以省略。im的类型是什么?作为不熟悉FMZM的人,我没有理由期望程序的行为与
n
声明为整数的行为相同。FMZM是否支持并使用用户定义的派生类型输入/输出?或者它应该打印什么?请查看上的示例代码,您将看到您应该首先使用IM\u格式('i200',n)
获取一个字符串,然后打印这个字符串。因为IM
是一个派生类型,所以=
,@Vladimir谢谢,这似乎起到了作用!