请问双双精度变量赋值时后缀加什么？

vvt

是的，就像这样
Real , parameter :: QP = Selected_real_kind( 18 )
Real(Kind=QP) :: a = 3.14_QP

vvt

的确是笔误，应该是
Integer , parameter :: QP = Selected_real_kind( 18 )
Real(Kind=QP) :: a = 3.14_QP

不写 dsqrt ，直接写 sqrt 即可，它适用于任何kind值。
详解 http://whyd.w.fcode.cn

vvt

向量，矩阵，是数学概念。
编程里只有一维数组，二维数组的概念。

通常我们用一维数组表示向量，二维数组表示矩阵。（但仅仅是表示，而非相等的概念）
其实 1*N 的二维数组，也可以表示向量。
matmul，必须对两个“二维数组”进行，并且结果也是二维数组（哪怕是1*1的二维数组呢）。显然，a(:,1) 是一维数组。（而 a(:,1:1）却是二维数组）

因此，你的选择可以是：
1. 改用 integer c(1,1) ; c = matmul( b , a(:,1:1) )
2. 改用 c = dot_product( b(1,:) , a(:,1) ) 此时是应用内积的算法