Mercurial > octave-nkf

comparison liboctave/DET.h @ 8335:64cf956a109c

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
templatize & fix DET
author Jaroslav Hajek <highegg@gmail.com>
date Wed, 19 Nov 2008 11:23:07 +0100
parents
children 9813c07ca946
comparison
equal deleted inserted replaced
8334:70dd33450061 8335:64cf956a109c
1 /*
2
3 Copyright (C) 2008 Jaroslav Hajek
4
5 This file is part of Octave.
6
7 Octave is free software; you can redistribute it and/or modify it
8 under the terms of the GNU General Public License as published by the
9 Free Software Foundation; either version 3 of the License, or (at your
10 option) any later version.
11
12 Octave is distributed in the hope that it will be useful, but WITHOUT
13 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
14 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
15 for more details.
16
17 You should have received a copy of the GNU General Public License
18 along with Octave; see the file COPYING. If not, see
19 <http://www.gnu.org/licenses/>.
20
21 */
22
23 #if !defined (octave_DET_h)
24 #define octave_DET_h 1
25
26 #include <cmath>
27 #include "oct-cmplx.h"
28 #include "lo-mappers.h"
29
30 template <class T>
31 class
32 OCTAVE_API
33 base_det
34 {
35 public:
36
37 base_det (T c = 0, int e = 0)
38 {
39 c2 = xlog2 (c, e2);
40 e2 += e;
41 }
42
43 base_det (T c, double e, double b)
44 {
45 e *= xlog2 (b);
46 e2 = e;
47 c *= xexp2 (e - e2);
48 int f;
49 c2 = xlog2 (c, f);
50 e2 += f;
51 }
52
53 base_det (const base_det& a) : c2(a.c2), e2(a.e2) { }
54
55 base_det& operator = (const base_det& a)
56 {
57 c2 = a.c2;
58 e2 = a.e2;
59 return *this;
60 }
61
62 T coef (void) const { return c2; }
63 int exp (void) const { return e2; }
64
65 T value () const { return c2 * static_cast<T> (std::ldexp (1.0, e2)); }
66 operator T () const { return value (); }
67
68 void operator *= (T t)
69 {
70 int e;
71 c2 *= xlog2 (t, e);
72 e2 += e;
73 }
74
75 private:
76
77 T c2;
78 int e2;
79 };
80
81 // Provide the old types by typedefs.
82 typedef base_det<double> DET;
83 typedef base_det<float> FloatDET;
84 typedef base_det<Complex> ComplexDET;
85 typedef base_det<FloatComplex> FloatComplexDET;
86
87 #endif