Question

求fortran大神来写个程序

已知一个1500*2000大小的矩阵 矩阵元素只有一个值为1 其余均为0
如何利用尽可能短的时间找到这个值为1的元素的坐标
http://pan.baidu.com/s/1pJM3AEf这个是那个矩阵的文件（格式为dat）

Answer 1

用maxloc命令，作用是找出数组中最大值的坐标。这里也就是1的坐标了。当数组是一维的时候，返回一个整数，当数组是n维时，返回一个n维数组，就是坐标值。

追问

回去试试

追答

举个简单的例子
program main
implicit none
real::a(3,2)
data a /1,2,6,4,11,90/
write(*,*)maxloc(a)
end
最后会输出3 ， 2。表示a(3,2)=90是最大的元素。
你这里只需要把dat文件导入a(1500,2000)的二维数组中即可。

追问

你看一下我这个代码，为何最后输出的位置始终是1 1呢，很明显那个非零元素不是在1 1那个位置
http://pan.baidu.com/s/1pJyNe8V 这个压缩包里面是我的f90文件和矩阵数据numbers.dat

追答

我试了下，也不知道问题出在什么地方，刚开始学习fortran，帮不了你什么忙，望见谅。如果找出了问题所在，有劳知会一声，共同学习！

在学习过程中，发现如果所有元素都是一样大，会输出1，1。不过这里问题应该不是这个。我查了一下书，书上有句话：“使用数组时超出范围是很危险的，绝对要避免。在编译过程中，通常都不会检查数组是否超过范围，而且很多情况下根本也无法发现这种错误。这个责任都是由程序员自行负担。”

追问

好的，我再研究一下

问题出在了read那里，应该用2000I1这个格式输入，不能用星号默认