vec_kern_unr_pref.h

Go to the documentation of this file.

#ifndef TBCI_VEC_KERN_UNR_PREF_H
#define TBCI_VEC_KERN_UNR_PREF_H  
#include "basics.h"

NAMESPACE_TBCI

#define ADD3(r,v1,v2,f1,f2)     r = v1 + v2
VKERN_TEMPL_3V(do_vec_vec_add, ADD3)

#define SUB3(r,v1,v2,f1,f2)     r = v1 - v2
VKERN_TEMPL_3V(do_vec_vec_sub, SUB3)

#define MUL3(r,v1,v2,f1,f2)     r = v1 * v2
VKERN_TEMPL_3V(do_vec_vec_mul, MUL3)

#define CMUL3(r,v1,v2,f1,f2)    r = CPLX__ conj(v1) * v2
VKERN_TEMPL_3V(do_vec_vec_cmul, CMUL3)

#define DIV3(r,v1,v2,f1,f2)     r = v1 / v2
VKERN_TEMPL_3V(do_vec_vec_div, DIV3)

#define CDIV3(r,v1,v2,f1,f2)    r = CPLX__ conj(v1) / v2
VKERN_TEMPL_3V(do_vec_vec_cdiv, CDIV3)


#define ADD2(r,v1,f1,f2)        r += v1
VKERN_TEMPL_2V(do_vec_add_vec, ADD2)

#define SUB2(r,v1,f1,f2)        r -= v1
VKERN_TEMPL_2V(do_vec_sub_vec, SUB2)

#define SUB2I(r,v1,f1,f2)       r  = v1 - r
VKERN_TEMPL_2V(do_vec_sub_vec_inv, SUB2I)

#define MUL2(r,v1,f1,f2)        r *= v1
VKERN_TEMPL_2V(do_vec_mul_vec, MUL2)

#define CMUL2(r,v1,f1,f2)       r  = CPLX__ conj(r) * v1
VKERN_TEMPL_2V(do_vec_cmul_vec, CMUL2)

#define CMUL2I(r,v1,f1,f2)      r *= CPLX__ conj(v1)
VKERN_TEMPL_2V(do_vec_cmul_vec_inv, CMUL2I)

#define DIV2(r,v1,f1,f2)        r /= v1
VKERN_TEMPL_2V(do_vec_div_vec, DIV2)

#define DIV2I(r,v1,f1,f2)       r = v1 / r
VKERN_TEMPL_2V(do_vec_div_vec_inv, DIV2I)

#define CDIV2(r,v1,f1,f2)       r  = CPLX__ conj(r) / v1
VKERN_TEMPL_2V(do_vec_cdiv_vec, CDIV2)

#define CDIV2I(r,v1,f1,f2)      r  = CPLX__ conj(v1) / r
VKERN_TEMPL_2V(do_vec_cdiv_vec_inv, CDIV2I)



#define ADD2NV(r,v1,f1,f2)      r = v1 + f2
VKERN_TEMPL_2V_C(do_vec_val_add, ADD2NV)

#define SUB2NV(r,v1,f1,f2)      r = v1 - f2
VKERN_TEMPL_2V_C(do_vec_val_sub, SUB2NV)

#define MUL2NV(r,v1,f1,f2)      r = v1 * f2
VKERN_TEMPL_2V_C(do_vec_val_mul, MUL2NV)


#define ADD2RV(r,v1,f1,f2)      r = f2 + v1
VKERN_TEMPL_2V_C(do_val_vec_add, ADD2RV)

#define SUB2RV(r,v1,f1,f2)      r = f2 - v1
VKERN_TEMPL_2V_C(do_val_vec_sub, SUB2RV)

#define MUL2RV(r,v1,f1,f2)      r = f2 * v1
VKERN_TEMPL_2V_C(do_val_vec_mul, MUL2RV)

#define DIV2RV(r,v1,f1,f2)      r = f2 / v1
VKERN_TEMPL_2V_C(do_val_vec_div, DIV2RV)


#define ADD1NV(r,f1,f2)         r += f2
VKERN_TEMPL_1V_C(do_vec_add_val, ADD1NV)

#define SUB1NV(r,f1,f2)         r -= f2
VKERN_TEMPL_1V_C(do_vec_sub_val, SUB1NV)

#define SUB1RV(r,f1,f2)         r  = f2 - r
VKERN_TEMPL_1V_C(do_val_sub_vec, SUB1RV)

#define MUL1NV(r,f1,f2)         r *= f2
VKERN_TEMPL_1V_C(do_vec_mul_val, MUL1NV)

#define DIV1NV(r,f1,f2)         r /= f2
VKERN_TEMPL_1V_C(do_vec_div_val, DIV1NV)

#define DIV1RV(r,f1,f2)         r = f2 / r
VKERN_TEMPL_1V_C(do_val_div_vec, DIV1RV)

#define ADD1RV(r,f1,f2)         r  = f2 + r
VKERN_TEMPL_1V_C(do_val_add_vec, ADD1RV)

// unused ...
#define MUL1RV(r,f1,f2)         r  = f2 * r;


/* TSVector stuff */

#define ADD2NS(r,v1,f1,f2)      r += f2*v1
VKERN_TEMPL_2V_C(do_vec_add_svc, ADD2NS)

#define SUB2NS(r,v1,f1,f2)      r -= f2*v1
VKERN_TEMPL_2V_C(do_vec_sub_svc, SUB2NS)

#define SUB2RS(r,v1,f1,f2)      r  = f2*v1 - r
VKERN_TEMPL_2V_C(do_vec_sub_svc_inv, SUB2RS)


#define ADD3NS(r,v1,v2,f1,f2)   r = v1 + f2*v2
VKERN_TEMPL_3V_C(do_vec_svc_add, ADD3NS)

#define SUB3NS(r,v1,v2,f1,f2)   r = v1 - f2*v2
VKERN_TEMPL_3V_C(do_vec_svc_sub, SUB3NS)


#define ADD3SN(r,v1,v2,f1,f2)   r = f2*v1 + v2
VKERN_TEMPL_3V_C(do_svc_vec_add, ADD3SN)

#define SUB3SN(r,v1,v2,f1,f2)   r = f2*v1 - v2
VKERN_TEMPL_3V_C(do_svc_vec_sub, SUB3SN)


#define ADD3SS(r,v1,v2,f1,f2)   r = f1*v1 + f2*v2
VKERN_TEMPL_3V_CC(do_svc_svc_add, ADD3SS)

#define SUB3SS(r,v1,v2,f1,f2)   r = f1*v1 - f2*v2
VKERN_TEMPL_3V_CC(do_svc_svc_sub, SUB3SS)


#define ADD2SN(r,v1,f1,f2)      r  = f2*r + v1
VKERN_TEMPL_2V_C(do_svc_add_vec, ADD2SN)

#define SUB2SN(r,v1,f1,f2)      r  = f2*r - v1
VKERN_TEMPL_2V_C(do_svc_sub_vec, SUB2SN)


#define ADD2SS(r,v1,f1,f2)      r  = f1*r + f2*v1
VKERN_TEMPL_2V_CC(do_svc_add_svc, ADD2SS)

#define SUB2SS(r,v1,f1,f2)      r  = f1*r - f2*v1
VKERN_TEMPL_2V_CC(do_svc_sub_svc, SUB2SS)


#define ADD2SV(r,v1,f1,f2)      r = f1*v1 + f2
VKERN_TEMPL_2V_CC(do_svc_val_add, ADD2SV)

#define SUB2SV(r,v1,f1,f2)      r = f1*v1 - f2
VKERN_TEMPL_2V_CC(do_svc_val_sub, SUB2SV)


#define ADD1SV(r,f1,f2)         r = f1*r + f2
VKERN_TEMPL_1V_CC(do_svc_add_val, ADD1SV)

#define SUB1SV(r,f1,f2)         r = f1*r - f2
VKERN_TEMPL_1V_CC(do_svc_sub_val, SUB1SV)


#define ADD2VS(r,v1,f1,f2)      r = f1 + f2*v1
VKERN_TEMPL_2V_CC(do_val_svc_add, ADD2VS)

#define SUB2VS(r,v1,f1,f2)      r = f1 - f2*v1
VKERN_TEMPL_2V_CC(do_val_svc_sub, SUB2VS)

#define DIV2VS(r,v1,f1,f2)      r = f1 / (f2*v1)
VKERN_TEMPL_2V_CC(do_val_svc_div, DIV2VS)


#define NEG2(r,v1,f1,f2)        r = -v1
VKERN_TEMPL_2V(do_vec_neg_vec, NEG2)

#define NEG1(r,f1,f2)           r = -r
VKERN_TEMPL_1V(do_vec_neg, NEG1)


#define DOT2(r,v1,f1,f2)        f2 += CPLX__ conj(r) * v1
VKERN_TEMPL_2V_T(do_vec_dot, DOT2, T)

#define MULT2(r,v1,f1,f2)       f2 += r * v1
VKERN_TEMPL_2V_T(do_vec_mult, MULT2, T)


#define FABS1(r,f1,f2)          f2 += fabssqr(r)
VKERN_TEMPL_1V_T(do_vec_fabssqr, FABS1, LONG_DOUBLE)

#define SQR1(r,f1,f2)           f2 += r*r
VKERN_TEMPL_1V_T(do_vec_sumsqr, SQR1, T)

#define SUM1(r,f1,f2)           f2 += r
VKERN_TEMPL_1V_T(do_vec_sum, SUM1, T)


#define COMP2(r,v1,f1,f2)       if (r != v1) ++f2
VKERN_TEMPL_2V_T(do_vv_comp, COMP2, volatile long)


// Used in do_bdmat_vec_mult
#define SUMMULT3(r,v1,v2,f1,f2) r += v1*v2
VKERN_TEMPL_3V(do_add_vec_vec_mul, SUMMULT3)
#define SUMCMULT3(r,v1,v2,f1,f2)        r += CPLX__ conj(v1)*v2
VKERN_TEMPL_3V(do_add_vec_vec_cmul, SUMCMULT3)


NAMESPACE_END

#endif  /* TBCI_VEC_KERN_UNR_PREF_H */

---