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

fixed 2 bugs as reported by Brendan Tracey

This commit is contained in:
wernsaar 2014-05-02 11:34:26 +02:00
parent 66198faab6
commit 025fc914cc
1 changed files with 14 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
interface/rotmg.c
View File

@ -26,7 +26,8 @@ USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
*****************************************************************************/ *****************************************************************************/
/************************************************************************************** /**************************************************************************************
* 2014/02/28 Saar * 2014/05/02 Saar
* fixed two bugs as reported by Brendan Tracey
* Test with lapack-3.5.0 : OK * Test with lapack-3.5.0 : OK
* *
**************************************************************************************/ **************************************************************************************/
@ -144,11 +145,14 @@ void CNAME(FLOAT *dd1, FLOAT *dd2, FLOAT *dx1, FLOAT dy1, FLOAT *dparam){
dflag = -ONE; dflag = -ONE;
} }
else else
{
if(dflag == ONE)
{ {
dh21 = -ONE; dh21 = -ONE;
dh12 = ONE; dh12 = ONE;
dflag = -ONE; dflag = -ONE;
} }
}
if( *dd1 <= RGAMSQ ) if( *dd1 <= RGAMSQ )
{ {
*dd1 = *dd1 * (GAM * GAM); *dd1 = *dd1 * (GAM * GAM);
@ -177,11 +181,14 @@ void CNAME(FLOAT *dd1, FLOAT *dd2, FLOAT *dx1, FLOAT dy1, FLOAT *dparam){
dflag = -ONE; dflag = -ONE;
} }
else else
{
if(dflag == ONE)
{ {
dh21 = -ONE; dh21 = -ONE;
dh12 = ONE; dh12 = ONE;
dflag = -ONE; dflag = -ONE;
} }
}
if( ABS(*dd2) <= RGAMSQ ) if( ABS(*dd2) <= RGAMSQ )
{ {
*dd2 = *dd2 * (GAM * GAM); *dd2 = *dd2 * (GAM * GAM);