diff -ur root-6.06.02.orig/misc/table/inc/TCernLib.h root-6.06.02/misc/table/inc/TCernLib.h --- root-6.06.02.orig/misc/table/inc/TCernLib.h 2016-03-03 10:36:03.000000000 +0100 +++ root-6.06.02/misc/table/inc/TCernLib.h 2016-03-07 19:15:34.726025891 +0100 @@ -15,7 +15,7 @@ #include "Rtypes.h" #include -// http://wwwinfo.cern.ch/asdoc/shortwrupsdir/f110/top.html +// https://cern-tex.web.cern.ch/cern-tex/shortwrupsdir/f110/top.html /////////////////////////////////////////////////////////////////////////////////////// // // @@ -25,8 +25,9 @@ // its negative, or a multiple of it, transpose a given matrix, build up a unit // // matrix, multiply a matrix by a diagonal (from left or from right) and may // // add the result to another matrix, add to square matrix the multiple of a diagonal // -// matrix, compute the products tex2html_wrap_inline191 (tex2html_wrap_inline193 denotes the transpose of row-wise without gaps end_html without gaps. // +// matrix, compute the products X=ABA' (A' denotes the transpose of A) and X=A'BA. // +// It is assumed that matrices are row-wise without gaps, contrary to the // +// Fortran convention. // // // /////////////////////////////////////////////////////////////////////////////////////// diff -ur root-6.06.02.orig/misc/table/src/TCernLib.cxx root-6.06.02/misc/table/src/TCernLib.cxx --- root-6.06.02.orig/misc/table/src/TCernLib.cxx 2016-03-03 10:36:03.000000000 +0100 +++ root-6.06.02/misc/table/src/TCernLib.cxx 2016-03-07 19:15:34.728025914 +0100 @@ -11,8 +11,8 @@ //////////////////////////////////////////////////////////////////////////////// // The set of methods to work with the plain matrix / vector -// "derived" from http://wwwinfo.cern.ch/asdoc/shortwrupsdir/f110/top.html -// "derived" from http://wwwinfo.cern.ch/asdoc/shortwrupsdir/f112/top.html +// "derived" from https://cern-tex.web.cern.ch/cern-tex/shortwrupsdir/f110/top.html +// "derived" from https://cern-tex.web.cern.ch/cern-tex/shortwrupsdir/f112/top.html // // Revision 1.7 2006/05/21 18:05:26 brun // Fix more coding conventions violations @@ -33,7 +33,7 @@ // // Revision 1.3 2003/04/03 17:39:39 fine // Make merge with ROOT 3.05.03 and add TR package -//122 +// // Revision 1.2 2003/02/04 23:35:20 fine // Clean up // @@ -154,27 +154,24 @@ //////////////////////////////////////////////////////////////////////////////// /// Matrix Multiplication +/// /// CERN PROGLIB# F110 MXMLRT .VERSION KERNFOR 2.00 720707 /// ORIG. 01/01/64 RKB -///BEGIN_HTML - see original documentation of CERNLIB package F110 - END_HTML - +// -- ENTRY MXMLRT +// -- C = A(I,J) X B(J,J) X A*(J,I) +// -- A* TANDS FOR A-TRANSPOSED +// mxmlrt (A,B,C,NI,NJ) IS EQUIVALENT TO +// CALL MXMPY (A,B,X,NI,NJ,NJ) +// CALL MXMPY1 (X,A,C,NI,NJ,NI) -// -- ENTRY MXMLRT */ -// -- C = A(I,J) X B(J,J) X A*(J,I) */ -// -- A* TANDS FOR A-TRANSPOSED */ -// mxmlrt (A,B,C,NI,NJ) IS EQUIVALENT TO */ -// CALL MXMPY (A,B,X,NI,NJ,NJ) */ -// CALL MXMPY1 (X,A,C,NI,NJ,NI) */ - -/* OR CALL MXMPY1 (B,A,Y,NJ,NJ,NI) */ -/* CALL MXMPY (A,Y,C,NI,NJ,NI) */ +// OR CALL MXMPY1 (B,A,Y,NJ,NJ,NI) +// CALL MXMPY (A,Y,C,NI,NJ,NI) // -- C = A*(I,J) X B(J,J) X A(J,I) @@ -191,6 +188,12 @@ //////////////////////////////////////////////////////////////////////////////// /// Matrix Multiplication (double precision) +/// +/// CERN PROGLIB# F110 MXMLRT .VERSION KERNFOR 2.00 720707 +/// ORIG. 01/01/64 RKB +/// +/// See original documentation of CERNLIB package +/// [F110](https://cern-tex.web.cern.ch/cern-tex/shortwrupsdir/f110/top.html) double *TCL::mxmlrt_0_(int n__, const double *a, const double *b, double *c, int ni,int nj) { @@ -215,36 +218,31 @@ for (int l = 1; l <= i; ++l,ia += j,++ib) b[ib] = a[ia]; } //////////////////////////////////////////////////////////////////////////////// +/// Matrix Transposition /// -/// Matrix Transposition /// CERN PROGLIB# F110 MXTRP .VERSION KERNFOR 1.0 650809 /// ORIG. 01/01/64 RKB -///BEGIN_HTML - see original documentation of CERNLIB package F110 - END_HTML - TCL_MXTRP(a, b, i, j) return b; } /* mxtrp */ //////////////////////////////////////////////////////////////////////////////// -/// Matrix Transposition (double precision) +/// Matrix Transposition (double precision) +/// /// CERN PROGLIB# F110 MXTRP .VERSION KERNFOR 1.0 650809 /// ORIG. 01/01/64 RKB -///BEGIN_HTML - see original documentation of CERNLIB package F110 - END_HTML - TCL_MXTRP(a, b, i, j) return b; @@ -283,38 +281,35 @@ //////////////////////////////////////////////////////////////////////////////// -/// /// Symmetric Multiplication of Rectangular Matrices +/// /// CERN PROGLIB# F112 TRAAT .VERSION KERNFOR 4.15 861204 -/// ORIG. 18/12/74 WH */ +/// ORIG. 18/12/74 WH +/// /// traat.F -- translated by f2c (version 19970219). /// -///BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML TCL_TRAAT(a, s, m, n) return s; } /* traat_ */ //////////////////////////////////////////////////////////////////////////////// -/// Symmetric Multiplication of Rectangular Matrices +/// Symmetric Multiplication of Rectangular Matrices +/// /// CERN PROGLIB# F112 TRAAT .VERSION KERNFOR 4.15 861204 -/// ORIG. 18/12/74 WH */ +/// ORIG. 18/12/74 WH +/// /// traat.F -- translated by f2c (version 19970219). /// -///BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML TCL_TRAAT(a, s, m, n) return s; } /* traat_ */ @@ -346,34 +341,34 @@ //////////////////////////////////////////////////////////////////////////////// /// Triangular - Rectangular Multiplication +/// /// CERN PROGLIB# F112 TRAL .VERSION KERNFOR 4.15 861204 /// ORIG. 18/12/74 WH +/// /// tral.F -- translated by f2c (version 19970219). -///BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML TCL_TRAL(a, u, b, m, n) return b; } /* tral_ */ //////////////////////////////////////////////////////////////////////////////// /// Triangular - Rectangular Multiplication +/// +/// CERN PROGLIB# F112 TRAL .VERSION KERNFOR 4.15 861204 +/// ORIG. 18/12/74 WH +/// /// tral.F -- translated by f2c (version 19970219). -/// CERN PROGLIB# F112 TRAL .VERSION KERNFOR 4.15 861204 */ -/// ORIG. 18/12/74 WH */ -///BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML TCL_TRAL(a, u, b, m, n) return b; } /* tral_ */ @@ -405,34 +400,34 @@ //////////////////////////////////////////////////////////////////////////////// /// Triangular - Rectangular Multiplication +/// /// CERN PROGLIB# F112 TRALT .VERSION KERNFOR 4.15 861204 /// ORIG. 18/12/74 WH +/// /// tralt.F -- translated by f2c (version 19970219). -///BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML TCL_TRALT(a, u, b, m, n) return b; } /* tralt_ */ //////////////////////////////////////////////////////////////////////////////// /// Triangular - Rectangular Multiplication +/// /// CERN PROGLIB# F112 TRALT .VERSION KERNFOR 4.15 861204 /// ORIG. 18/12/74 WH +/// /// tralt.F -- translated by f2c (version 19970219). -///BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML TCL_TRALT(a, u, b, m, n) return b; } /* tralt_ */ @@ -470,34 +465,34 @@ //////////////////////////////////////////////////////////////////////////////// /// Symmetric - Rectangular Multiplication -/// CERN PROGLIB# F112 TRAS .VERSION KERNFOR 4.15 861204 */ -/// ORIG. 18/12/74 WH */ +/// +/// CERN PROGLIB# F112 TRAS .VERSION KERNFOR 4.15 861204 +/// ORIG. 18/12/74 WH +/// /// tras.F -- translated by f2c (version 19970219). -///BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML TCL_TRAS(a, s, b, m, n) return b; } /* tras_ */ //////////////////////////////////////////////////////////////////////////////// /// Symmetric - Rectangular Multiplication -/// CERN PROGLIB# F112 TRAS .VERSION KERNFOR 4.15 861204 */ -/// ORIG. 18/12/74 WH */ +/// +/// CERN PROGLIB# F112 TRAS .VERSION KERNFOR 4.15 861204 +/// ORIG. 18/12/74 WH +/// /// tras.F -- translated by f2c (version 19970219). -///BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML TCL_TRAS(a, s, b, m, n) return b; } /* tras_ */ @@ -543,51 +538,51 @@ //////////////////////////////////////////////////////////////////////////////// /// Transformation of Symmetric Matrix -/// CERN PROGLIB# F112 TRASAT .VERSION KERNFOR 4.15 861204 */ -/// ORIG. 18/12/74 WH */ +/// +/// CERN PROGLIB# F112 TRASAT .VERSION KERNFOR 4.15 861204 +/// ORIG. 18/12/74 WH +/// /// trasat.F -- translated by f2c (version 19970219). -///BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML TCL_TRASAT(a, s, r__, m, n) return r__; } /* trasat_ */ //////////////////////////////////////////////////////////////////////////////// /// Transformation of Symmetric Matrix -/// CERN PROGLIB# F112 TRASAT .VERSION KERNFOR 4.15 861204 */ -/// ORIG. 18/12/74 WH */ +/// +/// CERN PROGLIB# F112 TRASAT .VERSION KERNFOR 4.15 861204 +/// ORIG. 18/12/74 WH +/// /// trasat.F -- translated by f2c (version 19970219). -///BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML TCL_TRASAT(a, s, r__, m, n) return r__; } /* trasat_ */ //////////////////////////////////////////////////////////////////////////////// /// Transformation of Symmetric Matrix -/// CERN PROGLIB# F112 TRASAT .VERSION KERNFOR 4.15 861204 */ -/// ORIG. 18/12/74 WH */ +/// +/// CERN PROGLIB# F112 TRASAT .VERSION KERNFOR 4.15 861204 +/// ORIG. 18/12/74 WH +/// /// trasat.F -- translated by f2c (version 19970219). -///BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML TCL_TRASAT(a, s, r__, m, n) return r__; } /* trasat_ */ @@ -596,17 +591,15 @@ //////////////////////////////////////////////////////////////////////////////// /// trata.F -- translated by f2c (version 19970219). -/// CERN PROGLIB# F112 TRATA .VERSION KERNFOR 4.15 861204 */ -/// ORIG. 18/12/74 WH */ -///BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML - /* Local variables */ int i__, j, ia, mn, ir, iat; double sum; @@ -636,22 +629,21 @@ return r__; } /* trata_ */ -//____________________________________________________________ -// trats.F -- translated by f2c (version 19970219). +//////////////////////////////////////////////////////////////////////////////// +/// trats.F -- translated by f2c (version 19970219). +/// +/// CERN PROGLIB# F112 TRATS .VERSION KERNFOR 4.15 861204 +/// ORIG. 18/12/74 WH +/// +/// See original documentation of CERNLIB package +/// [F112](https://cern-tex.web.cern.ch/cern-tex/shortwrupsdir/f112/top.html) + float *TCL::trats(const float *a, const float *s, float *b, int m, int n) { - //BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML /* Local variables */ int inds, i__, j, k, ia, ib, is; double sum; - /* CERN PROGLIB# F112 TRATS .VERSION KERNFOR 4.15 861204 */ - /* ORIG. 18/12/74 WH */ - /* Parameter adjuTments */ --b; --s; --a; @@ -684,25 +676,22 @@ return b; } /* trats_ */ -//____________________________________________________________ -// tratsa.F -- translated by f2c (version 19970219). -/* Subroutine */float *TCL::tratsa(const float *a, const float *s, float *r__, int m, int n) -{ - //BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML +//////////////////////////////////////////////////////////////////////////////// +/// tratsa.F -- translated by f2c (version 19970219). +/// +/// CERN PROGLIB# F112 TRATSA .VERSION KERNFOR 4.15 861204 +/// ORIG. 18/12/74 WH +/// +/// See original documentation of CERNLIB package +/// [F112](https://cern-tex.web.cern.ch/cern-tex/shortwrupsdir/f112/top.html) +float *TCL::tratsa(const float *a, const float *s, float *r__, int m, int n) +{ /* Local variables */ int imax, i__, j, k; int ia, ir, is, iaa, ind; double sum; - /* CERN PROGLIB# F112 TRATSA .VERSION KERNFOR 4.15 861204 */ - /* ORIG. 18/12/74 WH */ - - /* Parameter adjuTments */ --r__; --s; --a; @@ -743,26 +732,23 @@ return r__; } /* tratsa_ */ -//____________________________________________________________ -// trchlu.F -- translated by f2c (version 19970219). +//////////////////////////////////////////////////////////////////////////////// +/// trchlu.F -- translated by f2c (version 19970219). +/// +/// CERN PROGLIB# F112 TRCHLU .VERSION KERNFOR 4.16 870601 +/// ORIG. 18/12/74 W.HART +/// +/// See original documentation of CERNLIB package +/// [F112](https://cern-tex.web.cern.ch/cern-tex/shortwrupsdir/f112/top.html) + float *TCL::trchlu(const float *a, float *b, int n) { - //BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML /* Local variables */ int ipiv, kpiv, i__, j; double r__, dc; int id, kd; double sum; - - /* CERN PROGLIB# F112 TRCHLU .VERSION KERNFOR 4.16 870601 */ - /* ORIG. 18/12/74 W.HART */ - - /* Parameter adjuTments */ --b; --a; @@ -806,15 +792,17 @@ return b; } /* trchlu_ */ -//____________________________________________________________ -// trchul.F -- translated by f2c (version 19970219). -/* Subroutine */float *TCL::trchul(const float *a, float *b, int n) +//////////////////////////////////////////////////////////////////////////////// +/// trchul.F -- translated by f2c (version 19970219). +/// +/// CERN PROGLIB# F112 TRCHUL .VERSION KERNFOR 4.16 870601 +/// ORIG. 18/12/74 WH +/// +/// See original documentation of CERNLIB package +/// [F112](https://cern-tex.web.cern.ch/cern-tex/shortwrupsdir/f112/top.html) + +float *TCL::trchul(const float *a, float *b, int n) { - //BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML /* Local variables */ int ipiv, kpiv, i__; double r__; @@ -823,11 +811,6 @@ int id, kd; double sum; - - /* CERN PROGLIB# F112 TRCHUL .VERSION KERNFOR 4.16 870601 */ - /* ORIG. 18/12/74 WH */ - - /* Parameter adjuTments */ --b; --a; @@ -875,17 +858,15 @@ //////////////////////////////////////////////////////////////////////////////// /// trinv.F -- translated by f2c (version 19970219). -/// CERN PROGLIB# F112 TRINV .VERSION KERNFOR 4.15 861204 */ -/// ORIG. 18/12/74 WH */ -///BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML +/// +/// CERN PROGLIB# F112 TRINV .VERSION KERNFOR 4.15 861204 +/// ORIG. 18/12/74 WH +/// +/// See original documentation of CERNLIB package +/// [F112](https://cern-tex.web.cern.ch/cern-tex/shortwrupsdir/f112/top.html) +float *TCL::trinv(const float *t, float *s, int n) +{ int lhor, ipiv, lver, j; double sum = 0; double r__ = 0; @@ -934,22 +915,20 @@ return s; } /* trinv_ */ -//____________________________________________________________ -// trla.F -- translated by f2c (version 19970219). -/* Subroutine */float *TCL::trla(const float *u, const float *a, float *b, int m, int n) +//////////////////////////////////////////////////////////////////////////////// +/// trla.F -- translated by f2c (version 19970219). +/// +/// CERN PROGLIB# F112 TRLA .VERSION KERNFOR 4.15 861204 +/// ORIG. 18/12/74 WH +/// +/// See original documentation of CERNLIB package +/// [F112](https://cern-tex.web.cern.ch/cern-tex/shortwrupsdir/f112/top.html) + +float *TCL::trla(const float *u, const float *a, float *b, int m, int n) { int ipiv, ia, ib, iu; double sum; - /* CERN PROGLIB# F112 TRLA .VERSION KERNFOR 4.15 861204 */ - /* ORIG. 18/12/74 WH */ - //BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML - - /* Parameter adjuTments */ --b; --a; --u; @@ -981,22 +960,19 @@ } /* trla_ */ //////////////////////////////////////////////////////////////////////////////// +/// trlta.F -- translated by f2c (version 19970219). +/// +/// CERN PROGLIB# F112 TRLTA .VERSION KERNFOR 4.15 861204 +/// ORIG. 18/12/74 WH +/// +/// See original documentation of CERNLIB package +/// [F112](https://cern-tex.web.cern.ch/cern-tex/shortwrupsdir/f112/top.html) -/* trlta.F -- translated by f2c (version 19970219). -// Subroutine */float *TCL::trlta(const float *u, const float *a, float *b, int m, int n) +float *TCL::trlta(const float *u, const float *a, float *b, int m, int n) { int ipiv, mxpn, i__, nTep, ia, ib, iu, mx; double sum; - /* CERN PROGLIB# F112 TRLTA .VERSION KERNFOR 4.15 861204 */ - /* ORIG. 18/12/74 WH */ - //BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML - - /* Parameter adjuTments */ --b; --a; --u; @@ -1036,16 +1012,15 @@ //////////////////////////////////////////////////////////////////////////////// /// trpck.F -- translated by f2c (version 19970219). -/// CERN PROGLIB# F112 TRPCK .VERSION KERNFOR 2.08 741218 */ -/// ORIG. 18/12/74 WH */ -///BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML int i__, ia, ind, ipiv; /* Parameter adjuTments */ @@ -1072,17 +1047,15 @@ //////////////////////////////////////////////////////////////////////////////// /// trqsq.F -- translated by f2c (version 19970219). -/// CERN PROGLIB# F112 TRQSQ .VERSION KERNFOR 4.15 861204 */ -/// ORIG. 18/12/74 WH */ -///BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML - int indq, inds, imax, i__, j, k, l; int iq, ir, is, iqq; double sum; @@ -1140,16 +1113,15 @@ //////////////////////////////////////////////////////////////////////////////// /// trsa.F -- translated by f2c (version 19970219). -/// CERN PROGLIB# F112 TRSA .VERSION KERNFOR 4.15 861204 */ -/// ORIG. 18/12/74 WH */ -///BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML /* Local variables */ int inds, i__, j, k, ia, ib, is; double sum; @@ -1190,17 +1162,15 @@ //////////////////////////////////////////////////////////////////////////////// /// trsinv.F -- translated by f2c (version 19970219). +/// /// CERN PROGLIB# F112 TRSINV .VERSION KERNFOR 2.08 741218 -/// ORIG. 18/12/74 WH */ -///BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML - /* Function Body */ trchlu(g, gi, n); trinv(gi, gi, n); @@ -1209,17 +1179,15 @@ //////////////////////////////////////////////////////////////////////////////// /// trsmlu.F -- translated by f2c (version 19970219). -/// CERN PROGLIB# F112 TRSMLU .VERSION KERNFOR 4.15 861204 */ -/// ORIG. 18/12/74 WH */ -///BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML +/// +/// CERN PROGLIB# F112 TRSMLU .VERSION KERNFOR 4.15 861204 +/// ORIG. 18/12/74 WH +/// +/// See original documentation of CERNLIB package +/// [F112](https://cern-tex.web.cern.ch/cern-tex/shortwrupsdir/f112/top.html) +float *TCL::trsmlu(const float *u, float *s, int n) +{ /* Local variables */ int lhor, lver, i__, k, l, ind; double sum; @@ -1246,17 +1214,15 @@ //////////////////////////////////////////////////////////////////////////////// /// trsmul.F -- translated by f2c (version 19970219). -/// CERN PROGLIB# F112 TRSMUL .VERSION KERNFOR 4.15 861204 */ -/// ORIG. 18/12/74 WH */ -///BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML +/// +/// CERN PROGLIB# F112 TRSMUL .VERSION KERNFOR 4.15 861204 +/// ORIG. 18/12/74 WH +/// +/// See original documentation of CERNLIB package +/// [F112](https://cern-tex.web.cern.ch/cern-tex/shortwrupsdir/f112/top.html) +float *TCL::trsmul(const float *g, float *gi, int n) +{ /* Local variables */ int lhor, lver, lpiv, i__, j, k, ind; double sum; @@ -1284,18 +1250,15 @@ //////////////////////////////////////////////////////////////////////////////// /// trupck.F -- translated by f2c (version 19970219). +/// /// CERN PROGLIB# F112 TRUPCK .VERSION KERNFOR 2.08 741218 /// ORIG. 18/12/74 WH -///BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML - - int i__, im, is, iu, iv, ih, m2; /* Parameter adjuTments */ @@ -1336,25 +1299,20 @@ } /* trupck_ */ //////////////////////////////////////////////////////////////////////////////// +/// trsat.F -- translated by f2c (version 19970219). +/// +/// CERN PROGLIB# F112 TRSAT .VERSION KERNFOR 4.15 861204 +/// ORIG. 18/12/74 WH +/// +/// See original documentation of CERNLIB package +/// [F112](https://cern-tex.web.cern.ch/cern-tex/shortwrupsdir/f112/top.html) -/* trsat.F -- translated by f2c (version 19970219). -// Subroutine */ float *TCL::trsat(const float *s, const float *a, float *b, int m, int n) +float *TCL::trsat(const float *s, const float *a, float *b, int m, int n) { - //BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML - /* Local variables */ int inds, i__, j, k, ia, ib, is; double sum; - - /* CERN PROGLIB# F112 TRSAT .VERSION KERNFOR 4.15 861204 */ - /* ORIG. 18/12/74 WH */ - - /* Parameter adjuTments */ --b; --a; --s; @@ -1391,25 +1349,21 @@ // ------ double -//____________________________________________________________ -// trata.F -- translated by f2c (version 19970219). +//////////////////////////////////////////////////////////////////////////////// +/// trata.F -- translated by f2c (version 19970219). +/// +/// CERN PROGLIB# F112 TRATA .VERSION KERNFOR 4.15 861204 +/// ORIG. 18/12/74 WH +/// +/// See original documentation of CERNLIB package +/// [F112](https://cern-tex.web.cern.ch/cern-tex/shortwrupsdir/f112/top.html) + double *TCL::trata(const double *a, double *r__, int m, int n) { - //BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML - /* Local variables */ int i__, j, ia, mn, ir, iat; double sum; - - /* CERN PROGLIB# F112 TRATA .VERSION KERNFOR 4.15 861204 */ - /* ORIG. 18/12/74 WH */ - - /* Parameter adjuTments */ --r__; --a; @@ -1437,23 +1391,21 @@ return 0; } /* trata_ */ -//____________________________________________________________ -// trats.F -- translated by f2c (version 19970219). +//////////////////////////////////////////////////////////////////////////////// +/// trats.F -- translated by f2c (version 19970219). +/// +/// CERN PROGLIB# F112 TRATS .VERSION KERNFOR 4.15 861204 +/// ORIG. 18/12/74 WH +/// +/// See original documentation of CERNLIB package +/// [F112](https://cern-tex.web.cern.ch/cern-tex/shortwrupsdir/f112/top.html) + double *TCL::trats(const double *a, const double *s, double *b, int m, int n) { - //BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML /* Local variables */ int inds, i__, j, k, ia, ib, is; double sum; - - /* CERN PROGLIB# F112 TRATS .VERSION KERNFOR 4.15 861204 */ - /* ORIG. 18/12/74 WH */ - /* Parameter adjuTments */ --b; --s; --a; @@ -1487,24 +1439,22 @@ return 0; } /* trats_ */ -//____________________________________________________________ -// tratsa.F -- translated by f2c (version 19970219). -/* Subroutine */double *TCL::tratsa(const double *a, const double *s, double *r__, int m, int n) +//////////////////////////////////////////////////////////////////////////////// +/// tratsa.F -- translated by f2c (version 19970219). +/// +/// CERN PROGLIB# F112 TRATSA .VERSION KERNFOR 4.15 861204 +/// ORIG. 18/12/74 WH +/// +/// See original documentation of CERNLIB package +/// [F112](https://cern-tex.web.cern.ch/cern-tex/shortwrupsdir/f112/top.html) + +double *TCL::tratsa(const double *a, const double *s, double *r__, int m, int n) { - //BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML /* Local variables */ int imax, i__, j, k; int ia, ir, is, iaa, ind; double sum; - /* CERN PROGLIB# F112 TRATSA .VERSION KERNFOR 4.15 861204 */ - /* ORIG. 18/12/74 WH */ - - /* Parameter adjuTments */ --r__; --s; --a; @@ -1547,25 +1497,21 @@ //////////////////////////////////////////////////////////////////////////////// /// trchlu.F -- translated by f2c (version 19970219). -///BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML /* Local variables */ int ipiv, kpiv, i__, j; double r__, dc; int id, kd; double sum; - - /* CERN PROGLIB# F112 TRCHLU .VERSION KERNFOR 4.16 870601 */ - /* ORIG. 18/12/74 W.HART */ - - /* Parameter adjuTments */ --b; --a; @@ -1609,15 +1555,17 @@ return 0; } /* trchlu_ */ -//____________________________________________________________ -// trchul.F -- translated by f2c (version 19970219). +//////////////////////////////////////////////////////////////////////////////// +/// trchul.F -- translated by f2c (version 19970219). +/// +/// CERN PROGLIB# F112 TRCHUL .VERSION KERNFOR 4.16 870601 +/// ORIG. 18/12/74 WH +/// +/// See original documentation of CERNLIB package +/// [F112](https://cern-tex.web.cern.ch/cern-tex/shortwrupsdir/f112/top.html) + double *TCL::trchul(const double *a, double *b, int n) { - //BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML /* Local variables */ int ipiv, kpiv, i__; double r__; @@ -1626,11 +1574,6 @@ int id, kd; double sum; - - /* CERN PROGLIB# F112 TRCHUL .VERSION KERNFOR 4.16 870601 */ - /* ORIG. 18/12/74 WH */ - - /* Parameter adjuTments */ --b; --a; @@ -1677,17 +1620,15 @@ //////////////////////////////////////////////////////////////////////////////// /// trinv.F -- translated by f2c (version 19970219). -/// CERN PROGLIB# F112 TRINV .VERSION KERNFOR 4.15 861204 */ -/// ORIG. 18/12/74 WH */ /// -///BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML int lhor, ipiv, lver, j; double r__; int mx, ndTep, ind; @@ -1735,19 +1676,16 @@ } /* trinv_ */ //////////////////////////////////////////////////////////////////////////////// -/// /// trla.F -- translated by f2c (version 19970219). -/// CERN PROGLIB# F112 TRLA .VERSION KERNFOR 4.15 861204 */ -/// ORIG. 18/12/74 WH */ /// -///BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML int ipiv, ia, ib, iu; double sum; @@ -1782,17 +1720,15 @@ //////////////////////////////////////////////////////////////////////////////// /// trlta.F -- translated by f2c (version 19970219). +/// /// CERN PROGLIB# F112 TRLTA .VERSION KERNFOR 4.15 861204 /// ORIG. 18/12/74 WH -///BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML - int ipiv, mxpn, i__, nTep, ia, ib, iu, mx; double sum; @@ -1834,16 +1770,15 @@ //////////////////////////////////////////////////////////////////////////////// /// trpck.F -- translated by f2c (version 19970219). -/// CERN PROGLIB# F112 TRPCK .VERSION KERNFOR 2.08 741218 */ -/// ORIG. 18/12/74 WH */ -///BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML +/// +/// CERN PROGLIB# F112 TRPCK .VERSION KERNFOR 2.08 741218 +/// ORIG. 18/12/74 WH +/// +/// See original documentation of CERNLIB package +/// [F112](https://cern-tex.web.cern.ch/cern-tex/shortwrupsdir/f112/top.html) + +double *TCL::trpck(const double *s, double *u, int n) +{ int i__, ia, ind, ipiv; /* Parameter adjuTments */ @@ -1869,17 +1804,15 @@ //////////////////////////////////////////////////////////////////////////////// /// trqsq.F -- translated by f2c (version 19970219). -/// CERN PROGLIB# F112 TRQSQ .VERSION KERNFOR 4.15 861204 */ -/// ORIG. 18/12/74 WH */ -///BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML - int indq, inds, imax, i__, j, k, l; int iq, ir, is, iqq; double sum; @@ -1936,16 +1869,15 @@ //////////////////////////////////////////////////////////////////////////////// /// trsa.F -- translated by f2c (version 19970219). -/// CERN PROGLIB# F112 TRSA .VERSION KERNFOR 4.15 861204 */ -/// ORIG. 18/12/74 WH */ -///BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML /* Local variables */ int inds, i__, j, k, ia, ib, is; double sum; @@ -1985,17 +1917,15 @@ //////////////////////////////////////////////////////////////////////////////// /// trsinv.F -- translated by f2c (version 19970219). +/// /// CERN PROGLIB# F112 TRSINV .VERSION KERNFOR 2.08 741218 -/// ORIG. 18/12/74 WH */ -///BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML - /* Function Body */ trchlu(g, gi, n); trinv(gi, gi, n); @@ -2006,17 +1936,15 @@ //////////////////////////////////////////////////////////////////////////////// /// trsmlu.F -- translated by f2c (version 19970219). -/// CERN PROGLIB# F112 TRSMLU .VERSION KERNFOR 4.15 861204 */ -/// ORIG. 18/12/74 WH */ -///BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML +/// +/// CERN PROGLIB# F112 TRSMLU .VERSION KERNFOR 4.15 861204 +/// ORIG. 18/12/74 WH +/// +/// See original documentation of CERNLIB package +/// [F112](https://cern-tex.web.cern.ch/cern-tex/shortwrupsdir/f112/top.html) +double *TCL::trsmlu(const double *u, double *s, int n) +{ /* Local variables */ int lhor, lver, i__, k, l, ind; double sum; @@ -2043,17 +1971,15 @@ //////////////////////////////////////////////////////////////////////////////// /// trsmul.F -- translated by f2c (version 19970219). -/// CERN PROGLIB# F112 TRSMUL .VERSION KERNFOR 4.15 861204 */ -/// ORIG. 18/12/74 WH */ -///BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML +/// +/// CERN PROGLIB# F112 TRSMUL .VERSION KERNFOR 4.15 861204 +/// ORIG. 18/12/74 WH +/// +/// See original documentation of CERNLIB package +/// [F112](https://cern-tex.web.cern.ch/cern-tex/shortwrupsdir/f112/top.html) +double *TCL::trsmul(const double *g, double *gi, int n) +{ /* Local variables */ int lhor, lver, lpiv, i__, j, k, ind; double sum; @@ -2081,18 +2007,15 @@ //////////////////////////////////////////////////////////////////////////////// /// trupck.F -- translated by f2c (version 19970219). +/// /// CERN PROGLIB# F112 TRUPCK .VERSION KERNFOR 2.08 741218 /// ORIG. 18/12/74 WH -///BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML - - int i__, im, is, iu, iv, ih, m2; /* Parameter adjuTments */ @@ -2133,17 +2056,15 @@ //////////////////////////////////////////////////////////////////////////////// /// trsat.F -- translated by f2c (version 19970219) +/// /// CERN PROGLIB# F112 TRSAT .VERSION KERNFOR 4.15 861204 /// ORIG. 18/12/74 WH -///BEGIN_HTML - see original documentation of CERNLIB package F112 - END_HTML - /* Local variables */ int inds, i__, j, k, ia, ib, is; double sum;