Scribd
Upload
14 views1 page

Fortran Array

This document contains a Fortran program that demonstrates various matrix operations including initialization, addition, multiplication, and applying trigonometric functions. The program declares three 3x3 matrices, assigns initial values, and performs operations like summing, multiplying, and taking the sine of the matrices before outputting the results.

Uploaded by

Marcos A Santos
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as TXT, PDF, TXT or read online on Scribd
14 views1 page

Fortran Array

This document contains a Fortran program that demonstrates various matrix operations including initialization, addition, multiplication, and applying trigonometric functions. The program declares three 3x3 matrices, assigns initial values, and performs operations like summing, multiplying, and taking the sine of the matrices before outputting the results.

Uploaded by

Marcos A Santos
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as TXT, PDF, TXT or read online on Scribd
You are on page 1/ 1

program matrix

implicit none
real, parameter :: PI = 4.0 * atan( 1.0 ) ! The number pi
real, dimension(3,3) :: a, b, c ! Declare size of a, b and c
! real a(3,3), b(3,3), c(3,3) ! Alternative dimension statement
integer i, j ! Counters
character(len=*), parameter :: fmt = "( a, 3(/, 3(1x, f8.3)), / )"
! Format string for output

! Basic initialisation of matrices by assigning all values (very inefficient!)


a(1,1) = 1.0; a(1,2) = 2.0; a(1,3) = 3.0
a(2,1) = 4.0; a(2,2) = 5.0; a(2,3) = 6.0
a(3,1) = 7.0; a(3,2) = 8.0; a(3,3) = 9.0

b(1,1) = 10.0; b(1,2) = 20.0; b(1,3) = 30.0


b(2,1) = 40.0; b(2,2) = 50.0; b(2,3) = 60.0
b(3,1) = 70.0; b(3,2) = 80.0; b(3,3) = 90.0

! Alternative initialisation using data statements (note order)


! data a / 1.0, 4.0, 7.0, 2.0, 5.0, 8.0, 3.0, 6.0, 9.0 /
! data b / 10.0, 40.0, 70.0, 20.0, 50.0, 80.0, 30.0, 60.0, 90.0 /

! Alternative initialisation computing each element of a


! do j = 1, 3
! do i = 1, 3
! a(i,j) = ( i - 1 ) * 3 + j
! end do
! end do
! Then whole-array operation for b
! b = 10.0 * a

! Alternative initialisation using array constructor, implied do-loop and


reshape
! a = reshape( [ ( i, i = 1, 9 ) ], [ 3, 3 ] )
! a = transpose( a ) ! To get intended order!
! b = 10 * a

! Write out matrices (using implied do loops)


write( *, fmt ) "a", ( ( a(i,j), j = 1, 3 ), i = 1, 3 )
write( *, fmt ) "b", ( ( b(i,j), j = 1, 3 ), i = 1, 3 )

! Matrix sum
c = a + b
write( *, fmt ) "a+b", ( ( c(i,j), j = 1, 3 ), i = 1, 3 )

! "element-by-element" multiplication
c = a * b
write( *, fmt ) "a*b", ( ( c(i,j), j = 1, 3 ), i = 1, 3 )

! "Proper" matrix multiplication


c = matmul( a, b )
write( *, fmt ) "matmul(a,b)", ( ( c(i,j), j = 1, 3 ), i = 1, 3 )

! Some operation applied to all elements of a matrix


c = sin( b * PI / 180.0 )
write( *, fmt ) "sin(b)", ( ( c(i,j), j = 1, 3 ), i = 1, 3 )

end program matrix

You might also like