{"id":290757,"date":"2018-10-14T19:15:02","date_gmt":"2018-10-14T15:15:02","guid":{"rendered":"http:\/\/savepearlharbor.com\/?p=290757"},"modified":"-0001-11-30T00:00:00","modified_gmt":"-0001-11-29T21:00:00","slug":"","status":"publish","type":"post","link":"https:\/\/savepearlharbor.com\/?p=290757","title":{"rendered":"\u0424\u0443\u043d\u043a\u0446\u0438\u044f Math.Sin (double) \u0434\u043b\u044f GPU"},"content":{"rendered":"\n
\n

\u041f\u0440\u0435\u0434\u0438\u0441\u043b\u043e\u0432\u0438\u0435<\/h3>\n

\u041c\u043d\u0435 \u043f\u043e\u043d\u0430\u0434\u043e\u0431\u0438\u043b\u043e\u0441\u044c \u0432\u044b\u0447\u0438\u0441\u043b\u044f\u0442\u044c \u0434\u0443\u0433\u0443 \u0441 \u043f\u043e\u0432\u044b\u0448\u0435\u043d\u043d\u043e\u0439 \u0442\u043e\u0447\u043d\u043e\u0441\u0442\u044c\u044e \u043d\u0430 \u043f\u0440\u043e\u0446\u0435\u0441\u0441\u043e\u0440\u0435 \u0432\u0438\u0434\u0435\u043e\u043a\u0430\u0440\u0442\u044b \u0432 \u0440\u0435\u0436\u0438\u043c\u0435 \u0440\u0435\u0430\u043b\u044c\u043d\u043e\u0433\u043e \u0432\u0440\u0435\u043c\u0435\u043d\u0438. <\/p>\n

\u0410\u0432\u0442\u043e\u0440 \u043d\u0435 \u0441\u0442\u0430\u0432\u0438\u043b \u043f\u0435\u0440\u0435\u0434 \u0441\u043e\u0431\u043e\u0439 \u0446\u0435\u043b\u044c \u043f\u0440\u0435\u0432\u0437\u043e\u0439\u0442\u0438 \u0441\u0442\u0430\u043d\u0434\u0430\u0440\u0442\u043d\u0443\u044e \u0444\u0443\u043d\u043a\u0446\u0438\u044e System.Math.Sin() (C#) \u0438 \u0435\u0435 \u043d\u0435 \u0434\u043e\u0441\u0442\u0438\u0433.
<\/a>
\u0420\u0435\u0437\u0443\u043b\u044c\u0442\u0430\u0442 \u0440\u0430\u0431\u043e\u0442\u044b \u0438 \u043c\u043e\u0439 \u0432\u044b\u0431\u043e\u0440 (\u0434\u043b\u044f \u0442\u0435\u0445, \u043a\u043e\u043c\u0443 \u043d\u0435 \u0445\u043e\u0447\u0435\u0442\u0441\u044f \u0447\u0438\u0442\u0430\u0442\u044c):<\/i><\/p>\n

Sin_3(rad)<\/b><\/p>\n
\n
using System; class Math_d {     const double PI025      = Math.PI \/ 4;     \/\/\/ <summary> 2^17 = 131072 (1 \u043c\u0431), \u043e\u0448\u0438\u0431\u043a\u0430 \u043c\u0435\u043d\u044c\u0448\u0435 10000 (\u0441 \u043a\u043e\u043d\u0446\u0430), \u043f\u0440\u0438 2^21 = 22097152 (16 \u043c\u0431) \u043c\u0438\u043d\u0438\u043c\u0430\u043b\u044c\u043d\u0430\u044f \u043e\u0448\u0438\u0431\u043a\u0430 +-1 (\u0441 \u043a\u043e\u043d\u0446\u0430) (\u043f\u043e\u0447\u0442\u0438 \u043d\u0435\u0442 \u043e\u0448\u0438\u0431\u043a\u0438) <\/summary>     const int length_mem    = 22097152;     const double step_rad   = PI025 \/ length_mem;     \/\/\/ <summary> \u041c\u0430\u0441\u0441\u0438\u0432 \u043f\u0440\u0435\u0434\u0432\u0430\u0440\u0438\u0442\u0435\u043b\u044c\u043d\u043e \u0440\u0430\u0441\u0441\u0447\u0438\u0442\u0430\u043d\u043d\u043e\u0433\u043e sin, \u0438\u0441\u043f\u043e\u043b\u044c\u0437\u0443\u0435\u043c\u044b\u0439 \u0434\u043b\u044f \u043b\u0438\u043d\u0435\u0439\u043d\u043e\u0439 \u0438\u043d\u0442\u0435\u0440\u043f\u043e\u043b\u044f\u0446\u0438\u0438. <\/summary>     static double[] mem_sin;     \/\/\/ <summary> \u041c\u0430\u0441\u0441\u0438\u0432 \u043f\u0440\u0435\u0434\u0432\u0430\u0440\u0438\u0442\u0435\u043b\u044c\u043d\u043e \u0440\u0430\u0441\u0441\u0447\u0438\u0442\u0430\u043d\u043d\u043e\u0433\u043e cos, \u0438\u0441\u043f\u043e\u043b\u044c\u0437\u0443\u0435\u043c\u044b\u0439 \u0434\u043b\u044f \u043b\u0438\u043d\u0435\u0439\u043d\u043e\u0439 \u0438\u043d\u0442\u0435\u0440\u043f\u043e\u043b\u044f\u0446\u0438\u0438. <\/summary>     static double[] mem_cos;      \/\/\/ <summary> \u0418\u043d\u0438\u0446\u0438\u0438\u0440\u0443\u0435\u0442 \u0434\u0430\u043d\u043d\u044b\u0435, \u0432\u0440\u043e\u0434\u0435 \u043c\u0430\u0441\u0441\u0438\u0432\u043e\u0432 sin, \u043d\u0435\u043e\u0431\u0445\u043e\u0434\u0438\u043c\u044b\u0435 \u0432 \u0445\u043e\u0434\u0435 \u0440\u0430\u0441\u0447\u0435\u0442\u043e\u0432. <\/summary>     public static void Initialise()     {         Ini_Mem_Sin();         Ini_Mem_Cos();     }     public static double Sin_3(double rad)     {         double rad_025;         int i;          i = (int)(rad \/ PI025);     \/\/\u041d\u0430\u0445\u043e\u0434\u0438\u043c \u0438\u043d\u0434\u0435\u043a\u0441 \u0447\u0430\u0441\u0442\u0438         rad_025 = rad - PI025 * i;  \/\/\u041d\u0430\u0445\u043e\u0434\u0438\u043c \u043e\u0442\u0441\u0442\u0443\u043f \u043e\u0442 \u043d\u0430\u0447\u0430\u043b\u0430 \u0447\u0430\u0441\u0442\u0438 (\u0432 \u0440\u0430\u0434\u0438\u0430\u043d\u0430\u0445)         i = i & 7;                  \/\/\u041d\u0430\u0445\u043e\u0434\u0438\u043c \u043e\u0441\u0442\u0430\u0442\u043e\u043a \u043e\u0442 \u0434\u0435\u043b\u0435\u043d\u0438\u044f \u0438\u043d\u0434\u0435\u043a\u0441\u0430 \u043d\u0430 8          \/\/\u0418\u0449\u0435\u043c \u0438\u043d\u0434\u0435\u043a\u0441 \u043a\u0443\u0441\u043e\u0447\u043a\u0430 \u043a\u0440\u0443\u0433\u0430         switch (i)         {             case 0:                 return Sin_Lerp(rad_025);             case 1:                 return Cos_Lerp(PI025 - rad_025);             case 2:                 return Cos_Lerp(rad_025);             case 3:                 return Sin_Lerp(PI025 - rad_025);             case 4:                 return -Sin_Lerp(rad_025);             case 5:                 return -Cos_Lerp(PI025 - rad_025);             case 6:                 return -Cos_Lerp(rad_025);             case 7:                 return -Sin_Lerp(PI025 - rad_025);         }          return 0;     }      \/\/\/ <summary> \u041f\u043e\u0434\u0433\u043e\u0442\u0430\u0432\u043b\u0438\u0432\u0430\u0435\u0442 \u043c\u0430\u0441\u0441\u0438\u0432 sin \u0434\u043b\u044f \u043b\u0438\u043d\u0435\u0439\u043d\u043e\u0439 \u0438\u043d\u0442\u0435\u0440\u043f\u043e\u043b\u044f\u0446\u0438\u0438 <\/summary>     static void Ini_Mem_Sin()     {         double rad;         mem_sin = new double[length_mem];         for (int i = 0; i < length_mem; i++)         {             rad = (i * PI025) \/ length_mem;             mem_sin[i] = Math.Sin(rad);         }     }     \/\/\/ <summary> \u041f\u043e\u0434\u0433\u043e\u0442\u0430\u0432\u043b\u0438\u0432\u0430\u0435\u0442 \u043c\u0430\u0441\u0441\u0438\u0432 cos \u0434\u043b\u044f \u043b\u0438\u043d\u0435\u0439\u043d\u043e\u0439 \u0438\u043d\u0442\u0435\u0440\u043f\u043e\u043b\u044f\u0446\u0438\u0438 <\/summary>     static void Ini_Mem_Cos()     {         double rad;         mem_cos = new double[length_mem];         for (int i = 0; i < length_mem; i++)         {             rad = (i * PI025) \/ length_mem;             mem_cos[i] = Math.Cos(rad);         }     }     \/\/\/ <summary> \u0418\u0441\u043f\u043e\u043b\u044c\u0437\u0443\u0435\u0442 \u043b\u0438\u043d\u0435\u0439\u043d\u0443\u044e \u0438\u043d\u0442\u0435\u0440\u043f\u043e\u043b\u044f\u0446\u0438\u044e \u0434\u043b\u044f \u043f\u043e\u0438\u0441\u043a\u0430 sin \u043e\u0442 0 \u0434\u043e pi\/4. <\/summary>     \/\/\/ <param name="rad"> \u0423\u0433\u043e\u043b \u0432 \u0440\u0430\u0434\u0438\u0430\u043d\u0430\u0445 \u043e\u0442 0 \u0434\u043e pi\/4. <\/param>     static double Sin_Lerp(double rad)     {         int i_0;         int i_1;         double i_0d;         double percent;         double a;         double b;         double s;          percent = rad \/ PI025;         i_0d    = percent * length_mem;         i_0     = (int)i_0d;         i_1     = i_0 + 1;         a       = mem_sin[i_0];         b       = mem_sin[i_1];         s       = i_0d - i_0;          return Math_d.Lerp(a, b, s);     }     \/\/\/ <summary> \u0418\u0441\u043f\u043e\u043b\u044c\u0437\u0443\u0435\u0442 \u043b\u0438\u043d\u0435\u0439\u043d\u0443\u044e \u0438\u043d\u0442\u0435\u0440\u043f\u043e\u043b\u044f\u0446\u0438\u044e \u0434\u043b\u044f \u043f\u043e\u0438\u0441\u043a\u0430 cos \u043e\u0442 0 \u0434\u043e pi\/4. <\/summary>     \/\/\/ <param name="rad"> \u0423\u0433\u043e\u043b \u0432 \u0440\u0430\u0434\u0438\u0430\u043d\u0430\u0445 \u043e\u0442 0 \u0434\u043e pi\/4. <\/param>     static double Cos_Lerp(double rad)     {         int i_0;         int i_1;         double i_0d;         double percent;         double a;         double b;         double s;          percent = rad \/ PI025;         i_0d = percent * length_mem;         i_0 = (int)i_0d;         i_1 = i_0 + 1;         a = mem_cos[i_0];         b = mem_cos[i_1];         s = i_0d - i_0;          return Math_d.Lerp(a, b, s);     }     \/\/\/ <summary> \u041f\u0440\u043e\u0438\u0437\u0432\u043e\u0434\u0438\u0442 \u043b\u0438\u043d\u0435\u0439\u043d\u0443\u044e \u0438\u043d\u0442\u0435\u0440\u043f\u043e\u043b\u044f\u0446\u0438\u044e \u043c\u0435\u0436\u0434\u0443 \u0434\u0432\u0443\u043c\u044f \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u044f\u043c\u0438. (return a + s * (b - a)) <\/summary>     \/\/\/ <param name="a"> \u041d\u0430\u0447\u0430\u043b\u044c\u043d\u043e\u0435 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0435. <\/param>     \/\/\/ <param name="b"> \u041a\u043e\u043d\u0435\u0447\u043d\u043e\u0435 \u0437\u043d\u0430\u0447\u0435\u043d\u0438\u0435. <\/param>     \/\/\/ <param name="s"> \u041f\u0440\u043e\u0446\u0435\u043d\u0442 \u0438\u043d\u0442\u0435\u0440\u043f\u043e\u043b\u044f\u0446\u0438\u0438. 0 = a, 1 = b, 0.5 = \u0441\u0435\u0440\u0435\u0434\u0438\u043d\u0430 \u043c\u0435\u0436\u0434\u0443 a \u0438 b. <\/param>     public static double Lerp(double a, double b, double s)     {         return a + s * (b - a);     } } <\/code><\/pre>\n
  <\/div>\n<\/div>\n
  <\/p>\n
\u041f\u0440\u0438\u0447\u0438\u043d\u044b \u043f\u0443\u0431\u043b\u0438\u043a\u0430\u0446\u0438\u0438<\/h3>\n
  <\/p>\n