program main
implicit none
call input
call table
print *, 'The program has completed successfully'
stop
end
subroutine input
implicit none
integer ui, IERR
real X_min,X_max,dx,Y_min,Y_max,dy
real temp_x, temp_y
common /arguments/ X_min,X_max,dx,Y_min,Y_max,dy
parameter(ui=111)
open(ui,FILE='input.txt', status='old', IOSTAT=IERR, err=404)
read (ui,*) X_min,X_max,dx,Y_min,Y_max,dy
close(ui)
c Checking input values
if(dx .LE. 0 .OR. dy .LE. 0) then
write(6,*) 'Input values error: Wrong x/y step'
stop
elseif(X_min .GT. X_max .OR. Y_min .GT. Y_max) then
write(6,*) 'Input values error: Wrong min/max arguments'
stop
endif
!Checking if step is too low
!Example: x: 0.01; step: 0.000000001; - error
temp_x = X_min + dx
temp_y = Y_min + dy
if(temp_x .EQ. X_min .OR. temp_y .EQ. Y_min) then
write(6,*) 'Input values error: x/y step is too low'
stop
endif
return
c Input error
404 if(IERR .EQ. 29) then !FOR$IOS_FILNOTFOU
write(6,*) 'Error: File does not exist'
elseif(IERR .EQ. 43) then !FOR$IOS_FILNAMSPE
write(6,*) 'Error: File name specification error'
elseif(IERR .EQ. 9) then !FOR$IOS_PERACCFIL
write(6,*) 'Error: Permission to access file denied'
else
write(6,*) 'Error: Undefined error, code = ', IERR
endif
print *, 'Input file error', IERR
stop
end
subroutine cotg(value,eps_null,result,error_code,error_str)
implicit none
real result, value, eps_null
real temp
integer error_code
character*3 error_str
error_code = 0
error_str = ' '
temp = sin(value)
if(temp .EQ. 0) then
error_code = 1 !divided by zero
error_str = 'NaN'
return
elseif(abs(temp) .LE. eps_null) then
error_code = -1 !infinity
error_str = 'INF'
return
endif
result = cos(value)/temp
end
!proverit sverhu ili snizy razmeshenie functions
!0.00001 = a
!0.0001 = da
!log10(0.00009)
real function epsilon_power(a,da,epsilon_da)
implicit none
real a, da, epsilon_da
if(abs(a) .LT. 1) then
if(abs(a) .GT. epsilon_da) then
epsilon_power = aint(log10(abs(a))) - 4
else
epsilon_power = aint(log10(abs(a-da))) - 4
endif
else
epsilon_power = aint(log10(abs(a))) - 3
endif
end
subroutine find_steps(array,min,max,dx,n,eps)
implicit none
integer n
real array(*)
real epsilon_power
real temp,privious,min,max,dx
real power_temp, eps
n = 0
temp = min
!! 20 eto N dlya array
10 if(n .LT. 20) then
n = n + 1
!Check max
if(temp .GE. max .OR. abs(temp-max) .LE. eps) then
array(n) = max
return
endif
!End check max
if(abs(temp) .LT. eps) temp = 0
array(n) = temp
power_temp = epsilon_power(temp,dx,eps)
privious = anint(temp*10**(-power_temp))*10**power_temp
temp = temp + dx
!Skip invisible steps
20 power_temp = epsilon_power(temp,dx,eps)
if(abs(temp-privious) .LT. 0.9*10**power_temp) then
temp = temp + dx
goto 20
endif
!End skip invisible steps
temp = anint(temp*10**(-power_temp))*10**power_temp
goto 10
endif
min = temp
return
end
subroutine table
implicit none
integer uo, IERR, error_code
integer n,m,i,j
real array_x(20), array_y(20)
character*3 error_string
real X_min,X_max,dx,Y_min,Y_max,dY
real degreeToRadian, result
real eps_null, eps_dx, eps_dy
real temp
common /arguments/ X_min,X_max,dx,Y_min,Y_max,dY
common /const/ degreeToRadian
parameter(uo=112)
degreeToRadian = 3.1415926/180
!!!!!!!!!!!!!!!!!!!!!!!
! TASK: CHECK ERR=404
!!!!!!!!!!!1
open(uo,FILE='output.txt',status='unknown',IOSTAT=IERR,err=405)
eps_dx = 10**aint(log10(abs(dx)) - 3)
eps_dy = 10**aint(log10(abs(dy)) - 3)
if(eps_dx .LT. eps_dy) then
eps_null = 10*eps_dx*degreeToRadian !!!!!! 10 DELETE
else
eps_null = 10*eps_dy*degreeToRadian
endif
!x do while
10 continue
call find_steps(array_x,X_min,X_max,dx,n,eps_dx)
temp = Y_min
!y do while
20 continue
call find_steps(array_y,temp,Y_max,dy,m,eps_dy)
do i = 1, n + 1
write(uo,1000) !print ---
enddo
write(uo,1020)
write(uo,1030) 'y/x'
write(uo,1010) (array_x(i), i=1,n)
write(uo,1020)
do i = 1, n + 1
write(uo,1000) !print ---
enddo
write(uo,1020)
do i = 1,m
write(uo,1010) array_y(i)
do j = 1,n
call cotg((array_y(i)+array_x(j))*degreeToRadian,
* eps_null, result, error_code, error_string)
if(error_code .EQ. 0) then
write(uo,1010) result
else
write(uo,1030) error_string
endif
enddo
write(uo,1020)
enddo
do i = 1, n + 1
write(uo,1000) !print ---
enddo
write(uo,1020)
if(array_y(m) .NE. Y_max) goto 20 !end do while
if(array_x(n) .NE. X_max) goto 10 !end do while
close(uo)
return
c Formats
1000 format(15('-')\)
1010 format('| 'E11.4\' |')
1020 format(/\)
1030 format('| 'A\' |')
c Output file errors
405 if(IERR .EQ. 43) then !FOR$IOS_FILNAMSPE
write(6,*) 'Output error: File specification error'
elseif(IERR .EQ. 10) then !FOR$IOS_CANOVEEXI
write(6,*) 'Output error: Cannot overwrite existing file'
elseif(IERR .EQ. 21) then !FOR$IOS_DUPFILSPE
write(6,*) 'Output error: Duplicate file specifications'
else
write(6,*) 'Output error: Undefined error, code = ', IERR
endif
print *, 'Output file error ', IERR
stop
end
