Log of /trunk/math/matrix/inc/TVectorT.h
Parent Directory
Revision
30125 -
(
view)
(
download)
(
as text)
(
annotate)
-
[select for diffs]
Modified
Mon Sep 14 05:03:58 2009 UTC (5 years, 4 months ago) by
brun
File length: 17614 byte(s)
Diff to
previous 27658
From Eddy Offermann:
- added a Streamer function to TMatrixTSparse to remove a memory leak.
Updated LinkDef.h accordingly.
- added vector functionality as proposed by Peter D Barnes from LNL
(http://root.cern.ch/phpBB2/viewtopic.php?t=8351).
Mult function: double s = v1 * M * v2; // e.g., physics matrix element
OuterProduct : TMatrixD M = v1' * v2 // outer product of v1 and v2
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: 17278 byte(s)
Diff to
previous 22039
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
22039 -
(
view)
(
download)
(
as text)
(
annotate)
-
[select for diffs]
Modified
Thu Feb 7 05:48:31 2008 UTC (6 years, 11 months ago) by
brun
Original Path:
trunk/matrix/inc/TVectorT.h
File length: 17278 byte(s)
Diff to
previous 20882
From Eddy:
- remove bug in Decompose of all TDecompXXX classes : It is checked
whether
the matrix was already decomposed
- Added to TVector and the different TMatrix flavors the option to use
"Use" on constant arrys and vectors/matrices.
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/inc/TVectorT.h
File length: 15094 byte(s)
Diff to
previous 15181
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
15134 -
(
view)
(
download)
(
as text)
(
annotate)
-
[select for diffs]
Modified
Tue May 23 04:47:42 2006 UTC (8 years, 8 months ago) by
brun
Original Path:
trunk/matrix/inc/TVectorT.h
File length: 14998 byte(s)
Diff to
previous 15068
From Federico Carminati:
"I have implemented all copy and equal operators needed to silence all
warnings in AliRoot, as requested. I have implemented shallow copies as
would do the default operators synthetized by the compiler.
Most operators are protected. If users complain, you just have to move
them into the public area, but class derivation is of course supported.
It has been a terrible job, I have modified 278 files, but the changes
are backward compabile, and this goes a long way to permitting user to
use the effc++ flag with root headers."
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/inc/TVectorT.h
File length: 14857 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/inc/TVectorT.h
File length: 14845 byte(s)
Diff to
previous 14436
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
14436 -
(
view)
(
download)
(
as text)
(
annotate)
-
[select for diffs]
Modified
Tue Mar 28 10:00:17 2006 UTC (8 years, 9 months ago) by
brun
Original Path:
trunk/matrix/inc/TVectorT.h
File length: 14707 byte(s)
Diff to
previous 13939
Add two new functions
void Add(const TVectorT<Element> &v);
void Add(const TVectorT<Element> &v1, const TVectorT<Element> &v2);
The first function adds vector v to this.
The second function adds the 2 vectors v1 and v2 and store the result into this.
Revision
13939 -
(
view)
(
download)
(
as text)
(
annotate)
-
[select for diffs]
Modified
Thu Jan 26 16:31:01 2006 UTC (8 years, 11 months ago) by
brun
Original Path:
trunk/matrix/inc/TVectorT.h
File length: 14598 byte(s)
Diff to
previous 13929
From Eddy:
a patch that will hopefully resolve the Solaris issue .
From the error messages I deduce that the compiler does not like to
instantiate
a default templated argument in the function header . Therefore, I
replaced the default
value by the value 0.0 (which is an unusual value , the user would use
here "==")
and check in the routine and replace it with the appropriate epsilon.
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/inc/TVectorT.h
File length: 14670 byte(s)
Diff to
previous 13767
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.
