Last change
on this file since 1027 was
1027,
checked in by sam, 11 years ago

test: minor Remez algorithm tests.


Property svn:keywords set to
Id

File size:
967 bytes

Line  

1  // 

2  // Lol Engine  Sample math program: Chebyshev polynomials 

3  // 

4  // Copyright: (c) 20052011 Sam Hocevar <sam@hocevar.net> 

5  // This program is free software; you can redistribute it and/or 

6  // modify it under the terms of the Do What The Fuck You Want To 

7  // Public License, Version 2, as published by Sam Hocevar. See 

8  // http://sam.zoy.org/projects/COPYING.WTFPL for more details. 

9  // 

10  

11  #if defined HAVE_CONFIG_H 

12  # include "config.h" 

13  #endif 

14  

15  #include <cstring> 

16  #include <cstdio> 

17  

18  #include "core.h" 

19  

20  using namespace lol; 

21  using namespace std; 

22  

23  #include "remezmatrix.h" 

24  #include "remezsolver.h" 

25  

26  /* The function we want to approximate */ 

27  real myfun(real const &y) 

28  { 

29  real x = sqrt(y); 

30  return (sin(x)  x) / (x * y); 

31  } 

32  

33  real myerr(real const &y) 

34  { 

35  real x = sqrt(y); 

36  return sin(x) / (x * y); 

37  } 

38  

39  int main(void) 

40  { 

41  RemezSolver<6> solver; 

42  solver.Run(real::R_1 >> 400, real::R_PI_2 * real::R_PI_2, myfun, myerr, 40); 

43  

44  return EXIT_SUCCESS; 

45  } 

46  

Note: See
TracBrowser
for help on using the repository browser.