static COMMAND_FUNC( do_delvec )
{
Data_Obj *dp;
dp=pick_obj("");
if( dp==NULL ) return;
delvec(dp);
}
static void _delete_subobjects(QSP_ARG_DECL Data_Obj *dp) /** delete all subimages */
{
Node *np;
/*
* delvec() will remove the parent's child list
* after deleting last child.
*
* List elements are de-allocated by delvec()
*/
while( OBJ_CHILDREN( dp ) != NULL ){
np=QLIST_HEAD( OBJ_CHILDREN( dp ) );
delvec( (Data_Obj *) NODE_DATA(np) );
}
}
static void x_close_movie(QSP_ARG_DECL Movie *mvip)
{
delvec(QSP_ARG (Data_Obj *)mvip->mvi_data);
}
void inner(QSP_ARG_DECL Data_Obj *dpto,Data_Obj *dpfr1,Data_Obj *dpfr2)
{
//dimension_t _n; /* dot prod len */
dimension_t i,j;
Vec_Obj_Args oa1, *oap=&oa1;
//Dimension_Set sizes={{1,1,1,1,1}};
Dimension_Set *sizes;
index_t dst_indices[N_DIMENSIONS]={0,0,0,0,0};
index_t src1_indices[N_DIMENSIONS]={0,0,0,0,0};
index_t src2_indices[N_DIMENSIONS]={0,0,0,0,0};
Data_Obj *col_dp;
sizes=NEW_DIMSET;
for(i=0;i<N_DIMENSIONS;i++)
SET_DIMENSION(sizes,i,1);
#ifdef CAUTIOUS
clear_obj_args(oap);
#endif /* CAUTIOUS */
/* The types and precisions should be whatever is allowed by vdot,
* which is float, double, real and complex...
*/
if( ! is_good_for_inner(dpto,"inner") ) return;
if( ! is_good_for_inner(dpfr1,"inner") ) return;
if( ! is_good_for_inner(dpfr2,"inner") ) return;
/* we need to make sure that the types and precisions MATCH! */
if( ! prec_and_type_match(dpto,dpfr1,"inner") ) return;
if( ! prec_and_type_match(dpto,dpfr2,"inner") ) return;
if( OBJ_ROWS(dpto) != OBJ_ROWS(dpfr1) ){
sprintf(DEFAULT_ERROR_STRING,
"inner: dpto %s (%d) and first operand %s (%d) must have same # rows",
OBJ_NAME(dpto),OBJ_ROWS(dpto),OBJ_NAME(dpfr1),OBJ_ROWS(dpfr1));
NWARN(DEFAULT_ERROR_STRING);
return;
}
if( OBJ_COLS(dpto) != OBJ_COLS(dpfr2) ){
sprintf(DEFAULT_ERROR_STRING,
"inner: target %s (%d) and second operand %s (%d) must have same # columns",
OBJ_NAME(dpto),OBJ_COLS(dpto),OBJ_NAME(dpfr2),OBJ_COLS(dpfr2));
NWARN(DEFAULT_ERROR_STRING);
return;
}
if( OBJ_COLS(dpfr1) != OBJ_ROWS(dpfr2) ){
sprintf(DEFAULT_ERROR_STRING,
"inner: # cols of operand %s (%d) must match # rows of operand %s (%d)",
OBJ_NAME(dpfr1),OBJ_COLS(dpfr1),OBJ_NAME(dpfr2),OBJ_ROWS(dpfr2));
NWARN(DEFAULT_ERROR_STRING);
return;
}
//_n=OBJ_COLS(dpfr1); /* the length of each dot product we will compute */
if( IS_COMPLEX(dpto) ) SET_OA_ARGSTYPE(oap,COMPLEX_ARGS);
else SET_OA_ARGSTYPE(oap,REAL_ARGS);
/* vdot things it's inputs have the same shape, so if we are taking the inner
* product of a column vector with a row vector, we have to transpose one of
* the inputs...
*/
if( OBJ_ROWS(dpfr1) > 1 )
SET_OA_SRC1(oap,d_subscript(QSP_ARG dpfr1,0) ); /* subscript first row */
else
SET_OA_SRC1(oap,dpfr1); /* object is a row */
if( OBJ_COLS(dpto) > 1 )
col_dp=c_subscript(QSP_ARG dpfr2,0);
else
col_dp=dpfr2;
SET_OA_DEST(oap,mk_subimg(QSP_ARG dpto,0,0,"target pixel",1,1) );
//[sizes setDimensionAtIndex : 1 withValue : OBJ_ROWS(col_dp) ];
SET_DIMENSION(sizes,1,OBJ_ROWS(col_dp));
SET_DIMENSION(sizes,0,OBJ_COMPS(col_dp));
SET_OA_SRC2(oap,make_equivalence(QSP_ARG "_transposed_column",
col_dp,sizes,OBJ_PREC_PTR(col_dp)) );
for(i=0;i<OBJ_ROWS(dpto);i++){
src1_indices[2]=i;
SET_OBJ_DATA_PTR( OA_SRC1(oap), multiply_indexed_data(dpfr1,src1_indices) );
for(j=0;j<OBJ_COLS(dpto);j++){
dst_indices[2]=i; /* k_th component */
dst_indices[1]=j; /* k_th component */
SET_OBJ_DATA_PTR( OA_DEST(oap), multiply_indexed_data(dpto,dst_indices) );
src2_indices[1]=j;
SET_OBJ_DATA_PTR( OA_SRC2(oap), multiply_indexed_data(dpfr2,src2_indices) );
vdot(QSP_ARG oap);
}
}
delvec(QSP_ARG OA_SRC2(oap) ); /* "_transposed_column" */
if( OA_SRC1(oap) != dpfr1 )
delvec(QSP_ARG OA_SRC1(oap) );
if( col_dp != dpfr2 )
delvec(QSP_ARG col_dp);
delvec(QSP_ARG OA_DEST(oap) );
}
static int change_size(QSP_ARG_DECL Data_Obj *dst_dp,Data_Obj *src_dp )
{
Dimension_Set ef, *enlargement_factor=&ef;
Dimension_Set rf, *reduction_factor=&rf;
Vec_Obj_Args oa1, *oap=&oa1;
Dimension_Set size_ds, n_ds;
Dimension_Set *size_dsp=(&size_ds), *n_dsp=(&n_ds);
Data_Obj *src_ss_dp, *dst_ss_dp;
dimension_t i,j,k,l,m;
index_t offsets[N_DIMENSIONS]={0,0,0,0,0};
incr_t dst_incrs[N_DIMENSIONS], src_incrs[N_DIMENSIONS];
index_t dst_indices[N_DIMENSIONS]={0,0,0,0,0}, src_indices[N_DIMENSIONS]={0,0,0,0,0};
/* For simplicity, we don't allow size changes to be combined with conversions */
if( !dp_same_prec(QSP_ARG dst_dp,src_dp,"change_size") )
return(-1);
for(i=0;i<N_DIMENSIONS;i++){
if( OBJ_TYPE_DIM(dst_dp,i) > OBJ_TYPE_DIM(src_dp,i) ){
/* enlargement - subsample the destination */
SET_DIMENSION(enlargement_factor,i,
floor( OBJ_TYPE_DIM(dst_dp,i) / OBJ_TYPE_DIM(src_dp,i) ) );
SET_DIMENSION(reduction_factor,i, 0);
SET_DIMENSION(size_dsp,i, OBJ_TYPE_DIM(src_dp,i) );
SET_DIMENSION(n_dsp,i, DIMENSION(enlargement_factor,i) );
dst_incrs[i] = DIMENSION(n_dsp,i);
src_incrs[i] = 1;
} else {
/* reduction - subsample the source */
SET_DIMENSION(reduction_factor,i,
ceil( OBJ_TYPE_DIM(src_dp,i) / OBJ_TYPE_DIM(dst_dp,i) ) );
SET_DIMENSION(enlargement_factor,i, 0 );
SET_DIMENSION(size_dsp,i, floor( OBJ_TYPE_DIM(src_dp,i) /
DIMENSION(reduction_factor,i) ) );
/* We don't need to do this multiple times, just pick one and do it */
/*SET_DIMENSION(n_dsp,i, DIMENSION(reduction_factor,i) ); */
SET_DIMENSION(n_dsp,i, 1);
src_incrs[i] = DIMENSION(reduction_factor,i);
dst_incrs[i] = 1;
}
}
/* make the subsamples.
* the column increment is expressed in columns, etc.
*/
dst_ss_dp=make_subsamp(QSP_ARG "chngsize_dst_obj",dst_dp,size_dsp,offsets,dst_incrs);
src_ss_dp=make_subsamp(QSP_ARG "chngsize_src_obj",src_dp,size_dsp,offsets,src_incrs);
clear_obj_args(oap);
SET_OA_DEST(oap,dst_ss_dp);
SET_OA_SRC_OBJ(oap,0, src_ss_dp);
SET_OA_ARGSTYPE(oap, REAL_ARGS);
SET_OA_PFDEV(oap,OBJ_PFDEV(dst_dp));
for(i=0;i<DIMENSION(n_dsp,4);i++){ /* foreach sequence to copy */
if( dst_incrs[4] > 1 )
dst_indices[4]=i;
else src_indices[4]=i;
for(j=0;j<DIMENSION(n_dsp,3);j++){ /* foreach frame to copy */
if( dst_incrs[3] > 1 )
dst_indices[3]=j;
else src_indices[3]=j;
for(k=0;k<DIMENSION(n_dsp,2);k++){ /* foreach row */
if( dst_incrs[2] > 1 )
dst_indices[2]=k;
else src_indices[2]=k;
for(l=0;l<DIMENSION(n_dsp,1);l++){ /* foreach col */
if( dst_incrs[1] > 1 )
dst_indices[1]=l;
else src_indices[1]=l;
for(m=0;m<DIMENSION(n_dsp,0);m++){ /* foreach comp */
if( dst_incrs[0] > 1 )
dst_indices[0]=m;
else src_indices[0]=m;
/* relocate the appropriate subsample */
SET_OBJ_DATA_PTR(dst_ss_dp, multiply_indexed_data(dst_dp,dst_indices) );
SET_OBJ_DATA_PTR(src_ss_dp, multiply_indexed_data(src_dp,src_indices) );
// This doesn't check for cuda obj...
//vmov(oap);
perf_vfunc(QSP_ARG FVMOV, oap );
}
}
}
}
}
delvec(QSP_ARG dst_ss_dp);
delvec(QSP_ARG src_ss_dp);
SET_OBJ_FLAG_BITS(dst_dp, DT_ASSIGNED);
return(0);
}
static Data_Obj *insure_ram_obj(QSP_ARG_DECL Data_Obj *dp)
{
Data_Obj *tmp_dp;
char *tname;
Data_Area *save_ap;
Data_Obj *c_dp=NULL;
if( OBJ_IS_RAM(dp) ) return dp;
// This object lives on a different platform.
// We create a copy in RAM, and download the data
// using the platform download function.
save_ap = curr_ap;
curr_ap = ram_area_p;
tname = getbuf( strlen(OBJ_NAME(dp)) + strlen(DNAME_PREFIX) + 1 );
sprintf(tname,"%s%s",DNAME_PREFIX,OBJ_NAME(dp));
tmp_dp = dup_obj(QSP_ARG dp, tname);
givbuf(tname);
if( tmp_dp == NO_OBJ ){
// This can happen if the object is subscripted,
// as the bracket characters are illegal in names
return NO_OBJ;
}
curr_ap = save_ap;
// We can't download if the source data is not contiguous...
//
// We have a problem with bit precision, because the bits can
// be non-contiguous when the long words are - any time the number of columns
// is not evenly divided by the bits-per-word
if( (! IS_CONTIGUOUS(dp)) && ! HAS_CONTIGUOUS_DATA(dp) ){
Vec_Obj_Args oa1, *oap=&oa1;
advise("object is not contiguous, and does not have contiguous data...");
longlist(QSP_ARG dp);
save_ap = curr_ap;
curr_ap = OBJ_AREA( dp );
tname = getbuf( strlen(OBJ_NAME(dp)) + strlen(CNAME_PREFIX) + 1 );
sprintf(tname,"%s%s",CNAME_PREFIX,OBJ_NAME(dp));
c_dp = dup_obj(QSP_ARG dp, tname );
givbuf(tname);
curr_ap = save_ap;
// Now do the move...
setvarg2(oap,c_dp,dp);
if( IS_BITMAP(dp) ){
SET_OA_SBM(oap,dp);
SET_OA_SRC1(oap,NO_OBJ);
}
if( IS_REAL(dp) ) /* BUG case for QUAT too? */
OA_ARGSTYPE(oap) = REAL_ARGS;
else if( IS_COMPLEX(dp) ) /* BUG case for QUAT too? */
OA_ARGSTYPE(oap) = COMPLEX_ARGS;
else if( IS_QUAT(dp) ) /* BUG case for QUAT too? */
OA_ARGSTYPE(oap) = QUATERNION_ARGS;
else
//ERROR1("CAUTIOUS: insure_ram_obj: bad argset type!?");
assert( AERROR("insure_ram_obj: bad argset type!?") );
//fprintf(stderr,"insure_ram_obj: moving remote data to a contiguous object\n");
call_vfunc( QSP_ARG FIND_VEC_FUNC(FVMOV), oap );
//fprintf(stderr,"insure_ram_obj: DONE moving remote data to a contiguous object\n");
dp = c_dp;
}
gen_obj_dnload(QSP_ARG tmp_dp, dp);
if( c_dp != NO_OBJ )
delvec(QSP_ARG c_dp);
// BUG - when to delete?
// We try using the VOLATILE flag. This will work as long as
// the input object is not VOLATILE!?
SET_OBJ_FLAG_BITS(tmp_dp, DT_VOLATILE ) ;
return tmp_dp;
}
