ma->z = calloc(2 * ma->n + 1, sizeof(double));
ma->zswap = calloc(2 * ma->n + 1, sizeof(double));
ma->afs = calloc(2 * ma->n + 1, sizeof(double));
ma->afs1 = calloc(2 * ma->n + 1, sizeof(double));
ma->z = calloc(2 * ma->n + 1, sizeof(double));
ma->zswap = calloc(2 * ma->n + 1, sizeof(double));
ma->afs = calloc(2 * ma->n + 1, sizeof(double));
ma->afs1 = calloc(2 * ma->n + 1, sizeof(double));