回 帖 发 新 帖 刷新版面

主题:能否对矩阵进行这样的整体操作

一个矩阵a,如果它的每个元素都小于某个值eps=0.00001就执行某个操作,可以整体操作码?

if (a<eps) then

还需要对a里面的每个元素都操作吗?
if (a(i,j)<eps) then
?

回复列表 (共11个回复)

11 楼

[quote][quote]
用第二种或者forall。[/quote]
implicit none

do i=1,2*m
   do j=1,2*m
 hclc(i,j)=conjg(hl(j,i))
 hcrc(i,j)=conjg(hr(j,i))
   end do
 end do

这样表达有错吗?怎么提示出错:
Warning: This argument's data type is incompatible with this intrinsic procedure; procedure assumed EXTERNAL.   [HL]
 hclc(i,j)=conjg(hl(j,i))

这段话大体上是:这个数据类型与内部过程不一致;程序假定它为外部(变量)。

Error: This name does not have a type, and must have an explicit type.   [CONJG]
 hclc(i,j)=conjg(hl(j,i))
conjg是内部函数,怎么还需要定义类型呢?[em8][em8][em8][em10][em10][em10][em14][em14][em14][em14][/quote]


解决了吗?[em18]

我来回复

您尚未登录,请登录后再回复。点此登录或注册