mandelbrot.hplsql

declare size int = 100;
declare fac double = 2.0 / 100.0;

create procedure computeRow(y int)
begin
  declare Ci double = (y*fac - 1.0);

  for x in 0..size-1 loop
      declare Zr double = 0.0;
      declare Zi double = 0.0;
      declare Cr double = (x*fac - 1.5);
      declare i int = 50;
      declare ZrN double = 0;
      declare ZiN double = 0;
      loop
         Zi = 2.0 * Zr * Zi + Ci;
         Zr = ZrN - ZiN + Cr;
         ZiN = Zi * Zi;
         ZrN = Zr * Zr;
         i = i -1;
         exit when (ZiN + ZrN > 4.0 or i <= 0);
      end loop;
  end loop;
end;

declare start int = current_time_millis();

for y in 0..size-1 loop
  computeRow(y);
end loop;

print 'total: ' || (current_time_millis() - start);