52 #define max(a,b) ( (a) > (b) ? (a) : (b) )
53 #define min(a,b) ( (a) < (b) ? (a) : (b) )
57 #include "RTOp_obj_value_vtbl.h"
58 #include "RTOp_reduct_min_value.h"
64 static int RTOp_ROp_fraction_to_zero_boundary_apply_op(
68 , RTOp_ReductTarget reduct_obj )
75 RTOp_value_type *tau = (RTOp_value_type*)obj_data;
77 RTOp_value_type *alpha_max = (RTOp_value_type*)reduct_obj;
79 RTOp_index_type sub_dim;
81 const RTOp_value_type *v0_val;
84 const RTOp_value_type *v1_val;
88 register RTOp_index_type k;
90 RTOp_value_type alpha_max_ith;
95 if( num_vecs != 2 || ( num_vecs && vecs == NULL ) )
96 return RTOp_ERR_INVALID_NUM_VECS;
97 if( num_targ_vecs != 0 || ( num_targ_vecs && targ_vecs == NULL ) )
98 return RTOp_ERR_INVALID_NUM_TARG_VECS;
100 vecs[1].sub_dim != vecs[0].sub_dim
102 return RTOp_ERR_INCOMPATIBLE_VECS;
109 sub_dim = vecs[0].sub_dim;
111 v0_val = vecs[0].values;
112 v0_val_s = vecs[0].values_stride;
114 v1_val = vecs[1].values;
115 v1_val_s = vecs[1].values_stride;
120 for( k = 0; k < sub_dim; ++k, v0_val += v0_val_s, v1_val += v1_val_s )
123 alpha_max_ith = ((*v1_val) >= 0) ? 1.0 : (*tau)*(*v0_val)/(-(*v1_val));
125 (*alpha_max) = min( (*alpha_max), alpha_max_ith );
131 static int RTOp_ROp_fraction_to_zero_boundary_reduct_obj_reinit(
133 ,RTOp_ReductTarget reduct_obj
136 RTOp_value_type *alpha_max = (RTOp_value_type*)reduct_obj;
145 ,&RTOp_obj_value_vtbl
146 ,
"ROp_fraction_to_zero_boundary"
147 ,RTOp_ROp_fraction_to_zero_boundary_reduct_obj_reinit
148 ,RTOp_ROp_fraction_to_zero_boundary_apply_op
149 ,RTOp_reduct_min_value
150 ,RTOp_get_reduct_min_value_op
155 int RTOp_ROp_fraction_to_zero_boundary_construct( RTOp_value_type tau,
struct RTOp_RTOp* op )
161 op->vtbl = &RTOp_ROp_fraction_to_zero_boundary_vtbl;
162 op->vtbl->obj_data_vtbl->
obj_create(NULL,NULL,&op->obj_data);
163 return RTOp_ROp_fraction_to_zero_boundary_init(tau,op);
166 int RTOp_ROp_fraction_to_zero_boundary_destroy(
struct RTOp_RTOp* op )
168 op->vtbl->obj_data_vtbl->
obj_free(NULL,NULL,&op->obj_data);
174 int RTOp_ROp_fraction_to_zero_boundary_init( RTOp_value_type tau,
struct RTOp_RTOp* op )
176 RTOp_value_type *ptr_tau = (RTOp_value_type*)op->obj_data;
180 RTOp_value_type RTOp_ROp_fraction_to_zero_boundary_val(RTOp_ReductTarget reduct_obj)
182 return *((RTOp_value_type*)reduct_obj);
int(* obj_create)(const struct RTOp_obj_type_vtbl_t *vtbl, const void *instance_data, void **obj)
int(* obj_free)(const struct RTOp_obj_type_vtbl_t *vtbl, const void *instance_data, void **obj)