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.17
   111  */
   112  #if defined (__GNUC__)
   113  #pragma GCC diagnostic ignored "-Wimplicit-function-declaration"
   114  #pragma GCC diagnostic ignored "-Wincompatible-pointer-types"
   115  #if (__GNUC__ >= 14)
   116  #pragma GCC diagnostic ignored "-Wdeclaration-missing-parameter-type"
   117  #pragma GCC diagnostic ignored "-Wimplicit-int"
   118  #pragma GCC diagnostic ignored "-Wint-conversion"
   119  #pragma GCC diagnostic ignored "-Wreturn-mismatch"
   120  #endif
   121  #else
   122  #error VIF requires GCC
   123  #endif
   276  static CALLS __calls[__ncalls] = {
   277    {"dbdsdc", 0}, // subroutine
   278    {"dbdsqr", 0}, // subroutine
   279    {"dcombssq", 0}, // subroutine
   280    {"dgebd2", 0}, // subroutine
   281    {"dgebrd", 0}, // subroutine
   282    {"dgejsv", 0}, // subroutine
   283    {"dgelq2", 0}, // subroutine
   284    {"dgelqf", 0}, // subroutine
   285    {"dgelsd", 0}, // subroutine
   286    {"dgelss", 0}, // subroutine
   287    {"dgeqp3", 0}, // subroutine
   288    {"dgeqr2", 0}, // subroutine
   289    {"dgeqrf", 0}, // subroutine
   290    {"dgesdd", 0}, // subroutine
   291    {"dgesvd", 0}, // subroutine
   292    {"dgesvdq", 0}, // subroutine
   293    {"dgesvj", 0}, // subroutine
   294    {"dgsvj0", 0}, // subroutine
   295    {"dgsvj1", 0}, // subroutine
   296    {"disnan", 0}, // logical*4 function
   297    {"dlabad", 0}, // subroutine
   298    {"dlabrd", 0}, // subroutine
   299    {"dlacn2", 0}, // subroutine
   300    {"dlacpy", 0}, // subroutine
   301    {"dlaed6", 0}, // subroutine
   302    {"dlaisnan", 0}, // logical*4 function
   303    {"dlals0", 0}, // subroutine
   304    {"dlalsa", 0}, // subroutine
   305    {"dlalsd", 0}, // subroutine
   306    {"dlamrg", 0}, // subroutine
   307    {"dlange", 0}, // real*8 function
   308    {"dlanst", 0}, // real*8 function
   309    {"dlapmt", 0}, // subroutine
   310    {"dlapy2", 0}, // real*8 function
   311    {"dlaqp2", 0}, // subroutine
   312    {"dlaqps", 0}, // subroutine
   313    {"dlarfb", 0}, // subroutine
   314    {"dlarf", 0}, // subroutine
   315    {"dlarfg", 0}, // subroutine
   316    {"dlarft", 0}, // subroutine
   317    {"dlartg", 0}, // subroutine
   318    {"dlas2", 0}, // subroutine
   319    {"dlascl", 0}, // subroutine
   320    {"dlasd0", 0}, // subroutine
   321    {"dlasd1", 0}, // subroutine
   322    {"dlasd2", 0}, // subroutine
   323    {"dlasd3", 0}, // subroutine
   324    {"dlasd4", 0}, // subroutine
   325    {"dlasd5", 0}, // subroutine
   326    {"dlasd6", 0}, // subroutine
   327    {"dlasd7", 0}, // subroutine
   328    {"dlasd8", 0}, // subroutine
   329    {"dlasda", 0}, // subroutine
   330    {"dlasdq", 0}, // subroutine
   331    {"dlasdt", 0}, // subroutine
   332    {"dlaset", 0}, // subroutine
   333    {"dlasq1", 0}, // subroutine
   334    {"dlasq2", 0}, // subroutine
   335    {"dlasq3", 0}, // subroutine
   336    {"dlasq4", 0}, // subroutine
   337    {"dlasq5", 0}, // subroutine
   338    {"dlasq6", 0}, // subroutine
   339    {"dlasr", 0}, // subroutine
   340    {"dlasrt", 0}, // subroutine
   341    {"dlassq", 0}, // subroutine
   342    {"dlasv2", 0}, // subroutine
   343    {"dlaswp", 0}, // subroutine
   344    {"dlatrs", 0}, // subroutine
   345    {"dorg2r", 0}, // subroutine
   346    {"dorgbr", 0}, // subroutine
   347    {"dorgl2", 0}, // subroutine
   348    {"dorglq", 0}, // subroutine
   349    {"dorgqr", 0}, // subroutine
   350    {"dorm2r", 0}, // subroutine
   351    {"dormbr", 0}, // subroutine
   352    {"dorml2", 0}, // subroutine
   353    {"dormlq", 0}, // subroutine
   354    {"dormqr", 0}, // subroutine
   355    {"dpocon", 0}, // subroutine
   356    {"drscl", 0}, // subroutine
   357    {"ieeeck", 0}, // integer*4 function
   358    {"iladlc", 0}, // integer*4 function
   359    {"iladlr", 0}, // integer*4 function
   360    {"ilaenv", 0}, // integer*4 function
   361    {"iparmq", 0}, // integer*4 function
   362    {NULL, 0}
   363  };
   364  
   366  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_, 
   367  int_4 _p_ info_);
   369  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 
   370  _p_ work_, int_4 _p_ info_);
   373  real_8 _p_ tauq_, real_8 _p_ taup_, real_8 _p_ work_, int_4 _p_ info_);
   375  real_8 _p_ tauq_, real_8 _p_ taup_, real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ info_);
   377  , 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_ 
   378  v_, int_4 _p_ ldv_, real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ iwork_, int_4 _p_ info_);
   380  int_4 _p_ info_);
   382  int_4 _p_ lwork_, int_4 _p_ info_);
   384  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_
   385  , int_4 _p_ info_);
   387  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_)
   388  ;
   390  real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ info_);
   392  int_4 _p_ info_);
   394  int_4 _p_ lwork_, int_4 _p_ info_);
   396  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_, 
   397  int_4 _p_ info_);
   399  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_, 
   400  int_4 _p_ info_);
   402  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_ 
   403  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_, 
   404  int_4 _p_ lrwork_, int_4 _p_ info_);
   406  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_, 
   407  int_4 _p_ info_);
   409  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_, 
   410  int_4 _p_ nsweep_, real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ info_);
   412  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_, 
   413  real_8 _p_ tol_, int_4 _p_ nsweep_, real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ info_);
   417  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_);
   419  int_4 _p_ isave_);
   421  int_4 _p_ ldb_);
   423  real_8 _p_ finit_, real_8 _p_ tau_, int_4 _p_ info_);
   426  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_ 
   427  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_, 
   428  int_4 _p_ k_, real_8 _p_ c_, real_8 _p_ s_, real_8 _p_ work_, int_4 _p_ info_);
   430  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_, 
   431  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_ 
   432  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_);
   434  , 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_ 
   435  info_);
   437  index_);
   441  );
   444  real_8 _p_ tau_, real_8 _p_ vn1_, real_8 _p_ vn2_, real_8 _p_ work_);
   446  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_, 
   447  int_4 _p_ ldf_);
   449  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 
   450  _p_ work_, int_4 _p_ ldwork_);
   452  real_8 _p_ c_, int_4 _p_ ldc_, real_8 _p_ work_);
   455  real_8 _p_ tau_, real_8 _p_ t_, int_4 _p_ ldt_);
   459  , int_4 _p_ n_, real_8 _p_ a_, int_4 _p_ lda_, int_4 _p_ info_);
   461  real_8 _p_ vt_, int_4 _p_ ldvt_, int_4 _p_ smlsiz_, int_4 _p_ iwork_, real_8 _p_ work_, int_4 _p_ info_);
   463  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_ 
   464  work_, int_4 _p_ info_);
   466  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_
   467  , 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_, 
   468  int_4 _p_ idxq_, int_4 _p_ coltyp_, int_4 _p_ info_);
   470  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_, 
   471  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_);
   473  real_8 _p_ sigma_, real_8 _p_ work_, int_4 _p_ info_);
   475  dsigma_, real_8 _p_ work_);
   477  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 
   478  _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_ 
   479  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_);
   481  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_, 
   482  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_ 
   483  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_, 
   484  int_4 _p_ info_);
   486  real_8 _p_ difl_, real_8 _p_ difr_, int_4 _p_ lddifr_, real_8 _p_ dsigma_, real_8 _p_ work_, int_4 _p_ info_);
   488  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_, 
   489  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 
   490  _p_ c_, real_8 _p_ s_, real_8 _p_ work_, int_4 _p_ iwork_, int_4 _p_ info_);
   492  , 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 
   493  _p_ ldc_, real_8 _p_ work_, int_4 _p_ info_);
   495  ndimr_, int_4 _p_ msub_);
   497  , int_4 _p_ lda_);
   501  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_, 
   502  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_
   503  , real_8 _p_ tau_);
   505  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_ 
   506  ttype_, real_8 _p_ g_);
   508  , 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 
   509  _p_ ieee_, real_8 _p_ eps_);
   511  dmin1_, real_8 _p_ dmin2_, real_8 _p_ dn_, real_8 _p_ dnm1_, real_8 _p_ dnm2_);
   513  real_8 _p_ s_, real_8 _p_ a_, int_4 _p_ lda_);
   517  snr_, real_8 _p_ csr_, real_8 _p_ snl_, real_8 _p_ csl_);
   519  int_4 _p_ incx_);
   521  , int_4 _p_ lda_, real_8 _p_ x_, real_8 _p_ scale_, real_8 _p_ cnorm_, int_4 _p_ info_);
   523  real_8 _p_ work_, int_4 _p_ info_);
   525  real_8 _p_ tau_, real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ info_);
   527  real_8 _p_ work_, int_4 _p_ info_);
   529  real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ info_);
   531  real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ info_);
   533  int_4 _p_ lda_, real_8 _p_ tau_, real_8 _p_ c_, int_4 _p_ ldc_, real_8 _p_ work_, int_4 _p_ info_);
   535  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_, 
   536  int_4 _p_ info_);
   538  int_4 _p_ lda_, real_8 _p_ tau_, real_8 _p_ c_, int_4 _p_ ldc_, real_8 _p_ work_, int_4 _p_ info_);
   540  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_);
   542  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_);
   544  rcond_, real_8 _p_ work_, int_4 _p_ iwork_, int_4 _p_ info_);
   550  int_4 _p_ n4_);
   552  , int_4 _p_ lwork_);
   553  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_ 
   554  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_ 
   555  info_)
   556  { // ** body not listed **
   784  }
   785  
   786  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_ 
   787  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_, 
   788  int_4 _p_ info_)
   789  { // ** body not listed **
  1226  }
  1227  
  1228  int_4 _dcombssq (real_8 _p_ v1_, real_8 _p_ v2_)
  1229  { // ** body not listed **
  1245  }
  1246  
  1247  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_ 
  1248  tauq_, real_8 _p_ taup_, real_8 _p_ work_, int_4 _p_ info_)
  1249  { // ** body not listed **
  1320  }
  1321  
  1322  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_ 
  1323  tauq_, real_8 _p_ taup_, real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ info_)
  1324  { // ** body not listed **
  1411  }
  1412  
  1413  int_4 _dgejsv (char _p_ joba_, char _p_ jobu_, char _p_ jobv_, char _p_ jobr_, char _p_ jobt_, char _p_ jobp_, int_4 
  1414  _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_, 
  1415  int_4 _p_ ldv_, real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ iwork_, int_4 _p_ info_)
  1416  { // ** body not listed **
  2399  }
  2400  
  2401  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_ 
  2402  info_)
  2403  { // ** body not listed **
  2437  }
  2438  
  2439  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_ 
  2440  lwork_, int_4 _p_ info_)
  2441  { // ** body not listed **
  2514  }
  2515  
  2516  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_ 
  2517  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 
  2518  _p_ info_)
  2519  { // ** body not listed **
  2757  }
  2758  
  2759  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_ 
  2760  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_)
  2761  { // ** body not listed **
  3119  }
  3120  
  3121  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_ 
  3122  work_, int_4 _p_ lwork_, int_4 _p_ info_)
  3123  { // ** body not listed **
  3239  }
  3240  
  3241  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_ 
  3242  info_)
  3243  { // ** body not listed **
  3277  }
  3278  
  3279  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_ 
  3280  lwork_, int_4 _p_ info_)
  3281  { // ** body not listed **
  3354  }
  3355  
  3356  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_, 
  3357  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_)
  3358  { // ** body not listed **
  4001  }
  4002  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_ 
  4003  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_
  4004  )
  4005  { // ** body not listed **
  5600  }
  5601  
  5602  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_ 
  5603  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 
  5604  _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_ 
  5605  lrwork_, int_4 _p_ info_)
  5606  { // ** body not listed **
  6429  }
  6430  
  6431  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_ 
  6432  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_ 
  6433  info_)
  6434  { // ** body not listed **
  7327  }
  7328  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_ 
  7329  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_ 
  7330  nsweep_, real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ info_)
  7331  { // ** body not listed **
  7955  }
  7956  
  7957  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_, 
  7958  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_, 
  7959  int_4 _p_ nsweep_, real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ info_)
  7960  { // ** body not listed **
  8351  }
  8352  
  8353  logical_4 _disnan (real_8 _p_ din_)
  8354  { // ** body not listed **
  8359  }
  8360  
  8361  int_4 _dlabad (real_8 _p_ small_, real_8 _p_ large_)
  8362  { // ** body not listed **
  8369  }
  8370  
  8371  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_, 
  8372  real_8 _p_ tauq_, real_8 _p_ taup_, real_8 _p_ x_, int_4 _p_ ldx_, real_8 _p_ y_, int_4 _p_ ldy_)
  8373  { // ** body not listed **
  8497  }
  8498  
  8499  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_ 
  8500  isave_)
  8501  { // ** body not listed **
  8612  }
  8613  
  8614  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_
  8615  )
  8616  { // ** body not listed **
  8642  }
  8643  
  8644  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_ 
  8645  finit_, real_8 _p_ tau_, int_4 _p_ info_)
  8646  { // ** body not listed **
  8847  }
  8848  
  8849  logical_4 _dlaisnan (real_8 _p_ din1_, real_8 _p_ din2_)
  8850  { // ** body not listed **
  8855  }
  8856  
  8857  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 
  8858  _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_, 
  8859  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_ 
  8860  k_, real_8 _p_ c_, real_8 _p_ s_, real_8 _p_ work_, int_4 _p_ info_)
  8861  { // ** body not listed **
  9037  }
  9038  
  9039  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_, 
  9040  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 
  9041  _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_, 
  9042  real_8 _p_ givnum_, real_8 _p_ c_, real_8 _p_ s_, real_8 _p_ work_, int_4 _p_ iwork_, int_4 _p_ info_)
  9043  { // ** body not listed **
  9189  }
  9190  
  9191  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 
  9192  _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_)
  9193  { // ** body not listed **
  9428  }
  9429  
  9430  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_)
  9431  { // ** body not listed **
  9480  }
  9481  
  9482  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_)
  9483  { // ** body not listed **
  9545  }
  9546  
  9547  real_8 _dlanst (char _p_ norm_, int_4 _p_ n_, real_8 _p_ d_, real_8 _p_ e_)
  9548  { // ** body not listed **
  9601  }
  9602  
  9603  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_)
  9604  { // ** body not listed **
  9662  }
  9663  
  9664  real_8 _dlapy2 (real_8 _p_ x_, real_8 _p_ y_)
  9665  { // ** body not listed **
  9692  }
  9693  
  9694  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 
  9695  _p_ tau_, real_8 _p_ vn1_, real_8 _p_ vn2_, real_8 _p_ work_)
  9696  { // ** body not listed **
  9756  }
  9757  
  9758  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_ 
  9759  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_ 
  9760  ldf_)
  9761  { // ** body not listed **
  9859  }
  9860  
  9861  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 
  9862  _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_, 
  9863  int_4 _p_ ldwork_)
  9864  { // ** body not listed **
 10083  }
 10084  
 10085  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_ 
 10086  c_, int_4 _p_ ldc_, real_8 _p_ work_)
 10087  { // ** body not listed **
 10135  }
 10136  
 10137  int_4 _dlarfg (int_4 _p_ n_, real_8 _p_ alpha_, real_8 _p_ x_, int_4 _p_ incx_, real_8 _p_ tau_)
 10138  { // ** body not listed **
 10181  }
 10182  
 10183  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 
 10184  _p_ tau_, real_8 _p_ t_, int_4 _p_ ldt_)
 10185  { // ** body not listed **
 10298  }
 10299  
 10300  int_4 _dlartg (real_8 _p_ f_, real_8 _p_ g_, real_8 _p_ cs_, real_8 _p_ sn_, real_8 _p_ r_)
 10301  { // ** body not listed **
 10370  }
 10371  
 10372  int_4 _dlas2 (real_8 _p_ f_, real_8 _p_ g_, real_8 _p_ h_, real_8 _p_ ssmin_, real_8 _p_ ssmax_)
 10373  { // ** body not listed **
 10417  }
 10418  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 
 10419  _p_ n_, real_8 _p_ a_, int_4 _p_ lda_, int_4 _p_ info_)
 10420  { // ** body not listed **
 10586  }
 10587  
 10588  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_ 
 10589  vt_, int_4 _p_ ldvt_, int_4 _p_ smlsiz_, int_4 _p_ iwork_, real_8 _p_ work_, int_4 _p_ info_)
 10590  { // ** body not listed **
 10701  }
 10702  
 10703  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_, 
 10704  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_, 
 10705  int_4 _p_ info_)
 10706  { // ** body not listed **
 10766  }
 10767  
 10768  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_ 
 10769  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 
 10770  _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 
 10771  _p_ idxq_, int_4 _p_ coltyp_, int_4 _p_ info_)
 10772  { // ** body not listed **
 11000  }
 11001  
 11002  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_ 
 11003  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_ 
 11004  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_)
 11005  { // ** body not listed **
 11170  }
 11171  
 11172  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_ 
 11173  sigma_, real_8 _p_ work_, int_4 _p_ info_)
 11174  { // ** body not listed **
 11851  }
 11852  
 11853  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_, 
 11854  real_8 _p_ work_)
 11855  { // ** body not listed **
 11907  }
 11908  
 11909  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 
 11910  _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_, 
 11911  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 
 11912  _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_)
 11913  { // ** body not listed **
 11980  }
 11981  
 11982  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 
 11983  _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_ 
 11984  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 
 11985  _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_)
 11986  { // ** body not listed **
 12179  }
 12180  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 
 12181  _p_ difl_, real_8 _p_ difr_, int_4 _p_ lddifr_, real_8 _p_ dsigma_, real_8 _p_ work_, int_4 _p_ info_)
 12182  { // ** body not listed **
 12274  }
 12275  
 12276  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_, 
 12277  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 
 12278  _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_
 12279  , real_8 _p_ s_, real_8 _p_ work_, int_4 _p_ iwork_, int_4 _p_ info_)
 12280  { // ** body not listed **
 12457  }
 12458  
 12459  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 
 12460  _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_, 
 12461  real_8 _p_ work_, int_4 _p_ info_)
 12462  { // ** body not listed **
 12602  }
 12603  
 12604  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 
 12605  _p_ msub_)
 12606  { // ** body not listed **
 12637  }
 12638  
 12639  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 
 12640  _p_ lda_)
 12641  { // ** body not listed **
 12670  }
 12671  
 12672  int_4 _dlasq1 (int_4 _p_ n_, real_8 _p_ d_, real_8 _p_ e_, real_8 _p_ work_, int_4 _p_ info_)
 12673  { // ** body not listed **
 12741  }
 12742  
 12743  int_4 _dlasq2 (int_4 _p_ n_, real_8 _p_ z_, int_4 _p_ info_)
 12744  { // ** body not listed **
 13075  }
 13076  
 13077  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 
 13078  _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_
 13079  , 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_ 
 13080  tau_)
 13081  { // ** body not listed **
 13230  }
 13231  
 13232  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 
 13233  _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_, 
 13234  real_8 _p_ g_)
 13235  { // ** body not listed **
 13460  }
 13461  
 13462  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 
 13463  _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_ 
 13464  ieee_, real_8 _p_ eps_)
 13465  { // ** body not listed **
 13690  }
 13691  
 13692  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 
 13693  _p_ dmin2_, real_8 _p_ dn_, real_8 _p_ dnm1_, real_8 _p_ dnm2_)
 13694  { // ** body not listed **
 13796  }
 13797  
 13798  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_ 
 13799  s_, real_8 _p_ a_, int_4 _p_ lda_)
 13800  { // ** body not listed **
 14006  }
 14007  
 14008  int_4 _dlasrt (char _p_ id_, int_4 _p_ n_, real_8 _p_ d_, int_4 _p_ info_)
 14009  { // ** body not listed **
 14183  }
 14184  
 14185  int_4 _dlassq (int_4 _p_ n_, real_8 _p_ x_, int_4 _p_ incx_, real_8 _p_ scale_, real_8 _p_ sumsq_)
 14186  { // ** body not listed **
 14207  }
 14208  
 14209  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_, 
 14210  real_8 _p_ csr_, real_8 _p_ snl_, real_8 _p_ csl_)
 14211  { // ** body not listed **
 14326  }
 14327  
 14328  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_ 
 14329  incx_)
 14330  { // ** body not listed **
 14383  }
 14384  
 14385  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 
 14386  _p_ lda_, real_8 _p_ x_, real_8 _p_ scale_, real_8 _p_ cnorm_, int_4 _p_ info_)
 14387  { // ** body not listed **
 14735  }
 14736  
 14737  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_ 
 14738  work_, int_4 _p_ info_)
 14739  { // ** body not listed **
 14788  }
 14789  
 14790  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_
 14791  , real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ info_)
 14792  { // ** body not listed **
 14906  }
 14907  
 14908  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_ 
 14909  work_, int_4 _p_ info_)
 14910  { // ** body not listed **
 14963  }
 14964  
 14965  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_ 
 14966  work_, int_4 _p_ lwork_, int_4 _p_ info_)
 14967  { // ** body not listed **
 15056  }
 15057  
 15058  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_ 
 15059  work_, int_4 _p_ lwork_, int_4 _p_ info_)
 15060  { // ** body not listed **
 15149  }
 15150  
 15151  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_
 15152  , real_8 _p_ tau_, real_8 _p_ c_, int_4 _p_ ldc_, real_8 _p_ work_, int_4 _p_ info_)
 15153  { // ** body not listed **
 15233  }
 15234  
 15235  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_
 15236  , 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_)
 15237  { // ** body not listed **
 15372  }
 15373  
 15374  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_
 15375  , real_8 _p_ tau_, real_8 _p_ c_, int_4 _p_ ldc_, real_8 _p_ work_, int_4 _p_ info_)
 15376  { // ** body not listed **
 15456  }
 15457  
 15458  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_
 15459  , real_8 _p_ tau_, real_8 _p_ c_, int_4 _p_ ldc_, real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ info_)
 15460  { // ** body not listed **
 15578  }
 15579  
 15580  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_
 15581  , real_8 _p_ tau_, real_8 _p_ c_, int_4 _p_ ldc_, real_8 _p_ work_, int_4 _p_ lwork_, int_4 _p_ info_)
 15582  { // ** body not listed **
 15693  }
 15694  
 15695  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_, 
 15696  real_8 _p_ work_, int_4 _p_ iwork_, int_4 _p_ info_)
 15697  { // ** body not listed **
 15766  }
 15767  
 15768  int_4 _drscl (int_4 _p_ n_, real_8 _p_ sa_, real_8 _p_ sx_, int_4 _p_ incx_)
 15769  { // ** body not listed **
 15806  }
 15807  
 15808  int_4 _ieeeck (int_4 _p_ ispec_, real_4 _p_ zero_, real_4 _p_ one_)
 15809  { // ** body not listed **
 15889  }
 15890  
 15891  int_4 _iladlc (int_4 _p_ m_, int_4 _p_ n_, real_8 _p_ a_, int_4 _p_ lda_)
 15892  { // ** body not listed **
 15913  }
 15914  
 15915  int_4 _iladlr (int_4 _p_ m_, int_4 _p_ n_, real_8 _p_ a_, int_4 _p_ lda_)
 15916  { // ** body not listed **
 15937  }
 15938  
 15939  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_ 
 15940  n4_)
 15941  { // ** body not listed **
 16491  }
 16492  
 16493  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 
 16494  _p_ lwork_)
 16495  { // ** body not listed **
 16609  }
 16610  


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