Stilllegung des Forums
Das Forum wurde am 05.06.2023 nach über 20 Jahren stillgelegt (weitere Informationen und ein kleiner Rückblick).
Registrierungen, Anmeldungen und Postings sind nicht mehr möglich. Öffentliche Inhalte sind weiterhin zugänglich.
Das Team von spieleprogrammierer.de bedankt sich bei der Community für die vielen schönen Jahre.
Wenn du eine deutschsprachige Spieleentwickler-Community suchst, schau doch mal im Discord und auf ZFX vorbei!
Werbeanzeige
C-/C++-Quelltext |
|
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 |
#include <iostream> #include <ctime> #include <cmath> using namespace std; float FastSin0(const float& fAngle) { float fASqr = fAngle*fAngle; float fResult = 7.61e-03f; fResult *= fASqr; fResult -= 1.6605e-01f; fResult *= fASqr; fResult += 1.0f; fResult *= fAngle; return fResult; } float FastSin1(const float& fAngle) { float fASqr = fAngle*fAngle; float fResult = -2.39e-08; fResult *= fASqr; fResult += 2.7526e-06f; fResult *= fASqr; fResult -= 1.98409e-04f; fResult *= fASqr; fResult += 8.3333315e-03f; fResult *= fASqr; fResult -= 1.666666664e-01f; fResult *= fASqr; fResult += 1.0f; fResult *= fAngle; return fResult; } int fak(int n) { int x=1; while(n>1) { x*=n; n--; } return x; } float TnSin(float x, unsigned short n) { float Tn = 0; for (unsigned short k = 0; k<n; k++) { if(k%2) Tn += (-1.0f / fak(k)) * pow(x, k); else Tn += (1.0f / fak(k)) * pow(x, k); } return Tn; } int main(void) { long start, ende; float a1, a2, a3, a4, a5; float pi = atan(1.0f)*2.0f; //da nur halbes pi korrekt start = clock(); for(float i=0.0f; i<pi; i=i+0.000001f) FastSin0(i); ende = clock(); a1 = float(ende-start); cout<<"FastSin0 fertig"<<endl; start = clock(); for(float i=0.0f; i<pi; i=i+0.000001f) FastSin1(i); ende = clock(); a2 = float(ende-start); cout<<"FastSin1 fertig"<<endl; start = clock(); for(float i=0.0f; i<pi; i=i+0.000001f) sinf(i); ende = clock(); a3 = float(ende-start); cout<<"sinf fertig"<<endl; start = clock(); for(float i=0.0f; i<pi; i=i+0.000001f) sin(i); ende = clock(); a4 = float(ende-start); cout<<"sin fertig"<<endl; start = clock(); for(float i=0.0f; i<pi; i=i+0.000001f) TnSin(i, 5); ende = clock(); a5 = float(ende-start); cout<<"TnSin fertig"<<endl; cout<<"sinf: "<<a3/CLOCKS_PER_SEC<<" entspricht: "<<a3/a3*100<<"%"<<endl; cout<<"sin: "<<a4/CLOCKS_PER_SEC<<" entspricht: "<<a3/a4*100<<"%"<<endl; cout<<"FastSin0: "<<a2/CLOCKS_PER_SEC<<" entspricht: "<<a3/a1*100<<"%"<<endl; cout<<"FastSin1: "<<a1/CLOCKS_PER_SEC<<" entspricht: "<<a3/a2*100<<"%"<<endl; cout<<"TnSin(unoptimiert): "<<a5/CLOCKS_PER_SEC<<" entspricht: "<<a3/a5*100<<"%"<<endl; return 0; } |
Zitat
FastSin0 fertig
FastSin1 fertig
sinf fertig
sin fertig
TnSin fertig
sinf: 0.172 entspricht: 100%
sin: 0.219 entspricht: 78.5388%
FastSin0: 0.234 entspricht: 137.6%
FastSin1: 0.125 entspricht: 73.5043%
TnSin(unoptimiert): 1.531 entspricht: 11.2345%
Drücken Sie eine beliebige Taste . . .
C-/C++-Quelltext |
|
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 |
#include <iostream> #include <ctime> #include <cmath> float pi = atan(1.0f)*2.0f; //da nur halbes pi korrekt using namespace std; float FastSin0(const float& fAngle) { float fASqr = fAngle*fAngle; float fResult = 7.61e-03f; fResult *= fASqr; fResult -= 1.6605e-01f; fResult *= fASqr; fResult += 1.0f; fResult *= fAngle; return fResult; } float FastSin1(const float& fAngle) { float fASqr = fAngle*fAngle; float fResult = -2.39e-08; fResult *= fASqr; fResult += 2.7526e-06f; fResult *= fASqr; fResult -= 1.98409e-04f; fResult *= fASqr; fResult += 8.3333315e-03f; fResult *= fASqr; fResult -= 1.666666664e-01f; fResult *= fASqr; fResult += 1.0f; fResult *= fAngle; return fResult; } float TnSin(float x) { float tri = x*x*x; float Tn = 1.0f*x; Tn -= 1.0f/6.0f * tri; Tn += 1.0f/120.0f * tri * x * x; Tn -= 1.0f/5040.0f * tri * tri * x; return Tn; } int main(void) { long start, ende; float a1, a2, a3, a4, a5; start = clock(); for(int j=0; j<1e3; j++) for(float i=0.0f; i<pi; i=i+0.0001f) FastSin0(i); ende = clock(); a1 = float(ende-start); cout<<"FastSin0 fertig"<<endl; start = clock(); for(int j=0; j<1e3; j++) for(float i=0.0f; i<pi; i=i+0.0001f) FastSin1(i); ende = clock(); a2 = float(ende-start); cout<<"FastSin1 fertig"<<endl; start = clock(); for(int j=0; j<1e3; j++) for(float i=0.0f; i<pi; i=i+0.0001f) sinf(i); ende = clock(); a3 = float(ende-start); cout<<"sinf fertig"<<endl; start = clock(); for(int j=0; j<1e3; j++) for(float i=0.0f; i<pi; i=i+0.0001f) sin(i); ende = clock(); a4 = float(ende-start); cout<<"sin fertig"<<endl; start = clock(); for(int j=0; j<1e3; j++) for(float i=0.0f; i<pi; i=i+0.0001f) TnSin(i); ende = clock(); a5 = float(ende-start); cout<<"TnSin fertig"<<endl; start = clock(); cout<<"sinf: "<<a3/CLOCKS_PER_SEC<<" entspricht: "<<a3/a3*100<<"%"<<endl; cout<<"sin: "<<a4/CLOCKS_PER_SEC<<" entspricht: "<<a3/a4*100<<"%"<<endl; cout<<"FastSin0: "<<a2/CLOCKS_PER_SEC<<" entspricht: "<<a3/a1*100<<"%"<<endl; cout<<"FastSin1: "<<a1/CLOCKS_PER_SEC<<" entspricht: "<<a3/a2*100<<"%"<<endl; cout<<"TnSin: "<<a5/CLOCKS_PER_SEC<<" entspricht: "<<a3/a5*100<<"%"<<endl; return 0; } |
Zitat
FastSin0 fertig
FastSin1 fertig
sinf fertig
sin fertig
TnSin fertig
sinf: 1.656 entspricht: 100%
sin: 2.218 entspricht: 74.6619%
FastSin0: 2.078 entspricht: 121.765%
FastSin1: 1.36 entspricht: 79.692%
TnSin: 1.125 entspricht: 147.2%
Drücken Sie eine beliebige Taste . . .
C-/C++-Quelltext |
|
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 |
#include <iostream> #include <ctime> #include <cmath> float pi = atan(1.0f)*2.0f; //da nur halbes pi korrekt using namespace std; inline float FastSin0(const float& fAngle) { float fASqr = fAngle*fAngle; float fResult = 7.61e-03f; fResult *= fASqr; fResult -= 1.6605e-01f; fResult *= fASqr; fResult += 1.0f; fResult *= fAngle; return fResult; } inline float FastSin1(const float& fAngle) { float fASqr = fAngle*fAngle; float fResult = -2.39e-08; fResult *= fASqr; fResult += 2.7526e-06f; fResult *= fASqr; fResult -= 1.98409e-04f; fResult *= fASqr; fResult += 8.3333315e-03f; fResult *= fASqr; fResult -= 1.666666664e-01f; fResult *= fASqr; fResult += 1.0f; fResult *= fAngle; return fResult; } inline float TnSin(const float& x) { float tri = x*x*x; float Tn = 1.0f*x; Tn -= 1.0f/6.0f * tri; Tn += 1.0f/120.0f * tri * x * x; Tn -= 1.0f/5040.0f * tri * tri * x; Tn += 1.0f/362880.0f * tri * tri * tri; return Tn; } int main(void) { long start, ende; float a1, a2, a3, a4, a5; start = clock(); for(int j=0; j<1e3; j++) for(float i=0.0f; i<pi; i=i+0.0001f) FastSin0(i); ende = clock(); a1 = float(ende-start); cout<<"FastSin0 fertig"<<endl; start = clock(); for(int j=0; j<1e3; j++) for(float i=0.0f; i<pi; i=i+0.0001f) FastSin1(i); ende = clock(); a2 = float(ende-start); cout<<"FastSin1 fertig"<<endl; start = clock(); for(int j=0; j<1e3; j++) for(float i=0.0f; i<pi; i=i+0.0001f) sinf(i); ende = clock(); a3 = float(ende-start); cout<<"sinf fertig"<<endl; start = clock(); for(int j=0; j<1e3; j++) for(float i=0.0f; i<pi; i=i+0.0001f) sin(i); ende = clock(); a4 = float(ende-start); cout<<"sin fertig"<<endl; start = clock(); for(int j=0; j<1e3; j++) for(float i=0.0f; i<pi; i=i+0.0001f) TnSin(i); ende = clock(); a5 = float(ende-start); cout<<"TnSin fertig"<<endl; start = clock(); cout<<"sinf: "<<a3/CLOCKS_PER_SEC<<" entspricht: "<<a3/a3*100<<"%"<<endl; cout<<"sin: "<<a4/CLOCKS_PER_SEC<<" entspricht: "<<a3/a4*100<<"%"<<endl; cout<<"FastSin0: "<<a2/CLOCKS_PER_SEC<<" entspricht: "<<a3/a1*100<<"%"<<endl; cout<<"FastSin1: "<<a1/CLOCKS_PER_SEC<<" entspricht: "<<a3/a2*100<<"%"<<endl; cout<<"TnSin: "<<a5/CLOCKS_PER_SEC<<" entspricht: "<<a3/a5*100<<"%"<<endl; float p = 0.0f; for(float i=0.0f; i<pi; i+= 0.0001f) { float j = sin(i) - TnSin(i); if(p<j) p = j; } cout<<scientific<<"Genauigkeit TnSin: "<<p<<endl; p = 0.0f; for(float i=0.0f; i<pi; i+= 0.0001f) { float j = sin(i) - FastSin0(i); if(p<j) p = j; } cout<<"Genauigkeit FastSin0: "<<p<<endl; p = 0.0f; for(float i=0.0f; i<pi; i+= 0.0001f) { float j = sin(i) - FastSin1(i); if(p<j) p = j; } cout<<"Genauigkeit FastSin1: "<<p<<endl; return 0; } |
Zitat
FastSin0 fertig
FastSin1 fertig
sinf fertig
sin fertig
TnSin fertig
sinf: 1.578 entspricht: 100%
sin: 2.093 entspricht: 75.3942%
FastSin0: 2.047 entspricht: 118.825%
FastSin1: 1.328 entspricht: 77.0884%
TnSin: 1.313 entspricht: 120.183%
Genauigkeit TnSin: 1.192093e-007
Genauigkeit FastSin0: 1.642108e-004
Genauigkeit FastSin1: 1.192093e-007
Drücken Sie eine beliebige Taste . . .
Werbeanzeige