Fortran 等价
Fortran Equivalence
我想使用 Equivalence
。如何在 gfortran-5 中使用?
我正在使用:
Select Case (dst)
Case ("cm")
Equivalence :: (cm2Mm_r32, r32), (cm2Mm_r64, r64)
Call qcdiv (qb, qa, r32, r64, cm2Mm_r128)
Case ("m")
Equivalence :: (m2Mm_r32, r32), (m2Mm_r64, r64)
Call qcdiv (qb, qa, r32, r64, m2Mm_r128)
Case ("Km")
Equivalence :: (Km2Mm_r32, r32), (Km2Mm_r64, r64)
Call qcdiv (qb, qa, r32, r64, Km2Mm_r128)
End Select
等价语句必须位于声明部分。你不能用条件来控制它。
现代 Fortran 更喜欢 transfer()
转换位模式。如果不需要类型转换,使用指针。
我实际上在代码的一个地方使用了等价性,英特尔编译器无法优化 transfer()
,但这是一个例外,你在使用它之前应该三思而后行。
我想使用 Equivalence
。如何在 gfortran-5 中使用?
我正在使用:
Select Case (dst)
Case ("cm")
Equivalence :: (cm2Mm_r32, r32), (cm2Mm_r64, r64)
Call qcdiv (qb, qa, r32, r64, cm2Mm_r128)
Case ("m")
Equivalence :: (m2Mm_r32, r32), (m2Mm_r64, r64)
Call qcdiv (qb, qa, r32, r64, m2Mm_r128)
Case ("Km")
Equivalence :: (Km2Mm_r32, r32), (Km2Mm_r64, r64)
Call qcdiv (qb, qa, r32, r64, Km2Mm_r128)
End Select
等价语句必须位于声明部分。你不能用条件来控制它。
现代 Fortran 更喜欢 transfer()
转换位模式。如果不需要类型转换,使用指针。
我实际上在代码的一个地方使用了等价性,英特尔编译器无法优化 transfer()
,但这是一个例外,你在使用它之前应该三思而后行。