51 #define max(a,b) ( (a) > (b) ? (a) : (b) )
52 #define min(a,b) ( (a) < (b) ? (a) : (b) )
87 if( num_vecs != 2 || ( num_vecs && vecs == NULL ) )
89 if( num_targ_vecs != 0 || ( num_targ_vecs && targ_vecs == NULL ) )
92 vecs[1].sub_dim != vecs[0].sub_dim
113 for( k = 0; k < sub_dim; ++k, v0_val += v0_val_s, v1_val += v1_val_s )
116 gamma_i = fabs((*v1_val)) / ( 1.0 + fabs((*v0_val)) );
117 (*gamma) =
max( (*gamma), gamma_i );
const struct RTOp_RTOp_vtbl_t RTOp_ROp_max_rel_step_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
#define RTOp_ERR_INCOMPATIBLE_VECS
const RTOp_value_type * values
int RTOp_get_reduct_max_value_op(const struct RTOp_RTOp_vtbl_t *vtbl, const void *obj_data, RTOp_reduct_op_func_ptr_t *reduct_op_func_ptr)
#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
int RTOp_ROp_max_rel_step_construct(struct RTOp_RTOp *op)
int RTOp_ROp_max_rel_step_destroy(struct RTOp_RTOp *op)
const struct RTOp_RTOp_vtbl_t * vtbl
int RTOp_reduct_max_value(const struct RTOp_RTOp_vtbl_t *vtbl, const void *obj_data, RTOp_ReductTarget in_targ_obj, RTOp_ReductTarget inout_targ_obj)
static int RTOp_ROp_max_rel_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)
RTOp_value_type RTOp_ROp_max_rel_step_val(RTOp_ReductTarget reduct_obj)
#define RTOp_ERR_INVALID_NUM_TARG_VECS
const struct RTOp_obj_type_vtbl_t RTOp_obj_value_vtbl
Teuchos_Ordinal RTOp_index_type