void bli_cshiftdiag( conj_t conj, int offset, int m, int n, scomplex* sigma, scomplex* a, int a_rs, int a_cs ) { scomplex* alpha; scomplex sigma_conj; int i, j; bli_ccopys( conj, sigma, &sigma_conj ); i = j = 0; if ( offset < 0 ) i = -offset; else if ( offset > 0 ) j = offset; while ( i < m && j < n ) { alpha = a + i*a_rs + j*a_cs; alpha->real += sigma_conj.real; alpha->imag += sigma_conj.imag; ++i; ++j; } }
void bla_c_div(bla_scomplex *cp, const bla_scomplex *ap, const bla_scomplex *bp) { bli_ccopys( *ap, *cp ); bli_cinvscals( *bp, *cp ); }
void bla_c_div(singlecomplex *cp, singlecomplex *ap, singlecomplex *bp) { bli_ccopys( *ap, *cp ); bli_cinvscals( *bp, *cp ); }