17 lines
266 B
C
17 lines
266 B
C
|
#include <fixpoint/fp4816.h>
|
||
|
|
||
|
fp4816_t fp4816_taylor(fp4816_t v,int lev) {
|
||
|
int i;
|
||
|
fp4816_t result = v;
|
||
|
for (i=0;i<lev;i++) {
|
||
|
fp4816_t te = fp4816_taylor_element( v, i+1 );
|
||
|
if (i & 1) {
|
||
|
result += te;
|
||
|
} else {
|
||
|
result -= te;
|
||
|
};
|
||
|
};
|
||
|
return result;
|
||
|
};
|
||
|
|