Source

caml2geom / caml2geom.ml

(*
  Copyright (c) 2011 Barry Schwartz

  Permission is hereby granted, free of charge, to any person
  obtaining a copy of this software and associated documentation
  files (the "Software"), to deal in the Software without
  restriction, including without limitation the rights to use,
  copy, modify, merge, publish, distribute, sublicense, and/or sell
  copies of the Software, and to permit persons to whom the
  Software is furnished to do so, subject to the following
  conditions:

  The above copyright notice and this permission notice shall be
  included in all copies or substantial portions of the Software.

  THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
  EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES
  OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
  NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT
  HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY,
  WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
  FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR
  OTHER DEALINGS IN THE SOFTWARE.
*)

module type Fragment =
sig
  type t
  type output_type
  val default : unit -> t
  val is_zero : t -> bool
  val is_constant : t -> bool
  val is_finite : t -> bool
  val at0 : t -> output_type
  val at1 : t -> output_type
  val value_at : t -> float -> output_type
end

module D2 =
struct
  module type S =
  sig
    type elt
    type t
    val make : elt -> elt -> t
    val coord : t -> int -> elt
    val fst : t -> elt
    val snd : t -> elt
    val is_zero : t -> bool
    val is_constant : t -> bool
    val is_finite : t -> bool
  end

  module Make(Elt : Fragment) =
  struct
    type elt = Elt.t
    type t = elt * elt
    let make a b = (a,b)
    let fst (a,_) = a
    let snd (_,b) = b
    let is_zero (a,b) = Elt.is_zero a && Elt.is_zero b
    let is_constant (a,b) = Elt.is_constant a && Elt.is_constant b
    let is_finite (a,b) = Elt.is_finite a && Elt.is_finite b

    let coord (a,b) i =
      match i with
        | 0 -> a
        | 1 -> b
        | _ -> invalid_arg "Make(Elt).coord: index out of range"
  end
end

module Coord =
struct
  type t = float
  type dim2 = X | Y
  external _get_epsilon : unit -> t = "epsilon_wrapper"
  let _eps_value = _get_epsilon ()
  let epsilon = _eps_value
  external infinity : unit -> t = "infinity_wrapper"
  external are_near : t -> t -> bool = "are_near_wrapper"
  external are_near_eps : t -> t -> t -> bool = "are_near_eps_wrapper"
  external rel_error_bound : t -> t -> bool = "rel_error_bound_wrapper"
  external rel_error_bound_eps : t -> t -> t -> bool = "rel_error_bound_eps_wrapper"
end

module Int_coord =
struct
  type t = int
end

module Point =
struct
  type t
  external default : unit -> t = "new_default_point_wrapper"
  external make : Coord.t -> Coord.t -> t = "new_point_wrapper"
  external copy : t -> t = "copy_point_wrapper"
  external coord : t -> int -> Coord.t = "point_coord_wrapper"
  external set_coord : t -> int -> Coord.t -> unit = "set_point_coord_wrapper"
  external length : t -> Coord.t = "point_length_wrapper"
  external normalize : t -> unit = "point_normalize_wrapper"
  external ccw : t -> t = "point_ccw_wrapper"
  external cw : t -> t = "point_cw_wrapper"
  external neg : t -> t = "point_neg_wrapper"
  external set : t -> t -> unit = "point_set_wrapper"
  external add_to : t -> t -> unit = "point_add_to_wrapper"
  external sub_from : t -> t -> unit = "point_sub_from_wrapper"
  external mul_by : t -> Coord.t -> unit = "point_mul_by_wrapper"
  external div_by : t -> Coord.t -> unit = "point_div_by_wrapper"
end

module Interval =
struct
  type t
  external default : unit -> t = "new_default_interval_wrapper"
end

module Rect =
struct
  type t
  external default : unit -> t = "new_default_rect_wrapper"
end

module Bezier =
struct
  type t
  type output_type = Coord.t
  external default : unit -> t = "new_default_bezier_wrapper"
  external make_order0 : Coord.t -> t = "new_bezier_order0_wrapper"
  external make_order1 : Coord.t -> Coord.t -> t = "new_bezier_order1_wrapper"
  external make_order2 : Coord.t -> Coord.t -> Coord.t -> t = "new_bezier_order2_wrapper"
  external make_order3 : Coord.t -> Coord.t -> Coord.t -> Coord.t -> t = "new_bezier_order3_wrapper"
  external is_zero : t -> bool = "bezier_is_zero_wrapper"
  external is_constant : t -> bool = "bezier_is_constant_wrapper"
  external is_finite : t -> bool = "bezier_is_finite_wrapper"
  external at0 : t -> Coord.t = "bezier_at0_wrapper"
  external at1 : t -> Coord.t = "bezier_at1_wrapper"
  external value_at : t -> float -> Coord.t = "bezier_value_at_wrapper"
end

module Bezier_as_fragment : Fragment = Bezier

module type Curve =
sig
  type t
  val initial_point : t -> Point.t
  val final_point : t -> Point.t
  val is_degenerate : t -> bool
  val point_at : t -> Coord.t -> Point.t
  val value_at : t -> Coord.t -> Coord.dim2 -> Coord.t
  val point_and_derivatives : t -> Coord.t -> int -> Point.t array
  val set_initial : t -> Point.t -> unit
  val set_final : t -> Point.t -> unit
  val portion : t -> Coord.t -> Coord.t -> t
  val reverse : t -> t
  val derivative : t -> t
  val length : t -> Coord.t -> Coord.t
  val degrees_of_freedom : t -> int
end

module Bezier_curve =
struct
  type t
  external of_two_beziers : Bezier.t -> Bezier.t -> t = "bezier_curve_of_two_beziers_wrapper"
  external of_two_points : Point.t -> Point.t -> t = "bezier_curve_of_two_points_wrapper"
  external of_three_points : Point.t -> Point.t -> Point.t -> t = "bezier_curve_of_three_points_wrapper"
  external of_four_points : Point.t -> Point.t -> Point.t -> Point.t -> t = "bezier_curve_of_four_points_wrapper"
  external initial_point : t -> Point.t = "bezier_curve_initial_point_wrapper"
  external final_point : t -> Point.t = "bezier_curve_final_point_wrapper"
  external is_degenerate : t -> bool = "bezier_curve_is_degenerate_wrapper"
  external point_at : t -> Coord.t -> Point.t = "bezier_curve_point_at_wrapper"
  external value_at : t -> Coord.t -> Coord.dim2 -> Coord.t = "bezier_curve_value_at_wrapper"
  external point_and_derivatives : t -> Coord.t -> int -> Point.t array = "bezier_curve_point_and_derivatives_wrapper"
  external set_initial : t -> Point.t -> unit = "bezier_curve_set_initial_wrapper"
  external set_final : t -> Point.t -> unit = "bezier_curve_set_final_wrapper"
  external portion : t -> Coord.t -> Coord.t -> t = "bezier_curve_portion_wrapper"
  external reverse : t -> t = "bezier_curve_reverse_wrapper"
  external derivative : t -> t = "bezier_curve_derivative_wrapper"
  external length : t -> Coord.t -> Coord.t = "bezier_curve_length_wrapper"
  external degrees_of_freedom : t -> int = "bezier_curve_degrees_of_freedom_wrapper"
end

module Bezier_curve_as_curve : Curve = Bezier_curve