GravityField/.gitignore

@@ -13,7 +13,6 @@ WS/data/*
 *.gif
 *.log
 *.mp4
-*.dump
 
 *.stderr
 

GravityField/animation.f90

@@ -10,9 +10,9 @@ program animation
   implicit none
 
   ! some configuration constants
-  integer, parameter      :: S_WIDTH  = 1600
-  integer, parameter      :: S_HEIGHT = 1600
-  integer, parameter      :: NB_BODY  = 60
+  integer, parameter      :: S_WIDTH  = 1024
+  integer, parameter      :: S_HEIGHT = 1024
+  integer, parameter      :: NB_BODY  = 70
 
   !!! WARNING : global variable !!!
   type(massbody)          :: planets(NB_BODY)
@@ -71,7 +71,7 @@ subroutine deplace_les_planetes(moons, clipit)
   integer                 :: foo
   real                    :: depx, depy
 
-  integer, parameter      :: EE = 51
+  integer, parameter      :: EE = 45
   integer                 :: SW = S_WIDTH  - EE
   integer                 :: SH = S_HEIGHT - EE
 
@@ -100,7 +100,7 @@ subroutine deplace_les_planetes(moons, clipit)
       if (moons(foo)%posy .LT. EE)     moons(foo)%posy = SH
     endif
 
-    moons(foo)%heading = moons(foo)%heading + (0.73*(rand()-0.50))
+    moons(foo)%heading = moons(foo)%heading + (0.71*(rand()-0.50))
     if (moons(foo)%heading .GT. 6.2831853) moons(foo)%heading = 0.0
     if (moons(foo)%heading .LT. 0.0000001) moons(foo)%heading = 0.0
 

GravityField/essai.f90

@@ -9,71 +9,53 @@ program essai
 
   implicit none
 
-  call init_random()
+  ! some configuration constants
+  integer, parameter      :: S_WIDTH  = 800
+  integer, parameter      :: S_HEIGHT = 600
+  integer, parameter      :: NB_BODY  = 51
 
-  call essai_near_planet(9999, 4096)
+  type(massbody)          :: planets(NB_BODY)
+  integer                 :: foo
+  character(len=100)      :: filename
+
+  call init_random()
+  call create_some_planets(planets, 45e5, S_WIDTH, S_HEIGHT)
+
+  do foo=0, 1999
+    write (filename, "(a, i5.5, a)") 'WS/field/', foo, '.pgm'
+    call build_and_write_a_field(S_WIDTH, S_HEIGHT, planets, filename)
+    ! print *, trim(filename)
+    ! OMG! two magic numbers, wtf?
+    planets(1)%posx = planets(1)%posx + 3 + (4.5*rand())
+    planets(1)%posy = planets(1)%posy + 3 + (2.1*rand())
+    if (planets(1)%posx .GT. S_WIDTH)   planets(1)%posx = 0.0
+    if (planets(1)%posy .GT. S_HEIGHT)  planets(1)%posy = 0.0
+
+    call boulegue_les_astres(planets, 2.21)
+
+  enddo
 
   STOP 'BECAUSE YOLO'
 
 contains
 !-----------------------------------------------------------------------
 !-
-!     computation of thr nearest planet
+!    Et si on bougeait un peu tous ces corps planétaires ?
 !-
-subroutine essai_near_planet(nbplanets, szfield)
-  integer, intent(in)       :: nbplanets, szfield
+subroutine boulegue_les_astres(astres, factor)
+  type(massbody), intent(inout)     :: astres(:)
+  real, intent(in)                  :: factor
 
-  integer, dimension(:,:), allocatable     :: map
-  integer                                  :: ix, iy
-  real                                     :: fx, fy, dx, dy
-  integer                                  :: near, ipl, errcode
-  real                                     :: curdist, smalldist 
-  type(massbody)                           :: planets(nbplanets)
-
-  print *, "near planets test", nbplanets, szfield
-
-  allocate(map(szfield, szfield), stat=errcode)
-  map = -1
-
-  ! create some random bodies
-  do ipl=1, nbplanets
-    planets(ipl)%posx   = rand() * szfield
-    planets(ipl)%posy   = rand() * szfield
-    planets(ipl)%serial = ipl
-  end do
-  ! call save_bodies_to_txt_file(planets, "planets.txt")
-
-  ! loop over all the location of the field
-  do ix=1, szfield
-    fx = real(ix)
-    do iy=1, szfield
-      fy = real(iy)
-
-      near      = -1
-      smalldist = 1e37
-      ! loop over all the planet's bodies
-      do ipl=1, nbplanets
-        ! compute the "fake" distance
-        dx = fx - planets(ipl)%posx
-        dy = fy - planets(ipl)%posy
-        curdist = (dx*dx) + (dy*dy)
-        if (curdist .LT. smalldist) then
-          near = ipl
-          smalldist = curdist
-        endif
-      end do    ! loop on ipl  
-
-      map(ix, iy) = mod(near, 255)
-
-    enddo
-
-    write(0, *) "row", ix, "     on", szfield
+  integer                           :: foo
 
+  do foo = 2, ubound(astres, 1)
+      astres(foo)%posx = astres(foo)%posx + factor*(rand() - 0.5)
+      astres(foo)%posy = astres(foo)%posy + factor*(rand() - 0.5)
   enddo
 
-  call spit_as_pgm_8(map, "nearest.pgm")
 
 end subroutine
 !-----------------------------------------------------------------------
+!-----------------------------------------------------------------------
 
 end program

GravityField/plotation.sh

@@ -14,7 +14,7 @@ NUMERO=$1
 
 SRC=$(printf "WS/data/%05d.txt"    $NUMERO)
 DST=$(printf "WS/graph/%05d.png"   $NUMERO)
-TXT=$(printf "mb #%05d"            $NUMERO)
+TXT=$(printf "Mass bodies #%05d"   $NUMERO)
 
 echo "    " $SRC $DST "      " $TXT
 
@@ -23,20 +23,20 @@ set term png    size 512,512
 set xrange    [ 0.0 : 1024 ]
 set yrange    [ 0.0 : 1024 ]
 set output      "${DST}"
-set ytics       128
-set xtics       128
+set ytics       256
+set xtics       256
 set grid        front
 set title       "${TXT}"
-plot            "$SRC"  using 1:2 title "loc"
+plot            "$SRC"  using 1:2 title "location" 
 __EOC__
 }
 # ----------------------------------------------------------
 
-for foo in $(seq 0 5 1999)
+for foo in $(seq 0 8 1999)
 do
 	plot_a_map $foo
 done
 
-convert -delay 10 WS/graph/0????.png displacement.gif
+convert -delay 10 WS/graph/0????.png foo.gif
 
 # ----------------------------------------------------------

GravityField/raytrace.sh

@@ -24,14 +24,14 @@ outfile=$(printf "WS/troid/%05d.png" $clock)
 echo  $timestamp $texte $outfile
 
 convert	${TMPF}				\
-	-pointsize 22			\
+	-pointsize 24			\
 	-font Courier-Bold		\
 	-fill Orange			\
-	-annotate +15+22  "$timestamp"	\
-	-annotate +15+48  "$texte"	\
+	-annotate +15+16  "$timestamp"	\
+	-annotate +15+38  "$texte"	\
 	-pointsize 16			\
 	-gravity south-west		\
-	-annotate +15+9 "tTh & Konrad"	\
+	-annotate +15+6 "tTh & Konrad"	\
 	${outfile}
 }
 

GravityField/realfield.f90

@@ -72,7 +72,7 @@ subroutine create_some_planets(planets, coef, sx, sy)
       !-
       planets(1)%posx   = sx / 2
       planets(1)%posy   = sy / 2
-      planets(1)%mass   = 31e8
+      planets(1)%mass   = 29e8
       planets(1)%serial = 1337
       planets(1)%speed  = 6.666
     else
@@ -82,11 +82,10 @@ subroutine create_some_planets(planets, coef, sx, sy)
       planets(foo)%posx     = rand() * real(sx-1)
       planets(foo)%posy     = rand() * real(sy-1)
       planets(foo)%mass     = 7e6 + coef*foo
-      planets(foo)%heading  = 2 * 3.14159 * rand()
-      if (rand() .LT. 0.15) planets(foo)%speed = 3.14159
+      planets(foo)%heading  = 3.14159 * rand()
+      if (rand() .LT. 0.08) planets(foo)%speed = 3.14159
       planets(foo)%serial   = foo*2 + 120
     endif
-
     write (*, fmt) foo, planets(foo)
   enddo
 end subroutine
@@ -111,29 +110,11 @@ function compute_gravity(fx, fy, body)
   endif
 
 end function
-!-----------------------------------------------------------------------
-!-
-!   Export a massbody area to a text file. no error check, wtf ?
-!-
-subroutine save_bodies_to_txt_file (astres, fname)
-  type(massbody), intent(in)        :: astres(:)
-  character(len=*), intent(in)      :: fname
 
-  character(50)                     :: fmt
-  integer                           :: io, idx
-
-  write(0, "('saving planets to ', A20)") fname
-  fmt = "( 2(F9.3, ' ')   2(F9.3, '  '), F14.3,  I8)"
-  open(newunit=io, file=fname)
-  do idx = 1, ubound(astres, 1)
-    write(io, fmt) astres(idx)
-  enddo
-  close(io)
-end subroutine
 !-----------------------------------------------------------------------
 !-
 !   Compute the gravity field in a pre-allocated array relative
-!   to the massbody 'moon'. Nobody know where the magic numbers
+!   to the massbody 'moon'. Nobody know where the magic number
 !   come from, sorry.
 !-
 subroutine compute_a_field(field, moon)
@@ -177,7 +158,7 @@ subroutine build_and_write_a_field(szx, szy, moons, fname)
   field = 0.0
   do foo=1, ubound(moons, 1)
     call compute_a_field(tmpf, moons(foo))
-    tmpf = tmpf * 0.018
+    tmpf = tmpf * 0.019
     field = field + tmpf
   enddo
 
@@ -218,67 +199,9 @@ subroutine init_random()
   ! you MUST use it for initializing the initializer
   do t3=1, 4
     dummy = rand()
-    write(0, '("   dummy", I4, F9.6)') t3, dummy
+    write(0, *) 'dummy ', t3, dummy
   enddo
 
-end subroutine
-!-----------------------------------------------------------------------
-!-
-!    dump a field of reals numbers to disk - preliminary version
-!-
-subroutine dump_a_field_to_file(field, fname)
-  real, dimension(:,:), intent(in)  :: field
-  character(len=*), intent(in)      :: fname
-
-  integer                         :: header(8)
-  integer                         :: io
-
-  print *, "D) field size  ", ubound(field, 1), "W", ubound(field, 2), "H"
-  print *, "D) filename    ", fname
-
-  header     = 0
-  header(1)  = 574908040             ! magic number
-  header(2)  = 1                     ! this is a dump of real field
-  header(3)  = ubound(field, 1)
-  header(4)  = ubound(field, 2)
-  header(5)  = 666
-
-  open(newunit=io, file=fname, form='unformatted')
-  
-  write(io) header
-  write(io) field
-
-  close(io)
-
-end subroutine
-!-----------------------------------------------------------------------
-!-
-!    load a real field from file  - preliminary version
-!-
-subroutine load_a_field_from_file(field, fname)
-  real, dimension(:,:), intent(in)  :: field
-  character(len=*), intent(in)      :: fname
-
-  integer                         :: header(8)
-  integer                         :: io, foo
-
-  print *, "L) field size  ", ubound(field, 1), "W", ubound(field, 2), "H"
-
-  !-
-  ! how to check if the field array was valid ?
-  !-
-
-  open(newunit=io, file=fname, form='unformatted', status='old',     &
-                                                   action='read')
-  read(io) header
-  do foo=1, 8
-    print *, foo, header(foo)
-  enddo
-
-  STOP ' --- FUCKED UP BEYOND ALL REPAIR ---'
-
-  close(io)
-
 end subroutine
 !-----------------------------------------------------------------------
 !-----------------------------------------------------------------------

GravityField/vision.pov

@@ -31,8 +31,8 @@ height_field {
         translate <-0.5, 0, -0.5>
 	}
 texture {
-	pigment { color Gray70 }
-	finish { phong 0.55 }
+	pigment { color Gray80 }
+	finish { phong 0.50 }
 	}
 }
 
@@ -55,7 +55,7 @@ light_source { < -15, 2.3, 17>	color Green*0.25			}
 // ----------------------------------------------------------------------
 
 camera	{
-	location	<-8, 4-(1.5*NormClock), 1 + (5*NormClock)>
+	location	<-8, 4-NormClock, 1 + (5*NormClock)>
 	look_at		<0, 0, 0>
 	right		x*image_width/image_height
 	angle		33