/* A small test program for the new create / delete routines */
int main() {
pdl *bar;
pdl *foo = pdl_create(PDL_PERM);
int inds[2] = {1,1};
pdl_dump(foo);
pdl_reallocdims(foo,2);
foo->dims[0] = 5;
foo->dims[1] = 6;
pdl_reallocphysdata(foo);
pdl_dump(foo);
bar = pdl_createtrans(foo, pdl_affine_rectslice_transvtable);
pdl_dump(bar);
pdl_trans_affine_rectslice *trans =
((pdl_trans_affine_rectslice *)(foo->trans));
trans->starts[0] = 1;
trans->ends[0] = 3;
trans->starts[1] = 2;
trans->ends[1] = 4;
trans->steps[1] = 2;
pdl_transchanged(bar);
pdl_dump(bar);
pdl_make_physical_affine(bar);
pdl_dump(bar);
pdl_make_physical(bar);
pdl_dump(bar);
pdl_set(bar,2.0,inds);
pdl_changed(bar);
pdl_dump(foo);
pdl_make_physical_affine(foo);
}