매트릭스의 치수를 찾는 방법은 무엇입니까?

Question

나는 성격과 실수가 모두 포함 된 매트릭스를하고 난이 행렬을 읽고 (자체 크기를 찾습니다) 프로그램을 원한다. 내 행렬의 크기에 대한 내가 가진이 코드 (3,2)와

a b 13 15.5 13.2 
c d 16 16.75 19 
e f 19.2 12.2 18.2 

:

! A fortran95 program for G95 
Program Project2nd 
implicit none 
character(len=40), allocatable :: a(:,:) 
integer i,j,k,n,m,l,st 
character(len=40) d 
n=0; m=1; j=1; 

open(10,file=& 
'/Users/dariakowsari/Documents/Physics/Programming/Fortran95-Projects/Project2nd/input.txt', & 
    IOstat=st) 

do while (st == 0) 
    read(10,*,IOstat=st) d 
    n=n+1 
end do 

st=0 
do j=1,m 
    do while (st == 0) 
    allocate(a(1,m)) 
     read(10,*,IOstat=st) (a(1,j),j=1,m) 
      m=m+1 
     deallocate(a) 
end do 

print*, n,m 


end 

여기 내 매트릭스입니다 : 여기 내 코드입니다.

Answer 1

은 나를 위해 컴파일되지 않지만 몇 가지 변경 후 내가 당신에게 유사한 결과를 얻을 수 있었다 의미 귀하의 예제 코드에서 몇 가지 오류가 있습니다. * 업데이트 내 다른 (지금 삭제) 대답 코멘트에 @francescalus에서 언급 한 바와 같이, 그 접근 방식은 정의되지 않은 동작을 포함하고 같은 적절한 해결책이 아니다. 이것은 파일보다 더 많은 엘리먼트를 읽는 것으로부터 발생했다.)

다음은이 정의되지 않은 동작을하지 말아주세요 다른 방법,하지만 아마도 매우 비효율적이다. 당신이 줄 수를 카운트했다으로

Program Project2nd 
    implicit none 
    character(len=40), allocatable :: a(:) 
    integer, allocatable :: ind(:) 
    integer, parameter :: maxElements = 100 
    integer i,j,n,m,st 
    character(len=40) d 
    n=0; 

    open(10,file='mat.txt',IOstat=st) 
    !Find number of lines 
    do while (st == 0) 
    read(10,*,IOstat=st) d 
    if(st ==0) n=n+1 
    end do 
    !Move back to the start of the file 
    rewind(10) 

    !Read all of the data 
    do m=n,maxElements,n 
    allocate(a(m)) 
    read(10,*,IOstat=st) a 
    deallocate(a) 
    rewind(10) 
    if(st.ne.0) exit 
    enddo 
    m = m -n !Need to roll back m by one iteration to get the last which worked. 
    if(mod(m,n).ne.0) then 
    print*,"Error: Number of elements not divisible by number of rows." 
    stop 
    endif 
    !Number of columns = n_elements/nrow 
    m=m/n 
    print*, n,m 
end Program Project2nd 

는 기본적으로이 그러나 읽기 (즉, 일 == 0) 성공시에만 N을 증가 할주의, 동일한 코드를 사용합니다. st가 0이 아니 자마자 블록을 종료하지 않는다는 것을 유의하십시오. 일단 블록이 끝나면 끝이납니다. 그 다음에는 파일을 되 감아 파일 시작 부분에서 다음 읽기가 시작되도록해야합니다. 당신은 당신이 정말로 다음

여기

st = 0 ; m = n 
    do while (st==0) 
    allocate(a(m)) 
    read(10,*,IOstat=st) a 
    deallocate(a) 
    rewind(10) 
    if(st.ne.0) then 
     m = m - n !Go back to value of m that worked    
     exit 
    endif 
    m=m+n 
    enddo 

Answer 2

승 수행하는 방법처럼 뭔가 두 번째 do 루프를 교체이를 방지하려는 경우 오히려 maxElement을 지정할 필요가 없습니다 거라고 언급 이전 코멘트에서

/o 되감기.

implicit none 
    character(len=100) wholeline 
    character(len=20), allocatable :: c(:) 
    integer iline,io,ni,nums 
    open(20,file='testin.dat') 
    iline=0 
    do while(.true.) 
    read(20,'(a)',iostat=io)wholeline 
    if(io.ne.0)exit 
    iline=iline+1 
    ni=lineitems(wholeline) 
    allocate(c(ni)) 
    read(wholeline,*)c 
    nums=ctnums(c) 
    write(*,*)'line',iline,' contains ',ni,'items',nums, 
$  'are numbers' 
    deallocate(c) 
    enddo 
    write(*,*)'total lines is ',iline 
    contains 

    integer function ctnums(c) 
    ! count the number of items in a character array that are numbers 
    ! this is a template, 
    ! obviously you could assign the numbers to a real array here 
    character(len=*), allocatable :: c(:)  
    real f 
    integer i,io 
    ctnums=0 
    do i = 1,size(c) 
    read(c(i),*,iostat=io)f 
    if(io.eq.0)ctnums=ctnums+1 
    enddo 
    end function 

    integer function lineitems(line) 
    ! count the number of items in a space delimited string 
    integer,parameter ::maxitems=100 
    character(len=*) line 
    character(len=80) :: c(maxitems) 
    integer iline,io 
    lineitems=0 
    do iline=1,maxitems 
    read(line,*,iostat=io)c(:iline) 
    if(io.ne.0)return 
    lineitems=iline 
    enddo 
    if(lineitems.eq.maxitems)write(*,*)'warning maxitems reached' 
    end function 
    end 

출력

line 1 contains 5 items 3 are numbers 
line 2 contains 5 items 3 are numbers 
total lines is 2