asin durch inverse LUT?

[XBert]

Hi,
ich hab mich mal zur berechnung des cos/sin an das folgende tutorial gehalten und
würde jetzt gerne wissen wie ich mit der selben tabelle die inversen funktionen (arcus sinus reicht da ja acos=90°-asin) dazu berechne.
Hier das Tutorial: http://www.mikrocontroller.net/articles/AVR_Arithmetik/Sinus_und_Cosinus_(Lineare_Interpolation)

und hier mein Versuch:

int16_t linasin(int16_t input)
{
	uint8_t neg = input<0?1:0;

	input &= ~(1 << 15); //=abs()
	for (uint8_t i = 0; i < 91; i++) {
		if (input >= linsin_tab[i] && input < linsin_tab[i +1]) {
			uint16_t data0 = linsin_tab[i+1] - linsin_tab[i];
			uint16_t data1 = input - linsin_tab[i];
			int16_t res = i * 256 + 256 * data1 / data0;
			return neg ? -res :res;
		}
	}
}

MFG XBert

[XBert]

Hi,
ich hab das ganze doch noch hinbekommen; falls jemand ebenfalls darauf stoßen sollte, hier ist des Source-Code;
allerdings ist noch zu klären warum das Ergebnis an manchen stellen so weit vom soll-wert abweicht; (+-1 ist ok da ich durch 2^15 dividiere und nicht durch 2^15-1);
Ausgabe:

              std     |     diy
sin(30°):    16384    |    16385
asin(0.5):   7680     |    7547
asin(0.7):   11373    |    11373
asin(0.3):    4469    |    4468

und der Source-Code:

nt16_t linasin(int16_t input)
{
	uint8_t neg = input<0?1:0;
	uint16_t _input = abs(input);
	//input &= ~(1 << 15); //=abs()
	for (uint8_t i = 0; i < 91; i++) { 
		uint16_t data1 = pgm_read_word(&linsin_tab[i]);
		uint16_t data0 = pgm_read_word(&linsin_tab[i+1]);
		if (_input >= data1 && _input < data0) {
			data0 -= data1;
			data1 = _input - data1;
			int16_t res = i * 256 + 256 * data1 / data0;
			return neg ? -res :res;
		}
	}
	return 0;
}

MFG