Mercurial > op > rk_sakamoti

diff rkmain.f90 @ 3:2347cff808b0

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
chage a parameter of rkmod::calc
author "uncorrelated zombie" <uncorrelated@yahoo.co.jp>
date Sun, 21 Jul 2024 04:11:29 +0900
parents 206e90e91191
children 0ca8347e4c47
line wrap: on
line diff
--- a/rkmain.f90	Sun Jul 21 03:56:49 2024 +0900
+++ b/rkmain.f90	Sun Jul 21 04:11:29 2024 +0900
@@ -1,13 +1,12 @@
 program rkmain
   use rkmod
   implicit none
-  real(kind=dpkind),dimension(3) :: x
-  real(kind=dpkind) :: max_t
-  real(kind=dpkind),allocatable,dimension(:,:) :: r
+  real(kind=dpkind),dimension(3) :: x = [0.5, 0.5, 0.5]
+  integer :: max_i = 5000
+  real(kind=dpkind), allocatable, dimension(:, :) :: r
 
-  x = [0.5d0,0.5d0,0.5d0]
-  max_t = 5d1
-  call calc(x, max_t, r)
+  allocate(r(max_i, 4))
+  call calc(x, max_i, r)
   print *, r
 end program