FORTRAN 90 help- matrix-mul-add- input from user defined file name

[muzzy29]

Hi everyone

I am completely new to the Fortran. I have to make a program in fortran 90 to either add or multiply the matrix on user demand and ask user to enter the file name inorder to read the desired matrices.

I came up with the the following code which does not read the program and produces runtime error before the open statement.
And i am going to make subroutines for each case(add and multiply)
Can anyone help me identify the right way of making the source code in order to read the user defined input from input file and put the matrix on tothe screen as it is defined in the input file.
Following the formatting to input matrix

rowsa, colsa
2, 5, 6
8, 9, 7

Also fortran only does 10 * 10 matrix so how do i make sure it reads every entry of the program properly.

! Ask user for the desired operation
WRITE(*,*) "Enter the desired operation to perform as sum or multiply "
Read (*,*) operation

! Ask user to open file for the first matrix
WRITE(*,*) ""
WRITE(*,*) ""

WRITE(*,*) "Enter file name for the first matrix"
Read (*,*) inputA

OPEN (UNIT=20, file='inputA', STATUS='OLD')


! Read dimension of matrix A
READ (20,*) rowsA,colsA

! Read each matrix(array) elements
Do i= 1,rowsA

Do j=1,colsA

READ (11,*) A
WRITE (*,*) "The first matrix is ",A

End do

End do

close (20)


! Ask user to open file for the Second matrix
WRITE(*,*) ""
WRITE(*,*) ""

WRITE(*,*) "Enter file name for the second matrix"
Read (*,*) inputB

OPEN (UNIT=30, file='inputB', STATUS='OLD')


! Read dimension of matrix A
READ (30,*) rowsB,colsB

! Read each matrix(array) elements

Do i= 1,rowsB

Do j=1,colsB

READ (30,*) B
WRITE (*,*) "The second matrix is ",B

End do

End do

close (30)

END PROGRAM MAT_add_multiply

Thanks very much.

Muzzy

 

[gsal]

program mat_add_multiply
integer rowsA, colsA, rowsB, colsB, i, j
real A(100,100)
real B(100,100)
real C(100,100)
character operation, inputA*30, inputB*30

! Ask user for the desired operation
WRITE(*,*) "Enter the desired operation to perform as sum or multiply "
Read (*,*) operation

! Ask user to open file for the first matrix
WRITE(*,*) ""
WRITE(*,*) ""
WRITE(*,*) "Enter file name for the first matrix"
Read (*,*) inputA
OPEN (UNIT=20, file=inputA, STATUS='OLD') 

! Read dimension of matrix A
READ (20,*) rowsA,colsA
! Read each matrix(array) elements 
WRITE (*,*) "The first matrix is "
Do i= 1,rowsA
   READ (20,*) (A(i,j),j=1,colsA)
   write(*,"(100F7.2,$)") (A(i,j),j=1,colsA)
End do
close (20)

! Ask user to open file for the Second matrix
WRITE(*,*) ""
WRITE(*,*) ""
WRITE(*,*) "Enter file name for the second matrix"
Read (*,*) inputB
OPEN (UNIT=30, file=inputB, STATUS='OLD') 

! Read dimension of matrix A
READ (30,*) rowsB,colsB
! Read each matrix(array) elements 
WRITE (*,*) "The second matrix is "
Do i= 1,rowsB
   READ (30,*) (B(i,j),j=1,colsB)
   write(*,"(100F7.2,$)") (B(i,j),j=1,colsB)
End do
close (30)

WRITE(*,*) ""
WRITE(*,*) ""
WRITE(*,*) "The result is:"
WRITE(*,*) ""
If (operation == '*') then
   if (colsA.NE.rowsB) then 
      write(*,*) "Sorry, A*B cannot be performed;"
      write(*,*) "Please check matrices dimensions"
   else
      C = matmul(A(1:rowsA,1:colsA),B(1:rowsB,1:colsB))
      Do i= 1,rowsA
         WRITE(*,'(100F7.2)') (C(i,j),j=1,colsB)
      End do   
   end if
else if (operation == '+') then
   if (rowsA.NE.rowsB.or.colsA.NE.colsB) then  
      write(*,*) "Sorry, A+B cannot be performed;"
      write(*,*) "Please check matrices dimensions"
   else
      C = A(1:rowsA,1:colsA)+B(1:rowsB,1:colsB)
      Do i= 1,rowsA
         WRITE(*,'(100F7.2)') (C(i,j),j=1,colsB)
      End do   
   end if
end if

END PROGRAM MAT_add_multiply

 

[muzzy29]

THank you SO much it was really helpful :)

Just had a issue with multiplication part that seems like it just multiplies 1st row and 2nd row of A with first column of B and writing out third column of A as third column of result (matrix C) and 2nd column of matrix B in the 2nd column of result (matrix C). by adjusting the dimensions.

I am not completely aware of the loop that you coded up so I am not sure how would i fix that, if you can explain the formatting part while writing the result matrix and the looping part, I Will appreciate that.

Thank you again for your generous help :)

 

[gsal]

please post again and include (quote) the precise lines you want me to explain...and I'll be happy to do so.

 

[muzzy29]

I basically don't understand the loops which i included below.I do not understand that how fortran reads these loops and i want to write flowchart so if you can tell me how these loops works and the formatting statement that you used in write statement. Also there is some issue with 2nd loop (multiplication loop) and what is matmul? did not I need to make matrix C entries 0 before I add a*b entries in matrix C ?

"Do i= 1,rowsB
READ (30,*) (B(i,j),j=1,colsB)
write(*,"(100F7.2,$)") (B(i,j),j=1,colsB)
End do"



"C = matmul(A(1:rowsA,1:colsA),B(1:rowsB,1:colsB))
Do i= 1,rowsA
WRITE(*,'(100F7.2)') (C(i,j),j=1,colsB)
End do"




"C = A(1:rowsA,1:colsA)+B(1:rowsB,1:colsB)
Do i= 1,rowsA
WRITE(*,'(100F7.2)') (C(i,j),j=1,colsB)
End do "

Thank you

 

[gsal]

The DO-loops that you refer to take advantage of a Fortran feature called "implicit do-loop"...this do-loop does not need "do"/"end-do" enclosing statements; instead, it simply happens within a single set of parenthesis.

The construct "(B(i,j),j=1,colsB)" leaves i alone, but varies j from 1 to colsB.

The line
READ (30,*) (B(i,j),j=1,colsB)
reads an entire row of matrix B; more precisely, it reads the i-th row as the value for i remains constant but j is varied from 1 to colsB.

The loop
Do i= 1,rowsB
READ (30,*) (B(i,j),j=1,colsB)
End do
would seem equivalent to a double loop like this:
Do i= 1,rowsB
Do j= 1,colsB
READ (30,*) B(i,j)
End do
End do
but will not work as desired due to the way you want to format your files.

When reading sequential file like you are doing, the READ() command retrieves the next line in the file and reads from it as many values as requested; and so, the explicit double-do-loop above would only read the first value of every row and soon run out of file...for it to work, you would need to write your input file with one value per line, which looses it visual correspondance to a nice looking matrix.

The WRITE statements have also been written to take advantage of the implicit do-loop feature...because we do not know how many values per line we need to write, I have simply put a multiplier of 100 for the F7.2 specification...it will only use as many as needed.

The multiplication and addition operations as shown are taking advantage of Fortran-90 intrinsic features...just like the scalar addition c=a+b (c=3+7, for example) works in any other language, the array addition C=A+B works in Fortran-90 (a-la Matlab); the same goes for the intrinsic function 'matmul' (matrix multiplication)..it just does the double-loop for you, no need to do anything else.

The reason why I had to specify the dimensions of the matrices during addition and multiplication is because the matrices were declared larger than necessary since we did not know in advnace the necessary dimensions...it is possible to delay the creation of such matrices with "allocatable" but I did not want to introduce that at this time.

Don't be shy about "asking" questions to the comipler...by that, I mean, experiment with simple programs that test one feature at a time, compile then, and see how it behaves.

 

[YunzhE]

hey i tried the code when
"C = A(1:rowsA,1:colsA)+B(1:rowsB,1:colsB)
Do i= 1,rowsA
WRITE(*,'(100F7.2)') (C(i,j),j=1,colsB)
End do "

it count correctly in order when matrix 3,3 add together

but i can't find the correct answer for the multiplication with this, same matrix 3,3 multiply together

"C = matmul(A(1:rowsA,1:colsA),B(1:rowsB,1:colsB))
Do i= 1,rowsA
WRITE(*,'(100F7.2)') (C(i,j),j=1,colsB)
End do"

 

[gsal]

can you show me what matrices you are using and what results you are expecting and what results you are getting?

did you figure how to enter input to the program? the first answer is either the character '+' or the character '*', right? Then, the file names of the matrices...

Please list the files you are using, etc.

 

[YunzhE]

example:
file 1
3,3
2,1,2
3,1,0
1,1,1
file 2
3,3
2,2,2
2,2,2
2,2,2
yes i can do the addition part
the answer will be
4,3,4
5,3,2
3,3,3

but when the multiplication part it shows the answer is
4,2,4
6,2,0
2,2,2

which the answer should be
10,10,10
8,8,8
6,6,6

 

[gsal]

hhhmmm...my version of the program output the correct answer...