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.2.16
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
123 static struct {
124 int_4 l_, l2_, kmax_;
125 real_8 radix_, radixl_, rad2l_, dlg10r_;
126 } dxblk2_;
127 static struct {
128 int_4 nbitsf_;
129 } dxblk1_;
130 static struct {
131 int_4 nlg102_, mlg102_, lg102_[21];
132 } dxblk3_;
133 static struct {
134 int_4 l_, l2_, kmax_;
135 real_4 radix_, radixl_, rad2l_, dlg10r_;
136 } xblk2_;
137 static struct {
138 int_4 nbitsf_;
139 } xblk1_;
140 static struct {
141 int_4 nlg102_, mlg102_, lg102_[21];
142 } xblk3_;
144 _p_ ierror_);
147 _p_ id_, real_8 _p_ pqa_, int_4 _p_ ipqa_, int_4 _p_ ierror_);
149 dpn_, int_4 _p_ ipn_, int_4 _p_ isig_, int_4 _p_ ierror_);
151 x_, real_8 _p_ sx_, int_4 _p_ id_, real_8 _p_ pqa_, int_4 _p_ ipqa_, int_4 _p_ ierror_);
153 ipqa_, int_4 _p_ ierror_);
155 ipqa_, int_4 _p_ ierror_);
157 pqa_, int_4 _p_ ipqa_, int_4 _p_ ierror_);
159 x_, real_8 _p_ sx_, int_4 _p_ id_, real_8 _p_ pqa_, int_4 _p_ ipqa_, int_4 _p_ ierror_);
161 sx_, int_4 _p_ id_, real_8 _p_ pqa_, int_4 _p_ ipqa_, int_4 _p_ ierror_);
164 _p_ ierror_);
167 _p_ id_, real_4 _p_ pqa_, int_4 _p_ ipqa_, int_4 _p_ ierror_);
169 spn_, int_4 _p_ ipn_, int_4 _p_ isig_, int_4 _p_ ierror_);
171 x_, real_4 _p_ sx_, int_4 _p_ id_, real_4 _p_ pqa_, int_4 _p_ ipqa_, int_4 _p_ ierror_);
173 ipqa_, int_4 _p_ ierror_);
175 ipqa_, int_4 _p_ ierror_);
177 pqa_, int_4 _p_ ipqa_, int_4 _p_ ierror_);
179 x_, real_4 _p_ sx_, int_4 _p_ id_, real_4 _p_ pqa_, int_4 _p_ ipqa_, int_4 _p_ ierror_);
181 sx_, int_4 _p_ id_, real_4 _p_ pqa_, int_4 _p_ ipqa_, int_4 _p_ ierror_);
183 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_
184 ierror_)
185 { // ** body not listed **
360 }
361
362 int_4 _dxcon (real_8 _p_ x_, int_4 _p_ ix_, int_4 _p_ ierror_)
363 { // ** body not listed **
516 }
517
518 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_,
519 real_8 _p_ pqa_, int_4 _p_ ipqa_, int_4 _p_ ierror_)
520 { // ** body not listed **
632 }
633
634 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
635 _p_ ipn_, int_4 _p_ isig_, int_4 _p_ ierror_)
636 { // ** body not listed **
810 }
811
812 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_,
813 real_8 _p_ sx_, int_4 _p_ id_, real_8 _p_ pqa_, int_4 _p_ ipqa_, int_4 _p_ ierror_)
814 { // ** body not listed **
862 }
863
864 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_,
865 int_4 _p_ ierror_)
866 { // ** body not listed **
944 }
945
946 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_,
947 int_4 _p_ ierror_)
948 { // ** body not listed **
1035 }
1036
1037 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_,
1038 int_4 _p_ ipqa_, int_4 _p_ ierror_)
1039 { // ** body not listed **
1230 }
1231
1232 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_,
1233 real_8 _p_ sx_, int_4 _p_ id_, real_8 _p_ pqa_, int_4 _p_ ipqa_, int_4 _p_ ierror_)
1234 { // ** body not listed **
1306 }
1307
1308 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
1309 _p_ id_, real_8 _p_ pqa_, int_4 _p_ ipqa_, int_4 _p_ ierror_)
1310 { // ** body not listed **
1419 }
1420
1421 int_4 _dxset (int_4 _p_ irad_, int_4 _p_ nradpl_, real_8 _p_ dzero_, int_4 _p_ nbits_, int_4 _p_ ierror_)
1422 { // ** body not listed **
1607 }
1608
1609 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_
1610 ierror_)
1611 { // ** body not listed **
1786 }
1787
1788 int_4 _xcon (real_4 _p_ x_, int_4 _p_ ix_, int_4 _p_ ierror_)
1789 { // ** body not listed **
1942 }
1943
1944 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_,
1945 real_4 _p_ pqa_, int_4 _p_ ipqa_, int_4 _p_ ierror_)
1946 { // ** body not listed **
2058 }
2059
2060 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
2061 _p_ ipn_, int_4 _p_ isig_, int_4 _p_ ierror_)
2062 { // ** body not listed **
2236 }
2237
2238 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
2239 _p_ sx_, int_4 _p_ id_, real_4 _p_ pqa_, int_4 _p_ ipqa_, int_4 _p_ ierror_)
2240 { // ** body not listed **
2288 }
2289
2290 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
2291 _p_ ierror_)
2292 { // ** body not listed **
2370 }
2371
2372 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
2373 _p_ ierror_)
2374 { // ** body not listed **
2461 }
2462
2463 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
2464 _p_ ipqa_, int_4 _p_ ierror_)
2465 { // ** body not listed **
2656 }
2657
2658 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
2659 _p_ sx_, int_4 _p_ id_, real_4 _p_ pqa_, int_4 _p_ ipqa_, int_4 _p_ ierror_)
2660 { // ** body not listed **
2732 }
2733
2734 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
2735 _p_ id_, real_4 _p_ pqa_, int_4 _p_ ipqa_, int_4 _p_ ierror_)
2736 { // ** body not listed **
2845 }
2846
2847 int_4 _xset (int_4 _p_ irad_, int_4 _p_ nradpl_, real_4 _p_ dzero_, int_4 _p_ nbits_, int_4 _p_ ierror_)
2848 { // ** body not listed **
3033 }
3034
|
© 2002-2025 J.M. van der Veer (jmvdveer@xs4all.nl)
|
