floraachy
/
OpenBLAS
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Fix conditionals in the rescaling against GAMSQ

This commit is contained in:
Martin Kroeker 2018-02-18 12:54:52 +01:00 committed by GitHub
parent 53026dc63a
commit 72f14a0363
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
interface/rotmg.c
View File

@ -152,7 +152,7 @@ void CNAME(FLOAT *dd1, FLOAT *dd2, FLOAT *dx1, FLOAT dy1, FLOAT *dparam){
} }
if( *dd1 <= RGAMSQ ) if( *dd1 <= RGAMSQ )
{ {
while ((*dd1) <= RGAMSQ) { while (ABS(*dd1) <= RGAMSQ) {
*dd1 = *dd1 * (GAM * GAM); *dd1 = *dd1 * (GAM * GAM);
*dx1 = *dx1 / GAM; *dx1 = *dx1 / GAM;
dh11 = dh11 / GAM; dh11 = dh11 / GAM;
@ -161,7 +161,7 @@ void CNAME(FLOAT *dd1, FLOAT *dd2, FLOAT *dx1, FLOAT dy1, FLOAT *dparam){
} }
else else
{ {
while ((*dd1) >= GAMSQ) { while (ABS(*dd1) >= GAMSQ) {
*dd1 = *dd1 / (GAM * GAM); *dd1 = *dd1 / (GAM * GAM);
*dx1 = *dx1 * GAM; *dx1 = *dx1 * GAM;
dh11 = dh11 * GAM; dh11 = dh11 * GAM;