annotate base/LogRange.cpp @ 282:d9319859a4cf tip

(none)
author benoitrigolleau
date Fri, 31 Oct 2008 11:00:24 +0000
parents fc9323a41f5a
children
rev   line source
lbajardsilogic@0 1 /* -*- c-basic-offset: 4 indent-tabs-mode: nil -*- vi:set ts=8 sts=4 sw=4: */
lbajardsilogic@0 2
lbajardsilogic@0 3 /*
lbajardsilogic@0 4 Sonic Visualiser
lbajardsilogic@0 5 An audio file viewer and annotation editor.
lbajardsilogic@0 6 Centre for Digital Music, Queen Mary, University of London.
lbajardsilogic@0 7 This file copyright 2006 Chris Cannam.
lbajardsilogic@0 8
lbajardsilogic@0 9 This program is free software; you can redistribute it and/or
lbajardsilogic@0 10 modify it under the terms of the GNU General Public License as
lbajardsilogic@0 11 published by the Free Software Foundation; either version 2 of the
lbajardsilogic@0 12 License, or (at your option) any later version. See the file
lbajardsilogic@0 13 COPYING included with this distribution for more information.
lbajardsilogic@0 14 */
lbajardsilogic@0 15
lbajardsilogic@0 16 #include "LogRange.h"
lbajardsilogic@0 17
lbajardsilogic@0 18 #include <algorithm>
lbajardsilogic@0 19 #include <cmath>
lbajardsilogic@0 20
lbajardsilogic@0 21 void
lbajardsilogic@0 22 LogRange::mapRange(float &min, float &max, float logthresh)
lbajardsilogic@0 23 {
lbajardsilogic@0 24 if (min > max) std::swap(min, max);
lbajardsilogic@0 25 if (max == min) max = min + 1;
lbajardsilogic@0 26
lbajardsilogic@0 27 if (min >= 0.f) {
lbajardsilogic@0 28
lbajardsilogic@0 29 max = log10f(max); // we know max != 0
lbajardsilogic@0 30
lbajardsilogic@0 31 if (min == 0.f) min = std::min(logthresh, max);
lbajardsilogic@0 32 else min = log10f(min);
lbajardsilogic@0 33
lbajardsilogic@0 34 } else if (max <= 0.f) {
lbajardsilogic@0 35
lbajardsilogic@0 36 min = log10f(-min); // we know min != 0
lbajardsilogic@0 37
lbajardsilogic@0 38 if (max == 0.f) max = std::min(logthresh, min);
lbajardsilogic@0 39 else max = log10f(-max);
lbajardsilogic@0 40
lbajardsilogic@0 41 std::swap(min, max);
lbajardsilogic@0 42
lbajardsilogic@0 43 } else {
lbajardsilogic@0 44
lbajardsilogic@0 45 // min < 0 and max > 0
lbajardsilogic@0 46
lbajardsilogic@0 47 max = log10f(std::max(max, -min));
lbajardsilogic@0 48 min = std::min(logthresh, max);
lbajardsilogic@0 49 }
lbajardsilogic@0 50
lbajardsilogic@0 51 if (min == max) min = max - 1;
lbajardsilogic@0 52 }
lbajardsilogic@0 53
lbajardsilogic@0 54 float
lbajardsilogic@0 55 LogRange::map(float value, float thresh)
lbajardsilogic@0 56 {
lbajardsilogic@0 57 if (value == 0.f) return thresh;
lbajardsilogic@0 58 return log10f(fabsf(value));
lbajardsilogic@0 59 }
lbajardsilogic@0 60