Miroirs/newrpl
1
0
Fork 0
mirror of https://git.code.sf.net/p/newrpl/sources synced 2024-11-16 19:51:25 +01:00
Code Issues Projects Releases Packages Wiki Activity Actions

lib 64: TCHEBYCHEFF, LEGENDRE, HERMITE error handling

Browse source
This commit is contained in:
MPo 2016-10-26 16:42:15 +02:00
parent e0ee74694e
commit 914fe01061
1 changed files with 17 additions and 7 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

24
newrpl/lib-64-arithmetic.c
View file

@ -1850,6 +1850,9 @@ void LIB_HANDLER()
}
rplDropData(1);
WORDPTR *savestk=DSTop; // Drop arguments in case of error
if (n == 0) {
rplPushData((WORDPTR)(one_bint));
int elements = 1;
@ -1858,6 +1861,9 @@ void LIB_HANDLER()
rplDropData(elements);
rplPushData(newmat);
}
if(!newmat || Exceptions) {
if(DSTop>savestk) DSTop=savestk;
}
return;
}
else if (n == 1) {
@ -1871,8 +1877,6 @@ void LIB_HANDLER()
rplPushData((WORDPTR)(two_bint));
rplPushData((WORDPTR)(zero_bint));
break;
// default:
// break;
}
int elements = 2;
WORDPTR newmat=rplMatrixCompose(0,elements);
@ -1880,6 +1884,9 @@ void LIB_HANDLER()
rplDropData(elements);
rplPushData(newmat);
}
if(!newmat || Exceptions) {
if(DSTop>savestk) DSTop=savestk;
}
return;
}
else {
@ -1908,11 +1915,7 @@ void LIB_HANDLER()
rplOverwriteData(cur*(n+1)+1,(WORDPTR)(one_bint));
rplOverwriteData(oth*(n+1)+2,(WORDPTR)(two_bint));
break;
// default:
// break;
}
// rplOverwriteData(cur*(n+1)+1,(WORDPTR)(one_bint)); // [ ... 0 0 1]
// rplOverwriteData(oth*(n+1)+2,(WORDPTR)(one_bint)); // [ ... 0 1 0]
// recrsive formula
rplNumberToRReg(2, (WORDPTR)(two_bint));
@ -1953,7 +1956,10 @@ void LIB_HANDLER()
rplCheckResultAndError(&RReg[4]); // next
WORDPTR newnumber=rplNewReal(&RReg[4]);
if(!newnumber) return;
if(!newnumber || Exceptions) {
if(DSTop>savestk) DSTop=savestk;
return;
}
rplOverwriteData(cur*(n+1)+j+1,newnumber);
}
@ -1966,6 +1972,10 @@ void LIB_HANDLER()
rplDropData(elements);
rplPushData(newmat);
}
else {
if(DSTop>savestk) DSTop=savestk;
return;
}
}
return;