1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95
|
/* Ergo, version 3.8.2, a program for linear scaling electronic structure
* calculations.
* Copyright (C) 2023 Elias Rudberg, Emanuel H. Rubensson, Pawel Salek,
* and Anastasia Kruchinina.
*
* This program is free software: you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
* the Free Software Foundation, either version 3 of the License, or
* (at your option) any later version.
*
* This program is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with this program. If not, see <http://www.gnu.org/licenses/>.
*
* Primary academic reference:
* Ergo: An open-source program for linear-scaling electronic structure
* calculations,
* Elias Rudberg, Emanuel H. Rubensson, Pawel Salek, and Anastasia
* Kruchinina,
* SoftwareX 7, 107 (2018),
* <http://dx.doi.org/10.1016/j.softx.2018.03.005>
*
* For further information about Ergo, see <http://www.ergoscf.org>.
*/
/* This file belongs to the template_lapack part of the Ergo source
* code. The source files in the template_lapack directory are modified
* versions of files originally distributed as CLAPACK, see the
* Copyright/license notice in the file template_lapack/COPYING.
*/
#ifndef TEMPLATE_LAPACK_LAPY3_HEADER
#define TEMPLATE_LAPACK_LAPY3_HEADER
template<class Treal>
Treal template_lapack_lapy3(Treal *x, Treal *y, Treal *z__)
{
/* -- LAPACK auxiliary routine (version 3.0) --
Univ. of Tennessee, Univ. of California Berkeley, NAG Ltd.,
Courant Institute, Argonne National Lab, and Rice University
October 31, 1992
Purpose
=======
DLAPY3 returns sqrt(x**2+y**2+z**2), taking care not to cause
unnecessary overflow.
Arguments
=========
X (input) DOUBLE PRECISION
Y (input) DOUBLE PRECISION
Z (input) DOUBLE PRECISION
X, Y and Z specify the values x, y and z.
===================================================================== */
/* System generated locals */
Treal ret_val, d__1, d__2, d__3;
/* Local variables */
Treal xabs, yabs, zabs, w;
xabs = absMACRO(*x);
yabs = absMACRO(*y);
zabs = absMACRO(*z__);
/* Computing MAX */
d__1 = maxMACRO(xabs,yabs);
w = maxMACRO(d__1,zabs);
if (w == 0.) {
ret_val = 0.;
} else {
/* Computing 2nd power */
d__1 = xabs / w;
/* Computing 2nd power */
d__2 = yabs / w;
/* Computing 2nd power */
d__3 = zabs / w;
ret_val = w * template_blas_sqrt(d__1 * d__1 + d__2 * d__2 + d__3 * d__3);
}
return ret_val;
/* End of DLAPY3 */
} /* dlapy3_ */
#endif
|