Fortran内置函数SELECTED_REAL_KIND(p, r)，默认两个参数p是精度，r是范围。p是所需精度的十进制数值，r是以10^r次方表示的所需数值的范围。当执行该函数的时候，会返回达到或者超过指定精度或者范围的的实型数据的最小类别...

SOR迭代是在Gauss-Seidel迭代方法基础之上的进一步改进。其特征是取xk+1和xk的一个适当的加权平均来加快Gauss-Seidel收敛。对于方程组

当线性方程组的规模比较大时，采用高斯消元法需要太多时间。这时就要采用迭代法求解方程组了。高斯消元法是一个O(n^3)的浮点运算的有限序列，在经过有限步计算之后理论上得到的是精确解(无舍入误差时)。而迭代法在经过有...

与Jacobi迭代法密切相关的一种迭代方法叫做Gauss-Seidel迭代方法。Gauss-Seidel方法与Jacobi方法之间的差别是：在一个迭代步里，一旦未知变量值有更新，则立马投入使用。而不用像Jacobi方法那样下一个迭代步才使用。对于方...

这里的1为非测序个体, 2为测序个体. A11, A12, A21, A22可以由系谱构建的A矩阵提取. G为基因组构建的矩阵. H矩阵构建的相关代码见: 【GS专栏】全基因组选择中如何构建H矩阵....

当一个程序出现多层循环，对coder们来说这无疑是讨厌的，尤其是要跳出某一层循环就更加崩溃了，面对那么多的do-end do，怎么跳？到底怎么跳？Fortran语法允许在do循环前命名，这样一来就比...

一旦输出结果出现NaN，编译器不会给出任何错误提示，这个时候该如何调试程序呢？点击菜单栏的调试，最后一个为XXXX属性，打开对话框，左侧展开 Fortran分类，找到Floating Point选项...

Fortran中的tiny函数tiny(x)表示查询x的最小正值,x所能表示的最小的数，近似于0。tiny这个函数的参数，只与类型有关。两个数之间的差，不可能比这个数小。任意两个数 a 和 b，如果他们的差小于tiny(x)，则认为这两个数相等。在...

有这样一组数据：15 36 85 96 67 45 14 1428 23 45 56 45 45 34 1242 65 56 67 43 12 45 6073 85 96 93 82 65 16 29 这里想将其写入到一个A(4,8)的数组里面。 首先定义 ......

程序采用Fortran语言编写，在Intel编译器下调试通过。