1 /**************************************************************************************************
2 *                                                                                                 *
3 * This file is part of BLASFEO.                                                                   *
4 *                                                                                                 *
5 * BLASFEO -- BLAS for embedded optimization.                                                      *
6 * Copyright (C) 2019 by Gianluca Frison.                                                          *
7 * Developed at IMTEK (University of Freiburg) under the supervision of Moritz Diehl.              *
8 * All rights reserved.                                                                            *
9 *                                                                                                 *
10 * The 2-Clause BSD License                                                                        *
11 *                                                                                                 *
12 * Redistribution and use in source and binary forms, with or without                              *
13 * modification, are permitted provided that the following conditions are met:                     *
14 *                                                                                                 *
15 * 1. Redistributions of source code must retain the above copyright notice, this                  *
16 *    list of conditions and the following disclaimer.                                             *
17 * 2. Redistributions in binary form must reproduce the above copyright notice,                    *
18 *    this list of conditions and the following disclaimer in the documentation                    *
19 *    and/or other materials provided with the distribution.                                       *
20 *                                                                                                 *
21 * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND                 *
22 * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED                   *
23 * WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE                          *
24 * DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR                 *
25 * ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES                  *
26 * (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;                    *
27 * LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND                     *
28 * ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT                      *
29 * (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS                   *
30 * SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.                                    *
31 *                                                                                                 *
32 * Author: Gianluca Frison, gianluca.frison (at) imtek.uni-freiburg.de                             *
33 *                                                                                                 *
34 **************************************************************************************************/
35 
36 #include <stdlib.h>
37 #include <stdio.h>
38 
39 #if defined(LA_EXTERNAL_BLAS_WRAPPER)
40 #if defined(EXTERNAL_BLAS_BLIS)
41 #include "blis.h"
42 #elif defined(EXTERNAL_BLAS_MKL)
43 #include "mkl.h"
44 #else
45 #include "../include/s_blas.h"
46 #endif
47 #endif
48 
49 #include "../include/blasfeo_common.h"
50 #include "../include/blasfeo_s_aux.h"
51 
52 
53 
54 #define REAL float
55 
56 #define STRMAT blasfeo_smat
57 #define STRVEC blasfeo_svec
58 
59 #define GEMV_N_LIBSTR blasfeo_sgemv_n
60 #define GEMV_NT_LIBSTR blasfeo_sgemv_nt
61 #define GEMV_T_LIBSTR blasfeo_sgemv_t
62 #define SYMV_L_LIBSTR blasfeo_ssymv_l
63 #define TRMV_LNN_LIBSTR blasfeo_strmv_lnn
64 #define TRMV_LTN_LIBSTR blasfeo_strmv_ltn
65 #define TRMV_UNN_LIBSTR blasfeo_strmv_unn
66 #define TRMV_UTN_LIBSTR blasfeo_strmv_utn
67 #define TRSV_LNN_LIBSTR blasfeo_strsv_lnn
68 #define TRSV_LNN_MN_LIBSTR blasfeo_strsv_lnn_mn
69 #define TRSV_LNU_LIBSTR blasfeo_strsv_lnu
70 #define TRSV_LTN_LIBSTR blasfeo_strsv_ltn
71 #define TRSV_LTN_MN_LIBSTR blasfeo_strsv_ltn_mn
72 #define TRSV_LTU_LIBSTR blasfeo_strsv_ltu
73 #define TRSV_UNN_LIBSTR blasfeo_strsv_unn
74 #define TRSV_UTN_LIBSTR blasfeo_strsv_utn
75 
76 #define COPY scopy_
77 #define GEMV sgemv_
78 #define SYMV ssymv_
79 #define TRMV strmv_
80 #define TRSV strsv_
81 
82 
83 
84 #include "x_blas2_lib.c"
85 
86