Logic not trapping

Dev Shed Forums Sponsor:

August 30th, 2012, 10:38 AM

AstroTux

Contributing User

Join Date: Feb 2008

Posts: 593

Time spent in forums: 5 Days 21 m 29 sec
Reputation Power: 42

Logic not trapping

Hi,

This code is short and sweet. It is complete as posted. The problem is this: when the value of hdg_ins exceeds M_PI*2.0 (360 degrees) or goes negative (less than 0) the bounds check fails and doesn't execute.

1) Why
2) How to fix it?

Code:

{
	hdg_error = hdg_error + ((7.0/25920000.0)*M_PI);
	if(hdg_error >= (M_PI*2.0)) hdg_error = hdg_error - (M_PI*2.0);
	if(hdg_error < 0.0) hdg_error = hdg_error + M_PI*2;

	inertial_hdg = inertial_hdg + (rateyaw * cos(-roll)) + (ratepitch * sin(-roll));
	hdg_ins = (inertial_hdg * 0.025) + hdg_error;

	if(hdg_ins >= (M_PI*2.0)) hdg_ins = hdg_ins - (M_PI*2.0);
	if(hdg_ins < 0.0) hdg_ins = hdg_ins + (M_PI*2.0);
}

Many thanks in advance.

Best regards,
AstroTux.

August 30th, 2012, 01:41 PM

salem

Contributed User

Join Date: Jun 2005

Posts: 3,836

Folding Points: 153 Folding Title: Novice Folder

Time spent in forums: 2 Months 3 Weeks 2 Days 16 h 15 m 5 sec
Reputation Power: 1774

I think you would need to post a test case we can copy/paste/run for ourselves.

Compiling a bit of code and stepping through with gdb shows nothing unexpected here.

Code:

Breakpoint 1, main () at foo.c:18
18	    if(hdg_ins >= (M_PI*2.0)) hdg_ins = hdg_ins - (M_PI*2.0);
(gdb) list
13	    if(hdg_error < 0.0) hdg_error = hdg_error + M_PI*2;
14	
15	    inertial_hdg = inertial_hdg + (rateyaw * cos(-roll)) + (ratepitch * sin(-roll));
16	    hdg_ins = (inertial_hdg * 0.025) + hdg_error;
17	
18	    if(hdg_ins >= (M_PI*2.0)) hdg_ins = hdg_ins - (M_PI*2.0);
19	    if(hdg_ins < 0.0) hdg_ins = hdg_ins + (M_PI*2.0);
20	    }
21	    
22	    return 0;
(gdb) set hdg_ins = 7
(gdb) step
19	    if(hdg_ins < 0.0) hdg_ins = hdg_ins + (M_PI*2.0);
(gdb) step
22	    return 0;
(gdb) print hdg_ins
$5 = 0.71681469282041377

(gdb) run
Breakpoint 1, main () at foo.c:18
18	    if(hdg_ins >= (M_PI*2.0)) hdg_ins = hdg_ins - (M_PI*2.0);
(gdb) set hdg_ins = -7
(gdb) step
19	    if(hdg_ins < 0.0) hdg_ins = hdg_ins + (M_PI*2.0);
(gdb) step
22	    return 0;
(gdb) print hdg_ins
$6 = -0.71681469282041377
(gdb)

I set positive and negative out of range values, and it all seems fine.

The only thing I can think of is if hdg_ins is a float, and you're getting some weird rounding very close to the range limits, which results in the adjustments still being out of range.

__________________
If you dance barefoot on the broken glass of undefined behaviour, you've got to expect the occasional cut.
If at first you don't succeed, try writing your phone number on the exam paper

August 30th, 2012, 02:15 PM

dwise1_aol

Contributing User

Dev Shed God 2nd Plane (6000 - 6499 posts)

Join Date: Jan 2003

Location: USA

Posts: 6,125

Time spent in forums: 2 Months 2 Weeks 3 Days 17 h 54 sec
Reputation Power: 1949

It would also help if you gave us the test input values that you're using. Plus the variable declarations as salem alluded to.

August 30th, 2012, 03:03 PM

b49P23TIvg

Contributing User

Join Date: Aug 2011

Posts: 3,354

Time spent in forums: 1 Month 2 Weeks 3 Days 8 h 21 m 53 sec
Reputation Power: 383

use fmod, not if if if

Code:

#define PI2 (2*M_PI)
{
  hdg_error += (7.0/25920000.0)*M_PI;
  hdg_error = fmod(hdg_error+PI2,PI2);
  inertial_hdg = inertial_hdg + (rateyaw * cos(-roll)) + (ratepitch * sin(-roll));
  hdg_ins = (inertial_hdg * 0.025) + hdg_error;
  hdg_ins = fmod(hdg_ins+PI2,PI2);
}

__________________
[code]Code tags[/code] are essential for python code!

September 7th, 2012, 10:01 AM

AstroTux

Contributing User

Join Date: Feb 2008

Posts: 593

Time spent in forums: 5 Days 21 m 29 sec
Reputation Power: 42

Hi,

It was the float vs. integer thing.

e.g.

Code:

double result = 1/6;

1/6 treated like an integer. Destination type is double. How can I change this default behavior?

If I have:

Code:

double result = 1.0/6.0;

I get the correct result.

Best regards,
AstroTux.

September 7th, 2012, 12:19 PM

EEmaestro

Contributing User

Join Date: Apr 2012

Posts: 58

Time spent in forums: 1 Day 17 h 52 m 2 sec
Reputation Power: 2

Try this code & compare results.

Code:

     #include <stdio.h>

     int
     main (void)
     {
        double result;


        result = 1/6;
        printf ("\nresult = %f",result);

        result = 1./6.;
        printf ("\nresult = %f",result);

        result = (float) 1/6;
        printf ("\nresult = %f",result);

        result = (double) 1/6;
        printf ("\nresult = %f",result);

       return 0;
     }

Last edited by EEmaestro : September 7th, 2012 at 12:23 PM.

Logic not trapping

Developer Shed Advertisers and Affiliates