54 #define max(a,b) ( (a) > (b) ? (a) : (b) )
55 #define min(a,b) ( (a) < (b) ? (a) : (b) )
86 if( num_vecs != 0 || ( num_vecs && vecs == NULL ) )
88 if( num_targ_vecs != 1 || ( num_targ_vecs && targ_vecs == NULL ) )
97 z0_val = targ_vecs[0].
values;
104 for( k = 0; k < sub_dim; ++k, z0_val += z0_val_s )
107 assert(*z0_val >= 0);
108 (*z0_val) =
sqrt((*z0_val));
int(* obj_create)(const struct RTOp_obj_type_vtbl_t *vtbl, const void *instance_data, void **obj)
const struct RTOp_obj_type_vtbl_t * obj_data_vtbl
int RTOp_TOp_ele_wise_sqrt_destroy(struct RTOp_RTOp *op)
static int RTOp_TOp_ele_wise_sqrt_apply_op(const struct RTOp_RTOp_vtbl_t *vtbl, const void *obj_data, const int num_vecs, const struct RTOp_SubVector vecs[], const int num_targ_vecs, const struct RTOp_MutableSubVector targ_vecs[], RTOp_ReductTarget reduct_obj)
#define RTOp_ERR_INVALID_NUM_VECS
int(* obj_free)(const struct RTOp_obj_type_vtbl_t *vtbl, const void *instance_data, void **obj)
const struct RTOp_obj_type_vtbl_t RTOp_obj_null_vtbl
void sqrt(DVectorSlice *vs_lhs, const DVectorSlice &vs_rhs)
vs_lhs = sqrt(vs_rhs)
int RTOp_TOp_ele_wise_sqrt_construct(struct RTOp_RTOp *op)
const struct RTOp_RTOp_vtbl_t * vtbl
const struct RTOp_RTOp_vtbl_t RTOp_TOp_ele_wise_sqrt_vtbl
#define RTOp_ERR_INVALID_NUM_TARG_VECS
Teuchos_Ordinal RTOp_index_type