[dense-linear-algebra-libraries/linalg.git] / src / ti / linalg / blis / frame / 3 / her2k / attic / bli_her2k_int.c
diff --git a/src/ti/linalg/blis/frame/3/her2k/attic/bli_her2k_int.c b/src/ti/linalg/blis/frame/3/her2k/attic/bli_her2k_int.c
--- /dev/null
@@ -0,0 +1,163 @@
+/*
+
+ BLIS
+ An object-based framework for developing high-performance BLAS-like
+ libraries.
+
+ Copyright (C) 2014, The University of Texas at Austin
+
+ 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 at Austin 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 "blis.h"
+
+#define FUNCPTR_T her2k_fp
+
+typedef void (*FUNCPTR_T)( obj_t* a,
+ obj_t* bh,
+ obj_t* b,
+ obj_t* ah,
+ obj_t* c,
+ her2k_t* cntl );
+
+static FUNCPTR_T vars[2][4][3] =
+{
+ // lower
+ {
+ // unblocked optimized unblocked blocked
+ { NULL, NULL, bli_her2k_blk_var1f },
+ { NULL, bli_her2k_l_ker_var2, bli_her2k_blk_var2f },
+ { NULL, NULL, bli_her2k_blk_var3f },
+ { NULL, NULL, NULL },
+ },
+ // upper
+ {
+ // unblocked optimized unblocked blocked
+ { NULL, NULL, bli_her2k_blk_var1f },
+ { NULL, bli_her2k_u_ker_var2, bli_her2k_blk_var2f },
+ { NULL, NULL, bli_her2k_blk_var3f },
+ { NULL, NULL, NULL },
+ }
+};
+
+void bli_her2k_int( obj_t* alpha_abh,
+ obj_t* a,
+ obj_t* bh,
+ obj_t* alpha_bah,
+ obj_t* b,
+ obj_t* ah,
+ obj_t* beta,
+ obj_t* c,
+ her2k_t* cntl )
+{
+ obj_t a_local;
+ obj_t bh_local;
+ obj_t b_local;
+ obj_t ah_local;
+ obj_t c_local;
+ varnum_t n;
+ impl_t i;
+ bool_t uplo;
+ FUNCPTR_T f;
+
+ // Check parameters.
+ if ( bli_error_checking_is_enabled() )
+ bli_her2k_int_check( alpha_abh, a, bh, alpha_bah, b, ah, beta, c, cntl );
+
+ // If C has a zero dimension, return early.
+ if ( bli_obj_has_zero_dim( *c ) ) return;
+
+ // If A or B has a zero dimension, scale C by beta and return early.
+ if ( bli_obj_has_zero_dim( *a ) ||
+ bli_obj_has_zero_dim( *ah ) ||
+ bli_obj_has_zero_dim( *b ) ||
+ bli_obj_has_zero_dim( *bh ) )
+ {
+ bli_scalm( beta, c );
+ return;
+ }
+
+ // Alias A, B', B, and A' in case we need to update attached scalars.
+ bli_obj_alias_to( *a, a_local );
+ bli_obj_alias_to( *bh, bh_local );
+ bli_obj_alias_to( *b, b_local );
+ bli_obj_alias_to( *ah, ah_local );
+
+ // Alias C in case we need to induce a transposition.
+ bli_obj_alias_to( *c, c_local );
+
+ // If we are about to call a leaf-level implementation, and matrix C
+ // still needs a transposition, then we must induce one by swapping the
+ // strides and dimensions. Note that this transposition would normally
+ // be handled explicitly in the packing of C, but if C is not being
+ // packed, this is our last chance to handle the transposition.
+ if ( cntl_is_leaf( cntl ) && bli_obj_has_trans( *c ) )
+ {
+ bli_obj_induce_trans( c_local );
+ bli_obj_set_onlytrans( BLIS_NO_TRANSPOSE, c_local );
+ }
+
+ // If alpha_abh is non-unit, typecast and apply it to the scalar
+ // attached to B'.
+ if ( !bli_obj_equals( alpha_abh, &BLIS_ONE ) )
+ {
+ bli_obj_scalar_apply_scalar( alpha_abh, &bh_local );
+ }
+
+ // If alpha_bah is non-unit, typecast and apply it to the scalar
+ // attached to A'.
+ if ( !bli_obj_equals( alpha_bah, &BLIS_ONE ) )
+ {
+ bli_obj_scalar_apply_scalar( alpha_bah, &ah_local );
+ }
+
+ // If beta is non-unit, typecast and apply it to the scalar
+ // attached to C.
+ if ( !bli_obj_equals( beta, &BLIS_ONE ) )
+ {
+ bli_obj_scalar_apply_scalar( beta, &c_local );
+ }
+
+ // Set a bool based on the uplo field of c.
+ if ( bli_obj_root_is_lower( c_local ) ) uplo = 0;
+ else uplo = 1;
+
+ // Extract the variant number and implementation type.
+ n = cntl_var_num( cntl );
+ i = cntl_impl_type( cntl );
+
+ // Index into the variant array to extract the correct function pointer.
+ f = vars[uplo][n][i];
+
+ // Invoke the variant.
+ f( &a_local,
+ &bh_local,
+ &b_local,
+ &ah_local,
+ &c_local,
+ cntl );
+}
+