51 #define max(a,b) ( (a) > (b) ? (a) : (b) )
52 #define min(a,b) ( (a) < (b) ? (a) : (b) )
90 if( num_vecs != 1 || ( num_vecs && vecs == NULL ) )
92 if( num_targ_vecs != 1 || ( num_targ_vecs && targ_vecs == NULL ) )
95 targ_vecs[0].sub_dim != vecs[0].sub_dim
106 z0_val = targ_vecs[0].
values;
116 for( k = 0; k < sub_dim; ++k, v0_val += v0_val_s, z0_val += z0_val_s )
122 assert(*z0_val > -1e-50);
125 else if ( (*lower_or_upper) == 0 )
127 (*z0_val) = ((*v0_val) <= (*inf_bound_limit)) ? 0.0 : (*z0_val);
131 (*z0_val) = ((*v0_val) >= (*inf_bound_limit)) ? 0.0 : (*z0_val);
143 ,
"TOp_Correct_Multipliers"
176 *ptr_inf_bound_limit = inf_bound_limit;
177 *ptr_lower_or_upper = lower_or_upper;
const struct RTOp_obj_type_vtbl_t RTOp_obj_value_index_vtbl
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_Correct_Multipliers_destroy(struct RTOp_RTOp *op)
#define RTOp_ERR_INCOMPATIBLE_VECS
const RTOp_value_type * values
static int RTOp_TOp_Correct_Multipliers_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
const struct RTOp_RTOp_vtbl_t RTOp_TOp_Correct_Multipliers_vtbl
const struct RTOp_RTOp_vtbl_t * vtbl
int RTOp_TOp_Correct_Multipliers_construct(RTOp_value_type inf_bound_limit, RTOp_index_type lower_or_upper, struct RTOp_RTOp *op)
int RTOp_TOp_Correct_Multipliers_init(RTOp_value_type inf_bound_limit, RTOp_index_type lower_or_upper, struct RTOp_RTOp *op)
#define RTOp_ERR_INVALID_NUM_TARG_VECS
Teuchos_Ordinal RTOp_index_type