Simple approach to simulate water effect

Author

Message

skovron

20

Years of Service

User Offline

Joined: 14th Sep 2003

Location:

Posted: 26th Sep 2003 02:40 Edited at: 4th Oct 2003 00:00

Link

Just a simple Water Wave Effect
It generates waves on bitmap while you moving mouse and pressing left button.

The code is"ver unoptimised" so it rather run slowly but still it someone may find it useful. Experiment with some variables and remember to load rather small bitmaps, I used 320x240 becouse larger couse it run very very slow.

remstart
   Water Wave Effect
   by Mariusz Skowronski (2003)
remend

rem global vars
rem bitmap size and color depth
global g_bmpWidth as integer
global g_bmpHeight as integer
global g_bmpBpp as integer
global g_bytesPerPixel as integer

rem memblock pointers (for faster access)
global g_srcMemPtr as dword
global g_dstMemPtr as dword

rem screen dimensions
global g_scrWidth as integer
global g_scrHeight as integer

rem height array swap index (current and new)
global g_ch as integer
global g_nh as integer

rem properties
rem wave amplitude
global g_waveForce as integer
rem drop size
global g_dropRadius as integer
rem damping value (the smaller value the faster waves disappear) 1 - 100
global g_waveDamp as integer

rem init vars
g_scrWidth = 640
g_scrHeight = 480

g_waveForce = 300
g_dropRadius = 8
g_waveDamp = 20 : rem DON'T USE ZERO

rem begin
set display mode g_scrWidth, g_scrHeight, 16
sync rate 0: sync on

rem load bitmap (remember the smaller bitmap the faster it runs)
load bitmap "example1.jpg", 1
g_bmpWidth = bitmap width(1)
g_bmpHeight = bitmap height(1)
g_bmpBpp = bitmap depth(1)
g_bytesPerPixel = g_bmpBpp/8

rem create 2 bitmap memblocks
make memblock from bitmap 1, 1
make memblock from bitmap 2, 1

rem get memblock addresses
g_srcMemPtr = get memblock ptr(2)
g_dstMemPtr = get memblock ptr(1)

rem dont need this bitmap anymore
delete bitmap 1

rem allocate 2 height maps
dim g_heightMap(1, g_bmpWidth-1, g_bmpHeight-1) as integer
rem clear height maps
ClearHeights()

rem set current height map to 0
g_ch = 0: g_nh = 1

rem main loop
set current bitmap 0
do
   rem update heights
   UpdateHeights()
   rem update bitmap
   UpdateBitmap()
   rem uncomment following lie and rem this above to see it in greyscale
   rem DrawHeights()

rem create bitmap form memblock
   make bitmap from memblock 0, 1

sync

if mouseclick()=1
      rem make wave at point mx, my
      MakeWave(mousex(), mousey(), g_dropRadius, g_waveForce)
   endif

loop

delete memblock 1
delete memblock 2
rem that's all
end

rem this clears height maps
function ClearHeights()
   local x, y

for y=0 to g_bmpHeight-1
      for x=0 to g_bmpWidth-1
         g_heightMap(0, x, y) = 0
         g_heightMap(1, x, y) = 0
      next x
   next y
endfunction

rem this makes a new wave
function MakeWave(x, y, radius, mul)
   local u, v, sqru, sqrv, sqrr

sqrr = radius*radius
   rem check if not out of bitmap boundary
   if (x>radius) and (x<g_bmpWidth-1-radius) and (y>radius) and (y<g_bmpHeight-1-radius)
      for v=y-radius to y+radius
         sqrv = (v-y)*(v-y)
         for u=x-radius to x+radius
            sqru = (u-x)*(u-x)
            rem for speed reasons check (dx^2 + dy^2) <= r^2
            rem instead of real formula: sqrt(dx^2 + dy^) < r
            if (sqru+sqrv)<=sqrr:
               g_heightMap(g_ch,u,v) = mul*int(radius - sqrt(sqru+sqrv)+1)
            endif
         next u
      next v
   endif

endfunction

rem this is only for debug (shows waves in grayscale)
function DrawHeights()
   local x, y, height, color
   local dstAddr as dword
   rem for each height
   for y=0 to g_bmpHeight-1
      for x=0 to g_bmpWidth-1
         height = g_heightMap(g_ch, x, y)/10
         rem prepare gray color
         color = rgb(height+128, height+128, height+128)
         rem dont need this if useing one color depth
         select g_bytesPerPixel
            case 2: dstAddr = g_dstMemPtr +12+(y*g_bmpWidth+x)*2: endcase
            case 4: dstAddr = g_dstMemPtr +12+(y*g_bmpWidth+x)*4: endcase
         endselect
         *dstAddr = color
      next x
   next y
endfunction

rem the meat. this updates height map
function UpdateHeights()
   local x, y,  height, temp

rem we need skip bitmap edges couse we're making kind of area sampling
   for y=2 to g_bmpHeight-3
      for x=2 to g_bmpWidth-3

rem calculate height difference
         rem read 12 height values surrounding (x,y), divide by 6 so we've got twice as much as average value
         rem and then subtract current height at (x, y)
         height = g_heightMap(g_ch,x-1,y)
         height = height + g_heightMap(g_ch,x-2,y)
         height = height + g_heightMap(g_ch,x+1,y)
         height = height + g_heightMap(g_ch,x+2,y)
         height = height + g_heightMap(g_ch,x,y-1)
         height = height + g_heightMap(g_ch,x,y-2)
         height = height + g_heightMap(g_ch,x,y+1)
         height = height + g_heightMap(g_ch,x,y+2)
         height = height + g_heightMap(g_ch,x-1,y-1)
         height = height + g_heightMap(g_ch,x+1,y-1)
         height = height + g_heightMap(g_ch,x-1,y+1)
         height = height + g_heightMap(g_ch,x+1,y+1)
         height = int(height/6) - g_heightMap(g_nh,x,y)

rem add some dampig reducing amplitudes
         height = height - int(height/g_waveDamp + 1)

g_heightMap(g_nh, x, y) = height
      next x
   next y

rem swap height maps
   temp = g_ch
   g_ch = g_nh
   g_nh = temp
endfunction

rem this creates waved bitmap
function UpdateBitmap()
   local x, y, tx, ty, color
   local srcAddr, dstAddr as dword

rem skip bitmap edges
   for y=1 to g_bmpHeight-1
      for x=1 to g_bmpWidth-1

rem calculate pixel displacements based on height difference
         tx = x + (g_heightMap(g_ch, x-1, y) - g_heightMap(g_ch, x, y))/10: rem divide couse height values are too high
         ty = y + (g_heightMap(g_ch, x, y-1) - g_heightMap(g_ch, x, y))/10

rem if (tx, ty) still inside bitmap bounds
         if tx>=0 and tx<g_bmpWidth and ty>=0 and ty<g_bmpHeight

rem read color from coords we've calculated
            select g_bytesPerPixel
               case 2: srcAddr = g_srcMemPtr+12+(ty*g_bmpWidth+tx)*2: endcase
               case 4: srcAddr = g_srcMemPtr+12+(ty*g_bmpWidth+tx)*4: endcase
            endselect

color = *srcAddr
         else
            rem if out of bitmap bounds let the color be black
            color = 0
         endif

rem store this color at current position
         select g_bytesPerPixel
            case 2: dstAddr = g_dstMemPtr+12+(y*g_bmpWidth+x)*2: endcase
            case 4: dstAddr = g_dstMemPtr+12+(y*g_bmpWidth+x)*4: endcase
         endselect
         *dstAddr = color
      next x
   next y
endfunction

+ Code Snippet

remstart
   Water Wave Effect
   by Mariusz Skowronski (2003)
remend

rem global vars
rem bitmap size and color depth
global g_bmpWidth as integer
global g_bmpHeight as integer
global g_bmpBpp as integer
global g_bytesPerPixel as integer

rem memblock pointers (for faster access)
global g_srcMemPtr as dword
global g_dstMemPtr as dword

rem screen dimensions
global g_scrWidth as integer
global g_scrHeight as integer

rem height array swap index (current and new)
global g_ch as integer
global g_nh as integer


rem properties
rem wave amplitude
global g_waveForce as integer
rem drop size
global g_dropRadius as integer
rem damping value (the smaller value the faster waves disappear) 1 - 100
global g_waveDamp as integer


rem init vars
g_scrWidth = 640
g_scrHeight = 480

g_waveForce = 300
g_dropRadius = 8
g_waveDamp = 20 : rem DON'T USE ZERO

rem begin
set display mode g_scrWidth, g_scrHeight, 16
sync rate 0: sync on

rem load bitmap (remember the smaller bitmap the faster it runs)
load bitmap "example1.jpg", 1
g_bmpWidth = bitmap width(1)
g_bmpHeight = bitmap height(1)
g_bmpBpp = bitmap depth(1)
g_bytesPerPixel = g_bmpBpp/8


rem create 2 bitmap memblocks
make memblock from bitmap 1, 1
make memblock from bitmap 2, 1

rem get memblock addresses
g_srcMemPtr = get memblock ptr(2)
g_dstMemPtr = get memblock ptr(1)

rem dont need this bitmap anymore
delete bitmap 1

rem allocate 2 height maps
dim g_heightMap(1, g_bmpWidth-1, g_bmpHeight-1) as integer
rem clear height maps
ClearHeights()

rem set current height map to 0
g_ch = 0: g_nh = 1

rem main loop
set current bitmap 0
do
   rem update heights
   UpdateHeights()
   rem update bitmap
   UpdateBitmap()
   rem uncomment following lie and rem this above to see it in greyscale
   rem DrawHeights()

   rem create bitmap form memblock
   make bitmap from memblock 0, 1

   sync

   if mouseclick()=1
      rem make wave at point mx, my
      MakeWave(mousex(), mousey(), g_dropRadius, g_waveForce)
   endif

loop

delete memblock 1
delete memblock 2
rem that's all
end

rem this clears height maps
function ClearHeights()
   local x, y

   for y=0 to g_bmpHeight-1
      for x=0 to g_bmpWidth-1
         g_heightMap(0, x, y) = 0
         g_heightMap(1, x, y) = 0
      next x
   next y
endfunction


rem this makes a new wave
function MakeWave(x, y, radius, mul)
   local u, v, sqru, sqrv, sqrr

   sqrr = radius*radius
   rem check if not out of bitmap boundary
   if (x>radius) and (x<g_bmpWidth-1-radius) and (y>radius) and (y<g_bmpHeight-1-radius)
      for v=y-radius to y+radius
         sqrv = (v-y)*(v-y)
         for u=x-radius to x+radius
            sqru = (u-x)*(u-x)
            rem for speed reasons check (dx^2 + dy^2) <= r^2
            rem instead of real formula: sqrt(dx^2 + dy^) < r
            if (sqru+sqrv)<=sqrr:
               g_heightMap(g_ch,u,v) = mul*int(radius - sqrt(sqru+sqrv)+1)
            endif
         next u
      next v
   endif

endfunction

rem this is only for debug (shows waves in grayscale)
function DrawHeights()
   local x, y, height, color
   local dstAddr as dword
   rem for each height
   for y=0 to g_bmpHeight-1
      for x=0 to g_bmpWidth-1
         height = g_heightMap(g_ch, x, y)/10
         rem prepare gray color
         color = rgb(height+128, height+128, height+128)
         rem dont need this if useing one color depth
         select g_bytesPerPixel
            case 2: dstAddr = g_dstMemPtr +12+(y*g_bmpWidth+x)*2: endcase
            case 4: dstAddr = g_dstMemPtr +12+(y*g_bmpWidth+x)*4: endcase
         endselect
         *dstAddr = color
      next x
   next y
endfunction


rem the meat. this updates height map
function UpdateHeights()
   local x, y,  height, temp

   rem we need skip bitmap edges couse we're making kind of area sampling
   for y=2 to g_bmpHeight-3
      for x=2 to g_bmpWidth-3

         rem calculate height difference
         rem read 12 height values surrounding (x,y), divide by 6 so we've got twice as much as average value
         rem and then subtract current height at (x, y)
         height = g_heightMap(g_ch,x-1,y)
         height = height + g_heightMap(g_ch,x-2,y)
         height = height + g_heightMap(g_ch,x+1,y)
         height = height + g_heightMap(g_ch,x+2,y)
         height = height + g_heightMap(g_ch,x,y-1)
         height = height + g_heightMap(g_ch,x,y-2)
         height = height + g_heightMap(g_ch,x,y+1)
         height = height + g_heightMap(g_ch,x,y+2)
         height = height + g_heightMap(g_ch,x-1,y-1)
         height = height + g_heightMap(g_ch,x+1,y-1)
         height = height + g_heightMap(g_ch,x-1,y+1)
         height = height + g_heightMap(g_ch,x+1,y+1)
         height = int(height/6) - g_heightMap(g_nh,x,y)

         rem add some dampig reducing amplitudes
         height = height - int(height/g_waveDamp + 1)

         g_heightMap(g_nh, x, y) = height
      next x
   next y

   rem swap height maps
   temp = g_ch
   g_ch = g_nh
   g_nh = temp
endfunction


rem this creates waved bitmap
function UpdateBitmap()
   local x, y, tx, ty, color
   local srcAddr, dstAddr as dword

   rem skip bitmap edges
   for y=1 to g_bmpHeight-1
      for x=1 to g_bmpWidth-1

         rem calculate pixel displacements based on height difference
         tx = x + (g_heightMap(g_ch, x-1, y) - g_heightMap(g_ch, x, y))/10: rem divide couse height values are too high
         ty = y + (g_heightMap(g_ch, x, y-1) - g_heightMap(g_ch, x, y))/10

         rem if (tx, ty) still inside bitmap bounds
         if tx>=0 and tx<g_bmpWidth and ty>=0 and ty<g_bmpHeight

            rem read color from coords we've calculated
            select g_bytesPerPixel
               case 2: srcAddr = g_srcMemPtr+12+(ty*g_bmpWidth+tx)*2: endcase
               case 4: srcAddr = g_srcMemPtr+12+(ty*g_bmpWidth+tx)*4: endcase
            endselect

            color = *srcAddr
         else
            rem if out of bitmap bounds let the color be black
            color = 0
         endif

         rem store this color at current position
         select g_bytesPerPixel
            case 2: dstAddr = g_dstMemPtr+12+(y*g_bmpWidth+x)*2: endcase
            case 4: dstAddr = g_dstMemPtr+12+(y*g_bmpWidth+x)*4: endcase
         endselect
         *dstAddr = color
      next x
   next y
endfunction

Back to top

Profile PM

Hamish McHaggis

21

Years of Service

User Offline

Joined: 13th Dec 2002

Location: Modgnik Detinu

Posted: 27th Sep 2003 01:00

Link

Nout happens for me...

Brains are for idiots.

Athelon XP 1400 Plus - Nvidia Geforce MX400 - 256mb RAM

Back to top

Profile PM Email Website

Sedate Solution

21

Years of Service

User Offline

Joined: 15th Apr 2003

Location:

Posted: 27th Sep 2003 01:13

Link

You have posted a DEMO that requires an external file, can you incorporate the external file into the demo or give a link for it to be downloaded.

Back to top

Profile PM

skovron

20

Years of Service

User Offline

Joined: 14th Sep 2003

Location:

Posted: 27th Sep 2003 22:09

Link

Yea, sorry for that (you could use any picture)
but ofcourse, here is the link for example picture

http://www.arena.free-host.com/examples/example1.jpg

Back to top

Profile PM

D Man

21

Years of Service

User Offline

Joined: 3rd Oct 2002

Location: Germany

Posted: 28th Sep 2003 14:25

Link

This is very nice.
I think I'll try to optimize it a little bit.

God is real, unless declared integer.

Back to top

Profile PM Email

skovron

20

Years of Service

User Offline

Joined: 14th Sep 2003

Location:

Posted: 30th Sep 2003 01:43

Link

Just let me know when you finished. Id like to compare this 2 versions

Back to top

Profile PM

D Man

21

Years of Service

User Offline

Joined: 3rd Oct 2002

Location: Germany

Posted: 30th Sep 2003 13:27

Link

Well, at least I could double the frame rate.
I wished that I could do more, but it's really hard for me to see throw the code.
I think I'll try to understand this code fully

.

remstart
   Water Wave Effect
   by Mariusz Skowronski (2003)
remend

disable escapekey

rem global vars
rem bitmap size and color depth
global g_bmpWidth as integer
global g_bmpHeight as integer
global g_bmpBpp as integer
global g_bytesPerPixel as integer

rem memblock pointers (for faster access)
global g_srcMemPtr as dword
global g_dstMemPtr as dword

rem screen dimensions
global g_scrWidth as integer
global g_scrHeight as integer

rem height array swap index (current and new)
global g_ch as integer
global g_nh as integer

rem properties
rem wave amplitude
global g_waveForce as integer
rem drop size
global g_dropRadius as integer
rem damping value (the smaller value the faster waves disappear) 1 - 100
global g_waveDamp as integer

rem init vars
g_scrWidth = 640
g_scrHeight = 480

g_waveForce = 300
g_dropRadius = 8
g_waveDamp = 20 : rem DON'T USE ZERO

rem begin
set display mode g_scrWidth, g_scrHeight, 32
sync on : sync rate 0

rem load bitmap (remember the smaller bitmap the faster it runs)
load bitmap "example1.jpg", 1
g_bmpWidth = bitmap width(1)
g_bmpHeight = bitmap height(1)
g_bmpBpp = bitmap depth(1)
g_bytesPerPixel = g_bmpBpp/8

rem create 2 bitmap memblocks
make memblock from bitmap 2, 1

rem get memblock addresses
g_srcMemPtr = get memblock ptr(2)

rem dont need this bitmap anymore
delete bitmap 1

rem allocate 2 height maps
dim g_heightMap(1, g_bmpWidth-1, g_bmpHeight-1) as integer
rem clear height maps
ClearHeights()

rem set current height map to 0
g_ch = 0: g_nh = 1

set text opaque

do
   rem update heights
   UpdateHeights()
   rem update bitmap
   UpdateBitmap()
   rem uncomment following lie and rem this above to see it in greyscale
   rem DrawHeights()

text screen width()-text width("FPS: "+str$(screen fps())),screen height()-text size(),"FPS: "+str$(screen fps())
   sync

if mouseclick()=1
      rem make wave at point mx, my
      MakeWave(mousex(), mousey(), g_dropRadius, g_waveForce)
   endif

if escapekey() then exit

loop

delete memblock 2
rem that's all
end

rem this clears height maps
function ClearHeights()
   local x, y
   for y=0 to g_bmpHeight-1
      for x=0 to g_bmpWidth-1
         g_heightMap(0, x, y) = 0
         g_heightMap(1, x, y) = 0
      next x
   next y
endfunction

rem this makes a new wave
function MakeWave(x, y, radius, mul)
   local u, v, sqru, sqrv, sqrr

sqrr = radius*radius
   rem check if not out of bitmap boundary
   if (x>radius) and (x<g_bmpWidth-1-radius) and (y>radius) and (y<g_bmpHeight-1-radius)
      for v=y-radius to y+radius
         sqrv = (v-y)*(v-y)
         for u=x-radius to x+radius
            sqru = (u-x)*(u-x)
            rem for speed reasons check (dx^2 + dy^2) <= r^2
            rem instead of real formula: sqrt(dx^2 + dy^) < r
            if (sqru+sqrv)<=sqrr:
               g_heightMap(g_ch,u,v) = mul*int(radius - sqrt(sqru+sqrv)+1)
            endif
         next u
      next v
   endif

endfunction

rem this is only for debug (shows waves in grayscale)
function DrawHeights()
   local x, y, height, color
   local dstAddr as dword
   rem for each height
   for y=0 to g_bmpHeight-1
      for x=0 to g_bmpWidth-1
         height = g_heightMap(g_ch, x, y)/10
         rem prepare gray color
         color = rgb(height+128, height+128, height+128)
         rem dont need this if useing one color depth
         select g_bytesPerPixel
            case 2: dstAddr = g_dstMemPtr +12+(y*g_bmpWidth+x)*2: endcase
            case 4: dstAddr = g_dstMemPtr +12+(y*g_bmpWidth+x)*4: endcase
         endselect
         *dstAddr = color
      next x
   next y
endfunction

rem the meat. this updates height map
function UpdateHeights()
   local x, y,  height, temp

rem we need skip bitmap edges couse we're making kind of area sampling
   for y=2 to g_bmpHeight-3
      for x=2 to g_bmpWidth-3

rem calculate height difference
         rem read 12 height values surrounding (x,y), divide by 6 so we've got twice as much as average value
         rem and then subtract current height at (x, y)
         height = g_heightMap(g_ch,x-1,y)
         height = height + g_heightMap(g_ch,x-2,y)
         height = height + g_heightMap(g_ch,x+1,y)
         height = height + g_heightMap(g_ch,x+2,y)
         height = height + g_heightMap(g_ch,x,y-1)
         height = height + g_heightMap(g_ch,x,y-2)
         height = height + g_heightMap(g_ch,x,y+1)
         height = height + g_heightMap(g_ch,x,y+2)
         height = height + g_heightMap(g_ch,x-1,y-1)
         height = height + g_heightMap(g_ch,x+1,y-1)
         height = height + g_heightMap(g_ch,x-1,y+1)
         height = height + g_heightMap(g_ch,x+1,y+1)
         height = int(height/6) - g_heightMap(g_nh,x,y)

rem add some dampig reducing amplitudes
         height = height - int(height/g_waveDamp + 1)

g_heightMap(g_nh, x, y) = height
      next x
   next y

rem swap height maps
   temp = g_ch
   g_ch = g_nh
   g_nh = temp
endfunction

rem this creates waved bitmap
function UpdateBitmap()
   local x, y, tx, ty
   local srcAddr, dstAddr, color as dword

lock backbuffer
   backptr=get backbuffer ptr()
   rem skip bitmap edges
   for y=0 to g_bmpHeight-1
      for x=0 to g_bmpWidth-1

rem calculate pixel displacements based on height difference
         tx = x + (g_heightMap(g_ch, x-1, y) - g_heightMap(g_ch, x, y))/10: rem divide couse height values are too high
         ty = y + (g_heightMap(g_ch, x, y-1) - g_heightMap(g_ch, x, y))/10

rem if (tx, ty) still inside bitmap bounds
         if tx>=0 and tx<g_bmpWidth and ty>=0 and ty<g_bmpHeight
            rem read color from coords we've calculated
            select g_bytesPerPixel
               case 2 : copy memory backptr+12+(y*screen width()+x)*2,g_srcMemPtr+12+(ty*g_bmpWidth+tx)*2,2 : endcase
               case 4 : copy memory backptr+12+(y*screen width()+x)*4,g_srcMemPtr+12+(ty*g_bmpWidth+tx)*4,4 : endcase
            endselect
         endif
      next x
   next y
   unlock backbuffer
endfunction

+ Code Snippet

remstart
   Water Wave Effect
   by Mariusz Skowronski (2003)
remend

disable escapekey

rem global vars
rem bitmap size and color depth
global g_bmpWidth as integer
global g_bmpHeight as integer
global g_bmpBpp as integer
global g_bytesPerPixel as integer

rem memblock pointers (for faster access)
global g_srcMemPtr as dword
global g_dstMemPtr as dword

rem screen dimensions
global g_scrWidth as integer
global g_scrHeight as integer

rem height array swap index (current and new)
global g_ch as integer
global g_nh as integer


rem properties
rem wave amplitude
global g_waveForce as integer
rem drop size
global g_dropRadius as integer
rem damping value (the smaller value the faster waves disappear) 1 - 100
global g_waveDamp as integer


rem init vars
g_scrWidth = 640
g_scrHeight = 480

g_waveForce = 300
g_dropRadius = 8
g_waveDamp = 20 : rem DON'T USE ZERO

rem begin
set display mode g_scrWidth, g_scrHeight, 32
sync on : sync rate 0

rem load bitmap (remember the smaller bitmap the faster it runs)
load bitmap "example1.jpg", 1
g_bmpWidth = bitmap width(1)
g_bmpHeight = bitmap height(1)
g_bmpBpp = bitmap depth(1)
g_bytesPerPixel = g_bmpBpp/8


rem create 2 bitmap memblocks
make memblock from bitmap 2, 1

rem get memblock addresses
g_srcMemPtr = get memblock ptr(2)

rem dont need this bitmap anymore
delete bitmap 1

rem allocate 2 height maps
dim g_heightMap(1, g_bmpWidth-1, g_bmpHeight-1) as integer
rem clear height maps
ClearHeights()

rem set current height map to 0
g_ch = 0: g_nh = 1

set text opaque

do
   rem update heights
   UpdateHeights()
   rem update bitmap
   UpdateBitmap()
   rem uncomment following lie and rem this above to see it in greyscale
   rem DrawHeights()

   text screen width()-text width("FPS: "+str$(screen fps())),screen height()-text size(),"FPS: "+str$(screen fps())
   sync

   if mouseclick()=1
      rem make wave at point mx, my
      MakeWave(mousex(), mousey(), g_dropRadius, g_waveForce)
   endif

   if escapekey() then exit

loop

delete memblock 2
rem that's all
end

rem this clears height maps
function ClearHeights()
   local x, y
   for y=0 to g_bmpHeight-1
      for x=0 to g_bmpWidth-1
         g_heightMap(0, x, y) = 0
         g_heightMap(1, x, y) = 0
      next x
   next y
endfunction


rem this makes a new wave
function MakeWave(x, y, radius, mul)
   local u, v, sqru, sqrv, sqrr

   sqrr = radius*radius
   rem check if not out of bitmap boundary
   if (x>radius) and (x<g_bmpWidth-1-radius) and (y>radius) and (y<g_bmpHeight-1-radius)
      for v=y-radius to y+radius
         sqrv = (v-y)*(v-y)
         for u=x-radius to x+radius
            sqru = (u-x)*(u-x)
            rem for speed reasons check (dx^2 + dy^2) <= r^2
            rem instead of real formula: sqrt(dx^2 + dy^) < r
            if (sqru+sqrv)<=sqrr:
               g_heightMap(g_ch,u,v) = mul*int(radius - sqrt(sqru+sqrv)+1)
            endif
         next u
      next v
   endif

endfunction

rem this is only for debug (shows waves in grayscale)
function DrawHeights()
   local x, y, height, color
   local dstAddr as dword
   rem for each height
   for y=0 to g_bmpHeight-1
      for x=0 to g_bmpWidth-1
         height = g_heightMap(g_ch, x, y)/10
         rem prepare gray color
         color = rgb(height+128, height+128, height+128)
         rem dont need this if useing one color depth
         select g_bytesPerPixel
            case 2: dstAddr = g_dstMemPtr +12+(y*g_bmpWidth+x)*2: endcase
            case 4: dstAddr = g_dstMemPtr +12+(y*g_bmpWidth+x)*4: endcase
         endselect
         *dstAddr = color
      next x
   next y
endfunction


rem the meat. this updates height map
function UpdateHeights()
   local x, y,  height, temp

   rem we need skip bitmap edges couse we're making kind of area sampling
   for y=2 to g_bmpHeight-3
      for x=2 to g_bmpWidth-3

         rem calculate height difference
         rem read 12 height values surrounding (x,y), divide by 6 so we've got twice as much as average value
         rem and then subtract current height at (x, y)
         height = g_heightMap(g_ch,x-1,y)
         height = height + g_heightMap(g_ch,x-2,y)
         height = height + g_heightMap(g_ch,x+1,y)
         height = height + g_heightMap(g_ch,x+2,y)
         height = height + g_heightMap(g_ch,x,y-1)
         height = height + g_heightMap(g_ch,x,y-2)
         height = height + g_heightMap(g_ch,x,y+1)
         height = height + g_heightMap(g_ch,x,y+2)
         height = height + g_heightMap(g_ch,x-1,y-1)
         height = height + g_heightMap(g_ch,x+1,y-1)
         height = height + g_heightMap(g_ch,x-1,y+1)
         height = height + g_heightMap(g_ch,x+1,y+1)
         height = int(height/6) - g_heightMap(g_nh,x,y)

         rem add some dampig reducing amplitudes
         height = height - int(height/g_waveDamp + 1)

         g_heightMap(g_nh, x, y) = height
      next x
   next y

   rem swap height maps
   temp = g_ch
   g_ch = g_nh
   g_nh = temp
endfunction


rem this creates waved bitmap
function UpdateBitmap()
   local x, y, tx, ty
   local srcAddr, dstAddr, color as dword

   lock backbuffer
   backptr=get backbuffer ptr()
   rem skip bitmap edges
   for y=0 to g_bmpHeight-1
      for x=0 to g_bmpWidth-1

         rem calculate pixel displacements based on height difference
         tx = x + (g_heightMap(g_ch, x-1, y) - g_heightMap(g_ch, x, y))/10: rem divide couse height values are too high
         ty = y + (g_heightMap(g_ch, x, y-1) - g_heightMap(g_ch, x, y))/10

         rem if (tx, ty) still inside bitmap bounds
         if tx>=0 and tx<g_bmpWidth and ty>=0 and ty<g_bmpHeight
            rem read color from coords we've calculated
            select g_bytesPerPixel
               case 2 : copy memory backptr+12+(y*screen width()+x)*2,g_srcMemPtr+12+(ty*g_bmpWidth+tx)*2,2 : endcase
               case 4 : copy memory backptr+12+(y*screen width()+x)*4,g_srcMemPtr+12+(ty*g_bmpWidth+tx)*4,4 : endcase
            endselect
         endif
      next x
   next y
   unlock backbuffer
endfunction

God is real, unless declared integer.

Back to top

Profile PM Email

D Man

21

Years of Service

User Offline

Joined: 3rd Oct 2002

Location: Germany

Posted: 30th Sep 2003 14:41

Link

Ok...
I could increase the sync rate once more.
Now it runs with ~28 fps on my pc.

remstart
   Water Wave Effect
   by Mariusz Skowronski (2003)
remend

disable escapekey

rem global vars
rem bitmap size and color depth
global g_bmpWidth as integer
global g_bmpHeight as integer
global g_bmpBpp as integer
global g_bytesPerPixel as integer

rem memblock pointers (for faster access)
global g_srcMemPtr as dword
global g_dstMemPtr as dword

rem screen dimensions
global g_scrWidth as integer
global g_scrHeight as integer

rem height array swap index (current and new)
global g_ch as integer
global g_nh as integer

rem properties
rem wave amplitude
global g_waveForce as integer
rem drop size
global g_dropRadius as integer
rem damping value (the smaller value the faster waves disappear) 1 - 100
global g_waveDamp as integer

rem init vars
g_scrWidth = 640
g_scrHeight = 480

g_waveForce = 300
g_dropRadius = 8
g_waveDamp = 20 : rem DON'T USE ZERO

rem begin
set display mode g_scrWidth, g_scrHeight, 32
sync on : sync rate 0

rem load bitmap (remember the smaller bitmap the faster it runs)
load bitmap "example1.jpg", 1
g_bmpWidth = bitmap width(1)
g_bmpHeight = bitmap height(1)
g_bmpBpp = bitmap depth(1)
g_bytesPerPixel = g_bmpBpp/8

rem create bitmap memblock
make memblock from bitmap 2, 1

rem get memblock addresses
g_srcMemPtr = get memblock ptr(2)

rem dont need this bitmap anymore
delete bitmap 1

rem allocate 2 height maps
dim g_heightMap(1, g_bmpWidth-1, g_bmpHeight-1) as integer
rem clear height maps
ClearHeights()

rem set current height map to 0
g_ch = 0: g_nh = 1

make image from memblock 1, 2
paste image 1,0,0
sync
delete image 1

set text opaque

do
   UpdateHeights()

if mouseclick() then MakeWave(mousex(), mousey(), g_dropRadius, g_waveForce)
   if escapekey() then exit

text screen width()-text width("FPS: "+str$(screen fps())),screen height()-text size(),"FPS: "+str$(screen fps())
   sync
loop

delete memblock 2
rem that's all
end

rem this clears height maps
function ClearHeights()
   local x, y
   for y=0 to g_bmpHeight-1
      for x=0 to g_bmpWidth-1
         g_heightMap(0, x, y) = 0
         g_heightMap(1, x, y) = 0
      next x
   next y
endfunction

rem this makes a new wave
function MakeWave(x, y, radius, mul)
   local u, v, sqru, sqrv, sqrr

sqrr = radius*radius
   rem check if not out of bitmap boundary
   if (x>radius) and (x<g_bmpWidth-1-radius) and (y>radius) and (y<g_bmpHeight-1-radius)
      for v=y-radius to y+radius
         sqrv = (v-y)*(v-y)
         for u=x-radius to x+radius
            sqru = (u-x)*(u-x)
            rem for speed reasons check (dx^2 + dy^2) <= r^2
            rem instead of real formula: sqrt(dx^2 + dy^) < r
            if (sqru+sqrv)<=sqrr:
               g_heightMap(g_ch,u,v) = mul*int(radius - sqrt(sqru+sqrv)+1)
            endif
         next u
      next v
   endif

endfunction

rem this is only for debug (shows waves in grayscale)
function DrawHeights()
   local x, y, height, color
   local dstAddr as dword
   rem for each height
   for y=0 to g_bmpHeight-1
      for x=0 to g_bmpWidth-1
         height = g_heightMap(g_ch, x, y)/10
         rem prepare gray color
         color = rgb(height+128, height+128, height+128)
         rem dont need this if useing one color depth
         select g_bytesPerPixel
            case 2: dstAddr = g_dstMemPtr +12+(y*g_bmpWidth+x)*2: endcase
            case 4: dstAddr = g_dstMemPtr +12+(y*g_bmpWidth+x)*4: endcase
         endselect
         *dstAddr = color
      next x
   next y
endfunction

rem the meat. this updates height map
function UpdateHeights()
   local x, y,  height, temp

lock backbuffer

backptr=get backbuffer ptr()

rem we need skip bitmap edges couse we're making kind of area sampling
   for y=2 to g_bmpHeight-3
      for x=2 to g_bmpWidth-3

rem calculate height difference
         rem read 12 height values surrounding (x,y), divide by 6 so we've got twice as much as average value
         rem and then subtract current height at (x, y)
         height = g_heightMap(g_ch,x-1,y)
         inc height,g_heightMap(g_ch,x-2,y)+g_heightMap(g_ch,x+1,y)+g_heightMap(g_ch,x+2,y)+g_heightMap(g_ch,x,y-1)+g_heightMap(g_ch,x,y-2)+g_heightMap(g_ch,x,y+1)+g_heightMap(g_ch,x,y+2)+g_heightMap(g_ch,x-1,y-1)+g_heightMap(g_ch,x+1,y-1)+g_heightMap(g_ch,x-1,y+1)+g_heightMap(g_ch,x+1,y+1)
         height = int(height/6) - g_heightMap(g_nh,x,y)
         rem add some dampig reducing amplitudes
         dec height,height/g_waveDamp + 1
         g_heightMap(g_nh, x, y) = height

tx = x + (g_heightMap(g_ch, x-1, y) - height)/10: rem divide couse height values are too high
         ty = y + (g_heightMap(g_ch, x, y-1) - height)/10
         rem if (tx, ty) still inside bitmap bounds
         if tx>=0 and tx<g_bmpWidth and ty>=0 and ty<g_bmpHeight and (tx<>x or ty<>y)
            rem read color from coords we've calculated
            select g_bytesPerPixel
               case 2 : copy memory backptr+12+(y*screen width()+x)*2,g_srcMemPtr+12+(ty*g_bmpWidth+tx)*2,2 : endcase
               case 4 : copy memory backptr+12+(y*screen width()+x)*4,g_srcMemPtr+12+(ty*g_bmpWidth+tx)*4,4 : endcase
            endselect
         endif
      next x
   next y

unlock backbuffer

rem swap height maps
   temp = g_ch
   g_ch = g_nh
   g_nh = temp
endfunction

+ Code Snippet

remstart
   Water Wave Effect
   by Mariusz Skowronski (2003)
remend

disable escapekey

rem global vars
rem bitmap size and color depth
global g_bmpWidth as integer
global g_bmpHeight as integer
global g_bmpBpp as integer
global g_bytesPerPixel as integer

rem memblock pointers (for faster access)
global g_srcMemPtr as dword
global g_dstMemPtr as dword

rem screen dimensions
global g_scrWidth as integer
global g_scrHeight as integer

rem height array swap index (current and new)
global g_ch as integer
global g_nh as integer


rem properties
rem wave amplitude
global g_waveForce as integer
rem drop size
global g_dropRadius as integer
rem damping value (the smaller value the faster waves disappear) 1 - 100
global g_waveDamp as integer


rem init vars
g_scrWidth = 640
g_scrHeight = 480

g_waveForce = 300
g_dropRadius = 8
g_waveDamp = 20 : rem DON'T USE ZERO

rem begin
set display mode g_scrWidth, g_scrHeight, 32
sync on : sync rate 0

rem load bitmap (remember the smaller bitmap the faster it runs)
load bitmap "example1.jpg", 1
g_bmpWidth = bitmap width(1)
g_bmpHeight = bitmap height(1)
g_bmpBpp = bitmap depth(1)
g_bytesPerPixel = g_bmpBpp/8


rem create bitmap memblock
make memblock from bitmap 2, 1


rem get memblock addresses
g_srcMemPtr = get memblock ptr(2)

rem dont need this bitmap anymore
delete bitmap 1

rem allocate 2 height maps
dim g_heightMap(1, g_bmpWidth-1, g_bmpHeight-1) as integer
rem clear height maps
ClearHeights()

rem set current height map to 0
g_ch = 0: g_nh = 1

make image from memblock 1, 2
paste image 1,0,0
sync
delete image 1

set text opaque

do
   UpdateHeights()

   if mouseclick() then MakeWave(mousex(), mousey(), g_dropRadius, g_waveForce)
   if escapekey() then exit

   text screen width()-text width("FPS: "+str$(screen fps())),screen height()-text size(),"FPS: "+str$(screen fps())
   sync
loop

delete memblock 2
rem that's all
end

rem this clears height maps
function ClearHeights()
   local x, y
   for y=0 to g_bmpHeight-1
      for x=0 to g_bmpWidth-1
         g_heightMap(0, x, y) = 0
         g_heightMap(1, x, y) = 0
      next x
   next y
endfunction


rem this makes a new wave
function MakeWave(x, y, radius, mul)
   local u, v, sqru, sqrv, sqrr

   sqrr = radius*radius
   rem check if not out of bitmap boundary
   if (x>radius) and (x<g_bmpWidth-1-radius) and (y>radius) and (y<g_bmpHeight-1-radius)
      for v=y-radius to y+radius
         sqrv = (v-y)*(v-y)
         for u=x-radius to x+radius
            sqru = (u-x)*(u-x)
            rem for speed reasons check (dx^2 + dy^2) <= r^2
            rem instead of real formula: sqrt(dx^2 + dy^) < r
            if (sqru+sqrv)<=sqrr:
               g_heightMap(g_ch,u,v) = mul*int(radius - sqrt(sqru+sqrv)+1)
            endif
         next u
      next v
   endif

endfunction

rem this is only for debug (shows waves in grayscale)
function DrawHeights()
   local x, y, height, color
   local dstAddr as dword
   rem for each height
   for y=0 to g_bmpHeight-1
      for x=0 to g_bmpWidth-1
         height = g_heightMap(g_ch, x, y)/10
         rem prepare gray color
         color = rgb(height+128, height+128, height+128)
         rem dont need this if useing one color depth
         select g_bytesPerPixel
            case 2: dstAddr = g_dstMemPtr +12+(y*g_bmpWidth+x)*2: endcase
            case 4: dstAddr = g_dstMemPtr +12+(y*g_bmpWidth+x)*4: endcase
         endselect
         *dstAddr = color
      next x
   next y
endfunction

rem the meat. this updates height map
function UpdateHeights()
   local x, y,  height, temp

   lock backbuffer

   backptr=get backbuffer ptr()

   rem we need skip bitmap edges couse we're making kind of area sampling
   for y=2 to g_bmpHeight-3
      for x=2 to g_bmpWidth-3

         rem calculate height difference
         rem read 12 height values surrounding (x,y), divide by 6 so we've got twice as much as average value
         rem and then subtract current height at (x, y)
         height = g_heightMap(g_ch,x-1,y)
         inc height,g_heightMap(g_ch,x-2,y)+g_heightMap(g_ch,x+1,y)+g_heightMap(g_ch,x+2,y)+g_heightMap(g_ch,x,y-1)+g_heightMap(g_ch,x,y-2)+g_heightMap(g_ch,x,y+1)+g_heightMap(g_ch,x,y+2)+g_heightMap(g_ch,x-1,y-1)+g_heightMap(g_ch,x+1,y-1)+g_heightMap(g_ch,x-1,y+1)+g_heightMap(g_ch,x+1,y+1)
         height = int(height/6) - g_heightMap(g_nh,x,y)
         rem add some dampig reducing amplitudes
         dec height,height/g_waveDamp + 1
         g_heightMap(g_nh, x, y) = height

         tx = x + (g_heightMap(g_ch, x-1, y) - height)/10: rem divide couse height values are too high
         ty = y + (g_heightMap(g_ch, x, y-1) - height)/10
         rem if (tx, ty) still inside bitmap bounds
         if tx>=0 and tx<g_bmpWidth and ty>=0 and ty<g_bmpHeight and (tx<>x or ty<>y)
            rem read color from coords we've calculated
            select g_bytesPerPixel
               case 2 : copy memory backptr+12+(y*screen width()+x)*2,g_srcMemPtr+12+(ty*g_bmpWidth+tx)*2,2 : endcase
               case 4 : copy memory backptr+12+(y*screen width()+x)*4,g_srcMemPtr+12+(ty*g_bmpWidth+tx)*4,4 : endcase
            endselect
         endif
      next x
   next y

   unlock backbuffer

   rem swap height maps
   temp = g_ch
   g_ch = g_nh
   g_nh = temp
endfunction

God is real, unless declared integer.

Back to top

Profile PM Email

skovron

20

Years of Service

User Offline

Joined: 14th Sep 2003

Location:

Posted: 30th Sep 2003 18:11

Link

Wow! Great work!
I had ~26 with last improvement and now I've got ~36

by the way, didnt know that inc, dec are faster than standard +, - operators

Back to top

Profile PM

skovron

20

Years of Service

User Offline

Joined: 14th Sep 2003

Location:

Posted: 1st Oct 2003 13:04

Link

I saw that when the waves disapear the image stays a little "pixelized" (It looks nice too, like picture from some voxel engine). To solve this function UpdateHeights() should rather looks like this:

function UpdateHeights()
   local x, y,  height, temp

lock backbuffer

backptr=get backbuffer ptr()

rem we need skip bitmap edges couse we're making kind of area sampling
   for y=2 to g_bmpHeight-3
      for x=2 to g_bmpWidth-3

rem calculate height difference
         rem read 12 height values surrounding (x,y), divide by 6 so we've got twice as much as average value
         rem and then subtract current height at (x, y)
         height = g_heightMap(g_ch,x-1,y)
         inc height,g_heightMap(g_ch,x-2,y)+g_heightMap(g_ch,x+1,y)+g_heightMap(g_ch,x+2,y)+g_heightMap(g_ch,x,y-1)+g_heightMap(g_ch,x,y-2)+g_heightMap(g_ch,x,y+1)+g_heightMap(g_ch,x,y+2)+g_heightMap(g_ch,x-1,y-1)+g_heightMap(g_ch,x+1,y-1)+g_heightMap(g_ch,x-1,y+1)+g_heightMap(g_ch,x+1,y+1)
         height = int(height/6) - g_heightMap(g_nh,x,y)
         rem add some dampig reducing amplitudes
         dec height,height/g_waveDamp + 1
         g_heightMap(g_nh, x, y) = height

tx = x + (g_heightMap(g_ch, x-1, y) - height)/10: rem divide couse height values are too high
         ty = y + (g_heightMap(g_ch, x, y-1) - height)/10
         rem if (tx, ty) still inside bitmap bounds

rem HERE

if tx>=0 and tx<g_bmpWidth and ty>=0 and ty<g_bmpHeight 
            rem read color from coords we've calculated
            select g_bytesPerPixel
               case 2 : copy memory backptr+12+(y*screen width()+x)*2,g_srcMemPtr+12+(ty*g_bmpWidth+tx)*2,2 : endcase
               case 4 : copy memory backptr+12+(y*screen width()+x)*4,g_srcMemPtr+12+(ty*g_bmpWidth+tx)*4,4 : endcase
            endselect
         endif
      next x
   next y

unlock backbuffer

rem swap height maps
   temp = g_ch
   g_ch = g_nh
   g_nh = temp
endfunction

+ Code Snippet

function UpdateHeights()
   local x, y,  height, temp

   lock backbuffer

   backptr=get backbuffer ptr()

   rem we need skip bitmap edges couse we're making kind of area sampling
   for y=2 to g_bmpHeight-3
      for x=2 to g_bmpWidth-3

         rem calculate height difference
         rem read 12 height values surrounding (x,y), divide by 6 so we've got twice as much as average value
         rem and then subtract current height at (x, y)
         height = g_heightMap(g_ch,x-1,y)
         inc height,g_heightMap(g_ch,x-2,y)+g_heightMap(g_ch,x+1,y)+g_heightMap(g_ch,x+2,y)+g_heightMap(g_ch,x,y-1)+g_heightMap(g_ch,x,y-2)+g_heightMap(g_ch,x,y+1)+g_heightMap(g_ch,x,y+2)+g_heightMap(g_ch,x-1,y-1)+g_heightMap(g_ch,x+1,y-1)+g_heightMap(g_ch,x-1,y+1)+g_heightMap(g_ch,x+1,y+1)
         height = int(height/6) - g_heightMap(g_nh,x,y)
         rem add some dampig reducing amplitudes
         dec height,height/g_waveDamp + 1
         g_heightMap(g_nh, x, y) = height

         tx = x + (g_heightMap(g_ch, x-1, y) - height)/10: rem divide couse height values are too high
         ty = y + (g_heightMap(g_ch, x, y-1) - height)/10
         rem if (tx, ty) still inside bitmap bounds

rem HERE

         if tx>=0 and tx<g_bmpWidth and ty>=0 and ty<g_bmpHeight 
            rem read color from coords we've calculated
            select g_bytesPerPixel
               case 2 : copy memory backptr+12+(y*screen width()+x)*2,g_srcMemPtr+12+(ty*g_bmpWidth+tx)*2,2 : endcase
               case 4 : copy memory backptr+12+(y*screen width()+x)*4,g_srcMemPtr+12+(ty*g_bmpWidth+tx)*4,4 : endcase
            endselect
         endif
      next x
   next y

   unlock backbuffer

   rem swap height maps
   temp = g_ch
   g_ch = g_nh
   g_nh = temp
endfunction

I mean there souldnt be checked condition: and (tx<>x or ty<>y)

Back to top

Profile PM

D Man

21

Years of Service

User Offline

Joined: 3rd Oct 2002

Location: Germany

Posted: 1st Oct 2003 17:05

Link

Oh...you're right, but it is a little faster, when you don't syncronize what you don't have to.
It's better to use this function even if there is an unwanted effect at the start

.

function UpdateHeights()
   local x, y,  height, temp

lock backbuffer

backptr=get backbuffer ptr()

rem we need skip bitmap edges couse we're making kind of area sampling
   for y=2 to g_bmpHeight-3
      for x=2 to g_bmpWidth-3

rem calculate height difference
         rem read 12 height values surrounding (x,y), divide by 6 so we've got twice as much as average value
         rem and then subtract current height at (x, y)
         height = g_heightMap(g_ch,x-1,y)
         inc height,g_heightMap(g_ch,x-2,y)+g_heightMap(g_ch,x+1,y)+g_heightMap(g_ch,x+2,y)+g_heightMap(g_ch,x,y-1)+g_heightMap(g_ch,x,y-2)+g_heightMap(g_ch,x,y+1)+g_heightMap(g_ch,x,y+2)+g_heightMap(g_ch,x-1,y-1)+g_heightMap(g_ch,x+1,y-1)+g_heightMap(g_ch,x-1,y+1)+g_heightMap(g_ch,x+1,y+1)
         height = int(height/6) - g_heightMap(g_nh,x,y)
         rem add some dampig reducing amplitudes
         dec height,height/g_waveDamp + 1
         g_heightMap(g_nh, x, y) = height

tx = x + (g_heightMap(g_ch, x-1, y) - height)/10: rem divide couse height values are too high
         ty = y + (g_heightMap(g_ch, x, y-1) - height)/10
         rem if (tx, ty) still inside bitmap bounds
         if tx>=0 and tx<g_bmpWidth and ty>=0 and ty<g_bmpHeight and g_heightMap(g_ch, x-1, y)<>g_heightMap(g_ch, x, y-1)
         rem read color from coords we've calculated
            select g_bytesPerPixel
               case 2 : copy memory backptr+12+(y*screen width()+x)*2,g_srcMemPtr+12+(ty*g_bmpWidth+tx)*2,2 : endcase
               case 4 : copy memory backptr+12+(y*screen width()+x)*4,g_srcMemPtr+12+(ty*g_bmpWidth+tx)*4,4 : endcase
            endselect
         endif
      next x
   next y

unlock backbuffer

rem swap height maps
   temp = g_ch
   g_ch = g_nh
   g_nh = temp
endfunction

+ Code Snippet

function UpdateHeights()
   local x, y,  height, temp

   lock backbuffer

   backptr=get backbuffer ptr()

   rem we need skip bitmap edges couse we're making kind of area sampling
   for y=2 to g_bmpHeight-3
      for x=2 to g_bmpWidth-3

         rem calculate height difference
         rem read 12 height values surrounding (x,y), divide by 6 so we've got twice as much as average value
         rem and then subtract current height at (x, y)
         height = g_heightMap(g_ch,x-1,y)
         inc height,g_heightMap(g_ch,x-2,y)+g_heightMap(g_ch,x+1,y)+g_heightMap(g_ch,x+2,y)+g_heightMap(g_ch,x,y-1)+g_heightMap(g_ch,x,y-2)+g_heightMap(g_ch,x,y+1)+g_heightMap(g_ch,x,y+2)+g_heightMap(g_ch,x-1,y-1)+g_heightMap(g_ch,x+1,y-1)+g_heightMap(g_ch,x-1,y+1)+g_heightMap(g_ch,x+1,y+1)
         height = int(height/6) - g_heightMap(g_nh,x,y)
         rem add some dampig reducing amplitudes
         dec height,height/g_waveDamp + 1
         g_heightMap(g_nh, x, y) = height

         tx = x + (g_heightMap(g_ch, x-1, y) - height)/10: rem divide couse height values are too high
         ty = y + (g_heightMap(g_ch, x, y-1) - height)/10
         rem if (tx, ty) still inside bitmap bounds
         if tx>=0 and tx<g_bmpWidth and ty>=0 and ty<g_bmpHeight and g_heightMap(g_ch, x-1, y)<>g_heightMap(g_ch, x, y-1)
         rem read color from coords we've calculated
            select g_bytesPerPixel
               case 2 : copy memory backptr+12+(y*screen width()+x)*2,g_srcMemPtr+12+(ty*g_bmpWidth+tx)*2,2 : endcase
               case 4 : copy memory backptr+12+(y*screen width()+x)*4,g_srcMemPtr+12+(ty*g_bmpWidth+tx)*4,4 : endcase
            endselect
         endif
      next x
   next y

   unlock backbuffer

   rem swap height maps
   temp = g_ch
   g_ch = g_nh
   g_nh = temp
endfunction

God is real, unless declared integer.

Back to top

Profile PM Email

D Man

21

Years of Service

User Offline

Joined: 3rd Oct 2002

Location: Germany

Posted: 1st Oct 2003 17:12

Link

ok...forget about the last post.
There are more unwanted effects than wanted

.

God is real, unless declared integer.

Back to top

Profile PM Email

Hamish McHaggis

21

Years of Service

User Offline

Joined: 13th Dec 2002

Location: Modgnik Detinu

Posted: 1st Oct 2003 18:51

Link

OMG! That is an amazingly amazing effect. Looks so realistic. Great job

!

Brains are for idiots.

Athelon XP 1400 Plus - Nvidia Geforce MX400 - 256mb RAM

Back to top

Profile PM Email Website

skovron

20

Years of Service

User Offline

Joined: 14th Sep 2003

Location:

Posted: 1st Oct 2003 19:16

Link

Quote: "It's better to use this function even if there is an unwanted effect at the start "

I didnt mean the effect at the start but that when all waves stops, when image became "dissolved" (dontknow is it wright word), but now it dissapeared and it looks very good

and ofcourse compared to the first version it run almost like a rocket

Really good job D Man.

Back to top

Profile PM

D Man

21

Years of Service

User Offline

Joined: 3rd Oct 2002

Location: Germany

Posted: 1st Oct 2003 19:47

Link

Thank' you

, but you had the idea and I knew how to optimize it

.
Great job from you, too.

God is real, unless declared integer.

Back to top

Profile PM Email

skovron

20

Years of Service

User Offline

Joined: 14th Sep 2003

Location:

Posted: 2nd Oct 2003 10:31

Link

Thanks!
I've got this code i Code Base so I'll update it with this version and I'll give you credit (wright word I hope).
What do you think about?

Back to top

Profile PM

D Man

21

Years of Service

User Offline

Joined: 3rd Oct 2002

Location: Germany

Posted: 2nd Oct 2003 14:07

Link

I definately agree that you should update the entry in the Code Base but you don't have to give me credits, as this was your idea, but if you wan't I've nothing against it

.

God is real, unless declared integer.

Back to top

Profile PM Email

skovron

20

Years of Service

User Offline

Joined: 14th Sep 2003

Location:

Posted: 3rd Oct 2003 13:43

Link

Well I've updated the simulation with 2 more effects.

Now under keys 1, 2, 3 you've got three effects
1) old one (move mouse and press lmb)
2) something I called rain
3) similar to the 2nd but drops are bigger and slower

here's the code

remstart
   Water Wave Effect
   by Mariusz Skowronski (2003)
   optimized by D Man
remend

#constant EFFECT_PATH 1
#constant EFFECT_RAIN1 2
#constant EFFECT_RAIN2 3

disable escapekey

rem global vars
rem bitmap size and color depth
global g_bmpWidth as integer
global g_bmpHeight as integer
global g_bmpBpp as integer
global g_bytesPerPixel as integer

rem memblock pointers (for faster access)
global g_srcMemPtr as dword
global g_dstMemPtr as dword

rem screen dimensions
global g_scrWidth as integer
global g_scrHeight as integer

rem height array swap index (current and new)
global g_ch as integer
global g_nh as integer

global g_effect as integer

rem properties
rem wave amplitude
global g_waveForce as integer
rem drop size
global g_dropRadius as integer
rem damping value (the smaller value the faster waves disappear) 1 - 100
global g_waveDamp as integer

rem init vars
g_scrWidth = 640
g_scrHeight = 480

g_waveForce = 300
g_dropRadius = 8
g_waveDamp = 20 : rem DON'T USE ZERO

g_effect = EFFECT_PATH

rem begin
set display mode g_scrWidth, g_scrHeight, 32
sync on : sync rate 0

rem load bitmap (remember the smaller bitmap the faster it runs)
load bitmap "example1.jpg", 1
g_bmpWidth = bitmap width(1)
g_bmpHeight = bitmap height(1)
g_bmpBpp = bitmap depth(1)
g_bytesPerPixel = g_bmpBpp/8

rem create bitmap memblock
make memblock from bitmap 2, 1

rem get memblock addresses
g_srcMemPtr = get memblock ptr(2)

rem dont need this bitmap anymore
delete bitmap 1

rem allocate 2 height maps
dim g_heightMap(1, g_bmpWidth-1, g_bmpHeight-1) as integer
rem clear height maps
ClearHeights()

rem set current height map to 0
g_ch = 0: g_nh = 1

`make image from memblock 1, 2
`paste image 1,0,0
sync
`delete image 1

set text opaque
randomize 100
lastTime = 0
do
   UpdateHeights()

select g_effect
      case EFFECT_PATH:
         if mouseclick() then MakeWave(mousex(), mousey(), g_dropRadius, g_waveForce)
      endcase

case EFFECT_RAIN1:
            MakeWave(rnd(g_bmpWidth), rnd(g_bmpHeight), rnd(2)+1, g_waveForce)
            MakeWave(rnd(g_bmpWidth), rnd(g_bmpHeight), rnd(2)+1, g_waveForce)
      endcase

case EFFECT_RAIN2:
         if timer()-lastTime > 2000
            lastTime = timer()
            MakeWave(rnd(g_bmpWidth), rnd(g_bmpHeight), rnd(6)+8, g_waveForce)
         endif
      endcase
   endselect

char$ = inkey$()
   select char$
      case "1"
          g_effect = EFFECT_PATH
          g_waveForce = 300 : g_dropRadius = 8 : g_waveDamp = 20
      endcase
      case "2"
         g_effect = EFFECT_RAIN1
         g_waveForce = 800 : g_waveDamp = 6
      endcase
      case "3"
         g_effect = EFFECT_RAIN2
         lastTime = timer()
         g_waveForce = 200 : g_waveDamp = 30
      endcase
   endselect

if escapekey() then exit

text screen width()-text width("FPS: "+str$(screen fps())), 10,"FPS: "+str$(screen fps())
   sync
loop

delete memblock 2
rem that's all
end

rem this clears height maps
function ClearHeights()
   local x, y
   for y=0 to g_bmpHeight-1
      for x=0 to g_bmpWidth-1
         g_heightMap(0, x, y) = 0
         g_heightMap(1, x, y) = 0
      next x
   next y
endfunction

rem this makes a new wave
function MakeWave(x, y, radius, mul)
   local u, v, sqru, sqrv, sqrr

sqrr = radius*radius
   rem check if not out of bitmap boundary
   if (x>radius) and (x<g_bmpWidth-1-radius) and (y>radius) and (y<g_bmpHeight-1-radius)
      for v=y-radius to y+radius
         sqrv = (v-y)*(v-y)
         for u=x-radius to x+radius
            sqru = (u-x)*(u-x)
            rem for speed reasons check (dx^2 + dy^2) <= r^2
            rem instead of real formula: sqrt(dx^2 + dy^) < r
            if (sqru+sqrv)<=sqrr:
               g_heightMap(g_ch,u,v) = mul*int(radius - sqrt(sqru+sqrv)+1)
            endif
         next u
      next v
   endif

endfunction

rem this is only for debug (shows waves in grayscale)
function DrawHeights()
   local x, y, height, color
   local dstAddr as dword
   rem for each height
   for y=0 to g_bmpHeight-1
      for x=0 to g_bmpWidth-1
         height = g_heightMap(g_ch, x, y)/10
         rem prepare gray color
         color = rgb(height+128, height+128, height+128)
         rem dont need this if useing one color depth
         select g_bytesPerPixel
            case 2: dstAddr = g_dstMemPtr +12+(y*g_bmpWidth+x)*2: endcase
            case 4: dstAddr = g_dstMemPtr +12+(y*g_bmpWidth+x)*4: endcase
         endselect
         *dstAddr = color
      next x
   next y
endfunction

rem the meat. this updates height map
function UpdateHeights()
   local x, y,  height, temp

lock backbuffer

backptr=get backbuffer ptr()

rem we need skip bitmap edges couse we're making kind of area sampling
   for y=2 to g_bmpHeight-3
      for x=2 to g_bmpWidth-3

rem calculate height difference
         rem read 12 height values surrounding (x,y), divide by 6 so we've got twice as much as average value
         rem and then subtract current height at (x, y)
         height = g_heightMap(g_ch,x-1,y)
         inc height,g_heightMap(g_ch,x-2,y)+g_heightMap(g_ch,x+1,y)+g_heightMap(g_ch,x+2,y)+g_heightMap(g_ch,x,y-1)+g_heightMap(g_ch,x,y-2)+g_heightMap(g_ch,x,y+1)+g_heightMap(g_ch,x,y+2)+g_heightMap(g_ch,x-1,y-1)+g_heightMap(g_ch,x+1,y-1)+g_heightMap(g_ch,x-1,y+1)+g_heightMap(g_ch,x+1,y+1)
         height = int(height/6) - g_heightMap(g_nh,x,y)
         rem add some dampig reducing amplitudes
         dec height,height/g_waveDamp + 1
         g_heightMap(g_nh, x, y) = height

tx = x + (g_heightMap(g_ch, x-1, y) - height)/10: rem divide couse height values are too high
         ty = y + (g_heightMap(g_ch, x, y-1) - height)/10
         rem if (tx, ty) still inside bitmap bounds
         `if tx>=0 and tx<g_bmpWidth and ty>=0 and ty<g_bmpHeight and g_heightMap(g_ch, x-1, y)<>g_heightMap(g_ch, x, y-1)
         if tx>=0 and tx<g_bmpWidth and ty>=0 and ty<g_bmpHeight and g_heightMap(g_ch, x, y)<>g_heightMap(g_nh, x, y)
         rem read color from coords we've calculated
            select g_bytesPerPixel
               case 2 : copy memory backptr+12+(y*screen width()+x)*2,g_srcMemPtr+12+(ty*g_bmpWidth+tx)*2,2 : endcase
               case 4 : copy memory backptr+12+(y*screen width()+x)*4,g_srcMemPtr+12+(ty*g_bmpWidth+tx)*4,4 : endcase
            endselect
         endif
      next x
   next y

unlock backbuffer

rem swap height maps
   temp = g_ch
   g_ch = g_nh
   g_nh = temp
endfunction

+ Code Snippet

remstart
   Water Wave Effect
   by Mariusz Skowronski (2003)
   optimized by D Man
remend

#constant EFFECT_PATH 1
#constant EFFECT_RAIN1 2
#constant EFFECT_RAIN2 3

disable escapekey

rem global vars
rem bitmap size and color depth
global g_bmpWidth as integer
global g_bmpHeight as integer
global g_bmpBpp as integer
global g_bytesPerPixel as integer

rem memblock pointers (for faster access)
global g_srcMemPtr as dword
global g_dstMemPtr as dword

rem screen dimensions
global g_scrWidth as integer
global g_scrHeight as integer

rem height array swap index (current and new)
global g_ch as integer
global g_nh as integer

global g_effect as integer

rem properties
rem wave amplitude
global g_waveForce as integer
rem drop size
global g_dropRadius as integer
rem damping value (the smaller value the faster waves disappear) 1 - 100
global g_waveDamp as integer


rem init vars
g_scrWidth = 640
g_scrHeight = 480

g_waveForce = 300
g_dropRadius = 8
g_waveDamp = 20 : rem DON'T USE ZERO

g_effect = EFFECT_PATH

rem begin
set display mode g_scrWidth, g_scrHeight, 32
sync on : sync rate 0

rem load bitmap (remember the smaller bitmap the faster it runs)
load bitmap "example1.jpg", 1
g_bmpWidth = bitmap width(1)
g_bmpHeight = bitmap height(1)
g_bmpBpp = bitmap depth(1)
g_bytesPerPixel = g_bmpBpp/8


rem create bitmap memblock
make memblock from bitmap 2, 1


rem get memblock addresses
g_srcMemPtr = get memblock ptr(2)

rem dont need this bitmap anymore
delete bitmap 1

rem allocate 2 height maps
dim g_heightMap(1, g_bmpWidth-1, g_bmpHeight-1) as integer
rem clear height maps
ClearHeights()

rem set current height map to 0
g_ch = 0: g_nh = 1

`make image from memblock 1, 2
`paste image 1,0,0
sync
`delete image 1

set text opaque
randomize 100
lastTime = 0
do
   UpdateHeights()

   select g_effect
      case EFFECT_PATH:
         if mouseclick() then MakeWave(mousex(), mousey(), g_dropRadius, g_waveForce)
      endcase

      case EFFECT_RAIN1:
            MakeWave(rnd(g_bmpWidth), rnd(g_bmpHeight), rnd(2)+1, g_waveForce)
            MakeWave(rnd(g_bmpWidth), rnd(g_bmpHeight), rnd(2)+1, g_waveForce)
      endcase

      case EFFECT_RAIN2:
         if timer()-lastTime > 2000
            lastTime = timer()
            MakeWave(rnd(g_bmpWidth), rnd(g_bmpHeight), rnd(6)+8, g_waveForce)
         endif
      endcase
   endselect

   char$ = inkey$()
   select char$
      case "1"
          g_effect = EFFECT_PATH
          g_waveForce = 300 : g_dropRadius = 8 : g_waveDamp = 20
      endcase
      case "2"
         g_effect = EFFECT_RAIN1
         g_waveForce = 800 : g_waveDamp = 6
      endcase
      case "3"
         g_effect = EFFECT_RAIN2
         lastTime = timer()
         g_waveForce = 200 : g_waveDamp = 30
      endcase
   endselect

   if escapekey() then exit

   text screen width()-text width("FPS: "+str$(screen fps())), 10,"FPS: "+str$(screen fps())
   sync
loop

delete memblock 2
rem that's all
end

rem this clears height maps
function ClearHeights()
   local x, y
   for y=0 to g_bmpHeight-1
      for x=0 to g_bmpWidth-1
         g_heightMap(0, x, y) = 0
         g_heightMap(1, x, y) = 0
      next x
   next y
endfunction


rem this makes a new wave
function MakeWave(x, y, radius, mul)
   local u, v, sqru, sqrv, sqrr

   sqrr = radius*radius
   rem check if not out of bitmap boundary
   if (x>radius) and (x<g_bmpWidth-1-radius) and (y>radius) and (y<g_bmpHeight-1-radius)
      for v=y-radius to y+radius
         sqrv = (v-y)*(v-y)
         for u=x-radius to x+radius
            sqru = (u-x)*(u-x)
            rem for speed reasons check (dx^2 + dy^2) <= r^2
            rem instead of real formula: sqrt(dx^2 + dy^) < r
            if (sqru+sqrv)<=sqrr:
               g_heightMap(g_ch,u,v) = mul*int(radius - sqrt(sqru+sqrv)+1)
            endif
         next u
      next v
   endif

endfunction

rem this is only for debug (shows waves in grayscale)
function DrawHeights()
   local x, y, height, color
   local dstAddr as dword
   rem for each height
   for y=0 to g_bmpHeight-1
      for x=0 to g_bmpWidth-1
         height = g_heightMap(g_ch, x, y)/10
         rem prepare gray color
         color = rgb(height+128, height+128, height+128)
         rem dont need this if useing one color depth
         select g_bytesPerPixel
            case 2: dstAddr = g_dstMemPtr +12+(y*g_bmpWidth+x)*2: endcase
            case 4: dstAddr = g_dstMemPtr +12+(y*g_bmpWidth+x)*4: endcase
         endselect
         *dstAddr = color
      next x
   next y
endfunction

rem the meat. this updates height map
function UpdateHeights()
   local x, y,  height, temp

   lock backbuffer

   backptr=get backbuffer ptr()

   rem we need skip bitmap edges couse we're making kind of area sampling
   for y=2 to g_bmpHeight-3
      for x=2 to g_bmpWidth-3

         rem calculate height difference
         rem read 12 height values surrounding (x,y), divide by 6 so we've got twice as much as average value
         rem and then subtract current height at (x, y)
         height = g_heightMap(g_ch,x-1,y)
         inc height,g_heightMap(g_ch,x-2,y)+g_heightMap(g_ch,x+1,y)+g_heightMap(g_ch,x+2,y)+g_heightMap(g_ch,x,y-1)+g_heightMap(g_ch,x,y-2)+g_heightMap(g_ch,x,y+1)+g_heightMap(g_ch,x,y+2)+g_heightMap(g_ch,x-1,y-1)+g_heightMap(g_ch,x+1,y-1)+g_heightMap(g_ch,x-1,y+1)+g_heightMap(g_ch,x+1,y+1)
         height = int(height/6) - g_heightMap(g_nh,x,y)
         rem add some dampig reducing amplitudes
         dec height,height/g_waveDamp + 1
         g_heightMap(g_nh, x, y) = height

         tx = x + (g_heightMap(g_ch, x-1, y) - height)/10: rem divide couse height values are too high
         ty = y + (g_heightMap(g_ch, x, y-1) - height)/10
         rem if (tx, ty) still inside bitmap bounds
         `if tx>=0 and tx<g_bmpWidth and ty>=0 and ty<g_bmpHeight and g_heightMap(g_ch, x-1, y)<>g_heightMap(g_ch, x, y-1)
         if tx>=0 and tx<g_bmpWidth and ty>=0 and ty<g_bmpHeight and g_heightMap(g_ch, x, y)<>g_heightMap(g_nh, x, y)
         rem read color from coords we've calculated
            select g_bytesPerPixel
               case 2 : copy memory backptr+12+(y*screen width()+x)*2,g_srcMemPtr+12+(ty*g_bmpWidth+tx)*2,2 : endcase
               case 4 : copy memory backptr+12+(y*screen width()+x)*4,g_srcMemPtr+12+(ty*g_bmpWidth+tx)*4,4 : endcase
            endselect
         endif
      next x
   next y

   unlock backbuffer

   rem swap height maps
   temp = g_ch
   g_ch = g_nh
   g_nh = temp
endfunction

and the bitmap:
http://www.arena.free-host.com/examples/example1.jpg

again: enjoy

Back to top

Profile PM

D Man

21

Years of Service

User Offline

Joined: 3rd Oct 2002

Location: Germany

Posted: 4th Oct 2003 00:21

Link

Good idea and thanks for the credits

.

God is real, unless declared integer.

Back to top

Profile PM Email

John H

Retired Moderator

21

Years of Service

User Offline

Joined: 14th Oct 2002

Location: Burlington, VT

Posted: 5th Oct 2003 00:25

Link

Great!

Current Project: Eternal Destiny
We need a texture artist!
If your interested email me - Johnrocs@aol.com

Back to top

Profile PM

John H

Retired Moderator

21

Years of Service

User Offline

Joined: 14th Oct 2002

Location: Burlington, VT

Posted: 5th Oct 2003 00:29

Link

Hey I made a screenshot for all you guys

Great effect, this should encourage you to try it out

Effect: 'Rain Drops'

RPGamer

Current Project: Eternal Destiny
We need a texture artist!
If your interested email me - Johnrocs@aol.com

Back to top

Profile PM

skovron

20

Years of Service

User Offline

Joined: 14th Sep 2003

Location:

Posted: 5th Oct 2003 00:58

Link

Hey thanx for your promotion

I thought this topic will die soon but.... now

Back to top

Profile PM

John H

Retired Moderator

21

Years of Service

User Offline

Joined: 14th Oct 2002

Location: Burlington, VT

Posted: 5th Oct 2003 03:11

Link

No prob Skovron

Just promoting a great snippet

Current Project: Eternal Destiny
We need a texture artist!
If your interested email me - Johnrocs@aol.com

Back to top

Profile PM

Mike Inel

21

Years of Service

User Offline

Joined: 14th Feb 2003

Location: Sa upuan ko po...

Posted: 7th Oct 2003 13:12

Link

"Simple Water Wave Effect"... It looks VERY complicated to me.....

Uh.....duh.....whut?

Back to top

Profile PM Email Website

skovron

20

Years of Service

User Offline

Joined: 14th Sep 2003

Location:

Posted: 8th Oct 2003 00:35

Link

I'm not good at english.

But still it's simple not becouse the code is easy to understand without comment but simple becouse it isnt accurate. I mean the algorithm rely on "simple" approximation instead of ray-tracing.

I dont know your programming skil but belive me that many programmers which know how to generate waves using grids would say that it is very easy

On the other side I agree that it is always hard to understand code written by somebody else especially when its not comment well (as my code is)

Back to top

Profile PM

Imperial Darksith

20

Years of Service

User Offline

Joined: 29th Aug 2003

Location: The StarForge

Posted: 17th Oct 2003 12:55

Link

I noticed that when playing KOTOR that they use something rather interesting... I believe that it's a bitmap that's ghosted and scrolled diagonally creating a wave effect. Looks really realistic too. I suspect that your code runs something a bit like this. Am I correct?

I find your lack of faith...disturbing.

Back to top

Profile PM Email

D Man

21

Years of Service

User Offline

Joined: 3rd Oct 2002

Location: Germany

Posted: 17th Oct 2003 13:07

Link

No, you're wrong

.

God is real, unless declared integer.

Back to top

Profile PM Email

skovron

20

Years of Service

User Offline

Joined: 14th Sep 2003

Location:

Posted: 18th Oct 2003 15:18

Link

Youre wrong. No ghosted and scrolled bitmaps.
It just for every pixel checks height difference (along x and y axis) then this differences added to pixel coordinates gives new pixel position and that pixel is being copied from original bitmap

take a look at pseudo code:

+ Code Snippet

for x=1 to width
 for y=1 to height
   difference_x = height_map(x-1, y) - height_map(x, y)
   difference_y = height_map(x, y-1) - height_map(x, y)
   color = get pixel from bitmap(x+difference_x, y+difference_y)
   put pixel on screen(x, y, color)
 next y
next x

Back to top

Profile PM

Lenny Briscoe

21

Years of Service

User Offline

Joined: 15th Sep 2002

Location:

Posted: 19th Oct 2003 23:31

Link

That's a nice effect. Have to experiment with it.

Lenny

Back to top

Profile PM

Dreamora

21

Years of Service

User Offline

Joined: 20th Sep 2002

Location: Switzerland

Posted: 26th Oct 2003 12:26

Link

remembers the effect that you can play with while WaveRace GC is loading the track

Dell Inspiron 8200,P4-M 1,7ghz, GF4 440go (45.23 - DX9b), 1GB RAM

Back to top

Profile PM Email Website

DaemoN_v2

20

Years of Service

User Offline

Joined: 13th Oct 2003

Location: Wichita

Posted: 3rd Nov 2003 17:26

Link

niiiiiiiiiiiice

[tT]

Back to top

Profile PM Email

Sorry your browser is not supported!

Code Snippets / Simple approach to simulate water effect