annotate fft/fftw/fftw-3.3.4/genfft/assoctable.ml @ 40:223f770b5341 kissfft-double tip

Try a double-precision kissfft
author Chris Cannam
date Wed, 07 Sep 2016 10:40:32 +0100
parents 26056e866c29
children
rev   line source
Chris@19 1 (*
Chris@19 2 * Copyright (c) 1997-1999 Massachusetts Institute of Technology
Chris@19 3 * Copyright (c) 2003, 2007-14 Matteo Frigo
Chris@19 4 * Copyright (c) 2003, 2007-14 Massachusetts Institute of Technology
Chris@19 5 *
Chris@19 6 * This program is free software; you can redistribute it and/or modify
Chris@19 7 * it under the terms of the GNU General Public License as published by
Chris@19 8 * the Free Software Foundation; either version 2 of the License, or
Chris@19 9 * (at your option) any later version.
Chris@19 10 *
Chris@19 11 * This program is distributed in the hope that it will be useful,
Chris@19 12 * but WITHOUT ANY WARRANTY; without even the implied warranty of
Chris@19 13 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
Chris@19 14 * GNU General Public License for more details.
Chris@19 15 *
Chris@19 16 * You should have received a copy of the GNU General Public License
Chris@19 17 * along with this program; if not, write to the Free Software
Chris@19 18 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
Chris@19 19 *
Chris@19 20 *)
Chris@19 21
Chris@19 22 (*************************************************************
Chris@19 23 * Functional associative table
Chris@19 24 *************************************************************)
Chris@19 25
Chris@19 26 (*
Chris@19 27 * this module implements a functional associative table.
Chris@19 28 * The table is parametrized by an equality predicate and
Chris@19 29 * a hash function, with the restriction that (equal a b) ==>
Chris@19 30 * hash a == hash b.
Chris@19 31 * The table is purely functional and implemented using a binary
Chris@19 32 * search tree (not balanced for now)
Chris@19 33 *)
Chris@19 34
Chris@19 35 type ('a, 'b) elem =
Chris@19 36 Leaf
Chris@19 37 | Node of int * ('a, 'b) elem * ('a, 'b) elem * ('a * 'b) list
Chris@19 38
Chris@19 39 let empty = Leaf
Chris@19 40
Chris@19 41 let lookup hash equal key table =
Chris@19 42 let h = hash key in
Chris@19 43 let rec look = function
Chris@19 44 Leaf -> None
Chris@19 45 | Node (hash_key, left, right, this_list) ->
Chris@19 46 if (hash_key < h) then look left
Chris@19 47 else if (hash_key > h) then look right
Chris@19 48 else let rec loop = function
Chris@19 49 [] -> None
Chris@19 50 | (a, b) :: rest -> if (equal key a) then Some b else loop rest
Chris@19 51 in loop this_list
Chris@19 52 in look table
Chris@19 53
Chris@19 54 let insert hash key value table =
Chris@19 55 let h = hash key in
Chris@19 56 let rec ins = function
Chris@19 57 Leaf -> Node (h, Leaf, Leaf, [(key, value)])
Chris@19 58 | Node (hash_key, left, right, this_list) ->
Chris@19 59 if (hash_key < h) then
Chris@19 60 Node (hash_key, ins left, right, this_list)
Chris@19 61 else if (hash_key > h) then
Chris@19 62 Node (hash_key, left, ins right, this_list)
Chris@19 63 else
Chris@19 64 Node (hash_key, left, right, (key, value) :: this_list)
Chris@19 65 in ins table