Here are some more shots with the random walking origin points follow paramedic, ordered structures. For +Roger Bagula :
https://plus.google.com/101799841244447089430/posts/PjpUnX438pa
_
#Fractal #Math #DLA #Field #Vector #Random
https://plus.google.com/101799841244447089430/posts/PjpUnX438pa
_
#Fractal #Math #DLA #Field #Vector #Random
‹




›
2016-04-04
4 Photos - View album
Not exactly what I had in mind ,
but very good work!Apr 5, 2016
Chris “M” ThomassonOwnerApr 7, 2016
I explained the idea before, took me an hour to get the old drive going and find the files:
Programmed 15 years ago in true basic: starts with a Sierpinski set
iFS and then does DLA on that base:I didn't comment it because I don't think I published it: didn't work very well: used too much memory :
(True Basic)
SET MODE "color"
SET WINDOW 0,1026,0,750
SET COLOR MIX (1) 0,0, 0
SET BACKGROUND COLOR "white"
DIM a(0 to 1026,0 to 750)
LET k=255
FOR i=1 to 1025
FOR j=1 to 749
LET a(i,j)=k
SET COLOR a(i,j)
PLOT i,j
NEXT j
NEXT i
FOR i=0 to 1026
LET a(i,0)=k-1
SET COLOR a(i,0)
PLOT i,0
LET a(i,750)=k-1
SET COLOR a(i,750)
PLOT i,750
LET a(i,749)=k-1
SET COLOR a(i,749)
PLOT i,749
NEXT i
FOR j=0 to 750
LET a(0,j)=k-1
SET COLOR a(0,j)
PLOT 0,j
LET a(1026,j)=k-1
SET COLOR a(1026,j)
PLOT 1026,j
LET a(1025,j)=k-1
SET COLOR a(1025,j)
PLOT 1025,j
NEXT j
LET x=1
LET y=1
LET c=0
LET s1=250
LET s2 =(s1)*750/1026
DIM a1(25),b1(25)
let m=3
FOR i=1 to m
LET a1(i)=cos(2*pi*i/m)
LET b1(i)=sin(2*pi*i/m)
NEXT i
RANDOMIZE
(Rem IFS ; If m=4 Sierpinski square space fill start)
FOR n= 1 TO 20000
LET c =RND
LET l=1+int(c*m)
LET x1=x/2+a1(l)
LET y1=y/2+b1(l)
LET x=x1
LET y=y1
SET COLOR 255
IF n>10 THEN
LET a(int( 1026/2+s1*x),int(750/2+s2*y))=k-1
SET COLOR 2
PLOT Int(1026/2+s1*x),Int(750/2+s2*y)
END IF
NEXT n
FOR n= 1 TO 100000000
LET Iin=int(rnd*1026)
LET Jin=int(rnd*750)
IF a(Iin, Jin)< k then
LET a(Iin, Jin)=a(Iin, Jin)-1
SET COLOR a(Iin, Jin)
PLOT Iin,Jin
FOR I=-1 +Iin to 1+Iin
FOR J=-1 +Jin to 1+Jin
IF a(abs(I), abs(J))<=k AND (I<>Iin OR J<>Jin) and (I>=0 and J>=0) and (I<=1026 and J<=750) then
IF a(I,J)>0 then LET a(I,J)=a(I,J)-1
SET COLOR a(I,J)
PLOT I,J
END IF
NEXT J
NEXT I
END IF
NEXT n
END
(This one plotted a circle first: then did DLA:)
(* True Basic*)
SET MODE "color"
SET WINDOW 0,1026,0,750
SET COLOR MIX (1) 0,0, 0
SET BACKGROUND COLOR "white"
DIM a(0 to 1026,0 to 750)
LET k=255
FOR i=1 to 1025
FOR j=1 to 749
LET a(i,j)=k
SET COLOR a(i,j)
PLOT i,j
NEXT j
NEXT i
FOR i=0 to 1026
LET a(i,0)=k-1
SET COLOR a(i,0)
PLOT i,0
LET a(i,750)=k-1
SET COLOR a(i,750)
PLOT i,750
LET a(i,749)=k-1
SET COLOR a(i,749)
PLOT i,749
NEXT i
FOR j=0 to 750
LET a(0,j)=k-1
SET COLOR a(0,j)
PLOT 0,j
LET a(1026,j)=k-1
SET COLOR a(1026,j)
PLOT 1026,j
LET a(1025,j)=k-1
SET COLOR a(1025,j)
PLOT 1025,j
NEXT j
(* rem circle plot: 600 points dense: in this resolution nearly solid*)
FOR l=0 to 2*Pi step 2*Pi/600
LET I0=int(1026/2+100*cos(l))
LET J0=int(750/2+100*sin(l))
LET a(I0,J0)=k-1
NEXT l
FOR n= 1 TO 100000000
LET Iin=int(rnd*1026)
LET Jin=int(rnd*750)
IF a(Iin, Jin)< k then
LET a(Iin, Jin)=a(Iin, Jin)-1
SET COLOR a(Iin, Jin)
PLOT Iin,Jin
FOR I=-1 +Iin to 1+Iin
FOR J=-1 +Jin to 1+Jin
IF a(abs(I), abs(J))<=k AND (I<>Iin OR J<>Jin) and (I>=0 and J>=0) and (I<=1026 and J<=750) then
IF a(I,J)>0 then LET a(I,J)=a(I,J)-1
SET COLOR a(I,J)
PLOT I,J
END IF
NEXT J
NEXT I
END IF
NEXT n
END
(Here is my Levy flight or "Jump" algorithm, I think;)
(* True Basic*)
SET MODE "color"
SET WINDOW 0,1026,0,750
SET COLOR MIX (1) 0,0, 0
SET BACKGROUND COLOR "white"
DIM a(0 to 1026,0 to 750)
LET k=255
FOR i=0 to 1026
FOR j=0 to 750
LET a(i,j)=k
SET COLOR a(i,j)
PLOT i,j
NEXT j
NEXT i
FOR i=0 to 1026
LET a(i,0)=k-1
SET COLOR a(i,0)
PLOT i,0
LET a(i,750)=k-1
SET COLOR a(i,750)
PLOT i,750
LET a(i,749)=k-1
SET COLOR a(i,749)
PLOT i,749
NEXT i
FOR j=0 to 750
LET a(0,j)=k-1
SET COLOR a(0,j)
PLOT 0,j
LET a(1026,j)=k-1
SET COLOR a(1026,j)
PLOT 1026,j
LET a(1025,j)=k-1
SET COLOR a(1025,j)
PLOT 1025,j
NEXT j
FOR n= 1 TO 100000000
LET Iin=int(rnd*1026)
LET Jin=int(rnd*750)
IF a(Iin, Jin)< k then
LET a(Iin, Jin)=a(Iin, Jin)-1
SET COLOR a(Iin, Jin)
PLOT Iin,Jin
LET I=INT(-2+MOD(5*RND,5))+Iin
LET J=INT(-2+MOD(5*RND,5))+Jin
IF I<0 THEN LET I=0
IF I>1026 THEN LET I=1026
IF J<0 THEN LET J=0
IF J>750 THEN LET J=750
IF a(abs(I), abs(J))<=k and I<>Iin and J<>Jin and I>=0 and J>=0 and I<=1026 and J<=750 then
IF a(I,J)>0 then LET a(I,J)=a(I,J)-1
SET COLOR a(I,J)
PLOT I,J
END IF
END IF
NEXT n
ENDApr 7, 2016