integer ilong;
intptr_t iofwr;
doublereal* wrkar = 0;
+ doublereal* wrkar_off;
integer iszwr;
integer ibb, ier = 0;
integer isz1, isz2, isz3, isz4;
if (ier > 0) {
goto L9013;
}
- double* wrkar_offset = reinterpret_cast<double*>(iofwr + wrkar);
+ wrkar_off = reinterpret_cast<double*>(iofwr + wrkar);
ipt1 = isz1;
ipt2 = ipt1 + isz2;
ipt3 = ipt2 + isz3;
/* --- Return 2*(IORDRU+1) coeff of 2*(IORDRU+1) polynoms of Hermite
--- */
- AdvApp2Var_ApproxF2var::mma1her_(iordru, wrkar_offset, iercod);
+ AdvApp2Var_ApproxF2var::mma1her_(iordru, wrkar_off, iercod);
if (*iercod > 0) {
goto L9100;
}
---- */
mma2cd3_(ndimen, nbpntu, &urootl[1], nbpntv, iordru, &sotbu1[1],
- &sotbu2[1], &ditbu1[1], &ditbu2[1], &wrkar_offset[ipt2], wrkar_offset,
+ &sotbu2[1], &ditbu1[1], &ditbu2[1], &wrkar_off[ipt2], wrkar_off,
&sosotb[sosotb_offset], &soditb[soditb_offset],
&disotb[disotb_offset], &diditb[diditb_offset]);
}
/* --- Return 2*(IORDRV+1) coeff of 2*(IORDRV+1) polynoms of Hermite
--- */
- AdvApp2Var_ApproxF2var::mma1her_(iordrv, &wrkar_offset[ipt1], iercod);
+ AdvApp2Var_ApproxF2var::mma1her_(iordrv, &wrkar_off[ipt1], iercod);
if (*iercod > 0) {
goto L9100;
}
---- */
mma2cd2_(ndimen, nbpntu, nbpntv, &vrootl[1], iordrv, &sotbv1[1],
- &sotbv2[1], &ditbv1[1], &ditbv2[1], &wrkar_offset[ipt3], &wrkar_offset[ipt1],
+ &sotbv2[1], &ditbv1[1], &ditbv2[1], &wrkar_off[ipt3], &wrkar_off[ipt1],
&sosotb[sosotb_offset], &soditb[soditb_offset],
&disotb[disotb_offset], &diditb[diditb_offset]);
}
if (*iordru >= 0 && *iordrv >= 0) {
mma2cd1_(ndimen, nbpntu, &urootl[1], nbpntv, &vrootl[1], iordru,
iordrv, &contr1[contr1_offset], &contr2[contr2_offset],
- &contr3[contr3_offset], &contr4[contr4_offset], &wrkar_offset[ipt2],
- &wrkar_offset[ipt3], wrkar_offset, &wrkar_offset[ipt1],
+ &contr3[contr3_offset], &contr4[contr4_offset], &wrkar_off[ipt2],
+ &wrkar_off[ipt3], wrkar_off, &wrkar_off[ipt1],
&sosotb[sosotb_offset], &soditb[soditb_offset],
&disotb[disotb_offset], &diditb[diditb_offset]);
}
logical ldbg;
intptr_t iofwr;
doublereal* wrkar = 0;
+ doublereal* wrkar_off;
integer iszwr;
integer ier;
integer isz1, isz2, isz3, isz4, isz5, isz6, isz7;
iszwr = isz1 + isz2 + isz3 + isz4 + isz5 + isz6 + isz7;
AdvApp2Var_SysBase anAdvApp2Var_SysBase;
anAdvApp2Var_SysBase.mcrrqst_(&c__8, &iszwr, wrkar, &iofwr, &ier);
- double* wrkar_offset = reinterpret_cast<double*> (iofwr + wrkar);
+ wrkar_off = reinterpret_cast<double*> (iofwr + wrkar);
if (ier > 0) {
goto L9013;
}
/* ----------------- Return Gauss coefficients of integration ----------------
*/
- AdvApp2Var_ApproxF2var::mmapptt_(ndjacu, nbpntu, iordru, wrkar_offset, iercod);
+ AdvApp2Var_ApproxF2var::mmapptt_(ndjacu, nbpntu, iordru, wrkar_off, iercod);
if (*iercod > 0) {
goto L9999;
}
- AdvApp2Var_ApproxF2var::mmapptt_(ndjacv, nbpntv, iordrv, &wrkar_offset[ipt1], iercod);
+ AdvApp2Var_ApproxF2var::mmapptt_(ndjacv, nbpntv, iordrv, &wrkar_off[ipt1], iercod);
if (*iercod > 0) {
goto L9999;
}
/* ------------------- Return max polynoms of Jacobi ------------
*/
- AdvApp2Var_ApproxF2var::mma2jmx_(ndjacu, iordru, &wrkar_offset[ipt5]);
- AdvApp2Var_ApproxF2var::mma2jmx_(ndjacv, iordrv, &wrkar_offset[ipt5]);
+ AdvApp2Var_ApproxF2var::mma2jmx_(ndjacu, iordru, &wrkar_off[ipt5]);
+ AdvApp2Var_ApproxF2var::mma2jmx_(ndjacv, iordrv, &wrkar_off[ipt5]);
/* ------ Calculate the coefficients and their contribution to the error ----
*/
mma2ce2_(numdec, ndimen, nbsesp, &ndimse[1], ndminu, ndminv, ndguli,
ndgvli, ndjacu, ndjacv, iordru, iordrv, nbpntu, nbpntv, &epsapr[1],
&sosotb[sosotb_offset], &disotb[disotb_offset], &soditb[soditb_offset],
- &diditb[diditb_offset], wrkar_offset, &wrkar_offset[ipt1],
- &wrkar_offset[ipt4], &wrkar_offset[ipt5], &wrkar_offset[ipt6], &wrkar_offset[ipt2],
- &wrkar_offset[ipt3], &patjac[patjac_offset], &errmax[1], &errmoy[1], ndegpu,
+ &diditb[diditb_offset], wrkar_off, &wrkar_off[ipt1],
+ &wrkar_off[ipt4], &wrkar_off[ipt5], &wrkar_off[ipt6], &wrkar_off[ipt2],
+ &wrkar_off[ipt3], &patjac[patjac_offset], &errmax[1], &errmoy[1], ndegpu,
ndegpv, itydec, iercod);
if (*iercod > 0) {
goto L9999;
integer ideb1, ibid1, ibid2, ncfja, ndgre, ilong,
ndwrk;
doublereal* wrkar = 0;
+ doublereal* wrkar_off;
integer nupil;
intptr_t iofwr;
doublereal uvpav[4] /* was [2][2] */;
ndwrk = isz1 + isz2 + isz3 + isz4 + isz5;
anAdvApp2Var_SysBase.mcrrqst_(&c__8, &ndwrk, wrkar, &iofwr, &ier);
- double* wrkar_off = reinterpret_cast<double*>(iofwr + wrkar);
+ wrkar_off = reinterpret_cast<double*>(iofwr + wrkar);
if (ier > 0) {
goto L9013; }
/* --> For the parameters of discretisation (NBROOT+2 extremities). */