Log inRegister an accountBrowse CSDbHelp & documentationFacts & StatisticsThe forumsAvailable RSS-feeds on CSDbSupport CSDb Commodore 64 Scene Database
You are not logged in - nap
CSDb User Forums


Forums > C64 Coding > CL13 spring water- how ?
2015-08-04 13:10
Oswald

Registered: Apr 2002
Posts: 5020
CL13 spring water- how ?

anyone would care to share how the spring math works in there?

I've tried yesterday to make a simple spring, but it very rarely got close to resting, 8 out of 10 tweaks it "blows" up due to something google said is a problem with the euler equation. :)

so yes I'm trying to be lazy here and asking for the solution right away. Arok is less than 2 weeks away and I need to finish this fast :)

I've found verlet method, but dont understand it in the depth thats enough to turn it into 8 bit code.

for reference:

http://gamedevelopment.tutsplus.com/tutorials/make-a-splash-wit..

thanks,
2015-08-04 15:02
Axis/Oxyron

Registered: Apr 2007
Posts: 91
Hey Oswald,

the algorithm is relatively simple but alot of fiddling until it looks nice. The basic idea is to track height and velocity of every column. Every column is connected to its 4 neighbours to the left and the four neighbours to the right with dynamic springs (only for the y-axis). The stiffness of the 8 springs (factor on the applied force) is calculated with a gaussian weight table. The last step is mostly what makes it look smooth. I can give you a C-source from my prototype. But I guess it won´t help much. Its optimized to the brim and completely based on 8 and 16 bit values, because I planned to do the effect in realtime. So changing a single value, will blow everything up! ;o)
2015-08-04 15:52
Oswald

Registered: Apr 2002
Posts: 5020
thanks!

maybe with some tweakage I can get the euler method working, only I didnt expect this to take days, rather only some hours to get right.

this is what I do, the problem is that it almost never gets close to resting. maybe using a damp table where x is not divided with 2^x will help.

springtest	
        lda #40
	sta springhh
	sta springhl ;spring height
	lda #00
	sta springsl
	sta springsh ;spring speed
	
u	lda #$40
-	cmp $d012
	bne -

	ldx springhh ;spring height HI byte
	lda dampl,x  ;0 height is assumed resting
	clc          ;damp table is 16 bit, x/2^7 signed used to change velocity
	adc springsl 
	sta springsl ;spring speed LO
	lda damph,x
	adc springsh 
	sta springsh ;spring speed HI
	

	lda springhl
	clc
	adc springsl
	sta springhl ;spring height + spring speed 16 bit
	lda springhh
	adc springsh
	sta springhh
	sta $d000     ;to see how it moves
        jmp u
2015-08-04 16:43
Axis/Oxyron

Registered: Apr 2007
Posts: 91
Basically you are on the right way. The damping table is the place for tweaking. But I´m not sure what you want to reach. Only moving a rasterbar or sprite on a wavy motion or building a 2d-waveshape like in CL13? For reaching the latter, you definitely need to take the neighbours into account.
2015-08-04 16:56
Oswald

Registered: Apr 2002
Posts: 5020
well, first I need to get 1 spring right, only after that can I move on to spread the forces.

I hope the problem is that I only use a x/2^7 damping table, and dividing by some other value or making it non linear will help.
2015-08-04 23:12
Copyfault

Registered: Dec 2001
Posts: 466
Never did such a routine before, but sounds interesting for sure.

If I understand the article linked in the original posting _and_ your code snipplet correctly it seems that
i) the dampl/h-tables you already have in the code are actually acceleration tables
ii) you need to add real damping depending on the current velocity

something like
u	lda #$40
-	cmp $d012
	bne -

	ldy springhh ;spring height HI byte
	
lda acc_lo,x  ;0 height is assumed resting
	clc          ;damp table is 16 bit, x/2^7 signed used to change velocity
	adc springsl 
	sta springsl ;spring speed LO
	lda acc_hi,x
	adc springsh 
	sta springsh ;spring speed HI
	

	lda springhl
	clc
	adc springsl
	sta springhl ;spring height + spring speed 16 bit
	lda springhh
	adc springsh
	sta springhh
	sta $d000     ;to see how it moves
        jmp u

maybe already does the trick...
2015-08-04 23:35
Copyfault

Registered: Dec 2001
Posts: 466
sorry, I f**ked up my last posting, at least the "code part".

Another try ;))
u	 lda #$40
-	 cmp $d012
	 bne -

 	ldy springhh ;spring height HI byte
	ldx springsh ;speed value hi byte for damp calculation
        lda acc_lo,y  ;0 height is assumed resting
	sec          ;damp table is 16 bit, x/2^7 signed used to change velocity
        sbc damp_lo,x ;formula says "acc = k*x - d*speed"
	adc springsl
        sta springsl ;spring speed LO
	lda acc_hi,y
	adc springsh
	sbc damp_hi,x
        sta springsh ;spring speed HI
	
 	lda springhl
	clc
	adc springsl
	sta springhl ;spring height + spring speed 16 bit
	lda springhh
	adc springsh
	sta springhh
	sta $d000     ;to see how it moves
        jmp u

Maybe your damp-table is a combination of both acc_lo/hi and damp_lo/hi, but how can you get the dependencies right with only one table?

Guess Axis knows best, he's gone through it already.

Let's see where this leads to, I'm alread looking forward to this year's Arok party;)
2015-08-05 06:31
JackAsser

Registered: Jun 2002
Posts: 1989
Bah, just animate it using delta compressed motion on each column. :)
2015-08-05 06:50
Oswald

Registered: Apr 2002
Posts: 5020
copyfault,

formula says:

float x = Height - TargetHeight;
float acceleration = -k * x;

Position += Velocity;
Velocity += acceleration;

I took targetheight to be 0, so acceleration = damptable(height)



Position += Velocity; this is the last part in my code
Velocity += acceleration; this is the adc damp,x part

aha gotcha ! I mixed things up, I wasnt doing any dampening at all, that was just the normal spring part.

thanks! :)
2015-08-05 07:10
Axis/Oxyron

Registered: Apr 2007
Posts: 91
In C my code looks like:
velocity=velocities[x];
for (int dx=-4 ; dx<=4 ; dx++)
{
  delta=height[x+dx]-waterplane;
  velocity-=delta*weight[dx];
}
velocities[x]=velocity;
height[x]+=velocity*damping;

In my case weight is a gaussian distance table. If you only consider a single value its just a constant that can be used for tweaking. Perhaps it can even be optimized away, because my weight table sums up to 256.
Damping is a constant factor for tweaking motions.

So a simplified version considering only 1 column would look like this:
delta=height[x]-waterplane;
velocities[x]-=delta*weight;		;perhaps weight is 1 and can be removed
height[x]+=velocities[x]*damping;
2015-08-05 07:21
Oswald

Registered: Apr 2002
Posts: 5020
thanks that will nail it, basically I should've feed the speed into a table that always reduces it a "bit".
 
... 5 posts hidden. Click here to view all posts....
 
Previous - 1 | 2 - Next
RefreshSubscribe to this thread:

You need to be logged in to post in the forum.

Search the forum:
Search   for   in  
All times are CET.
Search CSDb
Advanced
Users Online
Mixer
Fred/Channel 4
Jazzcat/Onslaught
Sh0ckTr00per/Hokuto ..
Ghost/Quantum
Bieno/Commodore Plus
Guests online: 122
Top Demos
1 Next Level  (9.8)
2 Mojo  (9.7)
3 Coma Light 13  (9.7)
4 Edge of Disgrace  (9.6)
5 Comaland 100%  (9.6)
6 No Bounds  (9.6)
7 Uncensored  (9.6)
8 Wonderland XIV  (9.6)
9 Memento Mori  (9.6)
10 Bromance  (9.5)
Top onefile Demos
1 It's More Fun to Com..  (9.7)
2 Party Elk 2  (9.7)
3 Cubic Dream  (9.6)
4 Copper Booze  (9.5)
5 TRSAC, Gabber & Pebe..  (9.5)
6 Rainbow Connection  (9.5)
7 Wafer Demo  (9.5)
8 Dawnfall V1.1  (9.5)
9 Quadrants  (9.5)
10 Daah, Those Acid Pil..  (9.5)
Top Groups
1 Nostalgia  (9.3)
2 Oxyron  (9.3)
3 Booze Design  (9.3)
4 Censor Design  (9.3)
5 Crest  (9.3)
Top Musicians
1 Vincenzo  (9.8)
2 Rob Hubbard  (9.7)
3 Stinsen  (9.7)
4 Jeroen Tel  (9.6)
5 Linus  (9.6)

Home - Disclaimer
Copyright © No Name 2001-2024
Page generated in: 0.092 sec.