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

Fix condition in both second scaling loops

This commit is contained in:
Martin Kroeker 2018-02-18 12:37:09 +01:00 committed by GitHub
parent 6940c59a88
commit 798f1595d5
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
1 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
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 (ABS(*dd1) <= RGAMSQ) { while ((*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 (ABS(*dd1) <= GAMSQ) { while ((*dd1) >= GAMSQ) {
*dd1 = *dd1 / (GAM * GAM); *dd1 = *dd1 / (GAM * GAM);
*dx1 = *dx1 * GAM; *dx1 = *dx1 * GAM;
dh11 = dh11 * GAM; dh11 = dh11 * GAM;
@ -197,7 +197,7 @@ void CNAME(FLOAT *dd1, FLOAT *dd2, FLOAT *dx1, FLOAT dy1, FLOAT *dparam){
} }
else else
{ {
while (ABS(*dd2) <= GAMSQ) { while (ABS(*dd2) >= GAMSQ) {
*dd2 = *dd2 / (GAM * GAM); *dd2 = *dd2 / (GAM * GAM);
dh21 = dh21 * GAM; dh21 = dh21 * GAM;
dh22 = dh22 * GAM; dh22 = dh22 * GAM;