#lang typed/racket/base

;; Line renderers.

(require typed/racket/class racket/match racket/math racket/list racket/sequence
         plot/utils
         (only-in math/statistics stddev)
         "../common/type-doc.rkt"
         "../common/utils.rkt")

(provide (all-defined-out))

;; ===================================================================================================
;; Lines, parametric, polar

(: lines-render-proc (-> (Listof (Vectorof Real))
                         Plot-Color Nonnegative-Real Plot-Pen-Style
                         Nonnegative-Real
                         (U String #f)
                         2D-Render-Proc))
(define ((lines-render-proc vs color width style alpha label) area)
  (send area put-alpha alpha)
  (send area put-pen color width style)
  (send area put-lines vs)
  
  (cond [label  (line-legend-entry label color width style)]
        [else   empty]))

(:: lines
    (->* [(Sequenceof (Sequenceof Real))]
         [#:x-min (U Real #f) #:x-max (U Real #f)
          #:y-min (U Real #f) #:y-max (U Real #f)
          #:color Plot-Color
          #:width Nonnegative-Real
          #:style Plot-Pen-Style
          #:alpha Nonnegative-Real
          #:label (U String #f)]
         renderer2d))
(define (lines vs
               #:x-min [x-min #f] #:x-max [x-max #f]
               #:y-min [y-min #f] #:y-max [y-max #f]
               #:color [color (line-color)]
               #:width [width (line-width)]
               #:style [style (line-style)]
               #:alpha [alpha (line-alpha)]
               #:label [label #f])
  (define fail/kw (make-raise-keyword-error 'lines))
  (cond
    [(and x-min (not (rational? x-min)))  (fail/kw "#f or rational" '#:x-min x-min)]
    [(and x-max (not (rational? x-max)))  (fail/kw "#f or rational" '#:x-max x-max)]
    [(and y-min (not (rational? y-min)))  (fail/kw "#f or rational" '#:y-min y-min)]
    [(and y-max (not (rational? y-max)))  (fail/kw "#f or rational" '#:y-max y-max)]
    [(not (rational? width))  (fail/kw "rational?" '#:width width)]
    [(or (> alpha 1) (not (rational? alpha)))  (fail/kw "real in [0,1]" '#:alpha alpha)]
    [else
     (let ([vs  (sequence->listof-vector 'lines vs 2)])
       (define rvs (filter vrational? vs))
       (cond [(empty? rvs)  (renderer2d #f #f #f #f)]
             [else
              (match-define (list (vector #{rxs : (Listof Real)} #{rys : (Listof Real)}) ...) rvs)
              (let ([x-min  (if x-min x-min (apply min* rxs))]
                    [x-max  (if x-max x-max (apply max* rxs))]
                    [y-min  (if y-min y-min (apply min* rys))]
                    [y-max  (if y-max y-max (apply max* rys))])
                (renderer2d (vector (ivl x-min x-max) (ivl y-min y-max)) #f default-ticks-fun
                            (lines-render-proc vs color width style alpha label)))]))]))

(:: parametric
    (->* [(-> Real (Sequenceof Real)) Real Real]
         [#:x-min (U Real #f) #:x-max (U Real #f)
          #:y-min (U Real #f) #:y-max (U Real #f)
          #:samples Positive-Integer
          #:color Plot-Color
          #:width Nonnegative-Real
          #:style Plot-Pen-Style
          #:alpha Nonnegative-Real
          #:label (U String #f)]
         renderer2d))
(define (parametric f t-min t-max
                    #:x-min [x-min #f] #:x-max [x-max #f]
                    #:y-min [y-min #f] #:y-max [y-max #f]
                    #:samples [samples (line-samples)]
                    #:color [color (line-color)]
                    #:width [width (line-width)]
                    #:style [style (line-style)]
                    #:alpha [alpha (line-alpha)]
                    #:label [label #f])
  (define fail/pos (make-raise-argument-error 'parametric f t-min t-max))
  (define fail/kw (make-raise-keyword-error 'parametric))
  (cond
    [(not (rational? t-min))  (fail/pos "rational?" 1)]
    [(not (rational? t-max))  (fail/pos "rational?" 2)]
    [(and x-min (not (rational? x-min)))  (fail/kw "#f or rational" '#:x-min x-min)]
    [(and x-max (not (rational? x-max)))  (fail/kw "#f or rational" '#:x-max x-max)]
    [(and y-min (not (rational? y-min)))  (fail/kw "#f or rational" '#:y-min y-min)]
    [(and y-max (not (rational? y-max)))  (fail/kw "#f or rational" '#:y-max y-max)]
    [(< samples 2)  (fail/kw "Integer >= 2" '#:samples samples)]
    [(not (rational? width))  (fail/kw "rational?" '#:width width)]
    [(or (> alpha 1) (not (rational? alpha)))  (fail/kw "real in [0,1]" '#:alpha alpha)]
    [else
     (let ([f  (λ ([t : Real]) (sequence-head-vector 'parametric (f t) 2))])
       (lines (map f (linear-seq t-min t-max samples))
              #:x-min x-min #:x-max x-max #:y-min y-min #:y-max y-max
              #:color color #:width width #:style style #:alpha alpha
              #:label label))]))

(:: polar
    (->* [(-> Real Real)]
         [Real Real
          #:x-min (U Real #f) #:x-max (U Real #f)
          #:y-min (U Real #f) #:y-max (U Real #f)
          #:samples Positive-Integer
          #:color Plot-Color
          #:width Nonnegative-Real
          #:style Plot-Pen-Style
          #:alpha Nonnegative-Real
          #:label (U String #f)]
         renderer2d))
(define (polar f [θ-min 0] [θ-max (* 2 pi)]
               #:x-min [x-min #f] #:x-max [x-max #f]
               #:y-min [y-min #f] #:y-max [y-max #f]
               #:samples [samples (line-samples)]
               #:color [color (line-color)]
               #:width [width (line-width)]
               #:style [style (line-style)]
               #:alpha [alpha (line-alpha)]
               #:label [label #f])
  (define fail/pos (make-raise-argument-error 'polar f θ-min θ-max))
  (define fail/kw (make-raise-keyword-error 'polar))
  (cond
    [(not (rational? θ-min))  (fail/pos "rational?" 1)]
    [(not (rational? θ-max))  (fail/pos "rational?" 2)]
    [(and x-min (not (rational? x-min)))  (fail/kw "#f or rational" '#:x-min x-min)]
    [(and x-max (not (rational? x-max)))  (fail/kw "#f or rational" '#:x-max x-max)]
    [(and y-min (not (rational? y-min)))  (fail/kw "#f or rational" '#:y-min y-min)]
    [(and y-max (not (rational? y-max)))  (fail/kw "#f or rational" '#:y-max y-max)]
    [(< samples 2)  (fail/kw "Integer >= 2" '#:samples samples)]
    [(not (rational? width))  (fail/kw "rational?" '#:width width)]
    [(or (> alpha 1) (not (rational? alpha)))  (fail/kw "real in [0,1]" '#:alpha alpha)]
    [else
     (lines (let ([θs  (linear-seq θ-min θ-max samples)])
              (map polar->cartesian θs (map f θs)))
            #:x-min x-min #:x-max x-max #:y-min y-min #:y-max y-max
            #:color color #:width width #:style style #:alpha alpha
            #:label label)]))

;; ===================================================================================================
;; Function

(: function-render-proc (-> Sampler Positive-Integer
                            Plot-Color Nonnegative-Real Plot-Pen-Style
                            Nonnegative-Real
                            (U String #f)
                            2D-Render-Proc))
(define ((function-render-proc f samples color width style alpha label) area)
  (match-define (vector x-ivl y-ivl) (send area get-bounds-rect))
  (match-define (sample xs ys y-min y-max) (f x-ivl samples))
  
  (send area put-alpha alpha)
  (send area put-pen color width style)
  (send area put-lines (map (λ ([x : Real] [y : Real]) (vector x y)) xs ys))
  
  (cond [label  (line-legend-entry label color width style)]
        [else   empty]))

(:: function
    (->* [(-> Real Real)]
         [(U Real #f) (U Real #f)
          #:y-min (U Real #f) #:y-max (U Real #f)
          #:samples Positive-Integer
          #:color Plot-Color
          #:width Nonnegative-Real
          #:style Plot-Pen-Style
          #:alpha Nonnegative-Real
          #:label (U String #f)]
         renderer2d))
(define (function f [x-min #f] [x-max #f]
                  #:y-min [y-min #f] #:y-max [y-max #f]
                  #:samples [samples (line-samples)]
                  #:color [color (line-color)]
                  #:width [width (line-width)]
                  #:style [style (line-style)]
                  #:alpha [alpha (line-alpha)]
                  #:label [label #f])
  (define fail/pos (make-raise-argument-error 'function f x-min x-max))
  (define fail/kw (make-raise-keyword-error 'function))
  (cond
    [(and x-min (not (rational? x-min)))  (fail/pos "rational?" 1)]
    [(and x-max (not (rational? x-max)))  (fail/pos "rational?" 2)]
    [(and y-min (not (rational? y-min)))  (fail/kw "#f or rational" '#:y-min y-min)]
    [(and y-max (not (rational? y-max)))  (fail/kw "#f or rational" '#:y-max y-max)]
    [(< samples 2)  (fail/kw "Integer >= 2" '#:samples samples)]
    [(not (rational? width))  (fail/kw "rational?" '#:width width)]
    [(or (> alpha 1) (not (rational? alpha)))  (fail/kw "real in [0,1]" '#:alpha alpha)]
    [else
     (define x-ivl (ivl x-min x-max))
     (define y-ivl (ivl y-min y-max))
     (let ([f  (function->sampler f x-ivl)])
       (renderer2d (vector x-ivl y-ivl)
                   (function-bounds-fun f samples)
                   default-ticks-fun
                   (function-render-proc f samples color width style alpha label)))]))

;; ===================================================================================================
;; Inverse function

(: inverse-render-proc (-> Sampler Positive-Integer
                           Plot-Color Nonnegative-Real Plot-Pen-Style
                           Nonnegative-Real
                           (U String #f)
                           2D-Render-Proc))
(define ((inverse-render-proc f samples color width style alpha label) area)
  (match-define (vector x-ivl y-ivl) (send area get-bounds-rect))
  (match-define (sample ys xs x-min x-max) (f y-ivl samples))
  
  (send area put-alpha alpha)
  (send area put-pen color width style)
  (send area put-lines (map (λ ([x : Real] [y : Real]) (vector x y)) xs ys))
  
  (cond [label  (line-legend-entry label color width style)]
        [else   empty]))

(: inverse
   (->* [(-> Real Real)]
         [(U Real #f) (U Real #f)
          #:x-min (U Real #f) #:x-max (U Real #f)
          #:samples Positive-Integer
          #:color Plot-Color
          #:width Nonnegative-Real
          #:style Plot-Pen-Style
          #:alpha Nonnegative-Real
          #:label (U String #f)]
         renderer2d))
(define (inverse f [y-min #f] [y-max #f]
                 #:x-min [x-min #f] #:x-max [x-max #f]
                 #:samples [samples (line-samples)]
                 #:color [color (line-color)]
                 #:width [width (line-width)]
                 #:style [style (line-style)]
                 #:alpha [alpha (line-alpha)]
                 #:label [label #f])
  (define fail/pos (make-raise-argument-error 'inverse f y-min y-max))
  (define fail/kw (make-raise-keyword-error 'inverse))
  (cond
    [(and y-min (not (rational? y-min)))  (fail/pos "rational?" 1)]
    [(and y-max (not (rational? y-max)))  (fail/pos "rational?" 2)]
    [(and x-min (not (rational? x-min)))  (fail/kw "#f or rational" '#:x-min x-min)]
    [(and x-max (not (rational? x-max)))  (fail/kw "#f or rational" '#:x-max x-max)]
    [(< samples 2)  (fail/kw "Integer >= 2" '#:samples samples)]
    [(not (rational? width))  (fail/kw "rational?" '#:width width)]
    [(or (> alpha 1) (not (rational? alpha)))  (fail/kw "real in [0,1]" '#:alpha alpha)]
    [else
     (define x-ivl (ivl x-min x-max))
     (define y-ivl (ivl y-min y-max))
     (define g (inverse->sampler f y-ivl))
     (renderer2d (vector x-ivl y-ivl)
                 (inverse-bounds-fun g samples)
                 default-ticks-fun
                 (inverse-render-proc g samples color width style alpha label))]))

;; ===================================================================================================
;; Kernel density estimation

(:: density
    (->* [(Sequenceof Real)]
         [Nonnegative-Real
          (U (Sequenceof Real) #f)
          #:x-min (U Real #f) #:x-max (U Real #f)
          #:y-min (U Real #f) #:y-max (U Real #f)
          #:samples Positive-Integer
          #:color Plot-Color
          #:width Nonnegative-Real
          #:style Plot-Pen-Style
          #:alpha Nonnegative-Real
          #:label (U String #f)]
         renderer2d))
(define (density xs [bw-adjust 1] [orig-ws #f]
                 #:x-min [x-min #f] #:x-max [x-max #f]
                 #:y-min [y-min #f] #:y-max [y-max #f]
                 #:samples [samples (line-samples)]
                 #:color [color (line-color)]
                 #:width [width (line-width)]
                 #:style [style (line-style)]
                 #:alpha [alpha (line-alpha)]
                 #:label [label #f])
  (define fail/pos (make-raise-argument-error 'density xs bw-adjust orig-ws))
  (define fail/kw (make-raise-keyword-error 'density))
  (define ws (if orig-ws (sequence->list orig-ws) #f))
  (cond
    [(not (rational? bw-adjust))  (fail/pos "rational?" 1)]
    [(and ws (not (andmap (λ ([w : Real]) (and (not (negative? w)) (rational? w))) ws)))
     (fail/pos "sequence of nonnegative rationals" 2)]
    [(and x-min (not (rational? x-min)))  (fail/kw "#f or rational" '#:x-min x-min)]
    [(and x-max (not (rational? x-max)))  (fail/kw "#f or rational" '#:x-max x-max)]
    [(and y-min (not (rational? y-min)))  (fail/kw "#f or rational" '#:y-min y-min)]
    [(and y-max (not (rational? y-max)))  (fail/kw "#f or rational" '#:y-max y-max)]
    [(< samples 2)  (fail/kw "Integer >= 2" '#:samples samples)]
    [(not (rational? width))  (fail/kw "rational?" '#:width width)]
    [(or (> alpha 1) (not (rational? alpha)))  (fail/kw "real in [0,1]" '#:alpha alpha)]
    [else
     (let ([xs  (sequence->list xs)])
       (define n (length xs))
       (define sd (stddev xs ws))
       (define h (max 1e-308
                      (* 1e-14 (apply max (map abs (filter rational? xs))))
                      (* bw-adjust 1.06 sd (expt n -0.2))))
       (define-values (f fx-min fx-max) (kde xs h ws))
       (let ([x-min  (if x-min x-min fx-min)]
             [x-max  (if x-max x-max fx-max)])
         (function f x-min x-max #:y-min y-min #:y-max y-max #:samples samples
                   #:color color #:width width #:style style #:alpha alpha #:label label)))]))
