void bli_zshiftdiag( conj_t conj, int offset, int m, int n, dcomplex* sigma, dcomplex* a, int a_rs, int a_cs ) { dcomplex* alpha; dcomplex sigma_conj; int i, j; bli_zcopys( 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_z_div(bla_dcomplex *cp, const bla_dcomplex *ap, const bla_dcomplex *bp) { bli_zcopys( *ap, *cp ); bli_zinvscals( *bp, *cp ); }
void bla_z_div(doublecomplex *cp, doublecomplex *ap, doublecomplex *bp) { bli_zcopys( *ap, *cp ); bli_zinvscals( *bp, *cp ); }