主题:能否对矩阵进行这样的整体操作
tianhy2010
[专家分:60] 发布于 2010-10-11 09:17:00
一个矩阵a,如果它的每个元素都小于某个值eps=0.00001就执行某个操作,可以整体操作码?
if (a<eps) then
还需要对a里面的每个元素都操作吗?
if (a(i,j)<eps) then
?
回复列表 (共11个回复)
11 楼
jstzhurj [专家分:4680] 发布于 2010-10-11 15:41:00
[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]
我来回复