File: parameter_11.f90

package info (click to toggle)
lfortran 0.45.0-1
  • links: PTS, VCS
  • area: main
  • in suites: sid, trixie
  • size: 46,332 kB
  • sloc: cpp: 137,068; f90: 51,260; python: 6,444; ansic: 4,277; yacc: 2,285; fortran: 806; sh: 524; makefile: 30; javascript: 15
file content (22 lines) | stat: -rw-r--r-- 478 bytes parent folder | download 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
 
program parameter_11
integer, parameter :: n = 3
real:: v(n), a(n,n), b(n,n)
integer :: i,j

v = [(i,i=1,n)]
a = diag_rsp(v)
b = reshape([((merge(i,0,i==j), i=1,n), j=1,n)], [n,n])
print *, b
print *, "all(a==b)", all(a==b)
if (.not. all(a==b)) error stop
contains
function diag_rsp(v) result(res)
    real, intent(in) :: v(:)
    real :: res(size(v),size(v))
    integer :: i
    res = 0
    do i = 1, size(v)
      res(i,i) = v(i)
    end do
end function diag_rsp
end program