/*@null@*/ static PyObject *
PySip_new(
PyTypeObject* type,
/*@unused@*/ PyObject* args,
/*@unused@*/ PyObject* kwds) {
PySip* self;
self = (PySip*)type->tp_alloc(type, 0);
if (self != NULL) {
sip_clear(&self->x);
}
return (PyObject*)self;
}
int
sip_init(
sip_t* sip,
const unsigned int a_order, const double* a,
const unsigned int b_order, const double* b,
const unsigned int ap_order, const double* ap,
const unsigned int bp_order, const double* bp,
const double* crpix /* [2] */) {
unsigned int a_size = 0;
unsigned int b_size = 0;
unsigned int ap_size = 0;
unsigned int bp_size = 0;
unsigned int scratch_size = 0;
int status = 0;
struct wcserr** err = NULL;
static const char *function = "sip_init";
assert(sip != NULL);
sip_clear(sip);
err = &(sip->err);
/* We we have one of A/B or AP/BP, we must have both. */
if ((a == NULL) ^ (b == NULL)) {
return wcserr_set(
SIP_ERRMSG(WCSERR_BAD_COORD_TRANS),
"Both A and B SIP transform must be defined");
}
if ((ap == NULL) ^ (bp == NULL)) {
return wcserr_set(
SIP_ERRMSG(WCSERR_BAD_COORD_TRANS),
"Both AP and BP SIP transform must be defined");
}
if (a != NULL) {
sip->a_order = a_order;
a_size = (a_order + 1) * (a_order + 1) * sizeof(double);
sip->a = malloc(a_size);
if (sip->a == NULL) {
sip_free(sip);
status = wcserr_set(
SIP_ERRMSG(WCSERR_MEMORY), "Memory allocation failed");
goto exit;
}
memcpy(sip->a, a, a_size);
if (a_order > scratch_size) {
scratch_size = a_order;
}
sip->b_order = b_order;
b_size = (b_order + 1) * (b_order + 1) * sizeof(double);
sip->b = malloc(b_size);
if (sip->b == NULL) {
sip_free(sip);
status = wcserr_set(
SIP_ERRMSG(WCSERR_MEMORY), "Memory allocation failed");
goto exit;
}
memcpy(sip->b, b, b_size);
if (b_order > scratch_size) {
scratch_size = b_order;
}
}
if (ap != NULL) {
sip->ap_order = ap_order;
ap_size = (ap_order + 1) * (ap_order + 1) * sizeof(double);
sip->ap = malloc(ap_size);
if (sip->ap == NULL) {
sip_free(sip);
status = wcserr_set(
SIP_ERRMSG(WCSERR_MEMORY), "Memory allocation failed");
goto exit;
}
memcpy(sip->ap, ap, ap_size);
if (ap_order > scratch_size) {
scratch_size = ap_order;
}
sip->bp_order = bp_order;
bp_size = (bp_order + 1) * (bp_order + 1) * sizeof(double);
sip->bp = malloc(bp_size);
if (sip->bp == NULL) {
sip_free(sip);
status = wcserr_set(
SIP_ERRMSG(WCSERR_MEMORY), "Memory allocation failed");
goto exit;
}
memcpy(sip->bp, bp, bp_size);
if (bp_order > scratch_size) {
scratch_size = bp_order;
}
}
scratch_size = (scratch_size + 1) * sizeof(double);
sip->scratch = malloc(scratch_size);
if (sip->scratch == NULL) {
sip_free(sip);
status = wcserr_set(
SIP_ERRMSG(WCSERR_MEMORY), "Memory allocation failed");
goto exit;
}
sip->crpix[0] = crpix[0];
sip->crpix[1] = crpix[1];
exit:
return status;
}
