52 #define max(a,b) ( (a) > (b) ? (a) : (b) )
53 #define min(a,b) ( (a) < (b) ? (a) : (b) )
97 if( num_vecs != 3 || ( num_vecs && vecs == NULL ) )
99 if( num_targ_vecs != 1 || ( num_targ_vecs && targ_vecs == NULL ) )
102 vecs[1].sub_dim != vecs[0].sub_dim
103 || vecs[2].sub_dim != vecs[0].sub_dim
104 || targ_vecs[0].sub_dim != vecs[0].sub_dim
114 z0_val = targ_vecs[0].
values;
129 for( k = 0; k < sub_dim; ++k, v0_val += v0_val_s, v1_val += v1_val_s, v2_val += v2_val_s, z0_val += z0_val_s )
132 (*z0_val) = -(*v1_val) + (*mu)*(*v0_val) + (*alpha)*(*v0_val)*(*v1_val)*(*v2_val);
143 ,
"TOp_multiplier_step"
int RTOp_TOp_multiplier_step_init(RTOp_value_type mu, RTOp_value_type alpha, struct RTOp_RTOp *op)
static int RTOp_TOp_multiplier_step_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)
int(* obj_create)(const struct RTOp_obj_type_vtbl_t *vtbl, const void *instance_data, void **obj)
const struct RTOp_obj_type_vtbl_t RTOp_obj_value_value_vtbl
const struct RTOp_obj_type_vtbl_t * obj_data_vtbl
int RTOp_TOp_multiplier_step_destroy(struct RTOp_RTOp *op)
#define RTOp_ERR_INCOMPATIBLE_VECS
const RTOp_value_type * values
#define RTOp_ERR_INVALID_NUM_VECS
int RTOp_TOp_multiplier_step_construct(RTOp_value_type mu, RTOp_value_type alpha, struct RTOp_RTOp *op)
const struct RTOp_RTOp_vtbl_t RTOp_TOp_multiplier_step_vtbl
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
const struct RTOp_RTOp_vtbl_t * vtbl
#define RTOp_ERR_INVALID_NUM_TARG_VECS
Teuchos_Ordinal RTOp_index_type