LAPACK  3.10.0 LAPACK: Linear Algebra PACKage

## ◆ F77_sspr2()

 void F77_sspr2 ( int * layout, char * uplow, int * n, float * alpha, float * x, int * incx, float * y, int * incy, float * ap )

Definition at line 529 of file c_sblas2.c.

530  {
531  float *A, *AP;
532  int i,j,k,LDA;
533  CBLAS_UPLO uplo;
534
535  get_uplo_type(uplow,&uplo);
536
537  if (*layout == TEST_ROW_MJR) {
538  LDA = *n;
539  A = ( float* )malloc( LDA*LDA*sizeof( float ) );
540  AP = ( float* )malloc( (((LDA+1)*LDA)/2)*sizeof( float ) );
541  if (uplo == CblasUpper) {
542  for( j=0, k=0; j<*n; j++ )
543  for( i=0; i<j+1; i++, k++ )
544  A[ LDA*i+j ]=ap[ k ];
545  for( i=0, k=0; i<*n; i++ )
546  for( j=i; j<*n; j++, k++ )
547  AP[ k ]=A[ LDA*i+j ];
548  }
549  else {
550  for( j=0, k=0; j<*n; j++ )
551  for( i=j; i<*n; i++, k++ )
552  A[ LDA*i+j ]=ap[ k ];
553  for( i=0, k=0; i<*n; i++ )
554  for( j=0; j<i+1; j++, k++ )
555  AP[ k ]=A[ LDA*i+j ];
556  }
557  cblas_sspr2( CblasRowMajor, uplo, *n, *alpha, x, *incx, y, *incy, AP );
558  if (uplo == CblasUpper) {
559  for( i=0, k=0; i<*n; i++ )
560  for( j=i; j<*n; j++, k++ )
561  A[ LDA*i+j ]=AP[ k ];
562  for( j=0, k=0; j<*n; j++ )
563  for( i=0; i<j+1; i++, k++ )
564  ap[ k ]=A[ LDA*i+j ];
565  }
566  else {
567  for( i=0, k=0; i<*n; i++ )
568  for( j=0; j<i+1; j++, k++ )
569  A[ LDA*i+j ]=AP[ k ];
570  for( j=0, k=0; j<*n; j++ )
571  for( i=j; i<*n; i++, k++ )
572  ap[ k ]=A[ LDA*i+j ];
573  }
574  free(A);
575  free(AP);
576  }
577  else
578  cblas_sspr2( CblasColMajor, uplo, *n, *alpha, x, *incx, y, *incy, ap );
579 }
CBLAS_UPLO
Definition: cblas.h:21
@ CblasUpper
Definition: cblas.h:21
@ CblasColMajor
Definition: cblas.h:19
@ CblasRowMajor
Definition: cblas.h:19
void cblas_sspr2(CBLAS_LAYOUT layout, CBLAS_UPLO Uplo, const CBLAS_INDEX N, const float alpha, const float *X, const CBLAS_INDEX incX, const float *Y, const CBLAS_INDEX incY, float *A)
Definition: cblas_sspr2.c:12
#define TEST_ROW_MJR
Definition: cblas_test.h:12
void get_uplo_type(char *type, CBLAS_UPLO *uplo)
Definition: auxiliary.c:18
#define LDA
Definition: example_user.c:12
Here is the call graph for this function:
Here is the caller graph for this function: