here the images:

I’ve been working yesterday in a Circle Packing, and I started with matsys pseudo-code.  I found a lot bugs, finally I decided to make a mixture between his algorithm and mine.  Circle Packing can be a very usefull algorithm and I think until now we have not seen their real power applied to architecture.

This is my little contribuition (pseudo_Code)

1_ Input: maxradius, boundary, iterations (not necessary the number of circles)

2_Add in a random position within the boundary a circle with maximum radius

3_Add a Random Point

4_If this Point is inside of the Circle just created, add another one and delete the old one

5_Add Line to this point (out of the circle) to the center of the circle, (here I use the same idea of matsys) “If the distance is greated than the maximum radius add a circle at that point with the maximum radius. If the distance is less than the maximum radius, add a circle at that point with the measured distance.”

6_ add another point and test if this point is within of any circle and add a line to every center circle to find the closest circle (step _5)

7_ repeat the steps 5  to 6 the number of the iterations

soon more images.

share, modify and send me an Image

best

Option Explicit
‘Script written by <Carlos de la Barrera>
‘Script copyrighted by <Carlos de la Barrera>

Call Main()
Sub Main()
Dim Pt, R
Pt = Array(rnd * 100, rnd * 100, 0)
R = FuncCir(Pt, 10)

Call Pack(R)
End Sub

Function Pack(R)
Dim Pt, i, count

count = Rhino.GetInteger(” number of circles “, 50, 2, 1000)
If IsNull (count) Then Exit Function

Dim ArrData()

For i = 0 To count
ReDim Preserve ArrData(i)
ArrData(i) = R

Pt = Array(Rnd * 100, rnd * 100, 0)

Dim j
Dim t: t = 0

For j = 0 To UBound(ArrData)
If Rhino.IsPointOnSurface(ArrData(j)(1)(0), Pt)Then
Call Rhino.print (”point in surface”)
t = t + 1

End If
Next

If (t = 0) Then
ReDim ArrLine(UBound(ArrData))
Dim k, line

For K = 0 To UBound(ArrData)
line = Rhino.AddLine(Pt, Rhino.CircleCenterPoint(ArrData(k)(0)))
ArrLine(k) = Line
Next

Dim d, d2, l
d = funcDist(Pt, ArrData(0)(0), ArrLine(0))

For l = 0 To UBound(ArrLine)
d2 = funcDist(Pt, ArrData(l)(0), ArrLine(l))

If d < d2 Then
‘Call Rhino.Print (”nothing”)

Else
d = d2

End If
Next

R = FuncCir(Pt, d)
Call Rhino.DeleteObjects (ArrLine)

End If
Next
End Function

Function FuncCir(p, D)

If D >= 10 Then
D = 10

End If

Dim plane
Dim R(1)

plane = Rhino.PlaneFromFrame(p, Array(1.0,0.0,0.0), Array(0.0,1.0,0.0))
R(0)   = Rhino.AddCircle(plane, D)
R(1) = Rhino.AddPlanarSrf(Array(R(0)))

FuncCir = R
End Function

Function funcDist(Pt, circleP, Line)
Dim Ap, j, d
Ap = Rhino.CurveCurveIntersection (circleP, Line)

For j = 0 To UBound(Ap)
d = Rhino.Distance(Pt, Ap(j, 1))
Next

funcDist = d
End Function

the other day I listen: - talk it’s easy,  ¡ show me the code ! - ]]> http://www.designemergente.org/rhinoscripting/circle-packing/feed/ http://www.designemergente.org/rhinoscripting/circle-packing/ http://feeds.feedburner.com/~r/designemergente/~3/357358363/ http://www.designemergente.org/rhinoscripting/sensitive-analisys-prato-italy/#comments Wed, 06 Aug 2008 09:12:50 +0000 Grupo Diseño Emergente http://www.designemergente.org/?p=423 I was invited to paticipate in a contest to the Comune di Prato in Italy. I have to say thank you to this cool group of architects from Lisboa: Duarte Cardoso Ferreira, Nuno Lobo Guerra, Ricardo Boaventura, Rita Frias Silva from ADOC architects and all the team of 4ARQ encharged of the 3D images.

I was encharge to develope by programming the plans and the general features of the school. We start working identify every tree in the ground and the area of the differents rooms for the school. Every tree was identify by a point and a minum area of influence (in the fisrt step of the algorithm) and each area of the room of the school was enter in to the algorithm to know the best performance of the build with out touch any tree.

Then the algorithm calculate a voronoi diagram (only with the trees) and inform the resultant areas and we could see where can “grow” our building. the next step was develope a second algorithm to generate pseudo-random points, between the trees and responding obviously to the influence area of them.

Then the process consist basicaly in create generations and test the results, if the dispositions of points allow build the school (first and more important test), if the result produce nice spaces, if we like it, if the algorithm respect entrance determined by us, etc. etc. after every generation, we pass again this second algorithm to produce in some points a random movement to explore new configurations and results.

I take this idea from Mutsuro Sasaki, he use an strategy based in the approach in successive steps to get the best result moving the points in pseudo random way.

Thi first step was to find the correct data store to find the points. then I put all this points inside of an array and then create the curve. and for the other curves I did the same, searching the even and odd of every point.

I know there are a lot kind of patterns in the textile insdustri, and also my Grandmother know how to create with the crochet technics differents patterns, some of them with very complex estructures. If we thought that an algortihm is and a mathematic or aritmetic rule that is repeated, (Cecil Balmond, Informal, page 112), we will see that my grandmother and all the people who elaborate knitted fabric are actually programming an Algortihm.

These images represent the process of my recursion work that started few weeks ago, with the post recursion, I will continue working in this kind of programming, and exploring differents types of recursion. there are a lot information of recursion in the web. and also example codes in differents languages like C++, C#, Java, Rvb, etc. etc.

In relation with the architecture you can see the web page of Marc Fornes, he explore the process of recursion in deep, you can find very cool information about this kind of programming, and other interesting data type. THE VERY MANY also is reference of my work and reference of many other architects who explore the new technologies.

Simple rules configurate comples and beatifull Images. When you try to create complex results under complicate formules you have to be sure that your results will crash, and also the results will be not good enough. the result of the Complex Beauty appears after many iterations, and depend extreamly of your code and how many iterations you want.

Le modèle ci-dessous représente une série d’expériences avec différentes configurations de la fonction.

I learned recursion using C++, but I never used In Rvb, so this is my first experiments in recursion in Rvb. and this time I add The source of the original code, add comments, modify it, share, make differents Calls and send me Images.

Best.

Option Explicit
‘Script written by Carlos de la Barrera’

Call Main()
Sub Main()
Dim pt1, pt2, pt3, trip

pt1 = Rhino.GetPoint(”sel fisrt point”)
If IsNull (pt1) Then Exit Sub

pt2 = Rhino.GetPoint(”sel second point”)
If IsNull (pt2) Then Exit Sub

pt3 = Rhino.GetPoint(”sel third point”)
If IsNull (pt3) Then Exit Sub

trip = Rhino.AddPolyline(Array(pt1, pt2, pt3, pt1))

Call myFractal(trip, pt1, pt2, pt3)
End Sub

Function myFractal(trip, Opt1, Opt2, Opt3)
Dim pt1, pt2, pt3

pt1 = myMidFunction(Opt1, Opt2)
pt2 = myMidFunction(Opt2, Opt3)
pt3 = myMidFunction(Opt3, Opt1)

trip = Rhino.AddPolyline(Array(pt1, pt2, pt3, pt1))

If Rhino.Distance(pt1, pt2) > 0.01 Then
Call myFractal(trip, pt1, pt2, pt3)

Else
Exit Function

End If
End Function

Function myMidFunction(P1, P2)

Dim MidP
MidP = Array((P1(0) + P2(0)) / 2, (P1(1) + P2(1)) / 2, (P1(2) + P2(2)) / 2)

myMidFunction = MidP

End Function