또 포트란 질문 드립니다.. 벡터 계산 관련

예를 들어 20x1 크기의 벡터 fn0, fn1 있다고 치면 ..

이걸 계산을 통해서 20x20 인 벡터를 fn2 을 만드려고 하는데요.

=> fn2 = fn1 - fn0(특정위치의 값) 이런식으로.. fn1 벡터에서 fn0의 특정 값을 뺄셈해서..

써보자면..

fn1(20x1크기) - fn0(20x1 중 n번째 값) = fn2(20x첫번째)

1                                   -1                  0
2                                   -1                  1
3                                   -1                  2
.                                       .    =           .
.                                       .                 .
.                                       .                 .
20                                  -1                 19

이런식으로요.

fn2의 20x1(첫번째)의 값으로 지정하고 그 다음은 20x2, 20x3, 20x4 . . . 20x20 이런식으로

fn2를 만드려고 합니다.

매트랩에서는 for loop 를 사용해서

for i = 1:20
for j = 1:20
   fn2(i,j) = fn1 - fn0(j)
end for
end for

이런식으로 하면 fn2가 자동적으로 매 루프마다 20x1, 20x2.... 20x20 와 같이 커졌는데

포트란에서는 이걸 어떻게 하는지 모르겠네요...

도움좀 부탁드립니다.

잡담 | 1507명이 읽었어요. 3.135.221.0

비공개 손님 … 2016-01-22 10:07:21

ㅋㅋㅋㅋ 지난 번 게시물에 댓글에 쓴대로 매트랩 떠나셨으면 행렬 다루기가 X같아지신겁니다. 포트란을 선택하신 모양이네요. 매트랩은 행렬 사이즈가 넘치면 지가 알아서 동적할당을 해서 행렬 사이즈 맞춰주는데 포트란부터는 그런거 기대하지 마세요. 일단 fn2 를 20 x 20으로 선언을 해야될거예요. 짐 제가 핸드폰이고 컴을 할 수 없는 상황이어서 자세한 답변이 어렵네요. 좀 따 글 자세히 읽어보고 댓글 달게요.

비공개 손님 … 2016-01-22 10:17:32

1// 네 ㅎㅎㅎ. 연습겸 이리저리 부딛히면서 시작했습니다 ㅎㅎ. 일단 fn2 를 20x20로 allocate를 해놓은 상태입니다 ㅎㅎ.

비공개 손님 … 2016-01-22 10:43:55

저도 포트란 오랜만에 보는지라 잘 모르겠네요.
일단 원하시는 코드를 써 보자면

program programName

real, dimension(20):: fn0, fn1
real, dimension(20,20):: fn2
! 위엣 줄은 real, dimension(:,:) 로 동적할당해서 매트랩처럼 사용은 할 수 있는데
! 매트랩처럼 null 값을 0으로 자동 바꿈 안 해줄 겁니다.
integer::i,j

do i=0:19 ! 이제 인덱스는 0부터 시작합니다.
do j=0:19
fn2(i,j) = fn1(j) - fn0(i)
! 본문의 매트랩 수식이 이렇게 되어야 하는거 같은데...요? 벡터 - 상수 계산이 포트란에서도 매트랩처럼
! 벡터 - 상수*(1,1,....,1)로 자동 계산되기는 하는데, 이제부터 인덱스 따져 주는게 좋습니다.
! 벡터 - 상수 같은거 편하다고 막 쓰시다가 나중에 코드 복잡해지면 골 때립니다.
end do
end do

print *, fn2
! 이부분 잘 모르겠네요. 암튼 fn2 값을 출력해서 확인해보세요.

end program programName

대충 이 정도가 되지 않을까요? 아 그리고, 앞으로 질문하실 때, 작성하신 코드 복붙해주시면 좋습니다.
저도 심심한데 로렌츠 방정식 돌려볼까 싶네요.

x' = s(y-x)
y' = rx-y-xz
z' = bz+xy

비공개 손님 … 2016-01-22 10:55:08

코딩 작업중인게 무지 더러워서... 잠시만요 정리좀 하고 바로 올리겠습니다.

비공개 손님 … 2016-01-22 10:58:48

아...아니.. 지금 당장 올리라는게 아니고, 앞으로 질문하실 때.... 문제가 되는 코드와 컴파일시 또는 실행시 나타나는 에러 메시지를 붙여주시면 되요... 이게 보통 프로그래밍 질문하는 관례입니다.

비공개 손님 … 2016-01-22 11:00:22

program HW3
implicit none

real, allocatable :: wavleng(:), anfq(:), img(:), fn0(:), fn1(:), fn2(:,:), fn3(:)
integer :: i, j, s, length
real :: dummy
real, parameter :: c = 2.998e+8

! ===== Calculating array size by reading data from the source file ===== !
open (unit=1, file='HWdata.dat')

length = 0
do
read(1,*,IOSTAT=s) dummy

if(s /= 0) exit

length = length + 1
end do

close(unit=1, status='keep')

! ===== Allocating array size ===== !
open(unit=1, file='HWdata.dat')

allocate(wavleng(length))
allocate(img(length))
allocate(fn0(length))
allocate(fn1(length-1))
allocate(fn2(length-1,length-1))
allocate(fn3(length-1))

! ===== Reading data from source ===== !
do i = 1,length
read(1,*) wavleng(i), dummy, img(i)
end do

! ===== Converting wavelength to angular frequency ===== !
anfq = 2*4*atan(1.0)*(c/wavleng)

! ===== Calculating function inside the integral by parts ===== !

! ===== w'*n_i(w') part ===== !
fn0 = anfq * img

! ===== Assigning dummy matrix for w ===== !
fn2 = 1

! ================================= !
do i = 1,(length-1)
fn1(i) =((anfq(i+1) + anfq(i))/2)**2
end do

! ====== Matrix calculation ===== !
do i = 1,(length-1)
do j = 1,(length-1)
fn2(i,j) = fn1 - anfq(j)
end do
end do

! ===== Writing results to a asdf.dat file ===== !
open(unit=2, file='asdf.dat')
! do i = 1,(length-1)
! do j = 1,(length-1)
write(2, fmt='(i4,$)') fn2
! end do
! end do

close(unit=2, status='keep')

deallocate(wavleng)
deallocate(img)
deallocate(anfq)
deallocate(fn0)
deallocate(fn1)
deallocate(fn2)
deallocate(fn3)

close(unit=1, status='keep')
end program HW1

이런식으로 했답니다...

비공개 손님 … 2016-01-22 11:05:11

질문있습니다.

do i=0:19 ! 이제 인덱스는 0부터 시작합니다.
do j=0:19
fn2(i,j) = fn1(j) - fn0(i)
end do
end do

이렇게 써주셨는데요. 그러면 fn1의 첫번째 행의 수랑 fn0의 첫번째 행의 수가 연산 되는거 아닌가요?

저는 fn1 벡터전체의 모든 값에서 fn0의 n번째 값을 빼는 연산을 하려고 하거든요...

비공개 손님 … 2016-01-22 11:16:12

/7 저거 do 문을 일일히 풀어보면
====시작=====================
fn2(0,0) = fn1(0) - fn0(0)
fn2(0,1) = fn1(1) - fn0(0)
fn2(0,2) = fn1(2) - fn0(0)
.
.
.
fn2(0,19) = fn1(19) - fn0(0)
========요기까지가 fn2의 첫번째 행=================

fn2(1,0) = fn1(0) - fn0(1)
fn2(1,1) = fn1(1) - fn0(1)
fn2(1,2) = fn1(2) - fn0(1)
.
.
.
fn2(1,19) = fn1(19) - fn0(1)

========요기까지가 fn2의 두번째 행=================
.
.
.
======== 여기서부터 fn2의 마지막행 ================

fn2(19,0) = fn1(0) - fn0(19)
fn2(19,1) = fn1(1) - fn0(19)
fn2(19,2) = fn1(2) - fn0(19)
.
.
.
fn2(19,19) = fn1(19) - fn0(19)

======== fn2의 마지막행 계산 완료.================

그리고 일단

do i = 1:20

요렇게 i 가 1부터 시작하는거 같은데 i가 0부터 시작해야 맞을 겁니다. i를 1부터 시작하게 할려면 배열 선언할 때,

real, dimention(1:20,1:20)::fn2

이렇게 인덱스를 1에서 20까지 쓰겠다고 알려줘야 합니다. 아니면 기본적으로 인덱스는 0부터 시작합니다. 수치해석 프로그래밍은 인덱스와의 전쟁이기도 하죠.

비공개 손님 … 2016-01-22 11:17:57

아 죄송. do i=0:19 라고 썼네요. 맷랩이랑 햇깔... do i=1,19 로 읽어주세요...ㅜ.ㅜ

비공개 손님 … 2016-01-22 11:21:50

아하!! 제가 잘못 이해하고 있었네요 ㅎㅎㅎ 감사합니다.

그리고 인덱스 바꾸겠습니다 ㅎ, 매트랩에서는... 1부터 인덱싱이 되었어서 버릇되었나보네요..

그리고 질문 또 있습니다!

예를들어 20x1 행렬이 있다면 19x1 로 간단히 마지막 행이나 끝 행을 짤라버리는
간단한 방법이 혹시 있을까요?

매트랩에서는 A(:, end-1) 이런식으로 하면 됐었는데 포트란에서는 모르겠어서 지금
시도 한것은 do loop을 n-1 개로 돌려서 또다른 행렬변수를 선언하고 있거든요...

do i = 1,(n-1)
fn3(i) = fn2(i)
end do

이런식으루요...

비공개 손님 … 2016-01-22 11:27:24

아아아아아아아ㅏㅇ
ㅈㅅ
인덱스 1부터 시작합니다.... ㅈㅅㅈㅅㅈㅅㅈㅅ ㅠ.ㅠ

비공개 손님 … 2016-01-22 11:27:29

아.. 그리고 또 질문있는데요! (죄송합니다.... 구글링 해봐도 명쾌한 해결방법이 없어서...)

변수를 print나 write 할때

write(*,*) fn2 이런식으로 하면,,, 그냥 한 개의 행에 값들이 쫘라라라랄라라라라라라라라라

나열되는데.... 행x열 이런식으로 효과적으로 하는 방법 없나요?

제가 시도해봤던것은..

write(*, '(E10.4, E10.4, E10.4 . . . .)') fn 이렇게 하는 방법밖에 모르는데요...

만약 열의 수가 1000개 10000개 이렇게 커지면.... 10000번 E10.4 를 적을수도 없고해서요...

그래서 write(*, *'(E10.4)') (fn, i = 1,1000) 했다가... 다루는 데이터가 1260x1261인데....

이게 1260x1261x1000 이 되어서... 158700000 의 행이 표기되서.. 컴터 멎었답니다.......

비공개 손님 … 2016-01-22 11:35:29

행렬의 인덱스 시작 번호는 포트란 버전에 따라 다른 거 같아요. 포트란 90/95이면 1부터 시작하는게 맞습니다. 기본 시작 인덱스가 0인지 아닌지 한 번 실험해보시고 하세요.

그리고 부분 행렬 추출하는 건 매트랩처럼 하셔도 되는 걸로 알고 있습니다.

http://www.mathcs.emory.edu/~cheung/Courses/561/Syllabus/6-Fortran/array4.html

비공개 손님 … 2016-01-22 11:43:31

답변 감사합니다!

현재 *.f03 써서 2003 사용하고 있습니다!

비공개 손님 … 2016-01-22 11:45:42

'(E10.4, E10.4, E10.4 . . . .)'

요게 문자열이 잖아요. 그니까, 일단 문자열 변수를 선언하고

character ( len =10000) :: c ! 이거도 동적할당 가능할거에요.

do i=1, N ! N 은 (fn의 열의개수)
if (i != N) then
c = c // 'E10.4,' ! c 뒤에 'E10.4,'라는 문자열을 추가해줍니다.
else
c = c // 'E10.4'
end if
end do

이런 식으로 fn의 열의 갯수에 따라 E10.4를 문자열에 붙여주는 코드를 만들어서

write(*, c) ....

이렇게나 아니면

write(*, format=c) ....

이런 식으로 해주시면 될거에요.

제가 짐 포트란이 없어서 확실한 코드를 드릴 수 가 없지만, 요점은 write 함수에 들어가는 문자열을 직접 제조할 수 있다는 겁니다.

비공개 손님 … 2016-01-22 11:57:44

이해했습니다! 정말 감사합니다 ㅎㅎㅎ 조금조금씩 알아가네요 ㅎㅎㅎㅎ

비공개 손님 … 2016-01-22 12:03:59

저... 혹시

if (i != N) 에서 " != " 이 논리 연산자인가요..?

비공개 손님 … 2016-01-22 12:08:01

17/ 죄송합니다.

부정을 나타내는 논리 연산자로 썼는데 틀렸죠...그게.. 맷랩과 자꾸 햇깔려서리...ㅡ.ㅡ

i < N

으로 해주세요.

비공개 손님 … 2016-01-22 12:09:20

감사합니다!!!!

댓글을 작성하실 수 없습니다.
(권한이 없는 회원레벨)

목록으로

TOP arrow_upward