Log of /trunk/math/matrix/src/TMatrixTSym.cxx
Parent Directory
Revision
48472 -
(
view)
(
download)
(
as text)
(
annotate)
-
[select for diffs]
Modified
Tue Feb 5 15:54:36 2013 UTC (23 months, 2 weeks ago) by
rdm
File length: 72449 byte(s)
Diff to
previous 46214
From Tobias Schlueter:
fix for TMatrixDSym::TMatrixDSym(const TMatrixDSym&, EMatrixCreatorsOp2,
const TMatrixDSym&) doesn't size new matrix. Fixes #98605.
Revision
22885 -
(
view)
(
download)
(
as text)
(
annotate)
-
[select for diffs]
Modified
Fri Mar 28 13:57:25 2008 UTC (6 years, 9 months ago) by
rdm
File length: 71880 byte(s)
Diff to
previous 20882
move the following directories under the new "math" meta directory:
mathcore
mathmore
fftw
foam
fumili
genvector
matrix
minuit
minuit2
mlp
physics
smatrix
splot
unuran
quadp
Revision
17734 -
(
view)
(
download)
(
as text)
(
annotate)
-
[select for diffs]
Modified
Tue Feb 6 15:47:59 2007 UTC (7 years, 11 months ago) by
brun
Original Path:
trunk/matrix/src/TMatrixTSym.cxx
File length: 71944 byte(s)
Diff to
previous 17648
In the Streamer functions replace the lines like
TPad::Class()->ReadBuffer(b, this, v, R__s, R__c);
TPad::Class()->WriteBuffer(b,this);
by
b.ReadClassBuffer(TPad::Class(), this, v, R__s, R__c);
b.WriteClassBuffer(TPad::Class(),this);
Revision
16458 -
(
view)
(
download)
(
as text)
(
annotate)
-
[select for diffs]
Modified
Fri Oct 6 06:52:34 2006 UTC (8 years, 3 months ago) by
brun
Original Path:
trunk/matrix/src/TMatrixTSym.cxx
File length: 71917 byte(s)
Diff to
previous 16161
From Eddy Offermann:
Previously, the matrix package contained a large amount of ASSERT
statements which are a nuisance for programs analyzing a series of
independent events like in high-energy physics .
Assert's were issued when for instance a division by zero was requested
or a matrix was invalid . Most algorithms made matrices/vectors invalid
after an error occured in an operation , like inversion of a singular matrix .
Unfortunately, not all assert's were accompanied by error messages .
This situation has been completely overhauled :
- All error conditions in the algorithms are now accompanied by error
messages .
- In all algorithms it is still asserted that vectors/matrices are
valid BUT only in very few cases is a matrix/vector made invalid :
for instance if memory is allocated with incorrect parameters .
- In case of division by zero, the division is skipped . In case of a
singular matrix, the inversion routine returns the original matrix .
In the past the result of an inversion could be checked through the
value of the returned determinant or checking whether the inverted matrix
was valid .
Since from now on, we never make the matrix invalid in this operation, the
latter check will not indicate a singular matrix anymore .
The decompostion classes TDecomp... have a backward-compatible
change in the interface which makes detection of singularity easier :
old interface :
void Invert (TMatrixD &inv);
TMatrixD Invert ();
new interface :
Bool_t Invert (TMatrixD &inv);
TMatrixD Invert (Bool_t &status);
TMatrixD Invert () { Bool_t status; return Invert(status);
}
The returned status is kFALSE in case of singularity .
The old situation is easily reproduced by setting the ROOT variable
gErrorAbortLevel to kError . This cause an exception when there is an error
message (Error...) ., In the past the matrix would be made invalid which
would cause the next operation to throw an exception .
Revision
16161 -
(
view)
(
download)
(
as text)
(
annotate)
-
[select for diffs]
Modified
Wed Aug 30 12:54:13 2006 UTC (8 years, 4 months ago) by
brun
Original Path:
trunk/matrix/src/TMatrixTSym.cxx
File length: 72901 byte(s)
Diff to
previous 15123
From Eddy:
the promised update/patch for
TMatrixD A; ......
B.Use(A);
C.Use(A);
C += B;
Currently, It is checked in these kinds of operation whether the
objects
are the same . Obviously, that is not enough because objects could be
using the same data location . So I have to check if the data location
is identical .
Revision
14745 -
(
view)
(
download)
(
as text)
(
annotate)
-
[select for diffs]
Modified
Wed Apr 19 08:22:26 2006 UTC (8 years, 9 months ago) by
rdm
Original Path:
trunk/matrix/src/TMatrixTSym.cxx
File length: 70150 byte(s)
Diff to
previous 14537
Change the TError.h macros:
Assert -> R__ASSERT
Check -> R__CHECK
Change the TCollection.h macro:
ForEach -> R__FOR_EACH
This to avoid potential problems due too trivial macro names.
The old macros will be removed in the next release. Currently
they will print out warning messages with the advice to move
to the new macro names.
Revision
14537 -
(
view)
(
download)
(
as text)
(
annotate)
-
[select for diffs]
Modified
Tue Apr 4 05:51:06 2006 UTC (8 years, 9 months ago) by
brun
Original Path:
trunk/matrix/src/TMatrixTSym.cxx
File length: 70249 byte(s)
Diff to
previous 14460
From Eddy:
- 2 weeks ago AMultB.. functions were made public . However, unlike any
other
public function except for the *= calls, it could allocate an object
. In addition
already Mult function existed (which called AMultB with constr=0) .
This situation is rectified in the following way :
The public Mult, MultT and TMult are introduced .
AMultB , AtMultB and AMultBt are now functions that do not belong
to the
TMatrixT class anymore, they accept now pointers as arguments .
None of these functions will create a new object . In internal
calls to AMultB family ,
first a call to Allocate is made (if necessary) .
AMultB is used in one of the TMatrixTSym::Similarity calls,
removing the creation of a
TMatrixT object .
- smatrix/test function reflect the changes above
Revision
14460 -
(
view)
(
download)
(
as text)
(
annotate)
-
[select for diffs]
Modified
Wed Mar 29 05:16:49 2006 UTC (8 years, 9 months ago) by
brun
Original Path:
trunk/matrix/src/TMatrixTSym.cxx
File length: 70735 byte(s)
Diff to
previous 14401
From Eddy:
1) TVectorT TMatrixT TMatrixTSym TMatrixTSparse : made sure that
IsValid Assert's
are only done when gMatrixCheck is kTRUE
2) TMatrixTSym::Similarity and SimilarityT create a temporary matrix .
If storage
is <= kWorkMax use a local array . This way less temporary
elements are
created .
3) In TMatrixTSym::Similarity and SimilarityT some shape checking was
missing
Revision
14401 -
(
view)
(
download)
(
as text)
(
annotate)
-
[select for diffs]
Modified
Thu Mar 23 16:41:25 2006 UTC (8 years, 10 months ago) by
brun
Original Path:
trunk/matrix/src/TMatrixTSym.cxx
File length: 68876 byte(s)
Diff to
previous 14387
From Eddy:
a modification for TMatrixTSym::Invert :
Replaced the Bunch-Kaufman decomposition in the TMatrixTSym::Invert
by LU decomposition because it is signifivantly faster .
However, Bunch-Kaufman guarantees unlike LU that the inverted matrix
is symmetric . In case of numerically difficult matrices, the user
should use the TDecompBK class to invert.
Revision
14336 -
(
view)
(
download)
(
as text)
(
annotate)
-
[select for diffs]
Modified
Mon Mar 20 21:43:44 2006 UTC (8 years, 10 months ago) by
pcanal
Original Path:
trunk/matrix/src/TMatrixTSym.cxx
File length: 66267 byte(s)
Diff to
previous 13929
Reduce direct dependencies on TClass.h, TROOT.h and TStreamerInfo.h.
Warning: This means that some file that relied on the indirect
inclusion of these header file might now fail to compile with
an error message mention that gROOT is no known or that TClass,
TROOT or TStreamerInfo is incompletely defined. Simply add the
proper include directive.
Revision
13929 -
(
view)
(
download)
(
as text)
(
annotate)
-
[select for diffs]
Modified
Wed Jan 25 18:49:03 2006 UTC (9 years ago) by
brun
Original Path:
trunk/matrix/src/TMatrixTSym.cxx
File length: 66247 byte(s)
Diff to
previous 13770
From Eddy:
1) Added to LinkDef.h the NormalEqn function from TDecompChol . They
got
lost in this file during the transition to the templated version of
the matrix library .
2) Replaced in the templated code everywhere references to DB_EPSILON
and DBL_MIN
by the templated version through <numeric_limits> .
This form allows you to request diffs between any two revisions of this file.
For each of the two "sides" of the diff,
enter a numeric revision.
