ATL_f77syr.c (bench/ATL

      1
      2
      3
      4
      5
      6
      7
      8
      9
     10
     11
     12
     13
     14
     15
     16
     17
     18
     19
     20
     21
     22
     23
     24
     25
     26
     27
     28
     29
     30
     31
     32
     33
     34
     35
     36
     37
     38
     39
     40
     41
     42
     43
     44
     45
     46
     47
     48
     49
     50
     51
     52
     53
     54
     55
     56
     57
     58
     59
     60
     61
     62
     63
     64
     65
     66
     67
     68
     69
     70
     71
     72
     73
     74
     75
     76
     77
     78
     79
     80
     81
     82
     83
     84
     85
     86
     87
     88
     89
     90
     91
     92
     93
     94
     95

/*

 *             Automatically Tuned Linear Algebra Software v3.10.1

 *                   (C) Copyright 1999 Antoine P. Petitet

 *

 * Redistribution and use in source and binary forms, with or without

 * modification, are permitted provided that the following conditions

 * are met:

 *   1. Redistributions of source code must retain the above copyright

 *      notice, this list of conditions and the following disclaimer.

 *   2. Redistributions in binary form must reproduce the above copyright

 *      notice, this list of conditions, and the following disclaimer in the

 *      documentation and/or other materials provided with the distribution.

 *   3. The name of the ATLAS group or the names of its contributers may

 *      not be used to endorse or promote products derived from this

 *      software without specific written permission.

 *

 * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS

 * ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED

 * TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR

 * PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE ATLAS GROUP OR ITS CONTRIBUTORS

 * BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR

 * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF

 * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS

 * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN

 * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)

 * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE

 * POSSIBILITY OF SUCH DAMAGE.

 *

 */

#include "atlas_misc.h"

#include "atlas_tst.h"

#include "atlas_f77blas.h"

void Mjoin( PATL, f77syr )

(

   const enum ATLAS_UPLO     UPLO,

   const int                 N,

   const SCALAR              ALPHA,

   const TYPE                * X,

   const int                 INCX,

   TYPE                      * A,

   const int                 LDA

)

{

#if defined( StringSunStyle )

   #if defined( ATL_FunkyInts )

   F77_INTEGER               ONE = 1;

   #else

   int                       ONE = 1;

   #endif

#elif defined( StringStructVal ) || defined( StringStructPtr )

   F77_CHAR                  fuplo;

#elif defined( StringCrayStyle )

   F77_CHAR                  fuplo;

#endif

   char                      cuplo;

#ifdef ATL_FunkyInts

   const F77_INTEGER         F77N = N, F77lda  = LDA, F77incx = INCX;

#else

   #define F77N              N

   #define F77lda            LDA

   #define F77incx           INCX

#endif

#ifdef TCPLX

   TYPE                      alpha[2];

   *alpha = *ALPHA; alpha[1] = ALPHA[1];

#else

   TYPE                      alpha = ALPHA;

#endif

   if( UPLO == AtlasUpper ) cuplo = 'U';

   else                     cuplo = 'L';

   if( INCX < 0 ) X -= ( ( 1 - N ) * INCX ) SHIFT;

#if   defined( StringSunStyle  )

   F77syr( &cuplo, &F77N, SADD alpha, X, &F77incx, A, &F77lda, ONE );

#elif defined( StringCrayStyle )

   fuplo = ATL_C2F_TransChar( cuplo );

   F77syr( fuplo,  &F77N, SADD alpha, X, &F77incx, A, &F77lda );

#elif defined( StringStructVal )

   fuplo.len = 1; fuplo.cp = &cuplo;

   F77syr( fuplo,  &F77N, SADD alpha, X, &F77incx, A, &F77lda );

#elif defined( StringStructPtr )

   fuplo.len = 1; fuplo.cp = &cuplo;

   F77syr( &fuplo, &F77N, SADD alpha, X, &F77incx, A, &F77lda );

#else

   (void) fprintf( stderr, "\n\nF77/C interface not defined!!\n\n" );

   exit( -1 );

#endif

}