mathlib-slatec-extended.c
1 //! @file mathlib-slatec-extended.c
2 //
3 //! @section copyright
4 //
5 // This file is part of VIF - vintage fortran compiler.
6 // Copyright 2020-2025 J. Marcel van der Veer <algol68g@xs4all.nl>.
7 //
8 //! @section license
9 //
10 // This program is free software; you can redistribute it and/or modify it
11 // under the terms of the gnu general public license as published by the
12 // free software foundation; either version 3 of the license, or
13 // (at your option) any later version.
14 //
15 // This program is distributed in the hope that it will be useful, but
16 // without any warranty; without even the implied warranty of merchantability
17 // or fitness for a particular purpose. See the GNU general public license for
18 // more details. you should have received a copy of the GNU general public
19 // license along with this program. If not, see <http://www.gnu.org/licenses/>.
20 //
21 //! @Synopsis
22 //
23 //! 22 subprograms from SLATEC-EXTENDED
24
25 // Compiled from Fortran source code by VIF.
26 // Selected subprograms are:
27 //
28 // DXADD DXNRMP DXPNRM DXQMU DXSET XCON XNRMP XPMUP XPQNU XQNU
29 // DXCON DXPMU DXPQNU DXQNU XADD XLEGF XPMU XPNRM XQMU XSET
30 // DXLEGF DXPMUP
31
32 // The license for SLATEC Fortran source code is:
33 //
34 // The SLATEC Common Mathematical Library was developed at
35 // US government research laboratories and is in the public domain.
36 //
37 // Repository: http://www.netlib.org/slatec/
38 //
39 // The SLATEC common mathematical library is issued by the following
40 //
41 // Air Force Weapons Laboratory, Albuquerque
42 // Lawrence Livermore National Laboratory, Livermore
43 // Los Alamos National Laboratory, Los Alamos
44 // National Institute of Standards and Technology, Washington
45 // National Energy Research Supercomputer Center, Livermore
46 // Oak Ridge National Laboratory, Oak Ridge
47 // Sandia National Laboratories, Albuquerque
48 // Sandia National Laboratories, Livermore
49 //
50 // All questions concerning the distribution of the library should be
51 // directed to the NATIONAL ENERGY SOFTWARE CENTER, 9700 Cass Ave.,
52 // Argonne, Illinois 60439, and not to the authors of the subprograms.
53 //
54 // * * * * * Notice * * * * *
55 //
56 // This material was prepared as an account of work sponsored by the
57 // United States Government. Neither the United States, nor the
58 // Department of Energy, nor the Department of Defense, nor any of
59 // their employees, nor any of their contractors, subcontractors, or
60 // their employees, makes any warranty, expressed or implied, or
61 // assumes any legal liability or responsibility for the accuracy,
62 // completeness, or usefulness of any information, apparatus, product,
63 // or process disclosed, or represents that its use would not infringe
64 // upon privately owned rights.
65
66 /*
67 Generated by VIF - experimental VIntage Fortran compiler.
68 VIF release 1.1.7
69 */
97 static CALLS __calls[__ncalls] = {
98 {"dxadd", 0}, // subroutine
99 {"dxcon", 0}, // subroutine
100 {"dxlegf", 0}, // subroutine
101 {"dxnrmp", 0}, // subroutine
102 {"dxpmu", 0}, // subroutine
103 {"dxpmup", 0}, // subroutine
104 {"dxpnrm", 0}, // subroutine
105 {"dxpqnu", 0}, // subroutine
106 {"dxqmu", 0}, // subroutine
107 {"dxqnu", 0}, // subroutine
108 {"dxset", 0}, // subroutine
109 {"xadd", 0}, // subroutine
110 {"xcon", 0}, // subroutine
111 {"xlegf", 0}, // subroutine
112 {"xnrmp", 0}, // subroutine
113 {"xpmu", 0}, // subroutine
114 {"xpmup", 0}, // subroutine
115 {"xpnrm", 0}, // subroutine
116 {"xpqnu", 0}, // subroutine
117 {"xqmu", 0}, // subroutine
118 {"xqnu", 0}, // subroutine
119 {"xset", 0}, // subroutine
120 {NULL, 0}
121 };
122 static struct {
123 int_4 l_, l2_, kmax_;
124 real_8 radix_, radixl_, rad2l_, dlg10r_;
125 } dxblk2_;
126 static struct {
127 int_4 nbitsf_;
128 } dxblk1_;
129 static struct {
130 int_4 nlg102_, mlg102_, lg102_[21];
131 } dxblk3_;
132 static struct {
133 int_4 l_, l2_, kmax_;
134 real_4 radix_, radixl_, rad2l_, dlg10r_;
135 } xblk2_;
136 static struct {
137 int_4 nbitsf_;
138 } xblk1_;
139 static struct {
140 int_4 nlg102_, mlg102_, lg102_[21];
141 } xblk3_;
143 _p_ ierror_);
146 _p_ id_, real_8 _p_ pqa_, int_4 _p_ ipqa_, int_4 _p_ ierror_);
148 dpn_, int_4 _p_ ipn_, int_4 _p_ isig_, int_4 _p_ ierror_);
150 x_, real_8 _p_ sx_, int_4 _p_ id_, real_8 _p_ pqa_, int_4 _p_ ipqa_, int_4 _p_ ierror_);
152 ipqa_, int_4 _p_ ierror_);
154 ipqa_, int_4 _p_ ierror_);
156 pqa_, int_4 _p_ ipqa_, int_4 _p_ ierror_);
158 x_, real_8 _p_ sx_, int_4 _p_ id_, real_8 _p_ pqa_, int_4 _p_ ipqa_, int_4 _p_ ierror_);
160 sx_, int_4 _p_ id_, real_8 _p_ pqa_, int_4 _p_ ipqa_, int_4 _p_ ierror_);
163 _p_ ierror_);
166 _p_ id_, real_4 _p_ pqa_, int_4 _p_ ipqa_, int_4 _p_ ierror_);
168 spn_, int_4 _p_ ipn_, int_4 _p_ isig_, int_4 _p_ ierror_);
170 x_, real_4 _p_ sx_, int_4 _p_ id_, real_4 _p_ pqa_, int_4 _p_ ipqa_, int_4 _p_ ierror_);
172 ipqa_, int_4 _p_ ierror_);
174 ipqa_, int_4 _p_ ierror_);
176 pqa_, int_4 _p_ ipqa_, int_4 _p_ ierror_);
178 x_, real_4 _p_ sx_, int_4 _p_ id_, real_4 _p_ pqa_, int_4 _p_ ipqa_, int_4 _p_ ierror_);
180 sx_, int_4 _p_ id_, real_4 _p_ pqa_, int_4 _p_ ipqa_, int_4 _p_ ierror_);
182 int_4 _dxadd (real_8 _p_ x_, int_4 _p_ ix_, real_8 _p_ y_, int_4 _p_ iy_, real_8 _p_ z_, int_4 _p_ iz_, int_4 _p_
183 ierror_)
184 { // ** body not listed **
359 }
360 int_4 _dxcon (real_8 _p_ x_, int_4 _p_ ix_, int_4 _p_ ierror_)
361 { // ** body not listed **
514 }
515 int_4 _dxlegf (real_8 _p_ dnu1_, int_4 _p_ nudiff_, int_4 _p_ mu1_, int_4 _p_ mu2_, real_8 _p_ theta_, int_4 _p_ id_,
516 real_8 _p_ pqa_, int_4 _p_ ipqa_, int_4 _p_ ierror_)
517 { // ** body not listed **
629 }
630 int_4 _dxnrmp (int_4 _p_ nu_, int_4 _p_ mu1_, int_4 _p_ mu2_, real_8 _p_ darg_, int_4 _p_ mode_, real_8 _p_ dpn_, int_4
631 _p_ ipn_, int_4 _p_ isig_, int_4 _p_ ierror_)
632 { // ** body not listed **
806 }
807 int_4 _dxpmu (real_8 _p_ nu1_, real_8 _p_ nu2_, int_4 _p_ mu1_, int_4 _p_ mu2_, real_8 _p_ theta_, real_8 _p_ x_,
808 real_8 _p_ sx_, int_4 _p_ id_, real_8 _p_ pqa_, int_4 _p_ ipqa_, int_4 _p_ ierror_)
809 { // ** body not listed **
857 }
858 int_4 _dxpmup (real_8 _p_ nu1_, real_8 _p_ nu2_, int_4 _p_ mu1_, int_4 _p_ mu2_, real_8 _p_ pqa_, int_4 _p_ ipqa_,
859 int_4 _p_ ierror_)
860 { // ** body not listed **
938 }
939 int_4 _dxpnrm (real_8 _p_ nu1_, real_8 _p_ nu2_, int_4 _p_ mu1_, int_4 _p_ mu2_, real_8 _p_ pqa_, int_4 _p_ ipqa_,
940 int_4 _p_ ierror_)
941 { // ** body not listed **
1028 }
1029 int_4 _dxpqnu (real_8 _p_ nu1_, real_8 _p_ nu2_, int_4 _p_ mu_, real_8 _p_ theta_, int_4 _p_ id_, real_8 _p_ pqa_,
1030 int_4 _p_ ipqa_, int_4 _p_ ierror_)
1031 { // ** body not listed **
1222 }
1223 int_4 _dxqmu (real_8 _p_ nu1_, real_8 _p_ nu2_, int_4 _p_ mu1_, int_4 _p_ mu2_, real_8 _p_ theta_, real_8 _p_ x_,
1224 real_8 _p_ sx_, int_4 _p_ id_, real_8 _p_ pqa_, int_4 _p_ ipqa_, int_4 _p_ ierror_)
1225 { // ** body not listed **
1297 }
1298 int_4 _dxqnu (real_8 _p_ nu1_, real_8 _p_ nu2_, int_4 _p_ mu1_, real_8 _p_ theta_, real_8 _p_ x_, real_8 _p_ sx_, int_4
1299 _p_ id_, real_8 _p_ pqa_, int_4 _p_ ipqa_, int_4 _p_ ierror_)
1300 { // ** body not listed **
1409 }
1410 int_4 _dxset (int_4 _p_ irad_, int_4 _p_ nradpl_, real_8 _p_ dzero_, int_4 _p_ nbits_, int_4 _p_ ierror_)
1411 { // ** body not listed **
1596 }
1597 int_4 _xadd (real_4 _p_ x_, int_4 _p_ ix_, real_4 _p_ y_, int_4 _p_ iy_, real_4 _p_ z_, int_4 _p_ iz_, int_4 _p_
1598 ierror_)
1599 { // ** body not listed **
1774 }
1775 int_4 _xcon (real_4 _p_ x_, int_4 _p_ ix_, int_4 _p_ ierror_)
1776 { // ** body not listed **
1929 }
1930 int_4 _xlegf (real_4 _p_ dnu1_, int_4 _p_ nudiff_, int_4 _p_ mu1_, int_4 _p_ mu2_, real_4 _p_ theta_, int_4 _p_ id_,
1931 real_4 _p_ pqa_, int_4 _p_ ipqa_, int_4 _p_ ierror_)
1932 { // ** body not listed **
2044 }
2045 int_4 _xnrmp (int_4 _p_ nu_, int_4 _p_ mu1_, int_4 _p_ mu2_, real_4 _p_ sarg_, int_4 _p_ mode_, real_4 _p_ spn_, int_4
2046 _p_ ipn_, int_4 _p_ isig_, int_4 _p_ ierror_)
2047 { // ** body not listed **
2221 }
2222 int_4 _xpmu (real_4 _p_ nu1_, real_4 _p_ nu2_, int_4 _p_ mu1_, int_4 _p_ mu2_, real_4 _p_ theta_, real_4 _p_ x_, real_4
2223 _p_ sx_, int_4 _p_ id_, real_4 _p_ pqa_, int_4 _p_ ipqa_, int_4 _p_ ierror_)
2224 { // ** body not listed **
2272 }
2273 int_4 _xpmup (real_4 _p_ nu1_, real_4 _p_ nu2_, int_4 _p_ mu1_, int_4 _p_ mu2_, real_4 _p_ pqa_, int_4 _p_ ipqa_, int_4
2274 _p_ ierror_)
2275 { // ** body not listed **
2353 }
2354 int_4 _xpnrm (real_4 _p_ nu1_, real_4 _p_ nu2_, int_4 _p_ mu1_, int_4 _p_ mu2_, real_4 _p_ pqa_, int_4 _p_ ipqa_, int_4
2355 _p_ ierror_)
2356 { // ** body not listed **
2443 }
2444 int_4 _xpqnu (real_4 _p_ nu1_, real_4 _p_ nu2_, int_4 _p_ mu_, real_4 _p_ theta_, int_4 _p_ id_, real_4 _p_ pqa_, int_4
2445 _p_ ipqa_, int_4 _p_ ierror_)
2446 { // ** body not listed **
2637 }
2638 int_4 _xqmu (real_4 _p_ nu1_, real_4 _p_ nu2_, int_4 _p_ mu1_, int_4 _p_ mu2_, real_4 _p_ theta_, real_4 _p_ x_, real_4
2639 _p_ sx_, int_4 _p_ id_, real_4 _p_ pqa_, int_4 _p_ ipqa_, int_4 _p_ ierror_)
2640 { // ** body not listed **
2712 }
2713 int_4 _xqnu (real_4 _p_ nu1_, real_4 _p_ nu2_, int_4 _p_ mu1_, real_4 _p_ theta_, real_4 _p_ x_, real_4 _p_ sx_, int_4
2714 _p_ id_, real_4 _p_ pqa_, int_4 _p_ ipqa_, int_4 _p_ ierror_)
2715 { // ** body not listed **
2824 }
2825 int_4 _xset (int_4 _p_ irad_, int_4 _p_ nradpl_, real_4 _p_ dzero_, int_4 _p_ nbits_, int_4 _p_ ierror_)
2826 { // ** body not listed **
3011 }
|
© 2002-2025 J.M. van der Veer (jmvdveer@xs4all.nl)
|
