Mercurial > octave-nkf
annotate liboctave/DET.h @ 14200:64d9f33313cc stable rc-3-6-0-1
3.6.0-rc1 release candidate
* configure.ac (AC_INIT): Version is now 3.6.0-rc1.
(OCTAVE_RELEASE_DATE): Now 2012-01-12.
| author | John W. Eaton <jwe@octave.org> |
|---|---|
| date | Thu, 12 Jan 2012 14:31:50 -0500 |
| parents | 72c96de7a403 |
| children |
| rev | line source |
|---|---|
| 8335 | 1 /* |
| 2 | |
|
14138
72c96de7a403
maint: update copyright notices for 2012
John W. Eaton <jwe@octave.org>
parents:
11586
diff
changeset
|
3 Copyright (C) 2008-2012 Jaroslav Hajek |
| 8335 | 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 base_det | |
| 33 { | |
| 34 public: | |
| 35 | |
|
11505
9a308e96194e
more data member initialization fixes
John W. Eaton <jwe@octave.org>
parents:
9245
diff
changeset
|
36 base_det (T c = 1, int e = 0) |
|
9a308e96194e
more data member initialization fixes
John W. Eaton <jwe@octave.org>
parents:
9245
diff
changeset
|
37 : c2 (), e2 () |
|
11586
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
38 { |
|
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
39 c2 = xlog2 (c, e2); |
|
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
40 e2 += e; |
| 8335 | 41 } |
| 42 | |
|
11586
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
43 base_det (T c, double e, double b) |
|
11505
9a308e96194e
more data member initialization fixes
John W. Eaton <jwe@octave.org>
parents:
9245
diff
changeset
|
44 : c2 (), e2 () |
|
11586
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
45 { |
| 8335 | 46 e *= xlog2 (b); |
| 47 e2 = e; | |
| 48 c *= xexp2 (e - e2); | |
| 49 int f; | |
| 50 c2 = xlog2 (c, f); | |
| 51 e2 += f; | |
| 52 } | |
| 53 | |
|
11505
9a308e96194e
more data member initialization fixes
John W. Eaton <jwe@octave.org>
parents:
9245
diff
changeset
|
54 base_det (const base_det& a) : c2 (a.c2), e2 (a.e2) { } |
| 8335 | 55 |
| 56 base_det& operator = (const base_det& a) | |
| 57 { | |
| 58 c2 = a.c2; | |
| 59 e2 = a.e2; | |
| 60 return *this; | |
| 61 } | |
| 62 | |
| 63 T coef (void) const { return c2; } | |
| 64 int exp (void) const { return e2; } | |
| 65 | |
| 66 T value () const { return c2 * static_cast<T> (std::ldexp (1.0, e2)); } | |
| 67 operator T () const { return value (); } | |
| 68 | |
|
8336
9813c07ca946
make det take advantage of matrix type
Jaroslav Hajek <highegg@gmail.com>
parents:
8335
diff
changeset
|
69 base_det square () const { return base_det (c2*c2, e2+e2); } |
|
9813c07ca946
make det take advantage of matrix type
Jaroslav Hajek <highegg@gmail.com>
parents:
8335
diff
changeset
|
70 |
| 8335 | 71 void operator *= (T t) |
| 72 { | |
| 73 int e; | |
| 74 c2 *= xlog2 (t, e); | |
| 75 e2 += e; | |
| 76 } | |
| 77 | |
| 78 private: | |
| 79 | |
| 80 T c2; | |
| 81 int e2; | |
| 82 }; | |
| 83 | |
| 84 // Provide the old types by typedefs. | |
| 85 typedef base_det<double> DET; | |
| 86 typedef base_det<float> FloatDET; | |
| 87 typedef base_det<Complex> ComplexDET; | |
| 88 typedef base_det<FloatComplex> FloatComplexDET; | |
| 89 | |
| 90 #endif |