Mercurial > octave-nkf
comparison liboctave/dbleDET.cc @ 7789:82be108cc558
First attempt at single precision tyeps
* * *
corrections to qrupdate single precision routines
* * *
prefer demotion to single over promotion to double
* * *
Add single precision support to log2 function
* * *
Trivial PROJECT file update
* * *
Cache optimized hermitian/transpose methods
* * *
Add tests for tranpose/hermitian and ChangeLog entry for new transpose code
author | David Bateman <dbateman@free.fr> |
---|---|
date | Sun, 27 Apr 2008 22:34:17 +0200 |
parents | 2eb392d058bb |
children |
comparison
equal
deleted
inserted
replaced
7788:45f5faba05a2 | 7789:82be108cc558 |
---|---|
50 void | 50 void |
51 DET::initialize10 (void) | 51 DET::initialize10 (void) |
52 { | 52 { |
53 if (c2 != 0.0) | 53 if (c2 != 0.0) |
54 { | 54 { |
55 double etmp = e2 / xlog2 (10); | 55 double etmp = e2 / xlog2 (static_cast<double>(10)); |
56 e10 = static_cast<int> (xround (etmp)); | 56 e10 = static_cast<int> (xround (etmp)); |
57 etmp -= e10; | 57 etmp -= e10; |
58 c10 = c2 * pow (10.0, etmp); | 58 c10 = c2 * pow (10.0, etmp); |
59 } | 59 } |
60 } | 60 } |
72 } | 72 } |
73 | 73 |
74 double | 74 double |
75 DET::value (void) const | 75 DET::value (void) const |
76 { | 76 { |
77 return base2 ? c2 * xexp2 (e2) : c10 * pow (10.0, e10); | 77 return base2 ? c2 * xexp2 (static_cast<double>(e2)) : c10 * pow (10.0, e10); |
78 } | 78 } |
79 | 79 |
80 /* | 80 /* |
81 ;;; Local Variables: *** | 81 ;;; Local Variables: *** |
82 ;;; mode: C++ *** | 82 ;;; mode: C++ *** |