@@ -335,9 +335,10 @@ void Euler::doEulerBackward()
335
335
*T = new double [_dimSys*_dimSys], // Iterationsmatrix
336
336
*jac = new double [_dimSys*_dimSys], // Jacobimatrix
337
337
*yHelp = new double [_dimSys],
338
- *fHelp = new double [_dimSys],
339
- *pHelp = new double [_dimSys]; // Hilfsvariale Pivotisierun
340
- memset (pHelp,0 ,_dimSys*sizeof (double ));
338
+ *fHelp = new double [_dimSys];
339
+
340
+ long int *pHelp = new long int [_dimSys]; // Hilfsvariale Pivotisierun
341
+ memset (pHelp,0 ,_dimSys*sizeof (long int ));
341
342
342
343
while ( _idid == 0 && _solverStatus != USER_STOP )
343
344
{
@@ -539,9 +540,9 @@ void Euler::doMidpoint()
539
540
*deltaZ = new double [_dimSys], // Hilfsvariable für Stufe
540
541
*f0 = new double [_dimSys],
541
542
*LSErhs = new double [_dimSys],
542
- *pHelp = new double [_dimSys],
543
543
*fHelp = new double [_dimSys]; // Hilfsvariable für rechte Seite
544
544
545
+ long int *pHelp = new long int [_dimSys]; // Hilfsvariale Pivotisierun
545
546
// Rechte Seite
546
547
double * k1 = new double [_dimSys];
547
548
0 commit comments