Предположим, у меня есть матрицы A и B, обе имеют размеры n x n. Теперь я хочу рассматривать каждый столбец A как матрицу столбцов ColA(i) [i = 1, n], а каждую строку матрицы B как матрицу строк RowB(j) [j = 1, n].
Теперь можно использовать функцию matmul для построения новой матрицы с именем C, элементы которой находятся путем умножения ColA и RowB для всех возможных значений i и j? Что-то вроде
do i = 1, n
do j= 1, n
C(i,j) = matmul( ColA(i), RowB(j) )
end do
end do
Что ж, я вижу три возможных ответа на ваш вопрос, но не уверен, какой из них вы имеете в виду ...
i из A и строки j из B, которое будет равно умножению матрицы на матрицу B*A:C = matmul( B,A )
C=colA * rowB^T, но тогда каждая комбинация i,j приведет к отдельной матрице C.C = matmul( A(:,i), B(j,:) )
C_ij = A_ij * B_ij. В Фортране это довольно просто:C = A*B
Выбирайте!
Ваш код не имеет смысла. ColA(i) не является вектором, это скаляр, вы не можете передавать скаляры для умножения матриц. Точно так же C(i,j) не матрица, это скаляр. C - это матрица.
Или, если вы имеете в виду ColA как вектор векторов, это невозможно в Фортране. C тогда должно быть матрицей матриц, что тоже невозможно. Это возможно с производными типами, но вам придется кодировать это самостоятельно.
ColA(i) не может быть способом определить вектор-столбец в fortran, но это действительно соответствует моему намерению. Итак, позвольте мне перефразировать, я просто хочу умножить каждый столбец A на каждую строку B. Я думаю, что это матричное умножение законно, и мне было интересно, могу ли я использовать matmul процедуру для этого. Если да, то как? [Если я вам все еще непонятен, дайте мне знать.]
- person hbaromega; 21.05.2014
