Skip to content

Instantly share code, notes, and snippets.

@Hailong

Hailong/Mabdelbrot.java

Forked from dstogov/Mabdelbrot.java
Created November 21, 2020 05:47
Show Gist options

  • Select an option

    Learn more about clone URLs
  • Save Hailong/9e13e6353abfb8fdde2f7ebbaf12d749 to your computer and use it in GitHub Desktop.

Select an option

Learn more about clone URLs
Save Hailong/9e13e6353abfb8fdde2f7ebbaf12d749 to your computer and use it in GitHub Desktop.
Download ZIP

Revisions

  1. @dstogov dstogov created this gist Feb 27, 2015.
    49 changes: 49 additions & 0 deletions Mabdelbrot.java
    View file
    Original file line number Diff line number Diff line change
    @@ -0,0 +1,49 @@
    // by Erik Wrenholt
    import java.util.*;

    class Mandelbrot
    {
    static int BAILOUT = 16;
    static int MAX_ITERATIONS = 1000;

    private static int iterate(float x, float y)
    {
    float cr = y-0.5f;
    float ci = x;
    float zi = 0.0f;
    float zr = 0.0f;
    int i = 0;
    while (true) {
    i++;
    float temp = zr * zi;
    float zr2 = zr * zr;
    float zi2 = zi * zi;
    zr = zr2 - zi2 + cr;
    zi = temp + temp + ci;
    if (zi2 + zr2 > BAILOUT)
    return i;
    if (i > MAX_ITERATIONS)
    return 0;
    }
    }

    public static void main(String args[])
    {
    Date d1 = new Date();
    int x,y;
    for (y = -39; y < 39; y++) {
    System.out.print("\n");
    for (x = -39; x < 39; x++) {
    if (iterate(x/40.0f,y/40.0f) == 0)
    System.out.print("*");
    else
    System.out.print(" ");

    }
    }
    Date d2 = new Date();
    long diff = d2.getTime() - d1.getTime();
    System.out.println("\nJava Elapsed " + diff/1000.0f);

    }
    }
    53 changes: 53 additions & 0 deletions b.c
    View file
    Original file line number Diff line number Diff line change
    @@ -0,0 +1,53 @@
    #include <stdio.h>
    #include <sys/time.h>

    #define BAILOUT 16
    #define MAX_ITERATIONS 1000

    int mandelbrot(double x, double y)
    {
    double cr = y - 0.5;
    double ci = x;
    double zi = 0.0;
    double zr = 0.0;
    int i = 0;

    while(1) {
    i ++;
    double temp = zr * zi;
    double zr2 = zr * zr;
    double zi2 = zi * zi;
    zr = zr2 - zi2 + cr;
    zi = temp + temp + ci;
    if (zi2 + zr2 > BAILOUT)
    return i;
    if (i > MAX_ITERATIONS)
    return 0;
    }

    }

    int main (int argc, const char * argv[]) {
    struct timeval aTv;
    gettimeofday(&aTv, NULL);
    long init_time = aTv.tv_sec;
    long init_usec = aTv.tv_usec;

    int x,y;
    for (y = -39; y < 39; y++) {
    printf("\n");
    for (x = -39; x < 39; x++) {
    int i = mandelbrot(x/40.0, y/40.0);
    if (i==0)
    printf("*");
    else
    printf(" ");
    }
    }
    printf ("\n");

    gettimeofday(&aTv,NULL);
    double query_time = (aTv.tv_sec - init_time) + (double)(aTv.tv_usec - init_usec)/1000000.0;
    printf ("C Elapsed %0.3f\n", query_time);
    return 0;
    }
    47 changes: 47 additions & 0 deletions b.lua
    View file
    Original file line number Diff line number Diff line change
    @@ -0,0 +1,47 @@
    -- By Erik Wrenholt

    local BAILOUT = 16
    local MAX_ITERATIONS = 1000

    function iterate(x,y)

    local cr = y-0.5
    local ci = x
    local zi = 0.0
    local zr = 0.0
    local i = 0

    while 1 do
    i = i+1
    local temp = zr * zi
    local zr2 = zr*zr
    local zi2 = zi*zi
    zr = zr2-zi2+cr
    zi = temp+temp+ci
    if (zi2+zr2 > BAILOUT) then
    return i
    end

    if (i > MAX_ITERATIONS) then
    return 0
    end
    end

    end

    function mandelbrot()
    local t = os.time()
    for y = -39, 38 do
    for x = -39, 38 do
    if (iterate(x/40.0, y/40) == 0) then
    io.write("*")
    else
    io.write(" ")
    end
    end
    io.write("\n")
    end
    io.write(string.format("Time Elapsed %d\n", os.time() - t))
    end

    mandelbrot()
    54 changes: 54 additions & 0 deletions b.php
    View file
    Original file line number Diff line number Diff line change
    @@ -0,0 +1,54 @@
    <?php
    define("BAILOUT",16);
    define("MAX_ITERATIONS",1000);

    class Mandelbrot
    {
    function Mandelbrot()
    {
    $d1 = microtime(1);
    for ($y = -39; $y < 39; $y++) {
    for ($x = -39; $x < 39; $x++) {

    if ($this->iterate($x/40.0,$y/40.0) == 0)
    echo("*");
    else
    echo(" ");

    }
    echo("\n");
    }
    $d2 = microtime(1);
    $diff = $d2 - $d1;
    printf("\nPHP Elapsed %0.3f\n", $diff);
    }

    function iterate($x,$y)
    {
    $cr = $y-0.5;
    $ci = $x;
    $zr = 0.0;
    $zi = 0.0;
    $i = 0;
    while (true) {
    $i++;
    $temp = $zr * $zi;
    $zr2 = $zr * $zr;
    $zi2 = $zi * $zi;
    $zr = $zr2 - $zi2 + $cr;
    $zi = $temp + $temp + $ci;
    if ($zi2 + $zr2 > BAILOUT)
    return $i;
    if ($i > MAX_ITERATIONS)
    return 0;
    }

    }


    }

    ob_start();
    $m = new Mandelbrot();
    ob_end_flush();
    ?>
    55 changes: 55 additions & 0 deletions b.pl
    View file
    Original file line number Diff line number Diff line change
    @@ -0,0 +1,55 @@
    $BAILOUT=16;
    $MAX_ITERATIONS=1000;

    $begin = time();

    sub mandelbrot {
    local $x = $_[0];
    local $y = $_[1];

    local $cr = $y - 0.5;
    local $ci = $x;
    local $zi = 0.0;
    local $zr = 0.0;
    local $i = 0;

    while (1)
    {
    $i = $i + 1;
    local $temp = $zr * $zi;
    local $zr2 = $zr * $zr;
    local $zi2 = $zi * $zi;
    $zr = $zr2 - $zi2 + $cr;
    $zi = $temp + $temp + $ci;
    if ($zi2 + $zr2 > $BAILOUT)
    {
    return $i;
    }
    if ($i > $MAX_ITERATIONS)
    {
    return 0;
    }
    }
    }

    for ($y = -39; $y < 39; $y++)
    {
    print("\n");
    for ($x = -39; $x < 39; $x++)
    {
    $i = mandelbrot($x/40.0, $y/40.0);
    if ($i == 0)
    {
    print("*");
    }
    else
    {
    print(" ");
    }
    }
    }
    print("\n");

    $end = time() - $begin;

    print "Perl Elapsed $end\n";
    42 changes: 42 additions & 0 deletions b.py
    View file
    Original file line number Diff line number Diff line change
    @@ -0,0 +1,42 @@
    import sys, time
    stdout = sys.stdout

    BAILOUT = 16
    MAX_ITERATIONS = 1000

    class Iterator:
    def __init__(self):
    print 'Rendering...'
    for y in range(-39, 39):
    stdout.write('\n')
    for x in range(-39, 39):
    i = self.mandelbrot(x/40.0, y/40.0)

    if i == 0:
    stdout.write('*')
    else:
    stdout.write(' ')

    def mandelbrot(self, x, y):
    cr = y - 0.5
    ci = x
    zi = 0.0
    zr = 0.0
    i = 0

    while True:
    i += 1
    temp = zr * zi
    zr2 = zr * zr
    zi2 = zi * zi
    zr = zr2 - zi2 + cr
    zi = temp + temp + ci

    if zi2 + zr2 > BAILOUT:
    return i
    if i > MAX_ITERATIONS:
    return 0

    t = time.time()
    Iterator()
    print '\nPython Elapsed %.02f' % (time.time() - t)
    46 changes: 46 additions & 0 deletions b.ru
    View file
    Original file line number Diff line number Diff line change
    @@ -0,0 +1,46 @@
    BAILOUT = 16
    MAX_ITERATIONS = 1000

    class Mandelbrot

    def initialize
    puts "Rendering"
    for y in -39...39 do
    puts
    for x in -39...39 do
    i = iterate(x/40.0,y/40.0)
    if (i == 0)
    print "*"
    else
    print " "
    end
    end
    end
    end

    def iterate(x,y)
    cr = y-0.5
    ci = x
    zi = 0.0
    zr = 0.0
    i = 0

    while(1)
    i += 1
    temp = zr * zi
    zr2 = zr * zr
    zi2 = zi * zi
    zr = zr2 - zi2 + cr
    zi = temp + temp + ci
    return i if (zi2 + zr2 > BAILOUT)
    return 0 if (i > MAX_ITERATIONS)
    end

    end

    end

    time = Time.now
    Mandelbrot.new
    puts
    puts "Ruby Elapsed %f" % (Time.now - time)
    45 changes: 45 additions & 0 deletions b.sq
    View file
    Original file line number Diff line number Diff line change
    @@ -0,0 +1,45 @@
    function mandelbrot(x, y) {
    local cr = y - 0.5;
    local ci = x;
    local zi = 0.0;
    local zr = 0.0;
    local i = 0;
    local BAILOUT = 16;
    local MAX_ITERATIONS = 1000;

    while(1) {
    i++;
    local temp = zr * zi;
    local zr2 = zr * zr;
    local zi2 = zi * zi;
    zr = zr2 - zi2 + cr;
    zi = temp + temp + ci;
    if (zi2 + zr2 > BAILOUT) {
    return i;
    }
    if (i > MAX_ITERATIONS) {
    return 0;
    }
    }
    }

    function mandelbrot_run() {
    local x;
    local y;

    // local date = Date();
    for (y = -39; y < 39; y++) {
    for (x = -39; x < 39; x++) {
    if (mandelbrot(x/40.0, y/40.0) == 0) {
    print("*");
    } else {
    print(" ");
    }
    }
    print("\n");
    }
    // local date2 = Date();
    // output += "\nJavaScript Elapsed " + (date2.getTime() - date.getTime()) / 1000;
    }

    mandelbrot_run();
    25 changes: 25 additions & 0 deletions b.txt
    View file
    Original file line number Diff line number Diff line change
    @@ -0,0 +1,25 @@
    LapTop (12 Feb 2015)
    ====================
    PHP7-JIT (JIT=on) 0.011
    gcc -O2 (4.9.2) 0.013
    LuaJIT-2.0.3 (JIT=on) 0.014
    gcc -O0 (4.9.2) 0.022
    HHVM-3.5.0 (JIT=on) 0.030
    Java-1.8.0 (JIT=on) 0.059
    LuaJIT-2.0.3 (JIT=off) 0.073
    Java-1.8.0 (JIT=off) 0.251
    PHP-7 0.281
    squirrel-3.0.4 0.335
    Lua-5.2.2 0.339
    PHP-5.6 0.379
    PHP-5.5 0.383
    PHP-5.4 0.406
    ruby-2.1.5 0.684
    PHP-5.3 0.855
    HHVM-3.5.0 (JIT=off) 0.978
    PHP-5.2 1.096
    python-2.7.8 1.128
    PHP-5.1 1.217
    perl-5.18.4 2.083
    PHP-4.4 4.209
    PHP-5.0 4.434