mathlib-lapack-svd.c

     1  //! @file mathlib-lapack-svd.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  //! 85 subprograms from LAPACK-SVD
    24  
    25  // Compiled from Fortran source code by VIF.
    26  // Selected subprograms are:
    27  //
    28  // DBDSDC DGELSS DGSVJ1 DLALSA DLARF  DLASD2 DLASDQ DLASQ6 DORGBR DORMQR
    29  // DBDSQR DGEQP3 DISNAN DLALSD DLARFB DLASD3 DLASDT DLASR  DORGL2 DPOCON
    30  // DCOMBS DGEQR2 DLABAD DLAMRG DLARFG DLASD4 DLASET DLASRT DORGLQ DRSCL
    31  // DGEBD2 DGEQRF DLABRD DLANGE DLARFT DLASD5 DLASQ1 DLASSQ DORGQR IEEECK
    32  // DGEBRD DGESDD DLACN2 DLANST DLARTG DLASD6 DLASQ2 DLASV2 DORM2R ILADLC
    33  // DGEJSV DGESVD DLACPY DLAPMT DLAS2  DLASD7 DLASQ3 DLASWP DORMBR ILADLR
    34  // DGELQ2 DGESVD DLAED6 DLAPY2 DLASCL DLASD8 DLASQ4 DLATRS DORML2 ILAENV
    35  // DGELQF DGESVJ DLAISN DLAQP2 DLASD0 DLASDA DLASQ5 DORG2R DORMLQ IPARMQ
    36  // DGELSD DGSVJ0 DLALS0 DLAQPS DLASD1
    37  
    38  // Selected SVD functionality from LAPACK. 
    39  // 
    40  // The original code is F90 from 2023, but could be made to compile 
    41  // with VIF trivially.
    42  // 
    43  // Included LAPACK routines of interest to VIF are:
    44  // 
    45  //   DGEJSV
    46  //   DGELSD 
    47  //   DGELSS 
    48  //   DGESDD
    49  //   DGESVD
    50  //   DGESVDQ
    51  //   DGESVJ
    52  // 
    53  // Source: netlib.org/lapack
    54  // LAPACK 3.12.0
    55  
    56  // The license for LAPACK Fortran source code is:
    57  //
    58  // Copyright (c) 1992-2013 The University of Tennessee and The University
    59  //                         of Tennessee Research Foundation.  All rights
    60  //                         reserved.
    61  // Copyright (c) 2000-2013 The University of California Berkeley. All
    62  //                         rights reserved.
    63  // Copyright (c) 2006-2013 The University of Colorado Denver.  All rights
    64  //                         reserved.
    65  // 
    66  // $COPYRIGHT$
    67  // 
    68  // Additional copyrights may follow
    69  // 
    70  // $HEADER$
    71  // 
    72  // Redistribution and use in source and binary forms, with or without
    73  // modification, are permitted provided that the following conditions are
    74  // met:
    75  // 
    76  // - Redistributions of source code must retain the above copyright
    77  //   notice, this list of conditions and the following disclaimer.
    78  // 
    79  // - Redistributions in binary form must reproduce the above copyright
    80  //   notice, this list of conditions and the following disclaimer listed
    81  //   in this license in the documentation and/or other materials
    82  //   provided with the distribution.
    83  // 
    84  // - Neither the name of the copyright holders nor the names of its
    85  //   contributors may be used to endorse or promote products derived from
    86  //   this software without specific prior written permission.
    87  // 
    88  // The copyright holders provide no reassurances that the source code
    89  // provided does not infringe any patent, copyright, or any other
    90  // intellectual property rights of third parties.  The copyright holders
    91  // disclaim any liability to any recipient for claims brought against
    92  // recipient by any third party for infringement of that parties
    93  // intellectual property rights.
    94  // 
    95  // THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
    96  // "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
    97  // LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
    98  // A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT
    99  // OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
   100  // SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
   101  // LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
   102  // DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
   103  // THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
   104  // (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
   105  // OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
   106  // `
   107  
   108  /*
   109  Generated by VIF - experimental VIntage Fortran compiler.
   110  VIF release 1.2.12
   111  */
   264  static CALLS __calls[__ncalls] = {
   265    {"dbdsdc", 0}, // subroutine
   266    {"dbdsqr", 0}, // subroutine
   267    {"dcombssq", 0}, // subroutine
   268    {"dgebd2", 0}, // subroutine
   269    {"dgebrd", 0}, // subroutine
   270    {"dgejsv", 0}, // subroutine
   271    {"dgelq2", 0}, // subroutine
   272    {"dgelqf", 0}, // subroutine
   273    {"dgelsd", 0}, // subroutine
   274    {"dgelss", 0}, // subroutine
   275    {"dgeqp3", 0}, // subroutine
   276    {"dgeqr2", 0}, // subroutine
   277    {"dgeqrf", 0}, // subroutine
   278    {"dgesdd", 0}, // subroutine
   279    {"dgesvd", 0}, // subroutine
   280    {"dgesvdq", 0}, // subroutine
   281    {"dgesvj", 0}, // subroutine
   282    {"dgsvj0", 0}, // subroutine
   283    {"dgsvj1", 0}, // subroutine
   284    {"disnan", 0}, // logical*4 function
   285    {"dlabad", 0}, // subroutine
   286    {"dlabrd", 0}, // subroutine
   287    {"dlacn2", 0}, // subroutine
   288    {"dlacpy", 0}, // subroutine
   289    {"dlaed6", 0}, // subroutine
   290    {"dlaisnan", 0}, // logical*4 function
   291    {"dlals0", 0}, // subroutine
   292    {"dlalsa", 0}, // subroutine
   293    {"dlalsd", 0}, // subroutine
   294    {"dlamrg", 0}, // subroutine
   295    {"dlange", 0}, // real*8 function
   296    {"dlanst", 0}, // real*8 function
   297    {"dlapmt", 0}, // subroutine
   298    {"dlapy2", 0}, // real*8 function
   299    {"dlaqp2", 0}, // subroutine
   300    {"dlaqps", 0}, // subroutine
   301    {"dlarfb", 0}, // subroutine
   302    {"dlarf", 0}, // subroutine
   303    {"dlarfg", 0}, // subroutine
   304    {"dlarft", 0}, // subroutine
   305    {"dlartg", 0}, // subroutine
   306    {"dlas2", 0}, // subroutine
   307    {"dlascl", 0}, // subroutine
   308    {"dlasd0", 0}, // subroutine
   309    {"dlasd1", 0}, // subroutine
   310    {"dlasd2", 0}, // subroutine
   311    {"dlasd3", 0}, // subroutine
   312    {"dlasd4", 0}, // subroutine
   313    {"dlasd5", 0}, // subroutine
   314    {"dlasd6", 0}, // subroutine
   315    {"dlasd7", 0}, // subroutine
   316    {"dlasd8", 0}, // subroutine
   317    {"dlasda", 0}, // subroutine
   318    {"dlasdq", 0}, // subroutine
   319    {"dlasdt", 0}, // subroutine
   320    {"dlaset", 0}, // subroutine
   321    {"dlasq1", 0}, // subroutine
   322    {"dlasq2", 0}, // subroutine
   323    {"dlasq3", 0}, // subroutine
   324    {"dlasq4", 0}, // subroutine
   325    {"dlasq5", 0}, // subroutine
   326    {"dlasq6", 0}, // subroutine
   327    {"dlasr", 0}, // subroutine
   328    {"dlasrt", 0}, // subroutine
   329    {"dlassq", 0}, // subroutine
   330    {"dlasv2", 0}, // subroutine
   331    {"dlaswp", 0}, // subroutine
   332    {"dlatrs", 0}, // subroutine
   333    {"dorg2r", 0}, // subroutine
   334    {"dorgbr", 0}, // subroutine
   335    {"dorgl2", 0}, // subroutine
   336    {"dorglq", 0}, // subroutine
   337    {"dorgqr", 0}, // subroutine
   338    {"dorm2r", 0}, // subroutine
   339    {"dormbr", 0}, // subroutine
   340    {"dorml2", 0}, // subroutine
   341    {"dormlq", 0}, // subroutine
   342    {"dormqr", 0}, // subroutine
   343    {"dpocon", 0}, // subroutine
   344    {"drscl", 0}, // subroutine
   345    {"ieeeck", 0}, // integer*4 function
   346    {"iladlc", 0}, // integer*4 function
   347    {"iladlr", 0}, // integer*4 function
   348    {"ilaenv", 0}, // integer*4 function
   349    {"iparmq", 0}, // integer*4 function
   350    {NULL, 0}
   351  };
   352  
   354  int_4 _p_ ldu_, real_8 _p_ vt_, int_4 _p_ ldvt_, real_8 _p_ q_, int_4 _p_ iq_, real_8 _p_ work_, int_4 _p_ iwork_, 
   355  int_4 _p_ info_);
   357  real_8 _p_ e_, real_8 _p_ vt_, int_4 _p_ ldvt_, real_8 _p_ u_, int_4 _p_ ldu_, real_8 _p_ c_, int_4 _p_ ldc_, real_8 
   358  _p_ work_, int_4 _p_ info_);
   361  real_8 _p_ tauq_, real_8 _p_ taup_, real_8 _p_ work_, int_4 _p_ info_);
   363  real_8 _p_ tauq_, real_8 _p_ taup_, real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ info_);
   365  , int_4 _p_ m_, int_4 _p_ n_, real_8 _p_ a_, int_4 _p_ lda_, real_8 _p_ sva_, real_8 _p_ u_, int_4 _p_ ldu_, real_8 _p_ 
   366  v_, int_4 _p_ ldv_, real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ iwork_, int_4 _p_ info_);
   368  int_4 _p_ info_);
   370  int_4 _p_ lwork_, int_4 _p_ info_);
   372  int_4 _p_ ldb_, real_8 _p_ s_, real_8 _p_ rcond_, int_4 _p_ rank_, real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ iwork_
   373  , int_4 _p_ info_);
   375  int_4 _p_ ldb_, real_8 _p_ s_, real_8 _p_ rcond_, int_4 _p_ rank_, real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ info_)
   376  ;
   378  real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ info_);
   380  int_4 _p_ info_);
   382  int_4 _p_ lwork_, int_4 _p_ info_);
   384  real_8 _p_ u_, int_4 _p_ ldu_, real_8 _p_ vt_, int_4 _p_ ldvt_, real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ iwork_, 
   385  int_4 _p_ info_);
   387  real_8 _p_ s_, real_8 _p_ u_, int_4 _p_ ldu_, real_8 _p_ vt_, int_4 _p_ ldvt_, real_8 _p_ work_, int_4 _p_ lwork_, 
   388  int_4 _p_ info_);
   390  int_4 _p_ n_, real_8 _p_ a_, int_4 _p_ lda_, real_8 _p_ s_, real_8 _p_ u_, int_4 _p_ ldu_, real_8 _p_ v_, int_4 _p_ 
   391  ldv_, int_4 _p_ numrank_, int_4 _p_ iwork_, int_4 _p_ liwork_, real_8 _p_ work_, int_4 _p_ lwork_, real_8 _p_ rwork_, 
   392  int_4 _p_ lrwork_, int_4 _p_ info_);
   394  int_4 _p_ lda_, real_8 _p_ sva_, int_4 _p_ mv_, real_8 _p_ v_, int_4 _p_ ldv_, real_8 _p_ work_, int_4 _p_ lwork_, 
   395  int_4 _p_ info_);
   397  real_8 _p_ sva_, int_4 _p_ mv_, real_8 _p_ v_, int_4 _p_ ldv_, real_8 _p_ eps_, real_8 _p_ sfmin_, real_8 _p_ tol_, 
   398  int_4 _p_ nsweep_, real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ info_);
   400  real_8 _p_ d_, real_8 _p_ sva_, int_4 _p_ mv_, real_8 _p_ v_, int_4 _p_ ldv_, real_8 _p_ eps_, real_8 _p_ sfmin_, 
   401  real_8 _p_ tol_, int_4 _p_ nsweep_, real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ info_);
   405  real_8 _p_ e_, real_8 _p_ tauq_, real_8 _p_ taup_, real_8 _p_ x_, int_4 _p_ ldx_, real_8 _p_ y_, int_4 _p_ ldy_);
   407  int_4 _p_ isave_);
   409  int_4 _p_ ldb_);
   411  real_8 _p_ finit_, real_8 _p_ tau_, int_4 _p_ info_);
   414  b_, int_4 _p_ ldb_, real_8 _p_ bx_, int_4 _p_ ldbx_, int_4 _p_ perm_, int_4 _p_ givptr_, int_4 _p_ givcol_, int_4 _p_ 
   415  ldgcol_, real_8 _p_ givnum_, int_4 _p_ ldgnum_, real_8 _p_ poles_, real_8 _p_ difl_, real_8 _p_ difr_, real_8 _p_ z_, 
   416  int_4 _p_ k_, real_8 _p_ c_, real_8 _p_ s_, real_8 _p_ work_, int_4 _p_ info_);
   418  ldb_, real_8 _p_ bx_, int_4 _p_ ldbx_, real_8 _p_ u_, int_4 _p_ ldu_, real_8 _p_ vt_, int_4 _p_ k_, real_8 _p_ difl_, 
   419  real_8 _p_ difr_, real_8 _p_ z_, real_8 _p_ poles_, int_4 _p_ givptr_, int_4 _p_ givcol_, int_4 _p_ ldgcol_, int_4 _p_ 
   420  perm_, real_8 _p_ givnum_, real_8 _p_ c_, real_8 _p_ s_, real_8 _p_ work_, int_4 _p_ iwork_, int_4 _p_ info_);
   422  , real_8 _p_ b_, int_4 _p_ ldb_, real_8 _p_ rcond_, int_4 _p_ rank_, real_8 _p_ work_, int_4 _p_ iwork_, int_4 _p_ 
   423  info_);
   425  index_);
   429  );
   432  real_8 _p_ tau_, real_8 _p_ vn1_, real_8 _p_ vn2_, real_8 _p_ work_);
   434  int_4 _p_ lda_, int_4 _p_ jpvt_, real_8 _p_ tau_, real_8 _p_ vn1_, real_8 _p_ vn2_, real_8 _p_ auxv_, real_8 _p_ f_, 
   435  int_4 _p_ ldf_);
   437  n_, int_4 _p_ k_, real_8 _p_ v_, int_4 _p_ ldv_, real_8 _p_ t_, int_4 _p_ ldt_, real_8 _p_ c_, int_4 _p_ ldc_, real_8 
   438  _p_ work_, int_4 _p_ ldwork_);
   440  real_8 _p_ c_, int_4 _p_ ldc_, real_8 _p_ work_);
   443  real_8 _p_ tau_, real_8 _p_ t_, int_4 _p_ ldt_);
   447  , int_4 _p_ n_, real_8 _p_ a_, int_4 _p_ lda_, int_4 _p_ info_);
   449  real_8 _p_ vt_, int_4 _p_ ldvt_, int_4 _p_ smlsiz_, int_4 _p_ iwork_, real_8 _p_ work_, int_4 _p_ info_);
   451  beta_, real_8 _p_ u_, int_4 _p_ ldu_, real_8 _p_ vt_, int_4 _p_ ldvt_, int_4 _p_ idxq_, int_4 _p_ iwork_, real_8 _p_ 
   452  work_, int_4 _p_ info_);
   454  real_8 _p_ alpha_, real_8 _p_ beta_, real_8 _p_ u_, int_4 _p_ ldu_, real_8 _p_ vt_, int_4 _p_ ldvt_, real_8 _p_ dsigma_
   455  , real_8 _p_ u2_, int_4 _p_ ldu2_, real_8 _p_ vt2_, int_4 _p_ ldvt2_, int_4 _p_ idxp_, int_4 _p_ idx_, int_4 _p_ idxc_, 
   456  int_4 _p_ idxq_, int_4 _p_ coltyp_, int_4 _p_ info_);
   458  int_4 _p_ ldq_, real_8 _p_ dsigma_, real_8 _p_ u_, int_4 _p_ ldu_, real_8 _p_ u2_, int_4 _p_ ldu2_, real_8 _p_ vt_, 
   459  int_4 _p_ ldvt_, real_8 _p_ vt2_, int_4 _p_ ldvt2_, int_4 _p_ idxc_, int_4 _p_ ctot_, real_8 _p_ z_, int_4 _p_ info_);
   461  real_8 _p_ sigma_, real_8 _p_ work_, int_4 _p_ info_);
   463  dsigma_, real_8 _p_ work_);
   465  vf_, real_8 _p_ vl_, real_8 _p_ alpha_, real_8 _p_ beta_, int_4 _p_ idxq_, int_4 _p_ perm_, int_4 _p_ givptr_, int_4 
   466  _p_ givcol_, int_4 _p_ ldgcol_, real_8 _p_ givnum_, int_4 _p_ ldgnum_, real_8 _p_ poles_, real_8 _p_ difl_, real_8 _p_ 
   467  difr_, real_8 _p_ z_, int_4 _p_ k_, real_8 _p_ c_, real_8 _p_ s_, real_8 _p_ work_, int_4 _p_ iwork_, int_4 _p_ info_);
   469  real_8 _p_ z_, real_8 _p_ zw_, real_8 _p_ vf_, real_8 _p_ vfw_, real_8 _p_ vl_, real_8 _p_ vlw_, real_8 _p_ alpha_, 
   470  real_8 _p_ beta_, real_8 _p_ dsigma_, int_4 _p_ idx_, int_4 _p_ idxp_, int_4 _p_ idxq_, int_4 _p_ perm_, int_4 _p_ 
   471  givptr_, int_4 _p_ givcol_, int_4 _p_ ldgcol_, real_8 _p_ givnum_, int_4 _p_ ldgnum_, real_8 _p_ c_, real_8 _p_ s_, 
   472  int_4 _p_ info_);
   474  real_8 _p_ difl_, real_8 _p_ difr_, int_4 _p_ lddifr_, real_8 _p_ dsigma_, real_8 _p_ work_, int_4 _p_ info_);
   476  e_, real_8 _p_ u_, int_4 _p_ ldu_, real_8 _p_ vt_, int_4 _p_ k_, real_8 _p_ difl_, real_8 _p_ difr_, real_8 _p_ z_, 
   477  real_8 _p_ poles_, int_4 _p_ givptr_, int_4 _p_ givcol_, int_4 _p_ ldgcol_, int_4 _p_ perm_, real_8 _p_ givnum_, real_8 
   478  _p_ c_, real_8 _p_ s_, real_8 _p_ work_, int_4 _p_ iwork_, int_4 _p_ info_);
   480  , real_8 _p_ d_, real_8 _p_ e_, real_8 _p_ vt_, int_4 _p_ ldvt_, real_8 _p_ u_, int_4 _p_ ldu_, real_8 _p_ c_, int_4 
   481  _p_ ldc_, real_8 _p_ work_, int_4 _p_ info_);
   483  ndimr_, int_4 _p_ msub_);
   485  , int_4 _p_ lda_);
   489  sigma_, real_8 _p_ desig_, real_8 _p_ qmax_, int_4 _p_ nfail_, int_4 _p_ iter_, int_4 _p_ ndiv_, logical_4 _p_ ieee_, 
   490  int_4 _p_ ttype_, real_8 _p_ dmin1_, real_8 _p_ dmin2_, real_8 _p_ dn_, real_8 _p_ dn1_, real_8 _p_ dn2_, real_8 _p_ g_
   491  , real_8 _p_ tau_);
   493  real_8 _p_ dmin1_, real_8 _p_ dmin2_, real_8 _p_ dn_, real_8 _p_ dn1_, real_8 _p_ dn2_, real_8 _p_ tau_, int_4 _p_ 
   494  ttype_, real_8 _p_ g_);
   496  , real_8 _p_ dmin_, real_8 _p_ dmin1_, real_8 _p_ dmin2_, real_8 _p_ dn_, real_8 _p_ dnm1_, real_8 _p_ dnm2_, logical_4 
   497  _p_ ieee_, real_8 _p_ eps_);
   499  dmin1_, real_8 _p_ dmin2_, real_8 _p_ dn_, real_8 _p_ dnm1_, real_8 _p_ dnm2_);
   501  real_8 _p_ s_, real_8 _p_ a_, int_4 _p_ lda_);
   505  snr_, real_8 _p_ csr_, real_8 _p_ snl_, real_8 _p_ csl_);
   507  int_4 _p_ incx_);
   509  , int_4 _p_ lda_, real_8 _p_ x_, real_8 _p_ scale_, real_8 _p_ cnorm_, int_4 _p_ info_);
   511  real_8 _p_ work_, int_4 _p_ info_);
   513  real_8 _p_ tau_, real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ info_);
   515  real_8 _p_ work_, int_4 _p_ info_);
   517  real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ info_);
   519  real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ info_);
   521  int_4 _p_ lda_, real_8 _p_ tau_, real_8 _p_ c_, int_4 _p_ ldc_, real_8 _p_ work_, int_4 _p_ info_);
   523  real_8 _p_ a_, int_4 _p_ lda_, real_8 _p_ tau_, real_8 _p_ c_, int_4 _p_ ldc_, real_8 _p_ work_, int_4 _p_ lwork_, 
   524  int_4 _p_ info_);
   526  int_4 _p_ lda_, real_8 _p_ tau_, real_8 _p_ c_, int_4 _p_ ldc_, real_8 _p_ work_, int_4 _p_ info_);
   528  int_4 _p_ lda_, real_8 _p_ tau_, real_8 _p_ c_, int_4 _p_ ldc_, real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ info_);
   530  int_4 _p_ lda_, real_8 _p_ tau_, real_8 _p_ c_, int_4 _p_ ldc_, real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ info_);
   532  rcond_, real_8 _p_ work_, int_4 _p_ iwork_, int_4 _p_ info_);
   538  int_4 _p_ n4_);
   540  , int_4 _p_ lwork_);
   541  int_4 _dbdsdc (char _p_ uplo_, char _p_ compq_, int_4 _p_ n_, real_8 _p_ d_, real_8 _p_ e_, real_8 _p_ u_, int_4 _p_ 
   542  ldu_, real_8 _p_ vt_, int_4 _p_ ldvt_, real_8 _p_ q_, int_4 _p_ iq_, real_8 _p_ work_, int_4 _p_ iwork_, int_4 _p_ 
   543  info_)
   544  { // ** body not listed **
   772  }
   773  
   774  int_4 _dbdsqr (char _p_ uplo_, int_4 _p_ n_, int_4 _p_ ncvt_, int_4 _p_ nru_, int_4 _p_ ncc_, real_8 _p_ d_, real_8 _p_ 
   775  e_, real_8 _p_ vt_, int_4 _p_ ldvt_, real_8 _p_ u_, int_4 _p_ ldu_, real_8 _p_ c_, int_4 _p_ ldc_, real_8 _p_ work_, 
   776  int_4 _p_ info_)
   777  { // ** body not listed **
  1214  }
  1215  
  1216  int_4 _dcombssq (real_8 _p_ v1_, real_8 _p_ v2_)
  1217  { // ** body not listed **
  1233  }
  1234  
  1235  int_4 _dgebd2 (int_4 _p_ m_, int_4 _p_ n_, real_8 _p_ a_, int_4 _p_ lda_, real_8 _p_ d_, real_8 _p_ e_, real_8 _p_ 
  1236  tauq_, real_8 _p_ taup_, real_8 _p_ work_, int_4 _p_ info_)
  1237  { // ** body not listed **
  1308  }
  1309  
  1310  int_4 _dgebrd (int_4 _p_ m_, int_4 _p_ n_, real_8 _p_ a_, int_4 _p_ lda_, real_8 _p_ d_, real_8 _p_ e_, real_8 _p_ 
  1311  tauq_, real_8 _p_ taup_, real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ info_)
  1312  { // ** body not listed **
  1399  }
  1400  
  1401  int_4 _dgejsv (char _p_ joba_, char _p_ jobu_, char _p_ jobv_, char _p_ jobr_, char _p_ jobt_, char _p_ jobp_, int_4 
  1402  _p_ m_, int_4 _p_ n_, real_8 _p_ a_, int_4 _p_ lda_, real_8 _p_ sva_, real_8 _p_ u_, int_4 _p_ ldu_, real_8 _p_ v_, 
  1403  int_4 _p_ ldv_, real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ iwork_, int_4 _p_ info_)
  1404  { // ** body not listed **
  2387  }
  2388  
  2389  int_4 _dgelq2 (int_4 _p_ m_, int_4 _p_ n_, real_8 _p_ a_, int_4 _p_ lda_, real_8 _p_ tau_, real_8 _p_ work_, int_4 _p_ 
  2390  info_)
  2391  { // ** body not listed **
  2425  }
  2426  
  2427  int_4 _dgelqf (int_4 _p_ m_, int_4 _p_ n_, real_8 _p_ a_, int_4 _p_ lda_, real_8 _p_ tau_, real_8 _p_ work_, int_4 _p_ 
  2428  lwork_, int_4 _p_ info_)
  2429  { // ** body not listed **
  2502  }
  2503  
  2504  int_4 _dgelsd (int_4 _p_ m_, int_4 _p_ n_, int_4 _p_ nrhs_, real_8 _p_ a_, int_4 _p_ lda_, real_8 _p_ b_, int_4 _p_ 
  2505  ldb_, real_8 _p_ s_, real_8 _p_ rcond_, int_4 _p_ rank_, real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ iwork_, int_4 
  2506  _p_ info_)
  2507  { // ** body not listed **
  2745  }
  2746  
  2747  int_4 _dgelss (int_4 _p_ m_, int_4 _p_ n_, int_4 _p_ nrhs_, real_8 _p_ a_, int_4 _p_ lda_, real_8 _p_ b_, int_4 _p_ 
  2748  ldb_, real_8 _p_ s_, real_8 _p_ rcond_, int_4 _p_ rank_, real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ info_)
  2749  { // ** body not listed **
  3107  }
  3108  
  3109  int_4 _dgeqp3 (int_4 _p_ m_, int_4 _p_ n_, real_8 _p_ a_, int_4 _p_ lda_, int_4 _p_ jpvt_, real_8 _p_ tau_, real_8 _p_ 
  3110  work_, int_4 _p_ lwork_, int_4 _p_ info_)
  3111  { // ** body not listed **
  3227  }
  3228  
  3229  int_4 _dgeqr2 (int_4 _p_ m_, int_4 _p_ n_, real_8 _p_ a_, int_4 _p_ lda_, real_8 _p_ tau_, real_8 _p_ work_, int_4 _p_ 
  3230  info_)
  3231  { // ** body not listed **
  3265  }
  3266  
  3267  int_4 _dgeqrf (int_4 _p_ m_, int_4 _p_ n_, real_8 _p_ a_, int_4 _p_ lda_, real_8 _p_ tau_, real_8 _p_ work_, int_4 _p_ 
  3268  lwork_, int_4 _p_ info_)
  3269  { // ** body not listed **
  3342  }
  3343  
  3344  int_4 _dgesdd (char _p_ jobz_, int_4 _p_ m_, int_4 _p_ n_, real_8 _p_ a_, int_4 _p_ lda_, real_8 _p_ s_, real_8 _p_ u_, 
  3345  int_4 _p_ ldu_, real_8 _p_ vt_, int_4 _p_ ldvt_, real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ iwork_, int_4 _p_ info_)
  3346  { // ** body not listed **
  3989  }
  3990  int_4 _dgesvd (char _p_ jobu_, char _p_ jobvt_, int_4 _p_ m_, int_4 _p_ n_, real_8 _p_ a_, int_4 _p_ lda_, real_8 _p_ 
  3991  s_, real_8 _p_ u_, int_4 _p_ ldu_, real_8 _p_ vt_, int_4 _p_ ldvt_, real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ info_
  3992  )
  3993  { // ** body not listed **
  5588  }
  5589  
  5590  int_4 _dgesvdq (char _p_ joba_, char _p_ jobp_, char _p_ jobr_, char _p_ jobu_, char _p_ jobv_, int_4 _p_ m_, int_4 _p_ 
  5591  n_, real_8 _p_ a_, int_4 _p_ lda_, real_8 _p_ s_, real_8 _p_ u_, int_4 _p_ ldu_, real_8 _p_ v_, int_4 _p_ ldv_, int_4 
  5592  _p_ numrank_, int_4 _p_ iwork_, int_4 _p_ liwork_, real_8 _p_ work_, int_4 _p_ lwork_, real_8 _p_ rwork_, int_4 _p_ 
  5593  lrwork_, int_4 _p_ info_)
  5594  { // ** body not listed **
  6417  }
  6418  
  6419  int_4 _dgesvj (char _p_ joba_, char _p_ jobu_, char _p_ jobv_, int_4 _p_ m_, int_4 _p_ n_, real_8 _p_ a_, int_4 _p_ 
  6420  lda_, real_8 _p_ sva_, int_4 _p_ mv_, real_8 _p_ v_, int_4 _p_ ldv_, real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ 
  6421  info_)
  6422  { // ** body not listed **
  7315  }
  7316  int_4 _dgsvj0 (char _p_ jobv_, int_4 _p_ m_, int_4 _p_ n_, real_8 _p_ a_, int_4 _p_ lda_, real_8 _p_ d_, real_8 _p_ 
  7317  sva_, int_4 _p_ mv_, real_8 _p_ v_, int_4 _p_ ldv_, real_8 _p_ eps_, real_8 _p_ sfmin_, real_8 _p_ tol_, int_4 _p_ 
  7318  nsweep_, real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ info_)
  7319  { // ** body not listed **
  7943  }
  7944  
  7945  int_4 _dgsvj1 (char _p_ jobv_, int_4 _p_ m_, int_4 _p_ n_, int_4 _p_ n1_, real_8 _p_ a_, int_4 _p_ lda_, real_8 _p_ d_, 
  7946  real_8 _p_ sva_, int_4 _p_ mv_, real_8 _p_ v_, int_4 _p_ ldv_, real_8 _p_ eps_, real_8 _p_ sfmin_, real_8 _p_ tol_, 
  7947  int_4 _p_ nsweep_, real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ info_)
  7948  { // ** body not listed **
  8339  }
  8340  
  8341  logical_4 _disnan (real_8 _p_ din_)
  8342  { // ** body not listed **
  8347  }
  8348  
  8349  int_4 _dlabad (real_8 _p_ small_, real_8 _p_ large_)
  8350  { // ** body not listed **
  8357  }
  8358  
  8359  int_4 _dlabrd (int_4 _p_ m_, int_4 _p_ n_, int_4 _p_ nb_, real_8 _p_ a_, int_4 _p_ lda_, real_8 _p_ d_, real_8 _p_ e_, 
  8360  real_8 _p_ tauq_, real_8 _p_ taup_, real_8 _p_ x_, int_4 _p_ ldx_, real_8 _p_ y_, int_4 _p_ ldy_)
  8361  { // ** body not listed **
  8485  }
  8486  
  8487  int_4 _dlacn2 (int_4 _p_ n_, real_8 _p_ v_, real_8 _p_ x_, int_4 _p_ isgn_, real_8 _p_ est_, int_4 _p_ kase_, int_4 _p_ 
  8488  isave_)
  8489  { // ** body not listed **
  8600  }
  8601  
  8602  int_4 _dlacpy (char _p_ uplo_, int_4 _p_ m_, int_4 _p_ n_, real_8 _p_ a_, int_4 _p_ lda_, real_8 _p_ b_, int_4 _p_ ldb_
  8603  )
  8604  { // ** body not listed **
  8630  }
  8631  
  8632  int_4 _dlaed6 (int_4 _p_ kniter_, logical_4 _p_ orgati_, real_8 _p_ rho_, real_8 _p_ d_, real_8 _p_ z_, real_8 _p_ 
  8633  finit_, real_8 _p_ tau_, int_4 _p_ info_)
  8634  { // ** body not listed **
  8835  }
  8836  
  8837  logical_4 _dlaisnan (real_8 _p_ din1_, real_8 _p_ din2_)
  8838  { // ** body not listed **
  8843  }
  8844  
  8845  int_4 _dlals0 (int_4 _p_ icompq_, int_4 _p_ nl_, int_4 _p_ nr_, int_4 _p_ sqre_, int_4 _p_ nrhs_, real_8 _p_ b_, int_4 
  8846  _p_ ldb_, real_8 _p_ bx_, int_4 _p_ ldbx_, int_4 _p_ perm_, int_4 _p_ givptr_, int_4 _p_ givcol_, int_4 _p_ ldgcol_, 
  8847  real_8 _p_ givnum_, int_4 _p_ ldgnum_, real_8 _p_ poles_, real_8 _p_ difl_, real_8 _p_ difr_, real_8 _p_ z_, int_4 _p_ 
  8848  k_, real_8 _p_ c_, real_8 _p_ s_, real_8 _p_ work_, int_4 _p_ info_)
  8849  { // ** body not listed **
  9025  }
  9026  
  9027  int_4 _dlalsa (int_4 _p_ icompq_, int_4 _p_ smlsiz_, int_4 _p_ n_, int_4 _p_ nrhs_, real_8 _p_ b_, int_4 _p_ ldb_, 
  9028  real_8 _p_ bx_, int_4 _p_ ldbx_, real_8 _p_ u_, int_4 _p_ ldu_, real_8 _p_ vt_, int_4 _p_ k_, real_8 _p_ difl_, real_8 
  9029  _p_ difr_, real_8 _p_ z_, real_8 _p_ poles_, int_4 _p_ givptr_, int_4 _p_ givcol_, int_4 _p_ ldgcol_, int_4 _p_ perm_, 
  9030  real_8 _p_ givnum_, real_8 _p_ c_, real_8 _p_ s_, real_8 _p_ work_, int_4 _p_ iwork_, int_4 _p_ info_)
  9031  { // ** body not listed **
  9177  }
  9178  
  9179  int_4 _dlalsd (char _p_ uplo_, int_4 _p_ smlsiz_, int_4 _p_ n_, int_4 _p_ nrhs_, real_8 _p_ d_, real_8 _p_ e_, real_8 
  9180  _p_ b_, int_4 _p_ ldb_, real_8 _p_ rcond_, int_4 _p_ rank_, real_8 _p_ work_, int_4 _p_ iwork_, int_4 _p_ info_)
  9181  { // ** body not listed **
  9416  }
  9417  
  9418  int_4 _dlamrg (int_4 _p_ n1_, int_4 _p_ n2_, real_8 _p_ a_, int_4 _p_ dtrd1_, int_4 _p_ dtrd2_, int_4 _p_ index_)
  9419  { // ** body not listed **
  9468  }
  9469  
  9470  real_8 _dlange (char _p_ norm_, int_4 _p_ m_, int_4 _p_ n_, real_8 _p_ a_, int_4 _p_ lda_, real_8 _p_ work_)
  9471  { // ** body not listed **
  9533  }
  9534  
  9535  real_8 _dlanst (char _p_ norm_, int_4 _p_ n_, real_8 _p_ d_, real_8 _p_ e_)
  9536  { // ** body not listed **
  9589  }
  9590  
  9591  int_4 _dlapmt (logical_4 _p_ forwrd_, int_4 _p_ m_, int_4 _p_ n_, real_8 _p_ x_, int_4 _p_ ldx_, int_4 _p_ k_)
  9592  { // ** body not listed **
  9650  }
  9651  
  9652  real_8 _dlapy2 (real_8 _p_ x_, real_8 _p_ y_)
  9653  { // ** body not listed **
  9680  }
  9681  
  9682  int_4 _dlaqp2 (int_4 _p_ m_, int_4 _p_ n_, int_4 _p_ offset_, real_8 _p_ a_, int_4 _p_ lda_, int_4 _p_ jpvt_, real_8 
  9683  _p_ tau_, real_8 _p_ vn1_, real_8 _p_ vn2_, real_8 _p_ work_)
  9684  { // ** body not listed **
  9744  }
  9745  
  9746  int_4 _dlaqps (int_4 _p_ m_, int_4 _p_ n_, int_4 _p_ offset_, int_4 _p_ nb_, int_4 _p_ kb_, real_8 _p_ a_, int_4 _p_ 
  9747  lda_, int_4 _p_ jpvt_, real_8 _p_ tau_, real_8 _p_ vn1_, real_8 _p_ vn2_, real_8 _p_ auxv_, real_8 _p_ f_, int_4 _p_ 
  9748  ldf_)
  9749  { // ** body not listed **
  9847  }
  9848  
  9849  int_4 _dlarfb (char _p_ side_, char _p_ trans_, char _p_ direct_, char _p_ storev_, int_4 _p_ m_, int_4 _p_ n_, int_4 
  9850  _p_ k_, real_8 _p_ v_, int_4 _p_ ldv_, real_8 _p_ t_, int_4 _p_ ldt_, real_8 _p_ c_, int_4 _p_ ldc_, real_8 _p_ work_, 
  9851  int_4 _p_ ldwork_)
  9852  { // ** body not listed **
 10071  }
 10072  
 10073  int_4 _dlarf (char _p_ side_, int_4 _p_ m_, int_4 _p_ n_, real_8 _p_ v_, int_4 _p_ incv_, real_8 _p_ tau_, real_8 _p_ 
 10074  c_, int_4 _p_ ldc_, real_8 _p_ work_)
 10075  { // ** body not listed **
 10123  }
 10124  
 10125  int_4 _dlarfg (int_4 _p_ n_, real_8 _p_ alpha_, real_8 _p_ x_, int_4 _p_ incx_, real_8 _p_ tau_)
 10126  { // ** body not listed **
 10169  }
 10170  
 10171  int_4 _dlarft (char _p_ direct_, char _p_ storev_, int_4 _p_ n_, int_4 _p_ k_, real_8 _p_ v_, int_4 _p_ ldv_, real_8 
 10172  _p_ tau_, real_8 _p_ t_, int_4 _p_ ldt_)
 10173  { // ** body not listed **
 10286  }
 10287  
 10288  int_4 _dlartg (real_8 _p_ f_, real_8 _p_ g_, real_8 _p_ cs_, real_8 _p_ sn_, real_8 _p_ r_)
 10289  { // ** body not listed **
 10358  }
 10359  
 10360  int_4 _dlas2 (real_8 _p_ f_, real_8 _p_ g_, real_8 _p_ h_, real_8 _p_ ssmin_, real_8 _p_ ssmax_)
 10361  { // ** body not listed **
 10405  }
 10406  int_4 _dlascl (char _p_ type_, int_4 _p_ kl_, int_4 _p_ ku_, real_8 _p_ cfrom_, real_8 _p_ cto_, int_4 _p_ m_, int_4 
 10407  _p_ n_, real_8 _p_ a_, int_4 _p_ lda_, int_4 _p_ info_)
 10408  { // ** body not listed **
 10574  }
 10575  
 10576  int_4 _dlasd0 (int_4 _p_ n_, int_4 _p_ sqre_, real_8 _p_ d_, real_8 _p_ e_, real_8 _p_ u_, int_4 _p_ ldu_, real_8 _p_ 
 10577  vt_, int_4 _p_ ldvt_, int_4 _p_ smlsiz_, int_4 _p_ iwork_, real_8 _p_ work_, int_4 _p_ info_)
 10578  { // ** body not listed **
 10689  }
 10690  
 10691  int_4 _dlasd1 (int_4 _p_ nl_, int_4 _p_ nr_, int_4 _p_ sqre_, real_8 _p_ d_, real_8 _p_ alpha_, real_8 _p_ beta_, 
 10692  real_8 _p_ u_, int_4 _p_ ldu_, real_8 _p_ vt_, int_4 _p_ ldvt_, int_4 _p_ idxq_, int_4 _p_ iwork_, real_8 _p_ work_, 
 10693  int_4 _p_ info_)
 10694  { // ** body not listed **
 10754  }
 10755  
 10756  int_4 _dlasd2 (int_4 _p_ nl_, int_4 _p_ nr_, int_4 _p_ sqre_, int_4 _p_ k_, real_8 _p_ d_, real_8 _p_ z_, real_8 _p_ 
 10757  alpha_, real_8 _p_ beta_, real_8 _p_ u_, int_4 _p_ ldu_, real_8 _p_ vt_, int_4 _p_ ldvt_, real_8 _p_ dsigma_, real_8 
 10758  _p_ u2_, int_4 _p_ ldu2_, real_8 _p_ vt2_, int_4 _p_ ldvt2_, int_4 _p_ idxp_, int_4 _p_ idx_, int_4 _p_ idxc_, int_4 
 10759  _p_ idxq_, int_4 _p_ coltyp_, int_4 _p_ info_)
 10760  { // ** body not listed **
 10988  }
 10989  
 10990  int_4 _dlasd3 (int_4 _p_ nl_, int_4 _p_ nr_, int_4 _p_ sqre_, int_4 _p_ k_, real_8 _p_ d_, real_8 _p_ q_, int_4 _p_ 
 10991  ldq_, real_8 _p_ dsigma_, real_8 _p_ u_, int_4 _p_ ldu_, real_8 _p_ u2_, int_4 _p_ ldu2_, real_8 _p_ vt_, int_4 _p_ 
 10992  ldvt_, real_8 _p_ vt2_, int_4 _p_ ldvt2_, int_4 _p_ idxc_, int_4 _p_ ctot_, real_8 _p_ z_, int_4 _p_ info_)
 10993  { // ** body not listed **
 11158  }
 11159  
 11160  int_4 _dlasd4 (int_4 _p_ n_, int_4 _p_ i_, real_8 _p_ d_, real_8 _p_ z_, real_8 _p_ delta_, real_8 _p_ rho_, real_8 _p_ 
 11161  sigma_, real_8 _p_ work_, int_4 _p_ info_)
 11162  { // ** body not listed **
 11839  }
 11840  
 11841  int_4 _dlasd5 (int_4 _p_ i_, real_8 _p_ d_, real_8 _p_ z_, real_8 _p_ delta_, real_8 _p_ rho_, real_8 _p_ dsigma_, 
 11842  real_8 _p_ work_)
 11843  { // ** body not listed **
 11895  }
 11896  
 11897  int_4 _dlasd6 (int_4 _p_ icompq_, int_4 _p_ nl_, int_4 _p_ nr_, int_4 _p_ sqre_, real_8 _p_ d_, real_8 _p_ vf_, real_8 
 11898  _p_ vl_, real_8 _p_ alpha_, real_8 _p_ beta_, int_4 _p_ idxq_, int_4 _p_ perm_, int_4 _p_ givptr_, int_4 _p_ givcol_, 
 11899  int_4 _p_ ldgcol_, real_8 _p_ givnum_, int_4 _p_ ldgnum_, real_8 _p_ poles_, real_8 _p_ difl_, real_8 _p_ difr_, real_8 
 11900  _p_ z_, int_4 _p_ k_, real_8 _p_ c_, real_8 _p_ s_, real_8 _p_ work_, int_4 _p_ iwork_, int_4 _p_ info_)
 11901  { // ** body not listed **
 11968  }
 11969  
 11970  int_4 _dlasd7 (int_4 _p_ icompq_, int_4 _p_ nl_, int_4 _p_ nr_, int_4 _p_ sqre_, int_4 _p_ k_, real_8 _p_ d_, real_8 
 11971  _p_ z_, real_8 _p_ zw_, real_8 _p_ vf_, real_8 _p_ vfw_, real_8 _p_ vl_, real_8 _p_ vlw_, real_8 _p_ alpha_, real_8 _p_ 
 11972  beta_, real_8 _p_ dsigma_, int_4 _p_ idx_, int_4 _p_ idxp_, int_4 _p_ idxq_, int_4 _p_ perm_, int_4 _p_ givptr_, int_4 
 11973  _p_ givcol_, int_4 _p_ ldgcol_, real_8 _p_ givnum_, int_4 _p_ ldgnum_, real_8 _p_ c_, real_8 _p_ s_, int_4 _p_ info_)
 11974  { // ** body not listed **
 12167  }
 12168  int_4 _dlasd8 (int_4 _p_ icompq_, int_4 _p_ k_, real_8 _p_ d_, real_8 _p_ z_, real_8 _p_ vf_, real_8 _p_ vl_, real_8 
 12169  _p_ difl_, real_8 _p_ difr_, int_4 _p_ lddifr_, real_8 _p_ dsigma_, real_8 _p_ work_, int_4 _p_ info_)
 12170  { // ** body not listed **
 12262  }
 12263  
 12264  int_4 _dlasda (int_4 _p_ icompq_, int_4 _p_ smlsiz_, int_4 _p_ n_, int_4 _p_ sqre_, real_8 _p_ d_, real_8 _p_ e_, 
 12265  real_8 _p_ u_, int_4 _p_ ldu_, real_8 _p_ vt_, int_4 _p_ k_, real_8 _p_ difl_, real_8 _p_ difr_, real_8 _p_ z_, real_8 
 12266  _p_ poles_, int_4 _p_ givptr_, int_4 _p_ givcol_, int_4 _p_ ldgcol_, int_4 _p_ perm_, real_8 _p_ givnum_, real_8 _p_ c_
 12267  , real_8 _p_ s_, real_8 _p_ work_, int_4 _p_ iwork_, int_4 _p_ info_)
 12268  { // ** body not listed **
 12445  }
 12446  
 12447  int_4 _dlasdq (char _p_ uplo_, int_4 _p_ sqre_, int_4 _p_ n_, int_4 _p_ ncvt_, int_4 _p_ nru_, int_4 _p_ ncc_, real_8 
 12448  _p_ d_, real_8 _p_ e_, real_8 _p_ vt_, int_4 _p_ ldvt_, real_8 _p_ u_, int_4 _p_ ldu_, real_8 _p_ c_, int_4 _p_ ldc_, 
 12449  real_8 _p_ work_, int_4 _p_ info_)
 12450  { // ** body not listed **
 12590  }
 12591  
 12592  int_4 _dlasdt (int_4 _p_ n_, int_4 _p_ lvl_, int_4 _p_ nd_, int_4 _p_ inode_, int_4 _p_ ndiml_, int_4 _p_ ndimr_, int_4 
 12593  _p_ msub_)
 12594  { // ** body not listed **
 12625  }
 12626  
 12627  int_4 _dlaset (char _p_ uplo_, int_4 _p_ m_, int_4 _p_ n_, real_8 _p_ alpha_, real_8 _p_ beta_, real_8 _p_ a_, int_4 
 12628  _p_ lda_)
 12629  { // ** body not listed **
 12658  }
 12659  
 12660  int_4 _dlasq1 (int_4 _p_ n_, real_8 _p_ d_, real_8 _p_ e_, real_8 _p_ work_, int_4 _p_ info_)
 12661  { // ** body not listed **
 12729  }
 12730  
 12731  int_4 _dlasq2 (int_4 _p_ n_, real_8 _p_ z_, int_4 _p_ info_)
 12732  { // ** body not listed **
 13063  }
 13064  
 13065  int_4 _dlasq3 (int_4 _p_ i0_, int_4 _p_ n0_, real_8 _p_ z_, int_4 _p_ pp_, real_8 _p_ dmin_, real_8 _p_ sigma_, real_8 
 13066  _p_ desig_, real_8 _p_ qmax_, int_4 _p_ nfail_, int_4 _p_ iter_, int_4 _p_ ndiv_, logical_4 _p_ ieee_, int_4 _p_ ttype_
 13067  , real_8 _p_ dmin1_, real_8 _p_ dmin2_, real_8 _p_ dn_, real_8 _p_ dn1_, real_8 _p_ dn2_, real_8 _p_ g_, real_8 _p_ 
 13068  tau_)
 13069  { // ** body not listed **
 13218  }
 13219  
 13220  int_4 _dlasq4 (int_4 _p_ i0_, int_4 _p_ n0_, real_8 _p_ z_, int_4 _p_ pp_, int_4 _p_ n0in_, real_8 _p_ dmin_, real_8 
 13221  _p_ dmin1_, real_8 _p_ dmin2_, real_8 _p_ dn_, real_8 _p_ dn1_, real_8 _p_ dn2_, real_8 _p_ tau_, int_4 _p_ ttype_, 
 13222  real_8 _p_ g_)
 13223  { // ** body not listed **
 13448  }
 13449  
 13450  int_4 _dlasq5 (int_4 _p_ i0_, int_4 _p_ n0_, real_8 _p_ z_, int_4 _p_ pp_, real_8 _p_ tau_, real_8 _p_ sigma_, real_8 
 13451  _p_ dmin_, real_8 _p_ dmin1_, real_8 _p_ dmin2_, real_8 _p_ dn_, real_8 _p_ dnm1_, real_8 _p_ dnm2_, logical_4 _p_ 
 13452  ieee_, real_8 _p_ eps_)
 13453  { // ** body not listed **
 13678  }
 13679  
 13680  int_4 _dlasq6 (int_4 _p_ i0_, int_4 _p_ n0_, real_8 _p_ z_, int_4 _p_ pp_, real_8 _p_ dmin_, real_8 _p_ dmin1_, real_8 
 13681  _p_ dmin2_, real_8 _p_ dn_, real_8 _p_ dnm1_, real_8 _p_ dnm2_)
 13682  { // ** body not listed **
 13784  }
 13785  
 13786  int_4 _dlasr (char _p_ side_, char _p_ pivot_, char _p_ direct_, int_4 _p_ m_, int_4 _p_ n_, real_8 _p_ c_, real_8 _p_ 
 13787  s_, real_8 _p_ a_, int_4 _p_ lda_)
 13788  { // ** body not listed **
 13994  }
 13995  
 13996  int_4 _dlasrt (char _p_ id_, int_4 _p_ n_, real_8 _p_ d_, int_4 _p_ info_)
 13997  { // ** body not listed **
 14171  }
 14172  
 14173  int_4 _dlassq (int_4 _p_ n_, real_8 _p_ x_, int_4 _p_ incx_, real_8 _p_ scale_, real_8 _p_ sumsq_)
 14174  { // ** body not listed **
 14195  }
 14196  
 14197  int_4 _dlasv2 (real_8 _p_ f_, real_8 _p_ g_, real_8 _p_ h_, real_8 _p_ ssmin_, real_8 _p_ ssmax_, real_8 _p_ snr_, 
 14198  real_8 _p_ csr_, real_8 _p_ snl_, real_8 _p_ csl_)
 14199  { // ** body not listed **
 14314  }
 14315  
 14316  int_4 _dlaswp (int_4 _p_ n_, real_8 _p_ a_, int_4 _p_ lda_, int_4 _p_ k1_, int_4 _p_ k2_, int_4 _p_ ipiv_, int_4 _p_ 
 14317  incx_)
 14318  { // ** body not listed **
 14371  }
 14372  
 14373  int_4 _dlatrs (char _p_ uplo_, char _p_ trans_, char _p_ diag_, char _p_ normin_, int_4 _p_ n_, real_8 _p_ a_, int_4 
 14374  _p_ lda_, real_8 _p_ x_, real_8 _p_ scale_, real_8 _p_ cnorm_, int_4 _p_ info_)
 14375  { // ** body not listed **
 14723  }
 14724  
 14725  int_4 _dorg2r (int_4 _p_ m_, int_4 _p_ n_, int_4 _p_ k_, real_8 _p_ a_, int_4 _p_ lda_, real_8 _p_ tau_, real_8 _p_ 
 14726  work_, int_4 _p_ info_)
 14727  { // ** body not listed **
 14776  }
 14777  
 14778  int_4 _dorgbr (char _p_ vect_, int_4 _p_ m_, int_4 _p_ n_, int_4 _p_ k_, real_8 _p_ a_, int_4 _p_ lda_, real_8 _p_ tau_
 14779  , real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ info_)
 14780  { // ** body not listed **
 14894  }
 14895  
 14896  int_4 _dorgl2 (int_4 _p_ m_, int_4 _p_ n_, int_4 _p_ k_, real_8 _p_ a_, int_4 _p_ lda_, real_8 _p_ tau_, real_8 _p_ 
 14897  work_, int_4 _p_ info_)
 14898  { // ** body not listed **
 14951  }
 14952  
 14953  int_4 _dorglq (int_4 _p_ m_, int_4 _p_ n_, int_4 _p_ k_, real_8 _p_ a_, int_4 _p_ lda_, real_8 _p_ tau_, real_8 _p_ 
 14954  work_, int_4 _p_ lwork_, int_4 _p_ info_)
 14955  { // ** body not listed **
 15044  }
 15045  
 15046  int_4 _dorgqr (int_4 _p_ m_, int_4 _p_ n_, int_4 _p_ k_, real_8 _p_ a_, int_4 _p_ lda_, real_8 _p_ tau_, real_8 _p_ 
 15047  work_, int_4 _p_ lwork_, int_4 _p_ info_)
 15048  { // ** body not listed **
 15137  }
 15138  
 15139  int_4 _dorm2r (char _p_ side_, char _p_ trans_, int_4 _p_ m_, int_4 _p_ n_, int_4 _p_ k_, real_8 _p_ a_, int_4 _p_ lda_
 15140  , real_8 _p_ tau_, real_8 _p_ c_, int_4 _p_ ldc_, real_8 _p_ work_, int_4 _p_ info_)
 15141  { // ** body not listed **
 15221  }
 15222  
 15223  int_4 _dormbr (char _p_ vect_, char _p_ side_, char _p_ trans_, int_4 _p_ m_, int_4 _p_ n_, int_4 _p_ k_, real_8 _p_ a_
 15224  , int_4 _p_ lda_, real_8 _p_ tau_, real_8 _p_ c_, int_4 _p_ ldc_, real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ info_)
 15225  { // ** body not listed **
 15360  }
 15361  
 15362  int_4 _dorml2 (char _p_ side_, char _p_ trans_, int_4 _p_ m_, int_4 _p_ n_, int_4 _p_ k_, real_8 _p_ a_, int_4 _p_ lda_
 15363  , real_8 _p_ tau_, real_8 _p_ c_, int_4 _p_ ldc_, real_8 _p_ work_, int_4 _p_ info_)
 15364  { // ** body not listed **
 15444  }
 15445  
 15446  int_4 _dormlq (char _p_ side_, char _p_ trans_, int_4 _p_ m_, int_4 _p_ n_, int_4 _p_ k_, real_8 _p_ a_, int_4 _p_ lda_
 15447  , real_8 _p_ tau_, real_8 _p_ c_, int_4 _p_ ldc_, real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ info_)
 15448  { // ** body not listed **
 15566  }
 15567  
 15568  int_4 _dormqr (char _p_ side_, char _p_ trans_, int_4 _p_ m_, int_4 _p_ n_, int_4 _p_ k_, real_8 _p_ a_, int_4 _p_ lda_
 15569  , real_8 _p_ tau_, real_8 _p_ c_, int_4 _p_ ldc_, real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ info_)
 15570  { // ** body not listed **
 15681  }
 15682  
 15683  int_4 _dpocon (char _p_ uplo_, int_4 _p_ n_, real_8 _p_ a_, int_4 _p_ lda_, real_8 _p_ anorm_, real_8 _p_ rcond_, 
 15684  real_8 _p_ work_, int_4 _p_ iwork_, int_4 _p_ info_)
 15685  { // ** body not listed **
 15754  }
 15755  
 15756  int_4 _drscl (int_4 _p_ n_, real_8 _p_ sa_, real_8 _p_ sx_, int_4 _p_ incx_)
 15757  { // ** body not listed **
 15794  }
 15795  
 15796  int_4 _ieeeck (int_4 _p_ ispec_, real_4 _p_ zero_, real_4 _p_ one_)
 15797  { // ** body not listed **
 15877  }
 15878  
 15879  int_4 _iladlc (int_4 _p_ m_, int_4 _p_ n_, real_8 _p_ a_, int_4 _p_ lda_)
 15880  { // ** body not listed **
 15901  }
 15902  
 15903  int_4 _iladlr (int_4 _p_ m_, int_4 _p_ n_, real_8 _p_ a_, int_4 _p_ lda_)
 15904  { // ** body not listed **
 15925  }
 15926  
 15927  int_4 _ilaenv (int_4 _p_ ispec_, char _p_ name_, char _p_ opts_, int_4 _p_ n1_, int_4 _p_ n2_, int_4 _p_ n3_, int_4 _p_ 
 15928  n4_)
 15929  { // ** body not listed **
 16479  }
 16480  
 16481  int_4 _iparmq (int_4 _p_ ispec_, char _p_ name_, char _p_ opts_, int_4 _p_ n_, int_4 _p_ ilo_, int_4 _p_ ihi_, int_4 
 16482  _p_ lwork_)
 16483  { // ** body not listed **
 16597  }
 16598  


© 2002-2025 J.M. van der Veer (jmvdveer@xs4all.nl)