RDS Datum aus MJD errechnen?

Das würde ich auf einem kleinen Controller einfach so machen:

1

int16_t d, m, y;

2

3

//--------------------------------------

4

5

void MJDValueToDMY (int16_t mjdvalue)

6

7

// result correct for 01.01.2010 ... 31.12.2099

8

9

{

10

  int16_t nd[12] = {31, 28, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31};

11

  int16_t j;

12

13

  d = 0;

14

  m = 0;

15

  y = 10;

16

17

  j = mjdvalue-55198+1;  // 55198 = MJD of 01.01.2010

18

19

  while (j>0)

20

    {

21

    d = j;

22

23

    m++;

24

    if (m==13) {m = 1; y++;};

25

26

    j = j-nd[m];

27

    // take leap years into account

28

    if (y%4==0 && m==2) j--;

29

    }

30

}

Die Idee dahinter sollte klar sein. Nach dem Aufruf der Funktion steht 
das Ergebnis in d, m, y (Tag, Monat, Jahr. Jahreszahl zweistellig, also 
nicht 2011 sondern 11).

Kannst natürlich auch dieses Monstrum implementieren, das immer als 
Umrechnungsmethode angegeben wird:

y = (4*mjd-1) / 146097;
j = 4*mjd-1-146097*y;
d = mjd / 4;
j = (4*d+3) / 1461;
d = 4*d+3-1461*mjd;
...

und noch ein paar Zeilen mehr. Das funzt dann auch zurück bis ins 
siebzehnte Jahrhundert oder so, gibt aber sicher mehr Code, und wenn der 
Controller die darin vorkommenden 16-Bit-Integer-Multiplikationen und 
Divisionen nicht hardwaremäßig kann, noch mehr.

Aus meiner Sicht ist das Array nicht korrekt!
Meines sieht so aus:
const uint8_t nd[13] = {0,31, 28, 31, 30, 31, 30, 31, 31, 30, 31, 30, 
31};

Ansonsten der gleiche Code.