如何将连续列的数对应的最小值提取出

978142355 · 发表于 2015-11-11 22:56:32

本人刚刚接触Fortran代码，想将图片红框中所有1的那一部分对应E的最小值获得,然后所有2对应E的最小值获得,3对应E的最小值……以此类推，然后将对应最小值那一行输出到文件中，比如：3个红色框的部分，将它们最小值获得后写入文件中变成：
Step number  39  scan point    1  E=  -571.251770020
Step number 5  scan point    2  E=  -571.252807617
Step number 4  scan point    3  E=  -571.253234863
……
这个代码是我自己尝试写的，不过问题还没解决，只解决了红框中那一列中有多少个数。还望各位有什么思路或者代码修改不吝惜赐教，谢谢。

[Fortran] 纯文本查看 复制代码

program takeminenergy
    implicit none
	integer i,j,num,num1,num2,lines,k,num4
    integer,allocatable::x(:),y(:)
    real,allocatable::z(:),minz(:),minumz(:)
    open(12,file='result.txt')!输入文件
	open(13,file='a.txt')   !输出文件
    open(14,file='b.txt')  !输出文件
	num1=1
	num2=0
    write(*,*) 'How many lines:' !读取文件有多少行
    read(*,*) lines
    allocate (x(lines))
    allocate (y(lines))
    allocate (minz(lines))
    allocate (minumz(lines))
    do i=1,lines
        read(12,"(12x,I4,12x,I6)") x(i),y(i)
	    if (y(i)== num1) then 
	        num2=num2+1
	    else if (y(i) /= num1) then
	        y(num1)=num2
	        num1=num1+1
	        num2=1
        end if
    end do
    write(*,*) num1
    
        do j=1,num1 
            write(13,"(I6)") y(j) !获取数字相同那一列有多少个，比如全是数字1的共有39个
        end do
        close(13)
        close(12)
        
        open(13,file='a.txt')
        open(12,file='result.txt')
        open(15,file='c.txt') !输出文件
        
        do j=1,num1 
            read(13,"(I6)") num4
            write(*,*) num4
            do k=1,num4
             read(12,"(38x,f16.9)") z(k)
             write(15,*)z(k)
            end do
            
            minz=minval(z)
            write(14,"(f16.9)") minz
        deallocate(z)
        end do
        
end

楚香饭 · 发表于 2015-11-12 09:33:14

本帖最后由楚香饭于 2015-11-12 09:44 编辑

你想得太复杂了，其实问题可以很简单。

[Fortran] 纯文本查看 复制代码

Program www_fcode_cn
  Implicit None
  real(kind=kind(0.0d0)) :: rE , rMinE
  integer :: ierr , iRepeat , iPoint , pointerlast , iMinRepeat
  character(len=12) :: c(5)
  Open( 12 , File = 'result.txt' )
  Open( 13 , File = 'out.txt' )
  pointerlast = 1
  rMinE  = 1.0e30
  Do
    Read( 12 , * , ioStat=ierr ) c(1) , c(2) , iRepeat , c(3) , c(4) , iPoint , c(5) , rE
    if ( ( ierr== 0 ) .and. (iPoint == pointerlast) ) then !// 同一个Point
      if( rE <= rMinE ) then
        rMinE      = rE
        iMinRepeat = iRepeat
      end if
    else !//不同Point
      write(13,'(a,1x,a,1x,i4,2x,a,1x,a,i6,2x,a,f16.9)') trim(c(1)) , trim(c(2)) , iMinRepeat , &
        trim(c(3)) , trim(c(4)) , pointerlast , trim(c(5)) , rMinE
      pointerlast = iPoint
      iMinRepeat  = iRepeat
      rMinE = rE
    end if    
    if ( ierr /= 0 ) Exit
  End Do
  Close( 12 )
  Close( 13 )
End Program www_fcode_cn

978142355 · 发表于 2015-11-12 10:08:28

群主太无敌了，唉，我该如何才能达到你的这种境界。目前我也是一边看彭老师的书，一边练习编程，也关注论坛和qq群，但是到真正自己编程的时候发现思想总是很复杂，想问问群主，为什么你的编程思想可以这么简单，我看这代码里也没有什么特别特殊的东西。

楚香饭 · 发表于 2015-11-12 10:16:15

这个没有别的方法，多写多练吧。

很多时候同一个问题是可能有不同的思路的，你把这个代码看懂了，就理解他的思路了。以后遇到同类问题也不妨换个思路看看。

978142355 · 发表于 2015-11-12 11:32:29

懂了，多谢多练，思想自然就有了，再次谢谢群主。

[求助] 如何将连续列的数对应的最小值提取出

规矩勋章

新人勋章

爱心勋章

水王勋章

元老勋章

热心勋章