Vector Gravity Simulation in realtime (physics project)

Author

Message

Pincho Paxton

21

Years of Service

User Offline

Joined: 8th Dec 2002

Location:

Posted: 22nd Oct 2009 23:27 Edited at: 22nd Oct 2009 23:29

Link

I am currently working on a simulation of gravity. If you look at my picture, I need to create the hexagons from vectors so that when a planet is put in the middle of them the nearest one shrinks, and pulls all of the other hexagons towards it. Any ideas?

Back to top

Profile PM Email Website

Libervurto

17

Years of Service

User Offline

Joined: 30th Jun 2006

Location: On Toast

Posted: 23rd Oct 2009 02:24

Link

Just thinking out loud, if you use sin and cos to draw the hexagons maybe if you scale sin and cos depending on how far away they are from a planet then they will stretch like that.

TGC Forum - converting error messages into sarcasm since 2002.

Back to top

Profile PM Email

Pincho Paxton

21

Years of Service

User Offline

Joined: 8th Dec 2002

Location:

Posted: 23rd Oct 2009 02:51 Edited at: 23rd Oct 2009 02:54

Link

Will that work if I keep adding more, and more planets to the screen? So I get cross gravity forces, and combined gravity forces?

How about a moon in orbit animation?

Back to top

Profile PM Email Website

Pincho Paxton

21

Years of Service

User Offline

Joined: 8th Dec 2002

Location:

Posted: 23rd Oct 2009 02:59 Edited at: 23rd Oct 2009 03:16

Link

If you are interested.. take a look at this picture.. maybe you will know what I am doing.. maybe not. This is a genuine NASA image, that I have edited with a Hexagon grid.. I find it very interesting...

Wouldn't it be funny if Dark Basic found Dark Matter!

Back to top

Profile PM Email Website

Neuro Fuzzy

16

Years of Service

User Offline

Joined: 11th Jun 2007

Location:

Posted: 23rd Oct 2009 05:11

Link

I would start out just by making a hexagonal point cloud. Check out this code (it's in DBPro, i don't know anything about DarkBASIC, but what you're trying to do seems cool)

+ Code Snippet

type hex
	a as integer
	b as integer
	c as integer
	x as float
	y as float
	xvel as float
	yvel as float
endtype
dim hexagon(1000) as hex
global startx as float = 1.0
global starty as float = 1.0
global sidelength as float = 30.0
curx as float = startx-sidelength
cury as float = starty-sidelength
switch as integer = 0
switch2 as boolean = 0
for a=1 to 1000

	if switch=1
		switch=0
		inc curx,2*sidelength
	else
		switch=1
		inc curx, sidelength
	endif
	hexagon(a).x=curx
	hexagon(a).y=cury
	if curx>screen width()
		if switch2=0 
			switch2=1
			switch=1
			inc cury,sidelength*sqrt(3)/2
			curx=startx
		else
			switch=0
			switch2=0
			inc cury,sidelength*sqrt(3)/2
			curx=startx+sidelength/2
		endif
	endif
next a
for a=1 to 1000
	dot hexagon(a).x,hexagon(a).y
next a
sync
wait key

in the code when you see "sidelength/2" and "sidelength*sqrt(3)/2", that's really just "sidelength*cos(60)" and "sidelength*sin(60)".

I'm not sure how to figure out which points connect to which. I'll think about it more, but it probably has to do with division, modulus, and some 1d to 2d conversions (cuz every point is stored in an array). I haven't really started to solve for that yet, but this whole thing seems like it could end up looking really cool XD

Back to top

Profile PM Email

Pincho Paxton

21

Years of Service

User Offline

Joined: 8th Dec 2002

Location:

Posted: 23rd Oct 2009 14:32 Edited at: 23rd Oct 2009 14:53

Link

Thanks! I have DBPro as well. Do you think that I can do something like this?......

set vector2,Map,hexagon(a).x,hexagon(a).y

Back to top

Profile PM Email Website

Pincho Paxton

21

Years of Service

User Offline

Joined: 8th Dec 2002

Location:

Posted: 23rd Oct 2009 16:34 Edited at: 23rd Oct 2009 16:45

Link

Nearly have it, I just need to figure out the diagonals, which are an offset of X/Y. I'm not sure when the code drops down a line.

set display mode 1152,864,16
Map = 1
type hex
   a as integer
   b as integer
   c as integer
   x as float
   y as float
   xvel as float
   yvel as float
endtype
dim hexagon(1000) as hex
Dim map(1000,1000)
global startx as float = 1.0
global starty as float = 1.0
global sidelength as float = 30.0
curx as float = startx-sidelength
cury as float = starty-sidelength
switch as integer = 0
switch2 as boolean = 0
for a=1 to 1000

if switch=1
      switch=0
      inc curx,2*sidelength
   else
      switch=1
      inc curx, sidelength
   endif
   hexagon(a).x=curx
   hexagon(a).y=cury
   if curx>screen width()
      if switch2=0
         switch2=1
         switch=1
         inc cury,sidelength*sqrt(3)/2
         curx=startx
      else
         switch=0
         switch2=0
         inc cury,sidelength*sqrt(3)/2
         curx=startx+sidelength/2
      endif
   endif
next a
for a=1 to 999 step 2
   dot hexagon(a).x,hexagon(a).y
   `set vector2,Map,hexagon(a).x,hexagon(a).y
   line hexagon(a).x,hexagon(a).y,hexagon(a+1).x, hexagon(a+1).y
next a
sync
wait key

+ Code Snippet

set display mode 1152,864,16
Map = 1
type hex
   a as integer
   b as integer
   c as integer
   x as float
   y as float
   xvel as float
   yvel as float
endtype
dim hexagon(1000) as hex
Dim map(1000,1000)
global startx as float = 1.0
global starty as float = 1.0
global sidelength as float = 30.0
curx as float = startx-sidelength
cury as float = starty-sidelength
switch as integer = 0
switch2 as boolean = 0
for a=1 to 1000

   if switch=1
      switch=0
      inc curx,2*sidelength
   else
      switch=1
      inc curx, sidelength
   endif
   hexagon(a).x=curx
   hexagon(a).y=cury
   if curx>screen width()
      if switch2=0
         switch2=1
         switch=1
         inc cury,sidelength*sqrt(3)/2
         curx=startx
      else
         switch=0
         switch2=0
         inc cury,sidelength*sqrt(3)/2
         curx=startx+sidelength/2
      endif
   endif
next a
for a=1 to 999 step 2
   dot hexagon(a).x,hexagon(a).y
   `set vector2,Map,hexagon(a).x,hexagon(a).y
   line hexagon(a).x,hexagon(a).y,hexagon(a+1).x, hexagon(a+1).y
next a
sync
wait key

Ok I am moving this to DB Pro, because it is now a DB Pro project.

Back to top

Profile PM Email Website

Neuro Fuzzy

16

Years of Service

User Offline

Joined: 11th Jun 2007

Location:

Posted: 24th Oct 2009 09:41

Link

type hex
	a as integer
	b as integer
	c as integer
	x as float
	y as float
	xvel as float
	yvel as float
endtype
dim hexagon(1000) as hex
global startx as float = 1.0
global starty as float = 1.0
global sidelength as float = 30.0
curx as float = startx-sidelength
cury as float = starty-sidelength
switch as integer = 0
switch2 as boolean = 0
global exn as float = 700.0

do
cls
	exn=mousex()
	starty=mousey()
	curx = startx-sidelength*2
	cury = starty
	switch as integer = 1
	switch2 as boolean = 1
	for a=1 to 1000
	
		if switch=1
			switch=0
			inc curx,2*sidelength
		else
			switch=1
			inc curx, sidelength
		endif

if curx>exn
			if switch2=0 
				switch2=1
				switch=0
				inc cury,sidelength*sqrt(3)/2
				curx=startx
			else
				switch=1
				switch2=0
				inc cury,sidelength*sqrt(3)/2
				curx=startx-sidelength/2
			endif
		endif
		hexagon(a).x=curx
		hexagon(a).y=cury		
	next a
	for a=1 to 1000
		dx#=0
		dy#=0
		if spacekey()
			f#=1/sqrt((hexagon(a).x-mousex())^2+(hexagon(a).y-mousey())^2)`*rnd(1)
			dx#=(hexagon(a).x-mousex())*f#*100
			dy#=(hexagon(a).y-mousey())*f#*100
		endif
		dot hexagon(a).x+dx#,hexagon(a).y+dy#
		
	next a
	a=get_length(exn,sidelength)
	b=get_length(exn-sidelength*1.5,sidelength)
	text 0,0,"In row 1, there are "+str$(a)+" points"
	text 0,16,"In row 2, there are "+str$(b)+" points"
	sync
	
loop
	
function get_length(width,length)
 	a = width/(3.0*length)
	b# = (width/(3.0*length)-a)*3
	a=a*2
	if b#>=0 then inc a
	if b#>=1 then inc a
endfunction a

+ Code Snippet

type hex
	a as integer
	b as integer
	c as integer
	x as float
	y as float
	xvel as float
	yvel as float
endtype
dim hexagon(1000) as hex
global startx as float = 1.0
global starty as float = 1.0
global sidelength as float = 30.0
curx as float = startx-sidelength
cury as float = starty-sidelength
switch as integer = 0
switch2 as boolean = 0
global exn as float = 700.0

do
cls
	exn=mousex()
	starty=mousey()
	curx = startx-sidelength*2
	cury = starty
	switch as integer = 1
	switch2 as boolean = 1
	for a=1 to 1000
	
		if switch=1
			switch=0
			inc curx,2*sidelength
		else
			switch=1
			inc curx, sidelength
		endif

		if curx>exn
			if switch2=0 
				switch2=1
				switch=0
				inc cury,sidelength*sqrt(3)/2
				curx=startx
			else
				switch=1
				switch2=0
				inc cury,sidelength*sqrt(3)/2
				curx=startx-sidelength/2
			endif
		endif
		hexagon(a).x=curx
		hexagon(a).y=cury		
	next a
	for a=1 to 1000
		dx#=0
		dy#=0
		if spacekey()
			f#=1/sqrt((hexagon(a).x-mousex())^2+(hexagon(a).y-mousey())^2)`*rnd(1)
			dx#=(hexagon(a).x-mousex())*f#*100
			dy#=(hexagon(a).y-mousey())*f#*100
		endif
		dot hexagon(a).x+dx#,hexagon(a).y+dy#
		
	next a
	a=get_length(exn,sidelength)
	b=get_length(exn-sidelength*1.5,sidelength)
	text 0,0,"In row 1, there are "+str$(a)+" points"
	text 0,16,"In row 2, there are "+str$(b)+" points"
	sync
	
loop
	
function get_length(width,length)
 	a = width/(3.0*length)
	b# = (width/(3.0*length)-a)*3
	a=a*2
	if b#>=0 then inc a
	if b#>=1 then inc a
endfunction a

I've done a bit more, i smoothed out some kinks in the point cloud generator, and figured out how to get the number of points in one row. It was pretty hard for me ta figure out, i just guess I'm not used to geometry and integers :\

Again, the above code is in DBPro. I figure it doesn't really matter what language it's coded in. Once you know how to do it, you can port it to darkbasic, c++, java, flash, whatever.

Still, now that i know how to find the number of dots in a row and generate a point cloud, i can start... err... connecting the dots.

Back to top

Profile PM Email

Neuro Fuzzy

16

Years of Service

User Offline

Joined: 11th Jun 2007

Location:

Posted: 24th Oct 2009 10:37

Link

Wait, scratch that, not all kinks are worked out. There are actually points generated offscreen to the left of the screen

Well, i'm going ta bed, i think i should have the lines done by tomorrow, and that'll be cool!

Back to top

Profile PM Email

Pincho Paxton

21

Years of Service

User Offline

Joined: 8th Dec 2002

Location:

Posted: 24th Oct 2009 15:45 Edited at: 24th Oct 2009 15:53

Link

Thanks! I was just beginning to realise that I need the dots in rows. Starting to understand what I need now. I've never used dots, and lines before, and the new commands were making me forget what I am actually trying to achieve. On the far right I would need some drawing tools to add the atoms. I also want to see if snowflakes form naturally from gravity, so I need a temperature setting as well, and how this effects the distribution of particles in the atmosphere, and how this particle arrangement alters the hexagons around it. Basically, what I am saying is that I need a gap on the far right.

Back to top

Profile PM Email Website

Pincho Paxton

21

Years of Service

User Offline

Joined: 8th Dec 2002

Location:

Posted: 24th Oct 2009 17:30 Edited at: 24th Oct 2009 17:34

Link

Had a play with your new code. This is cool! Press space.

type hex
   a as integer
   b as integer
   c as integer
   x as float
   y as float
   xvel as float
   yvel as float
endtype
dim hexagon(1000) as hex
global startx as float = 1.0
global starty as float = 1.0
global sidelength as float = 30.0
curx as float = startx-sidelength
cury as float = starty-sidelength
switch as integer = 0
switch2 as boolean = 0
global exn as float = 700.0

do
cls
   exn=screen width()
   starty=Screen height()
   curx = startx-sidelength*2
   cury = 30
   switch as integer = 1
   switch2 as boolean = 1
   for a=1 to 1000

if switch=1
         switch=0
         inc curx,2*sidelength
      else
         switch=1
         inc curx, sidelength
      endif

if curx>exn
         if switch2=0
            switch2=1
            switch=0
            inc cury,sidelength*sqrt(3)/2
            curx=startx
         else
            switch=1
            switch2=0
            inc cury,sidelength*sqrt(3)/2
            curx=startx-sidelength/2
         endif
      endif
      hexagon(a).x=curx
      hexagon(a).y=cury
   next a
   for a=1 to 1000
      dx#=0
      dy#=0
      if spacekey()
         f#=1/sqrt((hexagon(a).x-mousex())^2+(hexagon(a).y-mousey())^2)`*rnd(1)
         dx#=(hexagon(a).x-mousex())*f#*100
         dy#=(hexagon(a).y-mousey())*f#*100
      endif
      dot hexagon(a).x+dx#,hexagon(a).y+dy#

next a
   `a=get_length(exn,sidelength)
   `b=get_length(exn-sidelength*1.5,sidelength)
   `text 0,0,"In row 1, there are "+str$(a)+" points"
   `text 0,16,"In row 2, there are "+str$(b)+" points"
   sync

loop

function get_length(width,length)
    a = width/(3.0*length)
   b# = (width/(3.0*length)-a)*3
   a=a*2
   if b#>=0 then inc a
   if b#>=1 then inc a
endfunction a

+ Code Snippet

type hex
   a as integer
   b as integer
   c as integer
   x as float
   y as float
   xvel as float
   yvel as float
endtype
dim hexagon(1000) as hex
global startx as float = 1.0
global starty as float = 1.0
global sidelength as float = 30.0
curx as float = startx-sidelength
cury as float = starty-sidelength
switch as integer = 0
switch2 as boolean = 0
global exn as float = 700.0

do
cls
   exn=screen width()
   starty=Screen height()
   curx = startx-sidelength*2
   cury = 30
   switch as integer = 1
   switch2 as boolean = 1
   for a=1 to 1000

      if switch=1
         switch=0
         inc curx,2*sidelength
      else
         switch=1
         inc curx, sidelength
      endif

      if curx>exn
         if switch2=0
            switch2=1
            switch=0
            inc cury,sidelength*sqrt(3)/2
            curx=startx
         else
            switch=1
            switch2=0
            inc cury,sidelength*sqrt(3)/2
            curx=startx-sidelength/2
         endif
      endif
      hexagon(a).x=curx
      hexagon(a).y=cury
   next a
   for a=1 to 1000
      dx#=0
      dy#=0
      if spacekey()
         f#=1/sqrt((hexagon(a).x-mousex())^2+(hexagon(a).y-mousey())^2)`*rnd(1)
         dx#=(hexagon(a).x-mousex())*f#*100
         dy#=(hexagon(a).y-mousey())*f#*100
      endif
      dot hexagon(a).x+dx#,hexagon(a).y+dy#

   next a
   `a=get_length(exn,sidelength)
   `b=get_length(exn-sidelength*1.5,sidelength)
   `text 0,0,"In row 1, there are "+str$(a)+" points"
   `text 0,16,"In row 2, there are "+str$(b)+" points"
   sync

loop

function get_length(width,length)
    a = width/(3.0*length)
   b# = (width/(3.0*length)-a)*3
   a=a*2
   if b#>=0 then inc a
   if b#>=1 then inc a
endfunction a

Although I don't see the tiny hexagons which I want to shrink down?

Back to top

Profile PM Email Website

Pincho Paxton

21

Years of Service

User Offline

Joined: 8th Dec 2002

Location:

Posted: 24th Oct 2009 17:45 Edited at: 24th Oct 2009 17:51

Link

Oh fixed it! And guess what.... I got a perfect snowflake in the middle of the screen! Yay.. this is great! Now it need to work all over the screen, and not just in the middle.

type hex
   a as integer
   b as integer
   c as integer
   x as float
   y as float
   xvel as float
   yvel as float
endtype
dim hexagon(1000) as hex
global startx as float = 1.0
global starty as float = 1.0
global sidelength as float = 30.0
curx as float = startx-sidelength
cury as float = starty-sidelength
switch as integer = 0
switch2 as boolean = 0
global exn as float = 700.0

do
cls
   exn=screen width()
   starty=Screen height()
   curx = startx-sidelength*2
   cury = 30
   switch as integer = 1
   switch2 as boolean = 1
   for a=1 to 1000

if switch=1
         switch=0
         inc curx,2*sidelength
      else
         switch=1
         inc curx, sidelength
      endif

if curx>exn
         if switch2=0
            switch2=1
            switch=0
            inc cury,sidelength*sqrt(3)/2
            curx=startx
         else
            switch=1
            switch2=0
            inc cury,sidelength*sqrt(3)/2
            curx=startx-sidelength/2
         endif
      endif
      hexagon(a).x=curx
      hexagon(a).y=cury
   next a
   for a=1 to 1000
      dx#=0
      dy#=0
      if spacekey()
         f#=1/sqrt((hexagon(a).x-mousex())^2+(hexagon(a).y-mousey())^2)`*rnd(1)
         dx#=(hexagon(a).x-mousex())*f#*-100
         dy#=(hexagon(a).y-mousey())*f#*-100
      endif
      dot hexagon(a).x+dx#,hexagon(a).y+dy#

next a
   `a=get_length(exn,sidelength)
   `b=get_length(exn-sidelength*1.5,sidelength)
   `text 0,0,"In row 1, there are "+str$(a)+" points"
   `text 0,16,"In row 2, there are "+str$(b)+" points"
   sync

loop

function get_length(width,length)
    a = width/(3.0*length)
   b# = (width/(3.0*length)-a)*3
   a=a*2
   if b#>=0 then inc a
   if b#>=1 then inc a
endfunction a

+ Code Snippet


type hex
   a as integer
   b as integer
   c as integer
   x as float
   y as float
   xvel as float
   yvel as float
endtype
dim hexagon(1000) as hex
global startx as float = 1.0
global starty as float = 1.0
global sidelength as float = 30.0
curx as float = startx-sidelength
cury as float = starty-sidelength
switch as integer = 0
switch2 as boolean = 0
global exn as float = 700.0

do
cls
   exn=screen width()
   starty=Screen height()
   curx = startx-sidelength*2
   cury = 30
   switch as integer = 1
   switch2 as boolean = 1
   for a=1 to 1000

      if switch=1
         switch=0
         inc curx,2*sidelength
      else
         switch=1
         inc curx, sidelength
      endif

      if curx>exn
         if switch2=0
            switch2=1
            switch=0
            inc cury,sidelength*sqrt(3)/2
            curx=startx
         else
            switch=1
            switch2=0
            inc cury,sidelength*sqrt(3)/2
            curx=startx-sidelength/2
         endif
      endif
      hexagon(a).x=curx
      hexagon(a).y=cury
   next a
   for a=1 to 1000
      dx#=0
      dy#=0
      if spacekey()
         f#=1/sqrt((hexagon(a).x-mousex())^2+(hexagon(a).y-mousey())^2)`*rnd(1)
         dx#=(hexagon(a).x-mousex())*f#*-100
         dy#=(hexagon(a).y-mousey())*f#*-100
      endif
      dot hexagon(a).x+dx#,hexagon(a).y+dy#

   next a
   `a=get_length(exn,sidelength)
   `b=get_length(exn-sidelength*1.5,sidelength)
   `text 0,0,"In row 1, there are "+str$(a)+" points"
   `text 0,16,"In row 2, there are "+str$(b)+" points"
   sync

loop

function get_length(width,length)
    a = width/(3.0*length)
   b# = (width/(3.0*length)-a)*3
   a=a*2
   if b#>=0 then inc a
   if b#>=1 then inc a
endfunction a

Edit: Oh I know what it is, I need the mouse to jump to grid coordinates.

Back to top

Profile PM Email Website

Pincho Paxton

21

Years of Service

User Offline

Joined: 8th Dec 2002

Location:

Posted: 24th Oct 2009 18:18 Edited at: 24th Oct 2009 18:29

Link

Hmmm.. A bit more complicated. The mouse cursor need to jump to the centre of a hexagon no matter what size it is, and the hexagons change positions. So that is the coding that I need.

I think this is the end of the Big Bang theory.

Back to top

Profile PM Email Website

Pincho Paxton

21

Years of Service

User Offline

Joined: 8th Dec 2002

Location:

Posted: 25th Oct 2009 06:26 Edited at: 25th Oct 2009 06:33

Link

Still struggling to lock the cursor. Anyway, I've added a temperature gauge, and slider. The colder the temperature the smaller the hexagons. Which reminds me of the Einstein-Bose condensate experiment. At zero the hexagons become just a dot (the atoms in Einsteins experiment also became just a dot), with total collapse of gravity. It seemed strange to Einstein, but works absolutely fine with this model.. in fact you have no other option. So far Quantum Physics are nothing more than everyday physics... which is exactly what you want.

Back to top

Profile PM Email Website

demons breath

20

Years of Service

User Offline

Joined: 4th Oct 2003

Location: Surrey, UK

Posted: 25th Oct 2009 14:00

Link

I think what you'd need is to find out which hexagon the mouse is over. Then take the x and y positions of 2 opposite dots and simply put it in between at (x1+x2)/2 (y1+y2)/2. No matter which points you used, as long as they were on opposite sides, the point half way in between would still be where you wanted the mouse cursor to jump to.

Unfortunately, if you had the mouse jumping to there then if you were reading the mousex() and mousey() positions then it would probably always jump the mouse to the center of the same hexagon unless you moved the mouse really fast. I don't know, I haven't tried it, but maybe storing where the mouse should be in an array (1,1) at (0,0) and (0,1) and then where the mouse would jump to in (1,0) and (1,1).

I'm not too sure how to tell where the mouse is though, and whether it's over a certain hexagon. Which is probably the part you weren't sure of. So this is probably a pretty useless post. Oh well...

"A West Texas girl, just like me"
-Bush

Back to top

Profile PM Website

Pincho Paxton

21

Years of Service

User Offline

Joined: 8th Dec 2002

Location:

Posted: 25th Oct 2009 17:47

Link

I'll just do it equal distances for now, it might still work anyway.

Back to top

Profile PM Email Website

Pincho Paxton

21

Years of Service

User Offline

Joined: 8th Dec 2002

Location:

Posted: 27th Oct 2009 19:21

Link

I have this at the moment, but it's not perfectly symmetrical. Is there anything I can do to get this more symmetrical?

set display mode 1152,864,16
sync on
type hex
   a as integer
   b as integer
   c as integer
   x as float
   y as float
   xvel as float
   yvel as float
endtype
hide mouse
dim hexagon(1000) as hex
Dim centre(700,2) as integer
global startx as float = 1.0
global starty as float = 1.0
global sidelength as float = 30.0

curx as float = startx-sidelength
cury as float = starty-sidelength
switch as integer = 0
switch2 as boolean = 0
global exn as float = 700.0
Load Image "Cursor.jpg",1
Load Image "Temperature.jpg",2
Lock = 1
mindist = 1000
set sprite 1,1,0
set sprite 2,1,0
do

sprite 2, screen width()-100,0,2
cls
   exn=screen width()
   starty=Screen height()
   curx = startx-sidelength*2
   cury = 30
   switch as integer = 1
   switch2 as boolean = 1
   Count = 1
   for a=1 to 1000

if switch=1
         switch=0
         inc curx,2*sidelength
Centre(Count,1) = curx - ((sidelength+(sprite width(1)/2)))
Centre(Count,2) = cury - (sprite Height(1)/2)
Inc count
      else
         switch=1
         inc curx, sidelength
      endif

if curx>exn
         if switch2=0
            switch2=1
            switch=0
            inc cury,sidelength*sqrt(3)/2
            curx=startx

else
            switch=1
            switch2=0
            inc cury,sidelength*sqrt(3)/2
            curx=startx-sidelength/2
         endif
      endif
      hexagon(a).x=curx
      hexagon(a).y=cury

If lock = 1
inc i
if i > 400 then i = 1

if abs(mousex()- centre(i,1)) < sidelength-6 and abs( mousey()- centre(i,2)) < sidelength-6

if mousex() + ((sidelength)-6) >= centre(i,1) then X3 = centre(i,1)
if mousey() + ((sidelength)-6) >= centre(i,2) then Y3 = centre(i,2)
if mousex() - ((sidelength)-6) <= centre(i,1) then X3 = centre(i,1)
if mousey() - ((sidelength)-6) <= centre(i,2) then Y3 = centre(i,2)

Endif
else
X3 = Mousex()
Y3 = Mousey()

Endif

dx#=0
      dy#=0
      if spacekey()
         f#=1/sqrt((hexagon(a).x-x3)^2+(hexagon(a).y-y3)^2)`*rnd(1)
         dx#=(hexagon(a).x-x3)*f#*-100
         dy#=(hexagon(a).y-y3)*f#*-100
      endif
      dot hexagon(a).x+dx#,hexagon(a).y+dy#

next a

sprite 1,x3,y3,1
   sync

loop

+ Code Snippet

set display mode 1152,864,16
sync on
type hex
   a as integer
   b as integer
   c as integer
   x as float
   y as float
   xvel as float
   yvel as float
endtype
hide mouse
dim hexagon(1000) as hex
Dim centre(700,2) as integer
global startx as float = 1.0
global starty as float = 1.0
global sidelength as float = 30.0

curx as float = startx-sidelength
cury as float = starty-sidelength
switch as integer = 0
switch2 as boolean = 0
global exn as float = 700.0
Load Image "Cursor.jpg",1
Load Image "Temperature.jpg",2
Lock = 1
mindist = 1000
set sprite 1,1,0
set sprite 2,1,0
do

sprite 2, screen width()-100,0,2
cls
   exn=screen width()
   starty=Screen height()
   curx = startx-sidelength*2
   cury = 30
   switch as integer = 1
   switch2 as boolean = 1
   Count = 1
   for a=1 to 1000

      if switch=1
         switch=0
         inc curx,2*sidelength
Centre(Count,1) = curx - ((sidelength+(sprite width(1)/2)))
Centre(Count,2) = cury - (sprite Height(1)/2)
Inc count
      else
         switch=1
         inc curx, sidelength
      endif

      if curx>exn
         if switch2=0
            switch2=1
            switch=0
            inc cury,sidelength*sqrt(3)/2
            curx=startx

         else
            switch=1
            switch2=0
            inc cury,sidelength*sqrt(3)/2
            curx=startx-sidelength/2
         endif
      endif
      hexagon(a).x=curx
      hexagon(a).y=cury






If lock = 1
inc i
if i > 400 then i = 1

if abs(mousex()- centre(i,1)) < sidelength-6 and abs( mousey()- centre(i,2)) < sidelength-6


if mousex() + ((sidelength)-6) >= centre(i,1) then X3 = centre(i,1)
if mousey() + ((sidelength)-6) >= centre(i,2) then Y3 = centre(i,2)
if mousex() - ((sidelength)-6) <= centre(i,1) then X3 = centre(i,1)
if mousey() - ((sidelength)-6) <= centre(i,2) then Y3 = centre(i,2)


Endif
else
X3 = Mousex()
Y3 = Mousey()

Endif

      dx#=0
      dy#=0
      if spacekey()
         f#=1/sqrt((hexagon(a).x-x3)^2+(hexagon(a).y-y3)^2)`*rnd(1)
         dx#=(hexagon(a).x-x3)*f#*-100
         dy#=(hexagon(a).y-y3)*f#*-100
      endif
      dot hexagon(a).x+dx#,hexagon(a).y+dy#

   next a


  sprite 1,x3,y3,1
   sync

loop

Back to top

Profile PM Email Website

Caleb1994

15

Years of Service

User Offline

Joined: 10th Oct 2008

Location: The Internet you idiot!

Posted: 28th Oct 2009 01:35

Link

Ummmm that's dbp code. this is the dbc board man haha

New Site! Check it out \/

Back to top

Profile PM Website

Neuro Fuzzy

16

Years of Service

User Offline

Joined: 11th Jun 2007

Location:

Posted: 28th Oct 2009 02:22

Link

still, it's about the solution to the problem. The algorithm to solve for it still holds true in any language.

Back to top

Profile PM Email

Sorry your browser is not supported!

DarkBASIC Discussion / Vector Gravity Simulation in realtime (physics project)