mirror of
https://github.com/carp-lang/Carp.git
synced 2024-10-11 12:37:32 +03:00
68 lines
1.8 KiB
Plaintext
68 lines
1.8 KiB
Plaintext
;;; Taken from
|
|
;;; http://benchmarksgame.alioth.debian.org/u64q/program.php?test=mandelbrot&lang=gcc&id=2
|
|
|
|
(add-cflag "-Wall -pipe -O3 -fomit-frame-pointer -march=native")
|
|
|
|
(use Double)
|
|
(use Int)
|
|
(use String)
|
|
(use IO)
|
|
|
|
(def w 16000.0)
|
|
(def h 16000.0)
|
|
|
|
(system-include "stdio.h")
|
|
(register-type FILE)
|
|
(register stdout FILE)
|
|
(register putc (Fn [Int FILE] ()))
|
|
|
|
(defn main []
|
|
(do
|
|
(print &(str* @"P4\n" (str w) @" " (str h) @"\n"))
|
|
(let [iter 50
|
|
limit 2.0
|
|
byte_acc 0
|
|
bit_num 0]
|
|
(for [y 0 (to-int h)]
|
|
(for [x 0 (to-int w)]
|
|
(let [Zr 0.0
|
|
Zi 0.0
|
|
Tr 0.0
|
|
Ti 0.0
|
|
Cr (- (/ (* 2.0 (from-int x)) w) 1.5)
|
|
Ci (- (/ (* 2.0 (from-int y)) w) 1.0)]
|
|
(do
|
|
(let [i 0]
|
|
(while (and (< i iter)
|
|
(<= (+ Tr Ti) (* limit limit)))
|
|
(do
|
|
(set! Zi (+ (* 2.0 (* Zr Zi)) Ci))
|
|
(set! Zr (+ (- Tr Ti) Cr))
|
|
(set! Tr (* Zr Zr))
|
|
(set! Ti (* Zi Zi))
|
|
(set! i (inc i)))))
|
|
|
|
(set! byte_acc (* 2 byte_acc))
|
|
|
|
(when (<= (+ Tr Ti) (* limit limit))
|
|
(set! byte_acc (bit-or byte_acc 1)))
|
|
|
|
(set! bit_num (inc bit_num))
|
|
|
|
(if (= bit_num 8)
|
|
(do
|
|
(putc byte_acc stdout)
|
|
(set! byte_acc 0)
|
|
(set! bit_num 0))
|
|
(if (= x (- (to-int w) 1))
|
|
(do
|
|
(set! byte_acc
|
|
(bit-shift-left byte_acc (- 8 (mod (to-int w) 8))))
|
|
(putc byte_acc stdout)
|
|
(set! byte_acc 0)
|
|
(set! bit_num 0))
|
|
())))))))))
|
|
|
|
|
|
;;; (build)
|