Smallest enclosing circle problem: Difference between revisions

 

 

=={{header|Go}}==

{{trans|Wren}}

 

import (

fmt.Println(welzl(test))

}

 

{{out}}

{(0.500000, 0.500000), 0.707107}

{(1.000000, 1.000000), 5.000000}

</pre>

 

=={{header|Julia}}==

N-dimensional solution using the Welzl algorithm. Derived from the BoundingSphere.jl module at https://github.com/JuliaFEM. See also Bernd Gärtner's paper at people.inf.ethz.ch/gaertner/subdir/texts/own_work/esa99_final.pdf.

using LinearAlgebra, Random

 

 

testwelzl()

<pre>

For points: [[0.0, 0.0], [0.0, 1.0], [1.0, 0.0]]

</pre>

 

 

 

 

 

 

 

 

 

 

 

 

 

 

{{out}}

<pre>

=={{header|Python}}==

{{trans|Julia}}

 

class ProjectorStack:

print(" Center is at: ", nsphere.center)

print(" Radius is: ", np.sqrt(nsphere.sqradius), "\n")

<pre>

For points: [[0. 0.]

=={{header|Raku}}==

{{trans|Go}}

method gist { "({$.x~", "~$.y})" }

}

 

class C { has (P $.c, Numeric $.r); # Circle

method gist { "Center = " ~$.c.gist ~ " and Radius = $.r\n" }

 

 

}

 

# returns the square of the distance between two points

 

sub getCenter (\bx, \by, \cx, \cy --> P) {

P.new: x => (cy*b - by*c) / (2 * d), y => (bx*c - cx*b) / (2 * d)

} # returns the center of a circle defined by 3 points

 

sub circleFrom3 (P \a, P \b, P \c --> C) {

} # returns a unique circle that intersects 3 points

 

sub circleFrom2 (P \a, P \b --> C ) {

my \center = P.new: x => ((a.x + b.x) / 2), y => ((a.y + b.y) / 2) ;

} # returns smallest circle that intersects 2 points

 

sub secTrivial( @rs --> C ) {

when * == 1 { return C.new: c => @rs[0], r => 0 }

}

for 0, 1 X 1, 2 -> ( \i, \j ) {

}

} # returns smallest enclosing circle for n ≤ 3

 

my \p = @ps.shift;

} # helper function for Welzl method

 

# applies the Welzl algorithm to find the SEC

 

my @tests = (

}

 

{{out}}

<pre>Solution for smallest circle enclosing ((0, 0) (0, 1) (1, 0)) :

 

It is based on Welzl's algorithm and follows closely the C++ code [https://www.geeksforgeeks.org/minimum-enclosing-circle-set-2-welzls-algorithm/?ref=rp here].

 

var Rand = Random.new()

]

 

 

{{out}}