diff --git a/interface/lapack/larf.c b/interface/lapack/larf.c deleted file mode 100644 index 3b538c4a2..000000000 --- a/interface/lapack/larf.c +++ /dev/null @@ -1,109 +0,0 @@ -/*********************************************************************/ -/* Copyright 2009, 2010 The University of Texas at Austin. */ -/* All rights reserved. */ -/* */ -/* 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. */ -/* */ -/* THIS SOFTWARE IS PROVIDED BY THE UNIVERSITY OF TEXAS AT */ -/* AUSTIN ``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 UNIVERSITY OF TEXAS AT */ -/* AUSTIN OR 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. */ -/* */ -/* The views and conclusions contained in the software and */ -/* documentation are those of the authors and should not be */ -/* interpreted as representing official policies, either expressed */ -/* or implied, of The University of Texas at Austin. */ -/*********************************************************************/ - -#include -#include "common.h" -#ifdef FUNCTION_PROFILE -#include "functable.h" -#endif - -static int (*larf[])(blas_arg_t *, BLASLONG *, BLASLONG *, FLOAT *, FLOAT *, BLASLONG) = { - LARF_L, LARF_R, -}; - -int NAME(char *SIDE, blasint *M, blasint *N, FLOAT *v, blasint *incV, FLOAT *tau, FLOAT *c, blasint *ldC, FLOAT *work){ - - blas_arg_t args; - - FLOAT *buffer; -#ifdef PPC440 - extern -#endif - FLOAT *sa, *sb; - - char side_arg = *SIDE; - int side; - - PRINT_DEBUG_NAME; - - TOUPPER(side_arg); - - args.m = *M; - args.n = *N; - args.a = (void *)v; - args.lda = *incV; - args.c = (void *)c; - args.ldc = *ldC; - - args.alpha = (void *)tau; - - side = -1; - if (side_arg == 'L') side = 0; - if (side_arg == 'R') side = 1; - - if (args.m == 0 || args.n == 0) return 0; - -#ifndef COMPLEX - if (*tau == ZERO) return 0; -#else - if ((*(tau + 0) == ZERO) && (*(tau + 1) == ZERO)) return 0; -#endif - - IDEBUG_START; - - FUNCTION_PROFILE_START(); - -#ifndef PPC440 - buffer = (FLOAT *)blas_memory_alloc(1); - - sa = (FLOAT *)((BLASLONG)buffer + GEMM_OFFSET_A); - sb = (FLOAT *)(((BLASLONG)sa + ((GEMM_P * GEMM_Q * COMPSIZE * SIZE + GEMM_ALIGN) & ~GEMM_ALIGN)) + GEMM_OFFSET_B); -#endif - - larf[side](&args, NULL, NULL, sa, sb, 0); - -#ifndef PPC440 - blas_memory_free(buffer); -#endif - - FUNCTION_PROFILE_END(COMPSIZE * COMPSIZE, args.m * args.n, 2. / 3. * args.m * args.n * args.n); - - IDEBUG_END; - - return 0; -}