amd/blis
1
0
Fork 0
mirror of https://github.com/amd/blis.git synced 2026-05-11 17:50:00 +00:00
Code Issues Packages Projects Releases Wiki Activity
Files
9d10d7dd9bc92a993fea7162bfa5983f75506f49
blis/test/test_blis2.c
Field G. Van Zee b6ef84fad1 Allow ldim of packed micro-panels != MR, NR. 
Details:
- Made substantial changes throughout the framework to decouple the leading
  dimension (row or column stride) used within each packed micro-panel from
  the corresponding register blocksize. It appears advantageous on some
  systems to use, for example, packed micro-panels of A where the column
  stride is greater than MR (whereas previously it was always equal to MR).
- Changes include:
  - Added BLIS_EXTEND_[MNK]R_? macros, which specify how much extra padding
    to use when packing micro-panels of A and B.
  - Adjusted all packing routines and macro-kernels to use PACKMR and PACKNR
    where appropriate, instead of MR and NR.
  - Added pd field (panel dimension) to obj_t.
  - New interface to bli_packm_cntl_obj_create().
  - Renamed bli_obj_packed_length()/_width() macros to
    bli_obj_padded_length()/_width().
  - Removed local #defines for cache/register blocksizes in level-3 *_cntl.c.
  - Print out new cache and register blocksize extensions in test suite.
- Also added new BLIS_EXTEND_[MNK]C_? macros for future use in using a larger
  blocksize for edge cases, which can improve performance at the margins.
2013-04-21 15:00:24 -05:00

2478 lines
63 KiB
C
Raw Blame History

/*
BLIS
An object-based framework for developing high-performance BLAS-like
libraries.
Copyright (C) 2013, The University of Texas
Redistribution and use in source and binary forms, with or without
modification, are permitted provided that the following conditions are
met:
- Redistributions of source code must retain the above copyright
notice, this list of conditions and the following disclaimer.
- 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.
- Neither the name of The University of Texas nor the names of its
contributors may be used to endorse or promote products derived
from this software without specific prior 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 COPYRIGHT
HOLDER 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.
*/
#include <unistd.h>
#include "blis.h"
double FLA_Clock( void );
extern gemm_t* gemm_cntl;
// trans m n alpha a lda x incx beta y incy
void dgemv_( char*, int*, int*, double*, double*, int*, double*, int*, double*, double*, int* );
void dger_( int*, int*, double*, double*, int*, double*, int*, double*, int* );
void dsymv_( char*, int*, double*, double*, int*, double*, int*, double*, double*, int* );
void dsyr_( char*, int*, double*, double*, int*, double*, int* );
void dsyr2_( char*, int*, double*, double*, int*, double*, int*, double*, int* );
void dtrmv_( char*, char*, char*, int*, double*, int*, double*, int* );
void dtrsv_( char*, char*, char*, int*, double*, int*, double*, int* );
// trans trans m n k alpha a lda b ldb beta c ldc
void dgemm_( char*, char*, int*, int*, int*, double*, double*, int*, double*, int*, double*, double*, int* );
int main( int argc, char** argv )
{
#if 1
{
obj_t a, b, c, w;
//obj_t a1, b1, c1;
//obj_t a11, b11, c11;
obj_t a_pack, b_pack, c_pack;
obj_t c_save;
obj_t alpha, beta;
dim_t m, n, k;
//dim_t b_part;
dim_t p;
dim_t p_begin, p_end, p_inc;
int m_input, n_input, k_input;
num_t dt_a, dt_b, dt_c;
num_t dt_alpha, dt_beta;
int r, n_repeats;
blksz_t* mr;
blksz_t* nr;
blksz_t* kr;
blksz_t* mc;
blksz_t* nc;
blksz_t* kc;
blksz_t* ni;
scalm_t* scalm_cntl;
packm_t* packm_cntl_a;
packm_t* packm_cntl_b;
packm_t* trmm_packm_cntl_a;
packm_t* trmm_packm_cntl_b;
packm_t* trsm_packm_cntl_a;
packm_t* trsm_packm_cntl_b;
unpackm_t* trsm_unpackm_cntl_b;
gemm_t* gemm_cntl_bp_ke;
gemm_t* gemm_cntl_op_bp;
gemm_t* gemm_cntl_mm_op;
gemm_t* gemm_cntl_vl_mm;
herk_t* herk_cntl_bp_ke;
herk_t* herk_cntl_op_bp;
herk_t* herk_cntl_mm_op;
herk_t* herk_cntl_vl_mm;
her2k_t* her2k_cntl_bp_ke;
her2k_t* her2k_cntl_op_bp;
her2k_t* her2k_cntl_mm_op;
her2k_t* her2k_cntl_vl_mm;
trmm_t* trmm_cntl_bp_ke;
trmm_t* trmm_cntl_op_bp;
trmm_t* trmm_cntl_mm_op;
trmm_t* trmm_cntl_vl_mm;
trmm_t* trmm3_cntl_mm_op;
trsm_t* trsm_cntl_bp_ke;
trsm_t* trsm_cntl_op_bp;
trsm_t* trsm_cntl_mm_op;
trsm_t* trsm_cntl_vl_mm;
double dtime;
double dtime_save;
double gflops;
bli_init();
n_repeats = 3;
#if 1
p_begin = 32;
p_end = 1600;
p_inc = 32;
#else
p_begin = 768;
p_end = 768;
p_inc = 1;
#endif
#if 1
m_input = -2;
//m_input = 128;
//m_input = 256;
n_input = -1;
//n_input = 1024;
//k_input = -1;
//k_input = 256;
k_input = -1;
#else
p_begin = 16;
p_end = 16;
p_inc = 1;
m_input = 15;
k_input = 13;
n_input = 15;
#endif
dt_a = BLIS_DOUBLE;
//dt_a = BLIS_DCOMPLEX;
dt_b = BLIS_DOUBLE;
dt_c = BLIS_DOUBLE;
//dt_c = BLIS_DCOMPLEX;
dt_alpha = BLIS_DOUBLE;
dt_beta = BLIS_DOUBLE;
//p = p_begin;
for ( p = p_begin; p <= p_end; p += p_inc )
{
if ( m_input < 0 ) m = p * ( dim_t )abs(m_input);
else m = ( dim_t ) m_input;
if ( n_input < 0 ) n = p * ( dim_t )abs(n_input);
else n = ( dim_t ) n_input;
if ( k_input < 0 ) k = p * ( dim_t )abs(k_input);
else k = ( dim_t ) k_input;
bli_obj_create( dt_alpha, 1, 1, 0, 0, &alpha );
bli_obj_create( dt_beta, 1, 1, 0, 0, &beta );
#if 1
bli_obj_create( dt_a, m, k, 0, 0, &a );
bli_obj_create( dt_b, k, n, 0, 0, &b );
bli_obj_create( dt_c, m, n, 0, 0, &c );
bli_obj_create( dt_c, m, n, 0, 0, &w );
bli_obj_create( dt_c, m, n, 0, 0, &c_save );
bli_randm( &a );
bli_randm( &b );
bli_randm( &c );
bli_randm( &w );
#elif 0
bli_obj_create( dt_a, m, k, 0, 0, &a );
bli_obj_create( dt_b, m, k, 0, 0, &b );
bli_obj_create( dt_c, m, m, 0, 0, &c );
bli_obj_create( dt_c, m, m, 0, 0, &c_save );
bli_randm( &a );
bli_randm( &b );
bli_randm( &c );
bli_setm( &BLIS_ZERO, &c );
#else
bli_obj_create( dt_a, m, m, 0, 0, &a );
bli_obj_create( dt_b, m, n, 0, 0, &b );
bli_obj_create( dt_c, m, n, 0, 0, &c );
bli_obj_create( dt_c, m, n, 0, 0, &c_save );
bli_obj_set_struc( BLIS_TRIANGULAR, a );
//bli_obj_set_uplo( BLIS_UPPER, a );
bli_obj_set_uplo( BLIS_LOWER, a );
//bli_obj_set_diag_offset( 4, a );
bli_randm( &a );
bli_randm( &c );
bli_randm( &b );
#endif
bli_obj_init_pack( &a_pack );
bli_obj_init_pack( &b_pack );
bli_obj_init_pack( &c_pack );
#if 0
bli_obj_set_struc( BLIS_TRIANGULAR, a );
//bli_obj_set_uplo( BLIS_LOWER, a );
bli_obj_set_uplo( BLIS_UPPER, a );
bli_obj_set_diag( BLIS_UNIT_DIAG, a );
bli_setm( &BLIS_ZERO, &a );
bli_obj_set_struc( BLIS_GENERAL, a );
bli_obj_set_uplo( BLIS_DENSE, a );
bli_obj_set_diag( BLIS_NONUNIT_DIAG, a );
#endif
bli_setsc( (2.0/1.0), 0.0, &alpha );
bli_setsc( -(1.0/1.0), 0.0, &beta );
mr = bli_blksz_obj_create( 2, 4, 2, 2 );
kr = bli_blksz_obj_create( 1, 1, 1, 1 );
nr = bli_blksz_obj_create( 1, 4, 1, 1 );
mc = bli_blksz_obj_create( 128, 384, 128, 128 );
kc = bli_blksz_obj_create( 256, 384, 256, 256 );
nc = bli_blksz_obj_create( 512, 512, 512, 512 );
ni = bli_blksz_obj_create( 16, 32, 16, 16 );
scalm_cntl =
bli_scalm_cntl_obj_create( BLIS_UNBLOCKED,
BLIS_VARIANT1 );
packm_cntl_a =
bli_packm_cntl_obj_create( BLIS_BLOCKED,
BLIS_VARIANT2,
mr, NULL,
kr, NULL,
TRUE, // scale?
TRUE, // densify?
FALSE, // invert diagonal?
FALSE, // reverse iteration if upper?
FALSE, // reverse iteration if lower?
BLIS_PACKED_ROW_PANELS );
packm_cntl_b =
bli_packm_cntl_obj_create( BLIS_BLOCKED,
BLIS_VARIANT2,
kr, NULL,
nr, NULL,
FALSE, // scale?
FALSE, // densify?
FALSE, // invert diagonal?
FALSE, // reverse iteration if upper?
FALSE, // reverse iteration if lower?
BLIS_PACKED_COL_PANELS );
gemm_cntl_bp_ke =
bli_gemm_cntl_obj_create( BLIS_UNB_OPT,
BLIS_VARIANT2,
NULL, NULL, NULL, NULL,
NULL, NULL, NULL, NULL );
gemm_cntl_op_bp =
bli_gemm_cntl_obj_create( BLIS_BLOCKED,
BLIS_VARIANT4,
//BLIS_VARIANT1,
mc,
ni,
NULL,
packm_cntl_a,
packm_cntl_b,
NULL,
gemm_cntl_bp_ke,
NULL );
gemm_cntl_mm_op =
bli_gemm_cntl_obj_create( BLIS_BLOCKED,
BLIS_VARIANT3,
kc,
NULL,
scalm_cntl,
NULL,
NULL,
NULL,
gemm_cntl_op_bp,
NULL );
gemm_cntl_vl_mm =
bli_gemm_cntl_obj_create( BLIS_BLOCKED,
BLIS_VARIANT2,
nc,
NULL,
NULL,
NULL,
NULL,
NULL,
gemm_cntl_mm_op,
NULL );
herk_cntl_bp_ke =
bli_herk_cntl_obj_create( BLIS_UNB_OPT,
BLIS_VARIANT2,
NULL, NULL, NULL, NULL,
NULL, NULL, NULL, NULL );
herk_cntl_op_bp =
bli_herk_cntl_obj_create( BLIS_BLOCKED,
//BLIS_VARIANT4,
BLIS_VARIANT1,
mc,
ni,
NULL,
packm_cntl_a,
packm_cntl_b,
NULL,
herk_cntl_bp_ke,
NULL );
herk_cntl_mm_op =
bli_herk_cntl_obj_create( BLIS_BLOCKED,
BLIS_VARIANT3,
kc,
NULL,
scalm_cntl,
NULL,
NULL,
NULL,
herk_cntl_op_bp,
NULL );
herk_cntl_vl_mm =
bli_herk_cntl_obj_create( BLIS_BLOCKED,
BLIS_VARIANT2,
nc,
NULL,
NULL,
NULL,
NULL,
NULL,
herk_cntl_mm_op,
NULL );
her2k_cntl_bp_ke =
bli_her2k_cntl_obj_create( BLIS_UNB_OPT,
BLIS_VARIANT2,
NULL, NULL, NULL, NULL, NULL,
NULL, NULL, NULL, NULL );
her2k_cntl_op_bp =
bli_her2k_cntl_obj_create( BLIS_BLOCKED,
//BLIS_VARIANT4,
BLIS_VARIANT1,
mc,
ni,
NULL,
packm_cntl_a,
packm_cntl_b,
NULL,
her2k_cntl_bp_ke,
herk_cntl_bp_ke,
NULL );
her2k_cntl_mm_op =
bli_her2k_cntl_obj_create( BLIS_BLOCKED,
BLIS_VARIANT3,
kc,
NULL,
scalm_cntl,
NULL,
NULL,
NULL,
her2k_cntl_op_bp,
NULL,
NULL );
her2k_cntl_vl_mm =
bli_her2k_cntl_obj_create( BLIS_BLOCKED,
BLIS_VARIANT2,
nc,
NULL,
NULL,
NULL,
NULL,
NULL,
her2k_cntl_mm_op,
NULL,
NULL );
trmm_packm_cntl_a =
bli_packm_cntl_obj_create( BLIS_BLOCKED,
BLIS_VARIANT3,
mr,
kr,
TRUE, // scale?
TRUE, // densify?
FALSE, // invert diagonal?
FALSE, // reverse iteration if upper?
FALSE, // reverse iteration if lower?
BLIS_PACKED_ROW_PANELS );
trmm_packm_cntl_b =
bli_packm_cntl_obj_create( BLIS_BLOCKED,
BLIS_VARIANT2,
kr,
nr,
FALSE, // scale?
FALSE, // densify?
FALSE, // invert diagonal?
FALSE, // reverse iteration if upper?
FALSE, // reverse iteration if lower?
BLIS_PACKED_COL_PANELS );
trmm_cntl_bp_ke =
bli_trmm_cntl_obj_create( BLIS_UNB_OPT,
BLIS_VARIANT2,
NULL, NULL, NULL, NULL, NULL,
NULL, NULL, NULL, NULL );
trmm_cntl_op_bp =
bli_trmm_cntl_obj_create( BLIS_BLOCKED,
//BLIS_VARIANT4,
BLIS_VARIANT1,
mc,
ni,
NULL,
trmm_packm_cntl_a,
trmm_packm_cntl_b,
NULL,
trmm_cntl_bp_ke,
gemm_cntl_bp_ke,
NULL );
trmm_cntl_mm_op =
bli_trmm_cntl_obj_create( BLIS_BLOCKED,
BLIS_VARIANT3,
kc,
NULL,
NULL, //scalm_cntl,
NULL,
NULL,
NULL,
trmm_cntl_op_bp,
NULL,
NULL );
trmm_cntl_vl_mm =
bli_trmm_cntl_obj_create( BLIS_BLOCKED,
BLIS_VARIANT2,
nc,
NULL,
NULL,
NULL,
NULL,
NULL,
trmm_cntl_mm_op,
NULL,
NULL );
trmm3_cntl_mm_op =
bli_trmm_cntl_obj_create( BLIS_BLOCKED,
BLIS_VARIANT3,
kc,
NULL,
scalm_cntl,
NULL,
NULL,
NULL,
trmm_cntl_op_bp,
NULL,
NULL );
trsm_packm_cntl_a =
bli_packm_cntl_obj_create( BLIS_BLOCKED,
BLIS_VARIANT3,
mr, // IMPORTANT: n dim multiple must be mr to
mr, // support right and bottom-right edge cases
FALSE, // scale?
TRUE, // densify?
TRUE, // invert diagonal?
TRUE, // reverse iteration if upper?
FALSE, // reverse iteration if lower?
BLIS_PACKED_ROW_PANELS );
trsm_packm_cntl_b =
bli_packm_cntl_obj_create( BLIS_BLOCKED,
BLIS_VARIANT2,
mr, // IMPORTANT: m dim multiple must be mr since
nr, // B_pack is updated (ie: serves as C) in trsm
TRUE, // scale?
FALSE, // densify?
FALSE, // invert diagonal?
FALSE, // reverse iteration if upper?
FALSE, // reverse iteration if lower?
BLIS_PACKED_COL_PANELS );
trsm_unpackm_cntl_b =
bli_unpackm_cntl_obj_create( BLIS_BLOCKED,
BLIS_VARIANT2,
NULL );
trsm_cntl_bp_ke =
bli_trsm_cntl_obj_create( BLIS_UNB_OPT,
//BLIS_VARIANT2,
BLIS_VARIANT3,
NULL, NULL, NULL, NULL, NULL,
NULL, NULL, NULL, NULL );
trsm_cntl_op_bp =
bli_trsm_cntl_obj_create( BLIS_BLOCKED,
BLIS_VARIANT4,
//BLIS_VARIANT1,
mc,
ni,
NULL,
trsm_packm_cntl_a,
trsm_packm_cntl_b,
NULL,
trsm_cntl_bp_ke,
gemm_cntl_bp_ke,
NULL );
trsm_cntl_mm_op =
bli_trsm_cntl_obj_create( BLIS_BLOCKED,
BLIS_VARIANT3,
kc,
NULL,
NULL, //scalm_cntl,
NULL,
NULL,
NULL,
trsm_cntl_op_bp,
NULL,
NULL );
trsm_cntl_vl_mm =
bli_trsm_cntl_obj_create( BLIS_BLOCKED,
BLIS_VARIANT2,
nc,
NULL,
NULL,
NULL,
NULL,
NULL,
trsm_cntl_mm_op,
NULL,
NULL );
//bli_printm( "a", &a, "%8.1e", "" );
//bli_printm( "b", &b, "%8.1e", "" );
//bli_printm( "c", &c, "%8.1e", "" );
//bli_printm( "alpha", &alpha, "%8.1e", "" );
//bli_printm( "beta", &beta, "%8.1e", "" );
bli_copym( &c, &c_save );
dtime_save = 1.0e9;
for ( r = 0; r < n_repeats; ++r )
{
bli_copym( &c_save, &c );
dtime = bli_clock();
#if 0
//bli_mm_clear_smem();
bli_packm_init( &a, &a_pack, packm_cntl_a );
bli_packm_int( &BLIS_ONE, &a, &a_pack, packm_cntl_a );
bli_packm_init( &b, &b_pack, packm_cntl_b );
bli_packm_int( &BLIS_ONE, &b, &b_pack, packm_cntl_b );
bli_gemm_ker_var2( &BLIS_ONE,
&a_pack,
&b_pack,
&BLIS_ONE,
&c,
NULL );
#endif
#if 0
bli_packm_init( &b, &b_pack, packm_cntl_b );
{
dim_t bn_alg = 32;
dim_t bn_use, j;
dim_t n_trans = bli_obj_width_after_trans( b );
obj_t c1_fuse, b1_fuse;
bli_packm_init( &a, &a_pack, packm_cntl_a );
bli_packm_int( &BLIS_ONE, &a, &a_pack, packm_cntl_a );
for ( j = 0; j < n_trans; j += bn_alg )
{
bn_use = bli_min( bn_alg, n_trans - j );
bli_acquire_mpart_l2r( BLIS_SUBPART1, j, bn_use, &c, &c1_fuse );
bli_acquire_mpart_l2r( BLIS_SUBPART1, j, bn_use, &b_pack, &b1_fuse );
bli_packm_int( &BLIS_ONE, &b, &b1_fuse, packm_cntl_b );
bli_gemm_ker_var2( &BLIS_ONE,
&a_pack,
&b1_fuse,
&BLIS_ONE,
&c1_fuse,
NULL );
}
}
#endif
#if 1
//bli_obj_set_struc( BLIS_HERMITIAN, a );
//bli_obj_set_uplo( BLIS_LOWER, a );
//bli_obj_set_uplo( BLIS_UPPER, a );
//bli_printm( "a", &a, "%4.1f", "" );
//bli_printm( "b", &b, "%4.1f", "" );
//bli_printm( "c", &c, "%4.1f", "" );
bli_gemm_int( &BLIS_ONE,
//bli_gemm_int( &alpha,
&a,
&b,
&BLIS_ONE,
&c,
//gemm_cntl_op_bp );
gemm_cntl_mm_op );
//bli_copym( &c, &w );
//bli_printm( "c after", &c, "%4.1f", "" );
//exit(1);
#endif
#if 0
bli_obj_set_struc( BLIS_HERMITIAN, a );
bli_obj_set_uplo( BLIS_LOWER, a );
//bli_obj_set_uplo( BLIS_UPPER, a );
bli_error_checking_level_set( BLIS_NO_ERROR_CHECKING );
bli_gemm_int( &BLIS_ONE,
&a,
&b,
&BLIS_ONE,
&c,
//gemm_cntl_op_bp );
gemm_cntl_mm_op );
#endif
#if 0
bli_obj_set_struc( BLIS_HERMITIAN, c );
bli_obj_set_uplo( BLIS_UPPER, c );
obj_t a1, c1, ar = a;
bli_obj_toggle_conj( ar );
bli_obj_toggle_trans( ar );
dim_t bm = bli_min( bli_obj_length( a ), 128 );
bli_acquire_mpart_t2b( BLIS_SUBPART1, 0, bm, &a, &a1 );
bli_acquire_mpart_t2b( BLIS_SUBPART1, 0, bm, &c, &c1 );
dtime = bli_clock();
bli_packm_init( &a1, &a_pack, packm_cntl_a );
bli_packm_int( &BLIS_ONE, &a1, &a_pack, packm_cntl_a );
bli_packm_init( &ar, &b_pack, packm_cntl_b );
bli_packm_int( &BLIS_ONE, &ar, &b_pack, packm_cntl_b );
//bli_printm( "a", &a, "%4.1f", "" );
//bli_printm( "c", &c, "%4.1f", "" );
bli_herk_u_ker_var2( &BLIS_ONE,
&a_pack,
&b_pack,
&BLIS_ONE,
&c1,
NULL );
//bli_printm( "c after", &c, "%4.1f", "" );
//exit(1);
#endif
#if 0
bli_obj_set_struc( BLIS_HERMITIAN, c );
//bli_obj_set_uplo( BLIS_LOWER, c );
bli_obj_set_uplo( BLIS_UPPER, c );
bli_error_checking_level_set( BLIS_NO_ERROR_CHECKING );
obj_t ah;
bli_obj_alias_with_trans( BLIS_CONJ_TRANSPOSE, a, ah );
//bli_printm( "a", &a, "%4.1f", "" );
//bli_printm( "c", &c, "%4.1f", "" );
bli_herk_int( &BLIS_ONE,
&a,
&ah,
&BLIS_ONE,
&c,
herk_cntl_op_bp );
//herk_cntl_mm_op );
//herk_cntl_vl_mm );
//bli_printm( "c after", &c, "%4.1f", "" );
//exit(1);
#endif
#if 0
bli_obj_set_struc( BLIS_HERMITIAN, c );
//bli_obj_set_uplo( BLIS_LOWER, c );
bli_obj_set_uplo( BLIS_UPPER, c );
bli_error_checking_level_set( BLIS_NO_ERROR_CHECKING );
obj_t ah, bh;
bli_obj_alias_with_trans( BLIS_CONJ_TRANSPOSE, a, ah );
bli_obj_alias_with_trans( BLIS_CONJ_TRANSPOSE, b, bh );
//bli_printm( "a", &a, "%4.1f", "" );
//bli_printm( "b", &b, "%4.1f", "" );
//bli_printm( "c", &c, "%4.1f", "" );
/*
bli_her2k_int( &BLIS_ONE,
&a,
&bh,
&BLIS_ONE,
&b,
&ah,
&BLIS_ONE,
&c,
her2k_cntl_op_bp );
//her2k_cntl_mm_op );
//her2k_cntl_vl_mm );
*/
bli_herk_int( &BLIS_ONE,
&a,
&bh,
&BLIS_ONE,
&c,
herk_cntl_op_bp );
bli_herk_int( &BLIS_ONE,
&b,
&ah,
&BLIS_ONE,
&c,
herk_cntl_op_bp );
//bli_printm( "c after", &c, "%4.1f", "" );
//exit(1);
#endif
#if 0
bli_obj_set_struc( BLIS_TRIANGULAR, a );
bli_obj_set_uplo( BLIS_LOWER, a );
//bli_obj_set_diag( BLIS_UNIT_DIAG, a );
bli_packm_init( &a, &a_pack, packm_cntl_a );
bli_packm_int( &BLIS_ONE, &a, &a_pack, packm_cntl_a );
bli_packm_init( &c, &b_pack, packm_cntl_b );
bli_packm_int( &BLIS_ONE, &c, &b_pack, packm_cntl_b );
bli_printm( "a", &a, "%4.1f", "" );
bli_printm( "b", &c, "%4.1f", "" );
bli_trmm_l_ker_var2( &BLIS_ONE,
&a_pack,
&b_pack,
&BLIS_ZERO,
&c,
NULL );
bli_printm( "b after", &c, "%4.1f", "" );
exit(1);
#endif
#if 0
bli_printm( "b", &c, "%4.1f", "" );
bli_printm( "a", &a, "%4.1f", "" );
obj_t a1, a2;
obj_t c1, c2;
dim_t off1 = 0;
dim_t off2 = 8;
dim_t m1 = off2;
dim_t m2 = m - m1;
bli_acquire_mpart_t2b( BLIS_SUBPART1, off1, m1, &a, &a1 );
bli_acquire_mpart_t2b( BLIS_SUBPART1, off2, m2, &a, &a2 );
bli_acquire_mpart_t2b( BLIS_SUBPART1, off1, m1, &c, &c1 );
bli_acquire_mpart_t2b( BLIS_SUBPART1, off2, m2, &c, &c2 );
bli_packm_init( &c, &b_pack, trmm_packm_cntl_b );
bli_packm_int( &BLIS_ONE, &c, &b_pack, trmm_packm_cntl_b );
bli_packm_init( &a1, &a_pack, trmm_packm_cntl_a );
bli_packm_int( &BLIS_ONE, &a1, &a_pack, trmm_packm_cntl_a );
bli_trmm_u_ker_var2( &BLIS_ONE,
&a_pack,
&b_pack,
&BLIS_ZERO,
&c1,
NULL );
bli_packm_init( &a2, &a_pack, trmm_packm_cntl_a );
bli_packm_int( &BLIS_ONE, &a2, &a_pack, trmm_packm_cntl_a );
bli_trmm_u_ker_var2( &BLIS_ONE,
&a_pack,
&b_pack,
&BLIS_ZERO,
&c2,
NULL );
bli_printm( "b after", &c, "%4.1f", "" );
exit(1);
#endif
#if 0
bli_obj_set_struc( BLIS_TRIANGULAR, a );
bli_obj_set_uplo( BLIS_LOWER, a );
//bli_obj_set_uplo( BLIS_UPPER, a );
bli_error_checking_level_set( BLIS_NO_ERROR_CHECKING );
//bli_printm( "a", &a, "%4.1f", "" );
//bli_printm( "b", &c, "%4.1f", "" );
bli_trmm_int( BLIS_LEFT,
&BLIS_ONE,
&a,
&c,
&BLIS_ZERO,
&c,
trmm_cntl_mm_op );
//trmm_cntl_vl_mm );
//bli_printm( "b after", &c, "%4.1f", "" );
//exit(1);
#endif
#if 0
bli_obj_set_struc( BLIS_TRIANGULAR, a );
//bli_obj_set_uplo( BLIS_LOWER, a );
bli_obj_set_uplo( BLIS_UPPER, a );
bli_error_checking_level_set( BLIS_NO_ERROR_CHECKING );
bli_printm( "a", &a, "%4.1f", "" );
bli_printm( "b", &b, "%4.1f", "" );
bli_printm( "c", &c, "%4.1f", "" );
bli_trmm_int( BLIS_LEFT,
&BLIS_ONE,
&a,
&b,
&BLIS_ONE,
&c,
trmm_cntl_mm_op );
//trmm_cntl_vl_mm );
bli_printm( "c after", &c, "%4.1f", "" );
exit(1);
#endif
#if 0
bli_printm( "b", &c, "%4.1f", "" );
bli_printm( "a", &a, "%4.1f", "" );
obj_t a1, a2;
obj_t c1, c2;
dim_t off1 = 0;
dim_t off2 = 8;
dim_t m1 = off2;
dim_t m2 = m - m1;
bli_acquire_mpart_t2b( BLIS_SUBPART1, off1, m1, &a, &a1 );
bli_acquire_mpart_t2b( BLIS_SUBPART1, off2, m2, &a, &a2 );
bli_acquire_mpart_t2b( BLIS_SUBPART1, off1, m1, &c, &c1 );
bli_acquire_mpart_t2b( BLIS_SUBPART1, off2, m2, &c, &c2 );
bli_packm_init( &c, &b_pack, trsm_packm_cntl_b );
bli_packm_int( &BLIS_ONE, &c, &b_pack, trsm_packm_cntl_b );
bli_packm_init( &a1, &a_pack, trsm_packm_cntl_a );
bli_packm_int( &BLIS_ONE, &a1, &a_pack, trsm_packm_cntl_a );
bli_trsm_l_ker_var2( &BLIS_ONE,
&a_pack,
&b_pack,
&BLIS_ZERO,
&c1,
NULL );
bli_packm_init( &a2, &a_pack, trsm_packm_cntl_a );
bli_packm_int( &BLIS_ONE, &a2, &a_pack, trsm_packm_cntl_a );
bli_trsm_l_ker_var2( &BLIS_ONE,
&a_pack,
&b_pack,
&BLIS_ZERO,
&c2,
NULL );
bli_printm( "b after", &c, "%4.1f", "" );
//bli_unpackm_int( &b_pack, &d, trsm_unpackm_cntl_b );
//bli_printm( "d (b after unpack2)", &d, "%4.1f", "" );
exit(1);
#endif
#if 0
bli_printm( "b", &c, "%4.1f", "" );
bli_printm( "a", &a, "%4.1f", "" );
obj_t a1, a2;
obj_t c1, c2;
dim_t off1 = 0;
dim_t off2 = 8;
dim_t m1 = off2;
dim_t m2 = m - m1;
bli_acquire_mpart_t2b( BLIS_SUBPART1, off1, m1, &a, &a1 );
bli_acquire_mpart_t2b( BLIS_SUBPART1, off2, m2, &a, &a2 );
bli_acquire_mpart_t2b( BLIS_SUBPART1, off1, m1, &c, &c1 );
bli_acquire_mpart_t2b( BLIS_SUBPART1, off2, m2, &c, &c2 );
bli_packm_init( &c, &b_pack, trsm_packm_cntl_b );
bli_packm_int( &BLIS_ONE, &c, &b_pack, trsm_packm_cntl_b );
bli_packm_init( &a2, &a_pack, trsm_packm_cntl_a );
bli_packm_int( &BLIS_ONE, &a2, &a_pack, trsm_packm_cntl_a );
bli_trsm_u_ker_var2( &BLIS_ONE,
&a_pack,
&b_pack,
&BLIS_ZERO,
&c2,
NULL );
bli_packm_init( &a1, &a_pack, trsm_packm_cntl_a );
bli_packm_int( &BLIS_ONE, &a1, &a_pack, trsm_packm_cntl_a );
bli_trsm_u_ker_var2( &BLIS_ONE,
&a_pack,
&b_pack,
&BLIS_ZERO,
&c1,
NULL );
bli_printm( "b after", &c, "%4.1f", "" );
//bli_printm( "b after", &c, "%7.4f", "" );
exit(1);
#endif
#if 0
bli_error_checking_level_set( BLIS_NO_ERROR_CHECKING );
/*
obj_t a1, a2;
obj_t c1, c2;
dim_t off1 = 0;
dim_t off2 = 128;
dim_t m1 = off2;
dim_t m2 = m - m1;
bli_acquire_mpart_t2b( BLIS_SUBPART1, off1, m1, &a, &a1 );
bli_acquire_mpart_t2b( BLIS_SUBPART1, off2, m2, &a, &a2 );
bli_acquire_mpart_t2b( BLIS_SUBPART1, off1, m1, &c, &c1 );
bli_acquire_mpart_t2b( BLIS_SUBPART1, off2, m2, &c, &c2 );
bli_packm_init( &c, &b_pack, trsm_packm_cntl_b );
bli_packm_int( &BLIS_ONE, &c, &b_pack, trsm_packm_cntl_b );
bli_packm_init( &a1, &a_pack, trsm_packm_cntl_a );
bli_packm_int( &BLIS_ONE, &a1, &a_pack, trsm_packm_cntl_a );
bli_trsm_l_ker_var2( &BLIS_ONE,
&a_pack,
&b_pack,
&BLIS_ZERO,
&c1,
NULL );
bli_packm_init( &a2, &a_pack, trsm_packm_cntl_a );
bli_packm_int( &BLIS_ONE, &a2, &a_pack, trsm_packm_cntl_a );
bli_trsm_l_ker_var2( &BLIS_ONE,
&a_pack,
&b_pack,
&BLIS_ZERO,
&c2,
NULL );
//bli_unpackm_int( &b_pack, &d, trsm_unpackm_cntl_b );
*/
/*
bli_packm_init( &c, &b_pack, trsm_packm_cntl_b );
bli_packm_int( &BLIS_ONE, &c, &b_pack, trsm_packm_cntl_b );
bli_packm_init( &a, &a_pack, trsm_packm_cntl_a );
bli_packm_int( &BLIS_ONE, &a, &a_pack, trsm_packm_cntl_a );
bli_trsm_l_ker_var2( &BLIS_ONE,
&a_pack,
&b_pack,
&BLIS_ZERO,
&c,
NULL );
bli_unpackm_int( &b_pack, &c, trsm_unpackm_cntl_b );
*/
//bli_printm( "d (b after unpack2)", &d, "%4.1f", "" );
//exit(1);
#endif
#if 0
bli_obj_set_struc( BLIS_TRIANGULAR, a );
bli_obj_set_uplo( BLIS_LOWER, a );
bli_error_checking_level_set( BLIS_NO_ERROR_CHECKING );
bli_trsm_int( BLIS_LEFT,
&BLIS_ONE,
&a,
&c,
&BLIS_ZERO,
&c,
trsm_cntl_op_bp );
#endif
#if 0
//bli_obj_set_struc( BLIS_TRIANGULAR, a );
//bli_obj_set_uplo( BLIS_LOWER, a );
bli_printm( "a", &a, "%4.1f", "" );
bli_printm( "b", &c, "%4.1f", "" );
bli_trsm_int( BLIS_LEFT,
&BLIS_ONE,
&a,
&c,
&BLIS_ZERO,
&c,
trsm_cntl_mm_op );
//trmm_cntl_vl_mm );
bli_printm( "c after", &c, "%4.1f", "" );
exit(1);
#endif
#if 0
//bli_obj_set_struc( BLIS_TRIANGULAR, a );
//bli_obj_set_uplo( BLIS_LOWER, a );
bli_error_checking_level_set( BLIS_NO_ERROR_CHECKING );
bli_trsm_int( BLIS_LEFT,
&BLIS_ONE,
&a,
&c,
&BLIS_ZERO,
&c,
//trsm_cntl_op_bp );
trsm_cntl_mm_op );
//trmm_cntl_vl_mm );
#endif
#if 0
char transa = 'N';
char transb = 'N';
int mm = bli_obj_length( c );
int kk = bli_obj_width_after_trans( a );
int nn = bli_obj_width( c );
int lda = bli_obj_col_stride( a );
int ldb = bli_obj_col_stride( b );
int ldc = bli_obj_col_stride( c );
dgemm_( &transa,
&transb,
&mm, &nn, &kk,
alpha.buffer,
a.buffer, &lda,
b.buffer, &ldb,
beta.buffer,
c.buffer, &ldc );
#endif
#if 0
bli_gemv( &alpha, &a, &b1, &beta, &c1 );
//bli_gemv_unf_var2( &alpha, &a, &b1, &beta, &c1, NULL );
//bli_gemv( &alpha, &a11, &b1, &beta, &c11 );
//bli_gemv_unf_var1( &alpha, &a, &b1, &beta, &c1, NULL );
//bli_gemv_unb_var2( &alpha, &a11, &b1, &beta, &c1, NULL );
//bli_gemv_unf_var2( &alpha, &a11, &b11, &beta, &c11, NULL );
//bli_gemv( &alpha, &a, &b, &beta, &c );
//bli_gemv_unb_var1( &alpha, &a, &b, &beta, &c, NULL );
//bli_gemv_unf_var1( &alpha, &a, &b, &beta, &c, NULL );
//bli_gemv_unb_var2( &alpha, &a, &b, &beta, &c, NULL );
//bli_gemv_unf_var2( &alpha, &a, &b, &beta, &c, NULL );
//#else
char trans = 'N';
//char trans = 'T';
//int mm = bli_obj_length( a11 );
//int kk = bli_obj_width( a11 );
//int lda = bli_obj_col_stride( a11 );
int mm = bli_obj_length( a );
int kk = bli_obj_width( a );
int lda = bli_obj_col_stride( a );
//int one = 1;
int incb = bli_obj_vector_inc( b1 );
int incc = bli_obj_vector_inc( c1 );
dgemv_( &trans,
&mm, &kk,
alpha.buffer,
a.buffer, &lda,
b.buffer, &incb,
beta.buffer,
c.buffer, &incc );
#endif
#if 0
//bli_ger_unb_var1( &alpha, &a, &b, &c, NULL );
//bli_ger_unb_var2( &alpha, &a, &b, &c, NULL );
//#else
int one = 1; int mm = m; int nn = n;
dger_( &mm, &nn,
alpha.buffer,
a.buffer, &one,
b.buffer, &one,
c.buffer, &mm );
#endif
#if 0
bli_obj_set_uplo( BLIS_LOWER, a );
//bli_hemv_unb_var1( &alpha, &a, &b, &beta, &c, NULL );
//bli_hemv_unb_var3( &alpha, &a, &b, &beta, &c, NULL );
//bli_hemv_unf_var1a( &alpha, &a, &b, &beta, &c, NULL );
//bli_hemv_unf_var3a( &alpha, &a, &b, &beta, &c, NULL );
bli_hemv_unf_var1( &alpha, &a, &b, &beta, &c, NULL );
//bli_hemv_unf_var3( &alpha, &a, &b, &beta, &c, NULL );
//#else
int one = 1; int mm = m; char uplo = 'L';
dsymv_( &uplo,
&mm,
alpha.buffer,
a.buffer, &mm,
b.buffer, &one,
beta.buffer,
c.buffer, &one );
#endif
#if 0
bli_obj_set_uplo( BLIS_LOWER, c );
bli_her_unb_var2( &alpha, &a, &c, NULL );
//#else
int one = 1; int mm = m; char uplo = 'L';
dsyr_( &uplo,
&mm,
alpha.buffer,
a.buffer, &one,
c.buffer, &mm );
#endif
#if 0
bli_obj_set_uplo( BLIS_LOWER, c );
//bli_her2( &alpha, &a, &b, &c );
//bli_her2_unb_var1( &alpha, &a, &b, &c, NULL );
//bli_her2_unb_var4( &alpha, &a, &b, &c, NULL );
//bli_her2_unf_var1( &alpha, &a, &b, &c, NULL );
//bli_her2_unf_var4( &alpha, &a, &b, &c, NULL );
//#else
int one = 1; int mm = bli_obj_col_stride( c ); char uplo = 'L';
dsyr2_( &uplo,
&mm,
alpha.buffer,
a.buffer, &one,
b.buffer, &one,
c.buffer, &mm );
#endif
#if 0
bli_setsc( (1.0/1.0), &alpha );
bli_obj_set_uplo( BLIS_LOWER, c );
bli_obj_set_conjtrans( BLIS_NO_TRANSPOSE, c );
//bli_trmv( &alpha, &c, &a );
//bli_trmv_unb_var1( &alpha, &c, &a, NULL );
//bli_trmv_unb_var2( &alpha, &c, &a, NULL );
bli_trmv_unf_var1( &alpha, &c, &a, NULL );
//bli_trmv_unf_var2( &alpha, &c, &a, NULL );
//#else
int one = 1; int mm = m; char uplo = 'U'; char trans = 'T'; char diag = 'N';
dtrmv_( &uplo,
&trans,
&diag,
&mm,
c.buffer, &mm,
a.buffer, &one );
#endif
#if 0
bli_setsc( (1.0/1.0), &alpha );
bli_obj_set_uplo( BLIS_UPPER, c );
bli_obj_set_conjtrans( BLIS_NO_TRANSPOSE, c );
//bli_trsv( &alpha, &c, &a );
//bli_trsv_unb_var1( &alpha, &c, &a, NULL );
//bli_trsv_unb_var2( &alpha, &c, &a, NULL );
bli_trsv_unf_var1( &alpha, &c, &a, NULL );
//bli_trsv_unf_var2( &alpha, &c, &a, NULL );
//#else
int one = 1; int mm = m; char uplo = 'L'; char trans = 'T'; char diag = 'N';
dtrsv_( &uplo,
&trans,
&diag,
&mm,
c.buffer, &mm,
a.buffer, &one );
#endif
//obj_t c_part;
//dim_t ij = 3;
//dim_t bb = 2;
//bli_acquire_mpart_tl2br( BLIS_SUBPART22,
// ij,
// bb,
// &c,
// &c_part );
//bli_printm( "c_part", &c_part, "%8.1e", "" );
dtime = bli_clock() - dtime;
dtime_save = bli_min( dtime, dtime_save );
}
//bli_printm( "c after", &c, "%8.1e", "" );
gflops = ( 2.0 * m * k * n ) / ( dtime_save * 1.0e9 );
//gflops = ( 2.0 * a_pack.m * k * n - a_pack.m * a_pack.m * k ) / ( dtime_save * 1.0e9 );
if ( bli_obj_is_hermitian( c ) || bli_obj_is_triangular( a ) ) gflops /= 2.0;
//gflops /= 4.0;
//gflops = ( 2.0 * bli_obj_length( a11 ) * bli_obj_width( a11 ) ) / ( dtime_save * 1.0e9 );
//gflops = ( 2.0 * bli_obj_length( a ) * bli_obj_width( a ) ) / ( dtime_save * 1.0e9 );
printf( "data_blis( %2ld, 1:5 ) = [ %4lu %4lu %4lu %10.3e %6.3f ];\n",
(p - p_begin + 1)/p_inc + 1, m, k, n, dtime_save, gflops );
bli_obj_release_pack( &a_pack );
bli_obj_release_pack( &b_pack );
bli_obj_release_pack( &c_pack );
bli_blksz_obj_free( mr );
bli_blksz_obj_free( nr );
bli_blksz_obj_free( kr );
bli_blksz_obj_free( mc );
bli_blksz_obj_free( nc );
bli_blksz_obj_free( kc );
bli_blksz_obj_free( ni );
bli_cntl_obj_free( scalm_cntl );
bli_cntl_obj_free( packm_cntl_a );
bli_cntl_obj_free( packm_cntl_b );
bli_cntl_obj_free( trmm_packm_cntl_a );
bli_cntl_obj_free( trmm_packm_cntl_b );
bli_cntl_obj_free( trsm_packm_cntl_a );
bli_cntl_obj_free( trsm_packm_cntl_b );
bli_cntl_obj_free( trsm_unpackm_cntl_b );
bli_cntl_obj_free( gemm_cntl_bp_ke );
bli_cntl_obj_free( gemm_cntl_op_bp );
bli_cntl_obj_free( gemm_cntl_mm_op );
bli_cntl_obj_free( gemm_cntl_vl_mm );
bli_cntl_obj_free( herk_cntl_bp_ke );
bli_cntl_obj_free( herk_cntl_op_bp );
bli_cntl_obj_free( herk_cntl_mm_op );
bli_cntl_obj_free( herk_cntl_vl_mm );
bli_cntl_obj_free( her2k_cntl_bp_ke );
bli_cntl_obj_free( her2k_cntl_op_bp );
bli_cntl_obj_free( her2k_cntl_mm_op );
bli_cntl_obj_free( her2k_cntl_vl_mm );
bli_cntl_obj_free( trmm_cntl_bp_ke );
bli_cntl_obj_free( trmm_cntl_op_bp );
bli_cntl_obj_free( trmm_cntl_mm_op );
bli_cntl_obj_free( trmm_cntl_vl_mm );
bli_cntl_obj_free( trmm3_cntl_mm_op );
bli_cntl_obj_free( trsm_cntl_bp_ke );
bli_cntl_obj_free( trsm_cntl_op_bp );
bli_cntl_obj_free( trsm_cntl_mm_op );
bli_cntl_obj_free( trsm_cntl_vl_mm );
bli_obj_free( &alpha );
bli_obj_free( &beta );
bli_obj_free( &a );
bli_obj_free( &b );
bli_obj_free( &c );
bli_obj_free( &w );
bli_obj_free( &c_save );
}
bli_finalize();
return 0;
}
#endif
// -- Level 1d -----------------------------------------------------------------
#if 0
{
num_t dt_alpha, dt_a, dt_b;
obj_t alpha, a, b;
dim_t m, n;
dt_alpha = BLIS_DOUBLE;
dt_a = BLIS_DOUBLE;
dt_b = BLIS_DOUBLE;
m = 6;
n = 11;
bli_init();
bli_obj_create( dt_alpha, 1, 1, 0, 0, &alpha );
bli_obj_create( dt_a, n, m, 0, 0, &a );
bli_obj_create( dt_b, m, n, 0, 0, &b );
bli_setsc( -(3.0/1.0), &alpha );
bli_setm( &BLIS_TWO, &a );
bli_setm( &BLIS_ONE, &b );
bli_printm( "alpha", &alpha, "%4.1f", "" );
bli_printm( "a before", &a, "%4.1f", "" );
bli_printm( "b before", &b, "%4.1f", "" );
//bli_obj_set_struc( BLIS_SYMMETRIC, a );
bli_obj_set_struc( BLIS_GENERAL, a );
bli_obj_set_uplo( BLIS_DENSE, a );
//bli_obj_set_uplo( BLIS_UPPER, a );
//bli_obj_set_uplo( BLIS_LOWER, a );
bli_obj_set_diag_offset( -7, a );
bli_obj_set_diag( BLIS_UNIT_DIAG, a );
bli_obj_set_trans( BLIS_TRANSPOSE, a );
//bli_obj_set_conj( BLIS_CONJUGATE, a );
bli_axpyd( &alpha, &a, &b );
bli_printm( "b after", &b, "%4.1f", "" );
bli_obj_free( &alpha );
bli_obj_free( &a );
bli_obj_free( &b );
bli_finalize();
return 0;
}
#endif
#if 0
{
num_t dt_a, dt_b;
obj_t a, b;
dim_t m, n;
dt_a = BLIS_DOUBLE;
dt_b = BLIS_DOUBLE;
m = 6;
n = 11;
bli_init();
bli_obj_create( dt_a, n, m, 0, 0, &a );
bli_obj_create( dt_b, m, n, 0, 0, &b );
bli_setm( &BLIS_TWO, &a );
bli_setm( &BLIS_ZERO, &b );
bli_printm( "a before", &a, "%4.1f", "" );
bli_printm( "b before", &b, "%4.1f", "" );
//bli_obj_set_struc( BLIS_SYMMETRIC, a );
bli_obj_set_struc( BLIS_GENERAL, a );
bli_obj_set_uplo( BLIS_DENSE, a );
//bli_obj_set_uplo( BLIS_UPPER, a );
//bli_obj_set_uplo( BLIS_LOWER, a );
bli_obj_set_diag_offset( -7, a );
bli_obj_set_diag( BLIS_NONUNIT_DIAG, a );
bli_obj_set_trans( BLIS_TRANSPOSE, a );
//bli_obj_set_conj( BLIS_CONJUGATE, a );
bli_copyd( &a, &b );
bli_printm( "b after", &b, "%4.1f", "" );
bli_obj_free( &a );
bli_obj_free( &b );
bli_finalize();
return 0;
}
#endif
#if 0
{
num_t dt_beta, dt_a, dt_b;
obj_t beta, a, b;
dim_t m, n;
dt_beta = BLIS_DOUBLE;
dt_a = BLIS_DOUBLE;
dt_b = BLIS_DOUBLE;
m = 6;
n = 11;
bli_init();
bli_obj_create( dt_beta, 1, 1, 0, 0, &beta );
bli_obj_create( dt_a, n, m, 0, 0, &a );
bli_obj_create( dt_b, m, n, 0, 0, &b );
bli_setsc( -(3.0/1.0), &beta );
bli_setm( &BLIS_MINUS_TWO, &a );
bli_setm( &BLIS_ONE, &b );
bli_printm( "beta", &beta, "%4.1f", "" );
bli_printm( "a before", &a, "%4.1f", "" );
bli_printm( "b before", &b, "%4.1f", "" );
//bli_obj_set_struc( BLIS_SYMMETRIC, a );
bli_obj_set_struc( BLIS_GENERAL, a );
bli_obj_set_uplo( BLIS_DENSE, a );
//bli_obj_set_uplo( BLIS_UPPER, a );
//bli_obj_set_uplo( BLIS_LOWER, a );
bli_obj_set_diag_offset( 2, a );
bli_obj_set_diag( BLIS_NONUNIT_DIAG, a );
bli_obj_set_trans( BLIS_TRANSPOSE, a );
//bli_obj_set_conj( BLIS_CONJUGATE, a );
bli_scal2d( &beta, &a, &b );
bli_printm( "b after", &b, "%4.1f", "" );
bli_obj_free( &beta );
bli_obj_free( &a );
bli_obj_free( &b );
bli_finalize();
return 0;
}
#endif
#if 0
{
num_t dt_a, dt_alpha;
obj_t alpha, a;
dim_t m, n;
dt_alpha = BLIS_DOUBLE;
dt_a = BLIS_DOUBLE;
m = 6;
n = 11;
bli_init();
bli_obj_create( dt_alpha, 1, 1, 0, 0, &alpha );
bli_obj_create( dt_a, m, n, 0, 0, &a );
bli_setsc( -(4.0/1.0), &alpha );
bli_setm( &BLIS_TWO, &a );
bli_printm( "alpha", &alpha, "%4.1f", "" );
bli_printm( "a before", &a, "%4.1f", "" );
//bli_obj_set_struc( BLIS_GENERAL, a );
//bli_obj_set_uplo( BLIS_DENSE, a );
//bli_obj_set_struc( BLIS_SYMMETRIC, a );
bli_obj_set_struc( BLIS_GENERAL, a );
bli_obj_set_uplo( BLIS_DENSE, a );
//bli_obj_set_uplo( BLIS_UPPER, a );
//bli_obj_set_uplo( BLIS_LOWER, a );
bli_obj_set_diag_offset( -2, a );
bli_obj_set_diag( BLIS_UNIT_DIAG, a );
bli_scald( &alpha, &a );
bli_printm( "a after", &a, "%4.1f", "" );
bli_obj_free( &alpha );
bli_obj_free( &a );
bli_finalize();
return 0;
}
#endif
#if 0
{
num_t dt_a, dt_alpha;
obj_t alpha, a;
dim_t m, n;
dt_alpha = BLIS_DOUBLE;
dt_a = BLIS_DOUBLE;
m = 6;
n = 11;
bli_init();
bli_obj_create( dt_alpha, 1, 1, 0, 0, &alpha );
bli_obj_create( dt_a, m, n, 0, 0, &a );
bli_setsc( -(4.0/1.0), &alpha );
bli_setm( &BLIS_TWO, &a );
bli_printm( "alpha", &alpha, "%4.1f", "" );
bli_printm( "a before", &a, "%4.1f", "" );
//bli_obj_set_struc( BLIS_GENERAL, a );
//bli_obj_set_uplo( BLIS_DENSE, a );
//bli_obj_set_struc( BLIS_SYMMETRIC, a );
bli_obj_set_struc( BLIS_GENERAL, a );
bli_obj_set_uplo( BLIS_DENSE, a );
//bli_obj_set_uplo( BLIS_UPPER, a );
//bli_obj_set_uplo( BLIS_LOWER, a );
bli_obj_set_diag_offset( 7, a );
bli_obj_set_diag( BLIS_NONUNIT_DIAG, a );
bli_setd( &alpha, &a );
bli_printm( "a after", &a, "%4.1f", "" );
bli_obj_free( &alpha );
bli_obj_free( &a );
bli_finalize();
return 0;
}
#endif
// -- Level 1m -----------------------------------------------------------------
#if 0
{
num_t dt_alpha, dt_a, dt_b;
obj_t alpha, a, b;
dim_t m, n;
dt_alpha = BLIS_DOUBLE;
dt_a = BLIS_DOUBLE;
dt_b = BLIS_DOUBLE;
m = 6;
n = 11;
bli_init();
bli_obj_create( dt_alpha, 1, 1, 0, 0, &alpha );
bli_obj_create( dt_a, n, m, 0, 0, &a );
bli_obj_create( dt_b, m, n, 0, 0, &b );
bli_setsc( -(3.0/1.0), &alpha );
bli_setm( &BLIS_TWO, &a );
bli_setm( &BLIS_ONE, &b );
bli_printm( "alpha", &alpha, "%4.1f", "" );
bli_printm( "a before", &a, "%4.1f", "" );
bli_printm( "b before", &b, "%4.1f", "" );
bli_obj_set_struc( BLIS_SYMMETRIC, a );
//bli_obj_set_struc( BLIS_GENERAL, a );
//bli_obj_set_uplo( BLIS_DENSE, a );
bli_obj_set_uplo( BLIS_UPPER, a );
//bli_obj_set_uplo( BLIS_LOWER, a );
bli_obj_set_diag_offset( 2, a );
bli_obj_set_diag( BLIS_UNIT_DIAG, a );
bli_obj_set_trans( BLIS_TRANSPOSE, a );
//bli_obj_set_conj( BLIS_CONJUGATE, a );
bli_axpym( &alpha, &a, &b );
bli_printm( "b after", &b, "%4.1f", "" );
bli_obj_free( &alpha );
bli_obj_free( &a );
bli_obj_free( &b );
bli_finalize();
return 0;
}
#endif
#if 0
{
num_t dt_a, dt_b;
obj_t a, b;
dim_t m, n;
dt_a = BLIS_DOUBLE;
dt_b = BLIS_DOUBLE;
m = 11;
n = 6;
bli_init();
bli_obj_create( dt_a, n, m, 0, 0, &a );
bli_obj_create( dt_b, m, n, 0, 0, &b );
//bli_setm( &BLIS_MINUS_ONE, &a );
bli_randm( &a );
bli_setm( &BLIS_ZERO, &b );
bli_printm( "a before", &a, "%4.1f", "" );
bli_printm( "b before", &b, "%4.1f", "" );
//bli_obj_set_struc( BLIS_GENERAL, a );
//bli_obj_set_uplo( BLIS_DENSE, a );
bli_obj_set_diag_offset( 2, a );
bli_obj_set_struc( BLIS_TRIANGULAR, a );
//bli_obj_set_uplo( BLIS_UPPER, a );
bli_obj_set_uplo( BLIS_LOWER, a );
bli_obj_set_trans( BLIS_TRANSPOSE, a );
bli_copym( &a, &b );
bli_printm( "b after", &b, "%4.1f", "" );
bli_obj_free( &a );
bli_obj_free( &b );
bli_finalize();
return 0;
}
#endif
#if 0
{
num_t dt_a, dt_alpha;
obj_t alpha, a;
dim_t m, n;
dt_alpha = BLIS_DOUBLE;
dt_a = BLIS_DOUBLE;
m = 6;
n = 1;
bli_init();
bli_obj_create( dt_a, m, n, 0, 0, &a );
bli_obj_create( dt_alpha, 1, 1, 0, 0, &alpha );
bli_setm( &BLIS_ONE, &a );
bli_setsc( -(4.0/1.0), &alpha );
bli_printm( "alpha", &alpha, "%4.1f", "" );
bli_printm( "a before", &a, "%4.1f", "" );
//bli_obj_set_struc( BLIS_GENERAL, a );
//bli_obj_set_uplo( BLIS_DENSE, a );
bli_obj_set_diag_offset( -7, a );
bli_obj_set_struc( BLIS_TRIANGULAR, a );
bli_obj_set_uplo( BLIS_UPPER, a );
//bli_obj_set_uplo( BLIS_LOWER, a );
bli_scalm( &alpha, &a );
bli_printm( "a after", &a, "%4.1f", "" );
bli_obj_free( &a );
bli_obj_free( &alpha );
bli_finalize();
return 0;
}
#endif
#if 0
{
num_t dt_beta, dt_a, dt_b;
obj_t beta, a, b;
dim_t m, n;
dt_beta = BLIS_DOUBLE;
dt_a = BLIS_DOUBLE;
dt_b = BLIS_DOUBLE;
m = 6;
n = 11;
bli_init();
bli_obj_create( dt_beta, 1, 1, 0, 0, &beta );
bli_obj_create( dt_a, n, m, 0, 0, &a );
bli_obj_create( dt_b, m, n, 0, 0, &b );
bli_setsc( -(3.0/1.0), &beta );
bli_setm( &BLIS_TWO, &a );
bli_setm( &BLIS_ONE, &b );
bli_printm( "beta", &beta, "%4.1f", "" );
bli_printm( "a before", &a, "%4.1f", "" );
bli_printm( "b before", &b, "%4.1f", "" );
bli_obj_set_struc( BLIS_SYMMETRIC, a );
//bli_obj_set_struc( BLIS_GENERAL, a );
//bli_obj_set_uplo( BLIS_DENSE, a );
bli_obj_set_uplo( BLIS_UPPER, a );
//bli_obj_set_uplo( BLIS_LOWER, a );
bli_obj_set_diag_offset( -2, a );
bli_obj_set_diag( BLIS_UNIT_DIAG, a );
bli_obj_set_trans( BLIS_TRANSPOSE, a );
//bli_obj_set_conj( BLIS_CONJUGATE, a );
bli_scal2m( &beta, &a, &b );
bli_printm( "b after", &b, "%4.1f", "" );
bli_obj_free( &beta );
bli_obj_free( &a );
bli_obj_free( &b );
bli_finalize();
return 0;
}
#endif
#if 0
{
num_t dt_beta, dt_a;
obj_t beta, a, p;
dim_t m, n;
blksz_t* mult_m;
blksz_t* mult_n;
packm_t* cntl;
dt_beta = BLIS_DOUBLE;
dt_a = BLIS_DOUBLE;
m = 6;
n = 11;
bli_init();
bli_obj_create( dt_beta, 1, 1, 0, 0, &beta );
bli_obj_create( dt_a, n, m, 0, 0, &a );
bli_setsc( -(1.0/1.0), &beta );
//bli_setm( &BLIS_TWO, &a );
bli_randm( &a );
bli_printm( "beta", &beta, "%4.1f", "" );
bli_printm( "a before", &a, "%4.1f", "" );
//bli_obj_set_struc( BLIS_SYMMETRIC, a );
bli_obj_set_struc( BLIS_TRIANGULAR, a );
//bli_obj_set_struc( BLIS_GENERAL, a );
//bli_obj_set_uplo( BLIS_DENSE, a );
bli_obj_set_uplo( BLIS_UPPER, a );
//bli_obj_set_uplo( BLIS_LOWER, a );
bli_obj_set_diag_offset( -2, a );
bli_obj_set_diag( BLIS_NONUNIT_DIAG, a );
bli_obj_set_trans( BLIS_TRANSPOSE, a );
//bli_obj_set_conj( BLIS_CONJUGATE, a );
mult_m = bli_blksz_obj_create( 2, 4, 2, 2 );
mult_n = bli_blksz_obj_create( 1, 2, 1, 1 );
cntl = bli_packm_cntl_obj_create( BLIS_UNBLOCKED,
BLIS_VARIANT1,
mult_m,
mult_n,
//TRUE, // scale?
FALSE, // scale?
//FALSE, // densify?
TRUE, // densify?
BLIS_PACKED_ROWS );
//BLIS_PACKED_COLUMNS );
bli_obj_init_pack( &p );
bli_packm_int( &beta, &a, &p, cntl );
bli_printm( "p after", &p, "%4.1f", "" );
p.m = p.pack_mem->m;
p.n = p.pack_mem->n;
bli_printm( "p mem", &p, "%4.1f", "" );
bli_obj_release_pack( &p );
bli_obj_free( &beta );
bli_obj_free( &a );
bli_blksz_obj_free( mult_m );
bli_blksz_obj_free( mult_n );
bli_cntl_obj_free( cntl );
bli_finalize();
return 0;
}
#endif
#if 0
{
num_t dt_beta, dt_a;
obj_t beta, a, p;
dim_t m, n;
blksz_t* mult_m;
blksz_t* mult_n;
packm_t* cntl;
dt_beta = BLIS_DOUBLE;
dt_a = BLIS_DOUBLE;
m = 7;
n = 11;
bli_init();
bli_obj_create( dt_beta, 1, 1, 0, 0, &beta );
bli_obj_create( dt_a, m, n, 0, 0, &a );
bli_setsc( -(1.0/1.0), &beta );
//bli_setm( &BLIS_TWO, &a );
bli_randm( &a );
bli_printm( "beta", &beta, "%4.1f", "" );
bli_printm( "a before", &a, "%4.1f", "" );
//bli_obj_set_struc( BLIS_SYMMETRIC, a );
//bli_obj_set_struc( BLIS_TRIANGULAR, a );
bli_obj_set_struc( BLIS_GENERAL, a );
bli_obj_set_uplo( BLIS_DENSE, a );
//bli_obj_set_uplo( BLIS_UPPER, a );
//bli_obj_set_uplo( BLIS_LOWER, a );
//bli_obj_set_diag_offset( -2, a );
bli_obj_set_diag_offset( 0, a );
//bli_obj_set_diag( BLIS_UNIT_DIAG, a );
//bli_obj_set_trans( BLIS_TRANSPOSE, a );
//bli_obj_set_conj( BLIS_CONJUGATE, a );
mult_m = bli_blksz_obj_create( 2, 4, 2, 2 );
mult_n = bli_blksz_obj_create( 1, 2, 1, 1 );
cntl = bli_packm_cntl_obj_create( BLIS_BLOCKED,
BLIS_VARIANT1,
mult_m,
mult_n,
//TRUE, // scale?
FALSE, // scale?
//FALSE, // densify?
TRUE, // densify?
BLIS_PACKED_COL_PANELS );
bli_obj_init_pack( &p );
bli_packm_int( &beta, &a, &p, cntl );
p.m = p.pack_mem->m;
p.n = p.pack_mem->n;
p.rs = 1;
p.cs = p.pack_mem->m;
bli_printm( "p mem", &p, "%4.1f", "" );
bli_obj_release_pack( &p );
bli_obj_free( &beta );
bli_obj_free( &a );
bli_blksz_obj_free( mult_m );
bli_blksz_obj_free( mult_n );
bli_cntl_obj_free( cntl );
bli_finalize();
return 0;
}
#endif
// -- Level 2 ------------------------------------------------------------------
#if 0
{
obj_t alpha, beta, a, x, y, xK, yK;
num_t dt_alpha, dt_beta, dt_a, dt_x, dt_y;
dim_t m, n, K;
//dt_alpha = BLIS_DCOMPLEX;
//dt_beta = BLIS_DCOMPLEX;
//dt_a = BLIS_DCOMPLEX;
//dt_x = BLIS_DCOMPLEX;
//dt_y = BLIS_DCOMPLEX;
dt_alpha = BLIS_DOUBLE;
dt_beta = BLIS_DOUBLE;
dt_a = BLIS_DOUBLE;
dt_x = BLIS_DOUBLE;
dt_y = BLIS_DOUBLE;
m = 11;
n = 13;
K = 4;
bli_init();
//bli_obj_create( dt_a, m, n, n, 1, &a );
//bli_obj_create( dt_a, m, n, 0, 0, &a );
bli_obj_create( dt_a, m, n, 0, 0, &a );
bli_obj_create( dt_x, K, n, 0, 0, &xK );
bli_obj_create( dt_y, K, m, 0, 0, &yK );
bli_acquire_mpart_t2b( BLIS_SUBPART1, 0, 1, &xK, &x );
bli_acquire_mpart_t2b( BLIS_SUBPART1, 0, 1, &yK, &y );
bli_obj_create( dt_alpha, 1, 1, 0, 0, &alpha );
bli_obj_create( dt_beta, 1, 1, 0, 0, &beta );
bli_setm( &BLIS_ZERO, &a );
bli_setm( &BLIS_ZERO, &x );
bli_setm( &BLIS_ZERO, &y );
bli_randm( &a );
bli_randm( &x );
bli_randm( &y );
bli_setsc( (2.0/1.0), &alpha );
bli_setsc( -(1.0/1.0), &beta );
//bli_setsc( (1.0/1.0), &alpha );
//bli_setsc( (0.0/1.0), &beta );
bli_printm( "a", &a, "%4.1f", "" );
bli_printm( "x", &x, "%4.1f", "" );
bli_printm( "y", &y, "%4.1f", "" );
//bli_obj_set_conjtrans( BLIS_CONJ_NO_TRANSPOSE, a );
//bli_obj_set_conj( BLIS_CONJUGATE, x );
bli_gemv( &alpha, &a, &x, &beta, &y );
//bli_gemv_unb_var1( &alpha, &a, &x, &beta, &y, NULL );
//bli_gemv_unf_var1( &alpha, &a, &x, &beta, &y, NULL );
//bli_gemv_unb_var2( &alpha, &a, &x, &beta, &y, NULL );
//bli_gemv_unf_var2( &alpha, &a, &x, &beta, &y, NULL );
bli_printm( "y after", &y, "%4.1f", "" );
bli_obj_free( &a );
bli_obj_free( &xK );
bli_obj_free( &yK );
bli_obj_free( &alpha );
bli_obj_free( &beta );
bli_finalize();
return 0;
}
#endif
#if 0
{
obj_t alpha, a, x, y, xK, yK;
num_t dt_alpha, dt_a, dt_x, dt_y;
dim_t m, n, K;
//dt_alpha = BLIS_DCOMPLEX;
//dt_a = BLIS_DCOMPLEX;
//dt_x = BLIS_DCOMPLEX;
//dt_y = BLIS_DCOMPLEX;
dt_alpha = BLIS_DOUBLE;
dt_a = BLIS_DOUBLE;
dt_x = BLIS_DOUBLE;
dt_y = BLIS_DOUBLE;
m = 11;
n = 13;
K = 4;
bli_init();
bli_obj_create( dt_a, m, n, 0, 0, &a );
//bli_obj_create( dt_a, m, n, n, 1, &a );
bli_obj_create( dt_x, K, m, 0, 0, &xK );
bli_obj_create( dt_y, K, n, 0, 0, &yK );
bli_acquire_mpart_t2b( BLIS_SUBPART1, 0, 1, &xK, &x );
bli_acquire_mpart_t2b( BLIS_SUBPART1, 0, 1, &yK, &y );
bli_obj_create( dt_alpha, 1, 1, 0, 0, &alpha );
bli_setm( &BLIS_ZERO, &a );
bli_setm( &BLIS_ZERO, &x );
bli_setm( &BLIS_ZERO, &y );
bli_randm( &a );
bli_randm( &x );
bli_randm( &y );
bli_setsc( (2.0/1.0), &alpha );
bli_printm( "a", &a, "%4.1f", "" );
bli_printm( "x", &x, "%4.1f", "" );
bli_printm( "y", &y, "%4.1f", "" );
bli_obj_set_conj( BLIS_NO_CONJUGATE, x );
bli_obj_set_conj( BLIS_NO_CONJUGATE, y );
bli_ger( &alpha, &x, &y, &a );
//bli_ger_unb_var1( &alpha, &x, &y, &a, NULL );
//bli_ger_unb_var2( &alpha, &x, &y, &a, NULL );
bli_printm( "a after", &a, "%4.1f", "" );
bli_obj_free( &a );
bli_obj_free( &xK );
bli_obj_free( &yK );
bli_obj_free( &alpha );
bli_finalize();
return 0;
}
#endif
#if 0
{
obj_t alpha, beta, a, x, y, xK, yK;
num_t dt_alpha, dt_beta, dt_a, dt_x, dt_y;
dim_t m, K;
//dt_alpha = BLIS_DCOMPLEX;
//dt_beta = BLIS_DCOMPLEX;
//dt_a = BLIS_DCOMPLEX;
//dt_x = BLIS_DCOMPLEX;
//dt_y = BLIS_DCOMPLEX;
dt_alpha = BLIS_DOUBLE;
dt_beta = BLIS_DOUBLE;
dt_a = BLIS_DOUBLE;
dt_x = BLIS_DOUBLE;
dt_y = BLIS_DOUBLE;
m = 11;
K = 4;
bli_init();
bli_obj_create( dt_a, m, m, 0, 0, &a );
//bli_obj_create( dt_a, m, m, m, 1, &a );
bli_obj_create( dt_x, K, m, 0, 0, &xK );
bli_obj_create( dt_y, K, m, 0, 0, &yK );
bli_acquire_mpart_t2b( BLIS_SUBPART1, 0, 1, &xK, &x );
bli_acquire_mpart_t2b( BLIS_SUBPART1, 0, 1, &yK, &y );
bli_obj_create( dt_alpha, 1, 1, 0, 0, &alpha );
bli_obj_create( dt_beta, 1, 1, 0, 0, &beta );
bli_setm( &BLIS_ZERO, &a );
bli_setm( &BLIS_ZERO, &x );
bli_setm( &BLIS_ZERO, &y );
bli_randm( &a );
bli_randm( &x );
bli_randm( &y );
bli_setsc( (2.0/1.0), &alpha );
bli_setsc( -(1.0/1.0), &beta );
//bli_setsc( (1.0/1.0), &alpha );
//bli_setsc( (1.0/1.0), &beta );
bli_printm( "a", &a, "%4.1f", "" );
bli_printm( "x", &x, "%4.1f", "" );
bli_printm( "y", &y, "%4.1f", "" );
//bli_obj_set_uplo( BLIS_LOWER, a );
bli_obj_set_uplo( BLIS_UPPER, a );
bli_obj_set_conj( BLIS_NO_CONJUGATE, a );
//bli_obj_set_conj( BLIS_CONJUGATE, x );
bli_hemv( &alpha, &a, &x, &beta, &y );
//bli_hemv_unb_var1( &alpha, &a, &x, &beta, &y, NULL );
//bli_hemv_unb_var2( &alpha, &a, &x, &beta, &y, NULL );
//bli_hemv_unb_var3( &alpha, &a, &x, &beta, &y, NULL );
//bli_hemv_unb_var4( &alpha, &a, &x, &beta, &y, NULL );
//bli_hemv_unf_var1a( &alpha, &a, &x, &beta, &y, NULL );
//bli_hemv_unf_var1( &alpha, &a, &x, &beta, &y, NULL );
//bli_hemv_unf_var3a( &alpha, &a, &x, &beta, &y, NULL );
//bli_hemv_unf_var3( &alpha, &a, &x, &beta, &y, NULL );
bli_printm( "y after", &y, "%4.1f", "" );
bli_obj_free( &a );
bli_obj_free( &xK );
bli_obj_free( &yK );
bli_obj_free( &alpha );
bli_obj_free( &beta );
bli_finalize();
return 0;
}
#endif
#if 0
{
obj_t alpha, a, x, xK;
num_t dt_alpha, dt_a, dt_x;
dim_t m, K;
//dt_x = BLIS_DCOMPLEX;
//dt_a = BLIS_DCOMPLEX;
dt_x = BLIS_DOUBLE;
dt_a = BLIS_DOUBLE;
dt_alpha = bli_datatype_proj_to_real( dt_x );
m = 13;
K = 4;
bli_init();
bli_obj_create( dt_x, K, m, 0, 0, &xK );
//bli_obj_create( dt_a, m, m, 0, 0, &a );
bli_obj_create( dt_a, m, m, m, 1, &a );
bli_acquire_mpart_t2b( BLIS_SUBPART1, 0, 1, &xK, &x );
bli_obj_create( dt_alpha, 1, 1, 0, 0, &alpha );
bli_setm( &BLIS_ZERO, &x );
bli_setm( &BLIS_ZERO, &a );
bli_randm( &x );
bli_randm( &a );
bli_setsc( (2.0/1.0), &alpha );
//bli_setsc( (1.0/1.0), &alpha );
bli_printm( "x", &x, "%4.1f", "" );
bli_printm( "a", &a, "%4.1f", "" );
//bli_obj_set_uplo( BLIS_LOWER, a );
bli_obj_set_uplo( BLIS_UPPER, a );
//bli_obj_set_conj( BLIS_CONJUGATE, x );
bli_her( &alpha, &x, &a );
//bli_her_unb_var1( &alpha, &x, &a, NULL );
//bli_her_unb_var2( &alpha, &x, &a, NULL );
bli_printm( "a after", &a, "%4.1f", "" );
bli_obj_free( &xK );
bli_obj_free( &a );
bli_obj_free( &alpha );
bli_finalize();
return 0;
}
#endif
#if 0
{
obj_t alpha, a, x, y, xK, yK;
num_t dt_alpha, dt_a, dt_x, dt_y;
dim_t m, K;
//dt_x = BLIS_DCOMPLEX;
//dt_y = BLIS_DCOMPLEX;
//dt_a = BLIS_DCOMPLEX;
//dt_alpha = BLIS_DCOMPLEX;
dt_x = BLIS_DOUBLE;
dt_y = BLIS_DOUBLE;
dt_a = BLIS_DOUBLE;
dt_alpha = BLIS_DOUBLE;
m = 13;
K = 4;
bli_init();
bli_obj_create( dt_x, K, m, 0, 0, &xK );
bli_obj_create( dt_y, K, m, 0, 0, &yK );
//bli_obj_create( dt_a, m, m, m, 1, &a );
bli_obj_create( dt_a, m, m, 0, 0, &a );
bli_acquire_mpart_t2b( BLIS_SUBPART1, 0, 1, &xK, &x );
bli_acquire_mpart_t2b( BLIS_SUBPART1, 0, 1, &yK, &y );
bli_obj_create( dt_alpha, 1, 1, 0, 0, &alpha );
bli_setm( &BLIS_ZERO, &x );
bli_setm( &BLIS_ZERO, &y );
bli_setm( &BLIS_ZERO, &a );
bli_randm( &x );
bli_randm( &y );
bli_randm( &a );
//bli_randv( &alpha );
bli_setsc( (2.0/1.0), &alpha );
//bli_setsc( (1.0/1.0), &alpha );
bli_printm( "x", &x, "%4.1f", "" );
bli_printm( "y", &y, "%4.1f", "" );
bli_printm( "a", &a, "%4.1f", "" );
//bli_printv( "alpha", &alpha, "%4.1f", "" );
//bli_obj_set_uplo( BLIS_LOWER, a );
bli_obj_set_uplo( BLIS_UPPER, a );
//bli_obj_set_conj( BLIS_CONJUGATE, x );
bli_her2( &alpha, &x, &y, &a );
//bli_her2_unb_var1( &alpha, &x, &y, &a, NULL );
//bli_her2_unb_var2( &alpha, &x, &y, &a, NULL );
//bli_her2_unb_var3( &alpha, &x, &y, &a, NULL );
//bli_her2_unb_var4( &alpha, &x, &y, &a, NULL );
//bli_her2_unf_var4( &alpha, &x, &y, &a, NULL );
bli_printm( "a after", &a, "%4.1f", "" );
bli_obj_free( &xK );
bli_obj_free( &yK );
bli_obj_free( &a );
bli_obj_free( &alpha );
bli_finalize();
return 0;
}
#endif
#if 0
{
obj_t alpha, a, x, xK;
num_t dt_alpha, dt_a, dt_x;
dim_t m, K;
//dt_alpha = BLIS_DCOMPLEX;
//dt_a = BLIS_DCOMPLEX;
//dt_x = BLIS_DCOMPLEX;
dt_alpha = BLIS_DOUBLE;
dt_a = BLIS_DOUBLE;
dt_x = BLIS_DOUBLE;
m = 13;
K = 4;
bli_init();
//bli_obj_create( dt_a, m, m, 0, 0, &a );
bli_obj_create( dt_a, m, m, m, 1, &a );
bli_obj_create( dt_x, K, m, 0, 0, &xK );
bli_obj_create( dt_x, 1, 1, 0, 0, &alpha );
bli_acquire_mpart_t2b( BLIS_SUBPART1, 0, 1, &xK, &x );
bli_setm( &BLIS_ZERO, &a );
bli_setm( &BLIS_ZERO, &x );
bli_randm( &a );
bli_randm( &x );
bli_setsc( (1.0/1.0), &alpha );
bli_obj_set_uplo( BLIS_UPPER, a );
//bli_obj_set_uplo( BLIS_LOWER, a );
bli_obj_set_conjtrans( BLIS_NO_TRANSPOSE, a );
//bli_obj_set_conjtrans( BLIS_TRANSPOSE, a );
bli_obj_set_diag( BLIS_NONUNIT_DIAG, a );
bli_printm( "a", &a, "%4.1f", "" );
bli_printm( "x", &x, "%4.1f", "" );
bli_trmv( &alpha, &a, &x );
//bli_trmv_unb_var2( &alpha, &a, &x, NULL );
//bli_trmv_unb_var1( &alpha, &a, &x, NULL );
//bli_trmv_unf_var1( &alpha, &a, &x, NULL );
//bli_trmv_unf_var2( &alpha, &a, &x, NULL );
bli_printm( "x after", &x, "%4.1f", "" );
bli_obj_free( &alpha );
bli_obj_free( &a );
bli_obj_free( &xK );
bli_finalize();
return 0;
}
#endif
#if 0
{
obj_t alpha, a, x, xK;
num_t dt_alpha, dt_a, dt_x;
dim_t m, K;
dt_alpha = BLIS_DOUBLE;
dt_a = BLIS_DOUBLE;
dt_x = BLIS_DOUBLE;
m = 13;
K = 4;
bli_init();
//bli_obj_create( dt_a, m, m, 0, 0, &a );
bli_obj_create( dt_a, m, m, m, 1, &a );
bli_obj_create( dt_x, K, m, 0, 0, &xK );
bli_obj_create( dt_x, 1, 1, 0, 0, &alpha );
bli_acquire_mpart_t2b( BLIS_SUBPART1, 0, 1, &xK, &x );
bli_setm( &BLIS_ZERO, &a );
bli_setm( &BLIS_ZERO, &x );
bli_randm( &a );
bli_randm( &x );
bli_setsc( (2.0/1.0), &alpha );
bli_obj_set_uplo( BLIS_UPPER, a );
//bli_obj_set_uplo( BLIS_LOWER, a );
bli_obj_set_conjtrans( BLIS_NO_TRANSPOSE, a );
bli_obj_set_diag( BLIS_NONUNIT_DIAG, a );
bli_printm( "a", &a, "%4.1f", "" );
bli_printm( "x", &x, "%4.1f", "" );
bli_trsv( &alpha, &a, &x );
//bli_trsv_unb_var1( &alpha, &a, &x, NULL );
//bli_trsv_unb_var2( &alpha, &a, &x, NULL );
//bli_trsv_unf_var1( &alpha, &a, &x, NULL );
//bli_trsv_unf_var2( &alpha, &a, &x, NULL );
bli_printm( "x after", &x, "%4.1f", "" );
bli_obj_free( &alpha );
bli_obj_free( &a );
bli_obj_free( &xK );
bli_finalize();
return 0;
}
#endif
// -- Level 3 ------------------------------------------------------------------
#if 0
{
obj_t a, b;
obj_t alpha;
num_t dt_a, dt_b;
num_t dt_al;
dim_t m, n;
dt_a = BLIS_DCOMPLEX;
dt_b = BLIS_DCOMPLEX;
dt_al = BLIS_DCOMPLEX;
m = 9;
n = 5;
bli_init();
bli_obj_create( dt_a, m, m, 0, 0, &a );
//bli_obj_create( dt_b, m, n, 0, 0, &b );
bli_obj_create( dt_b, n, m, 0, 0, &b );
bli_obj_create( dt_al, 1, 1, 0, 0, &alpha );
bli_setm( &BLIS_ZERO, &a );
bli_setm( &BLIS_ZERO, &b );
bli_randm( &a );
bli_randm( &b );
//bli_setsc( (2.0/1.0), &alpha );
bli_setsc( (1.0/1.0), &alpha );
bli_obj_set_uplo( BLIS_LOWER, a );
//bli_obj_set_uplo( BLIS_UPPER, a );
bli_obj_set_conjtrans( BLIS_TRANSPOSE, a );
bli_obj_set_diag( BLIS_NONUNIT_DIAG, a );
bli_printm( "a", &a, "%4.1f", "" );
bli_printm( "b", &b, "%4.1f", "" );
//bli_trmm_lu_unb_var3( &alpha, &a, &b, NULL );
//bli_trmm_lu_unb_var2( &alpha, &a, &b, NULL );
//bli_trmm_lu_unb_var1( &alpha, &a, &b, NULL );
//bli_trmm_ll_unb_var3( &alpha, &a, &b, NULL );
//bli_trmm_ll_unb_var2( &alpha, &a, &b, NULL );
//bli_trmm_ll_unb_var1( &alpha, &a, &b, NULL );
//bli_trmm( BLIS_LEFT, &alpha, &a, &b );
bli_trmm( BLIS_RIGHT, &alpha, &a, &b );
bli_printm( "b after", &b, "%4.1f", "" );
bli_obj_free( &a );
bli_obj_free( &b );
bli_finalize();
return 0;
}
#endif
}
Reference in New Issue View Git Blame Copy Permalink