• 设为首页
  • 点击收藏
  • 手机版
    手机扫一扫访问
    迪恩网络手机版
  • 关注官方公众号
    微信扫一扫关注
    公众号

C++ FMS函数代码示例

原作者: [db:作者] 来自: [db:来源] 收藏 邀请

本文整理汇总了C++中FMS函数的典型用法代码示例。如果您正苦于以下问题:C++ FMS函数的具体用法?C++ FMS怎么用?C++ FMS使用的例子?那么恭喜您, 这里精选的函数代码示例或许可以为您提供帮助。



在下文中一共展示了FMS函数的20个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于我们的系统推荐出更棒的C++代码示例。

示例1: r2cbIII_5

static void r2cbIII_5(R *R0, R *R1, R *Cr, R *Ci, stride rs, stride csr, stride csi, INT v, INT ivs, INT ovs)
{
     DK(KP1_902113032, +1.902113032590307144232878666758764286811397268);
     DK(KP1_118033988, +1.118033988749894848204586834365638117720309180);
     DK(KP500000000, +0.500000000000000000000000000000000000000000000);
     DK(KP2_000000000, +2.000000000000000000000000000000000000000000000);
     DK(KP618033988, +0.618033988749894848204586834365638117720309180);
     {
	  INT i;
	  for (i = v; i > 0; i = i - 1, R0 = R0 + ovs, R1 = R1 + ovs, Cr = Cr + ivs, Ci = Ci + ivs, MAKE_VOLATILE_STRIDE(rs), MAKE_VOLATILE_STRIDE(csr), MAKE_VOLATILE_STRIDE(csi)) {
	       E T1, T2, T3, Tc, Ta, T8, T9;
	       T8 = Ci[WS(csi, 1)];
	       T9 = Ci[0];
	       T1 = Cr[WS(csr, 2)];
	       T2 = Cr[WS(csr, 1)];
	       T3 = Cr[0];
	       Tc = FMS(KP618033988, T8, T9);
	       Ta = FMA(KP618033988, T9, T8);
	       {
		    E T6, T4, T5, T7, Tb;
		    T6 = T3 - T2;
		    T4 = T2 + T3;
		    R0[0] = FMA(KP2_000000000, T4, T1);
		    T5 = FNMS(KP500000000, T4, T1);
		    T7 = FNMS(KP1_118033988, T6, T5);
		    Tb = FMA(KP1_118033988, T6, T5);
		    R0[WS(rs, 2)] = FNMS(KP1_902113032, Ta, T7);
		    R1[0] = -(FMA(KP1_902113032, Ta, T7));
		    R1[WS(rs, 1)] = FMS(KP1_902113032, Tc, Tb);
		    R0[WS(rs, 1)] = FMA(KP1_902113032, Tc, Tb);
	       }
	  }
     }
}
开发者ID:BGCX261,项目名称:zpr-rozpoznawanie-dzwieku-git,代码行数:34,代码来源:r2cbIII_5.c


示例2: hc2rIII_5

static void hc2rIII_5(const R *ri, const R *ii, R *O, stride ris, stride iis, stride os, INT v, INT ivs, INT ovs)
{
     DK(KP1_902113032, +1.902113032590307144232878666758764286811397268);
     DK(KP1_118033988, +1.118033988749894848204586834365638117720309180);
     DK(KP500000000, +0.500000000000000000000000000000000000000000000);
     DK(KP2_000000000, +2.000000000000000000000000000000000000000000000);
     DK(KP618033988, +0.618033988749894848204586834365638117720309180);
     INT i;
     for (i = v; i > 0; i = i - 1, ri = ri + ivs, ii = ii + ivs, O = O + ovs, MAKE_VOLATILE_STRIDE(ris), MAKE_VOLATILE_STRIDE(iis), MAKE_VOLATILE_STRIDE(os)) {
	  E T1, T2, T3, Tc, Ta, T8, T9;
	  T8 = ii[WS(iis, 1)];
	  T9 = ii[0];
	  T1 = ri[WS(ris, 2)];
	  T2 = ri[WS(ris, 1)];
	  T3 = ri[0];
	  Tc = FMS(KP618033988, T8, T9);
	  Ta = FMA(KP618033988, T9, T8);
	  {
	       E T6, T4, T5, T7, Tb;
	       T6 = T3 - T2;
	       T4 = T2 + T3;
	       O[0] = FMA(KP2_000000000, T4, T1);
	       T5 = FNMS(KP500000000, T4, T1);
	       T7 = FNMS(KP1_118033988, T6, T5);
	       Tb = FMA(KP1_118033988, T6, T5);
	       O[WS(os, 4)] = FNMS(KP1_902113032, Ta, T7);
	       O[WS(os, 1)] = -(FMA(KP1_902113032, Ta, T7));
	       O[WS(os, 3)] = FMS(KP1_902113032, Tc, Tb);
	       O[WS(os, 2)] = FMA(KP1_902113032, Tc, Tb);
	  }
     }
}
开发者ID:exic,项目名称:last.fm-dbus,代码行数:32,代码来源:hc2rIII_5.c


示例3: findMedianSortedArrays

    double findMedianSortedArrays(int A[], int m, int B[], int n) {
        int total = m + n ;
        if(total%2 == 1)
        	return FMS(A,m,B,n,total/2+1);
        else
        	return (FMS(A,m,B,n,total/2+1) + FMS(A,m,B,n,total/2))/2;

    }
开发者ID:qingchen021,项目名称:Leetcode,代码行数:8,代码来源:Median+of+Two+Sorted+Arrays.cpp


示例4: r2cfII_8

static void r2cfII_8(R *R0, R *R1, R *Cr, R *Ci, stride rs, stride csr, stride csi, INT v, INT ivs, INT ovs)
{
     DK(KP923879532, +0.923879532511286756128183189396788286822416626);
     DK(KP707106781, +0.707106781186547524400844362104849039284835938);
     DK(KP414213562, +0.414213562373095048801688724209698078569671875);
     {
	  INT i;
	  for (i = v; i > 0; i = i - 1, R0 = R0 + ivs, R1 = R1 + ivs, Cr = Cr + ovs, Ci = Ci + ovs, MAKE_VOLATILE_STRIDE(rs), MAKE_VOLATILE_STRIDE(csr), MAKE_VOLATILE_STRIDE(csi)) {
	       E Te, T8, Td, T5, Tj, Tl, Tf, Tb;
	       {
		    E T1, Th, T9, Ti, T4, Ta;
		    T1 = R0[0];
		    Th = R0[WS(rs, 2)];
		    {
			 E T2, T3, T6, T7;
			 T2 = R0[WS(rs, 1)];
			 T3 = R0[WS(rs, 3)];
			 T6 = R1[0];
			 T7 = R1[WS(rs, 2)];
			 T9 = R1[WS(rs, 3)];
			 Ti = T2 + T3;
			 T4 = T2 - T3;
			 Te = FMA(KP414213562, T6, T7);
			 T8 = FNMS(KP414213562, T7, T6);
			 Ta = R1[WS(rs, 1)];
		    }
		    Td = FNMS(KP707106781, T4, T1);
		    T5 = FMA(KP707106781, T4, T1);
		    Tj = FMA(KP707106781, Ti, Th);
		    Tl = FNMS(KP707106781, Ti, Th);
		    Tf = FMA(KP414213562, T9, Ta);
		    Tb = FMS(KP414213562, Ta, T9);
	       }
	       {
		    E Tk, Tg, Tc, Tm;
		    Tk = Te + Tf;
		    Tg = Te - Tf;
		    Tc = T8 + Tb;
		    Tm = Tb - T8;
		    Cr[WS(csr, 1)] = FMA(KP923879532, Tg, Td);
		    Cr[WS(csr, 2)] = FNMS(KP923879532, Tg, Td);
		    Ci[WS(csi, 3)] = FNMS(KP923879532, Tk, Tj);
		    Ci[0] = -(FMA(KP923879532, Tk, Tj));
		    Ci[WS(csi, 1)] = FMA(KP923879532, Tm, Tl);
		    Ci[WS(csi, 2)] = FMS(KP923879532, Tm, Tl);
		    Cr[0] = FMA(KP923879532, Tc, T5);
		    Cr[WS(csr, 3)] = FNMS(KP923879532, Tc, T5);
	       }
	  }
     }
}
开发者ID:AlexanderB1290,项目名称:scafacos,代码行数:51,代码来源:r2cfII_8.c


示例5: FMS

 double FMS(int a[], int m , int b[], int n, const int k)
 {
 	if(m > n )
 		return FMS(b,n,a,m,k);
 	if(m == 0)
 		return b[k-1];
 	if(k == 1)
 		return min(a[0],b[0]);
 	int pa = min(m,k/2);
 	int pb = k - pa;
 	if(a[pa-1] < b[pb-1])
 		return FMS(a+pa,m-pa,b,n,k-pa);
 	else
 		return FMS(a,m,b+pb,n-pb,k-pb);
 }
开发者ID:qingchen021,项目名称:Leetcode,代码行数:15,代码来源:Median+of+Two+Sorted+Arrays.cpp


示例6: hc2rIII_5

static void hc2rIII_5(const R *ri, const R *ii, R *O, stride ris, stride iis, stride os, int v, int ivs, int ovs)
{
     DK(KP2_000000000, +2.000000000000000000000000000000000000000000000);
     DK(KP1_118033988, +1.118033988749894848204586834365638117720309180);
     DK(KP500000000, +0.500000000000000000000000000000000000000000000);
     DK(KP1_175570504, +1.175570504584946258337411909278145537195304875);
     DK(KP1_902113032, +1.902113032590307144232878666758764286811397268);
     int i;
     for (i = v; i > 0; i = i - 1, ri = ri + ivs, ii = ii + ivs, O = O + ovs) {
	  E Ta, Tc, T1, T4, T5, T6, Tb, T7;
	  {
	       E T8, T9, T2, T3;
	       T8 = ii[WS(iis, 1)];
	       T9 = ii[0];
	       Ta = FMA(KP1_902113032, T8, KP1_175570504 * T9);
	       Tc = FNMS(KP1_902113032, T9, KP1_175570504 * T8);
	       T1 = ri[WS(ris, 2)];
	       T2 = ri[WS(ris, 1)];
	       T3 = ri[0];
	       T4 = T2 + T3;
	       T5 = FMS(KP500000000, T4, T1);
	       T6 = KP1_118033988 * (T3 - T2);
	  }
	  O[0] = FMA(KP2_000000000, T4, T1);
	  Tb = T6 - T5;
	  O[WS(os, 2)] = Tb + Tc;
	  O[WS(os, 3)] = Tc - Tb;
	  T7 = T5 + T6;
	  O[WS(os, 1)] = T7 - Ta;
	  O[WS(os, 4)] = -(T7 + Ta);
     }
}
开发者ID:OS2World,项目名称:MM-SOUND-PM123,代码行数:32,代码来源:hc2rIII_5.c


示例7: r2cfII_6

static void r2cfII_6(R *R0, R *R1, R *Cr, R *Ci, stride rs, stride csr, stride csi, INT v, INT ivs, INT ovs)
{
     DK(KP866025403, +0.866025403784438646763723170752936183471402627);
     DK(KP500000000, +0.500000000000000000000000000000000000000000000);
     {
	  INT i;
	  for (i = v; i > 0; i = i - 1, R0 = R0 + ivs, R1 = R1 + ivs, Cr = Cr + ovs, Ci = Ci + ovs, MAKE_VOLATILE_STRIDE(24, rs), MAKE_VOLATILE_STRIDE(24, csr), MAKE_VOLATILE_STRIDE(24, csi)) {
	       E T1, T9, T2, T3, T6, T7;
	       T1 = R0[0];
	       T9 = R1[WS(rs, 1)];
	       T2 = R0[WS(rs, 2)];
	       T3 = R0[WS(rs, 1)];
	       T6 = R1[WS(rs, 2)];
	       T7 = R1[0];
	       {
		    E Tc, T4, Ta, T8, T5, Tb;
		    Cr[WS(csr, 1)] = T1 + T2 - T3;
		    Tc = T2 + T3;
		    T4 = T3 - T2;
		    Ta = T6 + T7;
		    T8 = T6 - T7;
		    T5 = FMA(KP500000000, T4, T1);
		    Tb = FMA(KP500000000, Ta, T9);
		    Ci[WS(csi, 1)] = T9 - Ta;
		    Cr[WS(csr, 2)] = FMA(KP866025403, T8, T5);
		    Cr[0] = FNMS(KP866025403, T8, T5);
		    Ci[WS(csi, 2)] = FMS(KP866025403, Tc, Tb);
		    Ci[0] = -(FMA(KP866025403, Tc, Tb));
	       }
	  }
     }
}
开发者ID:SKA-ScienceDataProcessor,项目名称:FastImaging,代码行数:32,代码来源:r2cfII_6.c


示例8: r2hcII_6

static void r2hcII_6(const R *I, R *ro, R *io, stride is, stride ros, stride ios, INT v, INT ivs, INT ovs)
{
     DK(KP866025403, +0.866025403784438646763723170752936183471402627);
     DK(KP500000000, +0.500000000000000000000000000000000000000000000);
     INT i;
     for (i = v; i > 0; i = i - 1, I = I + ivs, ro = ro + ovs, io = io + ovs, MAKE_VOLATILE_STRIDE(is), MAKE_VOLATILE_STRIDE(ros), MAKE_VOLATILE_STRIDE(ios)) {
	  E T1, T9, T2, T3, T6, T7;
	  T1 = I[0];
	  T9 = I[WS(is, 3)];
	  T2 = I[WS(is, 4)];
	  T3 = I[WS(is, 2)];
	  T6 = I[WS(is, 5)];
	  T7 = I[WS(is, 1)];
	  {
	       E Tc, T4, Ta, T8, T5, Tb;
	       ro[WS(ros, 1)] = T1 + T2 - T3;
	       Tc = T2 + T3;
	       T4 = T3 - T2;
	       Ta = T6 + T7;
	       T8 = T6 - T7;
	       T5 = FMA(KP500000000, T4, T1);
	       Tb = FMA(KP500000000, Ta, T9);
	       io[WS(ios, 1)] = T9 - Ta;
	       ro[WS(ros, 2)] = FMA(KP866025403, T8, T5);
	       ro[0] = FNMS(KP866025403, T8, T5);
	       io[WS(ios, 2)] = FMS(KP866025403, Tc, Tb);
	       io[0] = -(FMA(KP866025403, Tc, Tb));
	  }
     }
}
开发者ID:exic,项目名称:last.fm-dbus,代码行数:30,代码来源:r2hcII_6.c


示例9: r2cbIII_7

static void r2cbIII_7(R *R0, R *R1, R *Cr, R *Ci, stride rs, stride csr, stride csi, INT v, INT ivs, INT ovs)
{
     DK(KP1_949855824, +1.949855824363647214036263365987862434465571601);
     DK(KP1_801937735, +1.801937735804838252472204639014890102331838324);
     DK(KP2_000000000, +2.000000000000000000000000000000000000000000000);
     DK(KP692021471, +0.692021471630095869627814897002069140197260599);
     DK(KP801937735, +0.801937735804838252472204639014890102331838324);
     DK(KP356895867, +0.356895867892209443894399510021300583399127187);
     DK(KP554958132, +0.554958132087371191422194871006410481067288862);
     {
	  INT i;
	  for (i = v; i > 0; i = i - 1, R0 = R0 + ovs, R1 = R1 + ovs, Cr = Cr + ivs, Ci = Ci + ivs, MAKE_VOLATILE_STRIDE(rs), MAKE_VOLATILE_STRIDE(csr), MAKE_VOLATILE_STRIDE(csi)) {
	       E Tn, Td, Tg, Ti, Tl, T8;
	       {
		    E T1, T9, Tb, Ta, T2, T4, Th, Tm, Tc, T3, Te;
		    T1 = Cr[WS(csr, 3)];
		    T9 = Ci[WS(csi, 1)];
		    Tb = Ci[0];
		    Ta = Ci[WS(csi, 2)];
		    T2 = Cr[WS(csr, 2)];
		    T4 = Cr[0];
		    Th = FMA(KP554958132, T9, Tb);
		    Tm = FNMS(KP554958132, Ta, T9);
		    Tc = FMA(KP554958132, Tb, Ta);
		    T3 = Cr[WS(csr, 1)];
		    Te = FNMS(KP356895867, T2, T4);
		    Tn = FNMS(KP801937735, Tm, Tb);
		    {
			 E Tf, Tk, T7, T5, Tj, T6;
			 Td = FMA(KP801937735, Tc, T9);
			 T5 = T2 + T3 + T4;
			 Tj = FNMS(KP356895867, T4, T3);
			 T6 = FNMS(KP356895867, T3, T2);
			 Tf = FNMS(KP692021471, Te, T3);
			 R0[0] = FMA(KP2_000000000, T5, T1);
			 Tk = FNMS(KP692021471, Tj, T2);
			 T7 = FNMS(KP692021471, T6, T4);
			 Tg = FNMS(KP1_801937735, Tf, T1);
			 Ti = FNMS(KP801937735, Th, Ta);
			 Tl = FNMS(KP1_801937735, Tk, T1);
			 T8 = FNMS(KP1_801937735, T7, T1);
		    }
	       }
	       R1[WS(rs, 2)] = FMS(KP1_949855824, Ti, Tg);
	       R0[WS(rs, 1)] = FMA(KP1_949855824, Ti, Tg);
	       R0[WS(rs, 2)] = FNMS(KP1_949855824, Tn, Tl);
	       R1[WS(rs, 1)] = -(FMA(KP1_949855824, Tn, Tl));
	       R0[WS(rs, 3)] = FNMS(KP1_949855824, Td, T8);
	       R1[0] = -(FMA(KP1_949855824, Td, T8));
	  }
     }
}
开发者ID:BGCX261,项目名称:zpr-rozpoznawanie-dzwieku-git,代码行数:52,代码来源:r2cbIII_7.c


示例10: r2cf_8

static void r2cf_8(R *R0, R *R1, R *Cr, R *Ci, stride rs, stride csr, stride csi, INT v, INT ivs, INT ovs)
{
     DK(KP707106781, +0.707106781186547524400844362104849039284835938);
     {
	  INT i;
	  for (i = v; i > 0; i = i - 1, R0 = R0 + ivs, R1 = R1 + ivs, Cr = Cr + ovs, Ci = Ci + ovs, MAKE_VOLATILE_STRIDE(rs), MAKE_VOLATILE_STRIDE(csr), MAKE_VOLATILE_STRIDE(csi)) {
	       E T4, T7, T3, Tj, Td, T5, T8, T9;
	       {
		    E T1, T2, Tb, Tc;
		    T1 = R0[0];
		    T2 = R0[WS(rs, 2)];
		    Tb = R1[WS(rs, 3)];
		    Tc = R1[WS(rs, 1)];
		    T4 = R0[WS(rs, 1)];
		    T7 = T1 - T2;
		    T3 = T1 + T2;
		    Tj = Tb + Tc;
		    Td = Tb - Tc;
		    T5 = R0[WS(rs, 3)];
		    T8 = R1[0];
		    T9 = R1[WS(rs, 2)];
	       }
	       {
		    E T6, Tf, Ta, Ti;
		    T6 = T4 + T5;
		    Tf = T4 - T5;
		    Ta = T8 - T9;
		    Ti = T8 + T9;
		    {
			 E Th, Tk, Te, Tg;
			 Th = T3 + T6;
			 Cr[WS(csr, 2)] = T3 - T6;
			 Tk = Ti + Tj;
			 Ci[WS(csi, 2)] = Tj - Ti;
			 Te = Ta + Td;
			 Tg = Td - Ta;
			 Cr[0] = Th + Tk;
			 Cr[WS(csr, 4)] = Th - Tk;
			 Ci[WS(csi, 3)] = FMA(KP707106781, Tg, Tf);
			 Ci[WS(csi, 1)] = FMS(KP707106781, Tg, Tf);
			 Cr[WS(csr, 1)] = FMA(KP707106781, Te, T7);
			 Cr[WS(csr, 3)] = FNMS(KP707106781, Te, T7);
		    }
	       }
	  }
     }
}
开发者ID:dstuck,项目名称:tinker_integrated_PIMC,代码行数:47,代码来源:r2cf_8.c


示例11: hc2cf2_20


//.........这里部分代码省略.........
					T4u = FNMS(KP559016994, T4t, T4s);
					T4C = FMA(KP559016994, T4t, T4s);
				   }
				   {
					E T2x, T4Q, T4p, T4r, T4R, T2A, T51, T55;
					{
					     E T4h, T50, T4X, T4o, T4Z;
					     T4E = T4d + T4g;
					     T4h = T4d - T4g;
					     Rm[WS(rs, 1)] = FMA(KP951056516, T4B, T4u);
					     Rp[WS(rs, 2)] = FNMS(KP951056516, T4B, T4u);
					     Rp[WS(rs, 6)] = FMA(KP951056516, T4D, T4C);
					     Rm[WS(rs, 5)] = FNMS(KP951056516, T4D, T4C);
					     T50 = T4W - T4V;
					     T4X = T4V + T4W;
					     T4o = T4k - T4n;
					     T4F = T4k + T4n;
					     T2x = T2v + T2w;
					     T4Q = T2v - T2w;
					     Im[WS(rs, 9)] = T4X - T4Y;
					     T4Z = FMA(KP250000000, T4X, T4Y);
					     T4p = FMA(KP618033988, T4o, T4h);
					     T4r = FNMS(KP618033988, T4h, T4o);
					     T4R = T2z - T2y;
					     T2A = T2y + T2z;
					     T51 = FNMS(KP559016994, T50, T4Z);
					     T55 = FMA(KP559016994, T50, T4Z);
					}
					{
					     E T49, T48, T2B, T4a, T4q;
					     T2B = T2x + T2A;
					     T49 = T2x - T2A;
					     Ip[WS(rs, 2)] = FMA(KP951056516, T54, T51);
					     Im[WS(rs, 1)] = FMS(KP951056516, T54, T51);
					     Ip[WS(rs, 6)] = FMA(KP951056516, T56, T55);
					     Im[WS(rs, 5)] = FMS(KP951056516, T56, T55);
					     Rp[0] = T2u + T2B;
					     T48 = FNMS(KP250000000, T2B, T2u);
					     T4a = FMA(KP559016994, T49, T48);
					     T4q = FNMS(KP559016994, T49, T48);
					     T4U = FMA(KP618033988, T4Q, T4R);
					     T4S = FNMS(KP618033988, T4R, T4Q);
					     Rm[WS(rs, 3)] = FMA(KP951056516, T4p, T4a);
					     Rp[WS(rs, 4)] = FNMS(KP951056516, T4p, T4a);
					     Rp[WS(rs, 8)] = FMA(KP951056516, T4r, T4q);
					     Rm[WS(rs, 7)] = FNMS(KP951056516, T4r, T4q);
					}
				   }
			      }
			      {
				   E T3O, T5u, T5w, T5o, T5q, T5n;
				   {
					E T5m, T5l, T2I, T4O, T3N, T3L, T2X, T5s, T4N, T5t, T3c, T3v, T3K, T4G;
					T5m = T3u + T3n;
					T3v = T3n - T3u;
					T3K = T3C - T3J;
					T5l = T3J + T3C;
					T3O = T2C + T2H;
					T2I = T2C - T2H;
					T4O = T4E - T4F;
					T4G = T4E + T4F;
					T3N = FMA(KP618033988, T3v, T3K);
					T3L = FNMS(KP618033988, T3K, T3v);
					T2X = T2P + T2W;
					T5s = T2P - T2W;
					Ip[0] = T4G + T4M;
开发者ID:Aegisub,项目名称:fftw3,代码行数:67,代码来源:hc2cf2_20.c


示例12: DK


//.........这里部分代码省略.........
					T1V = FNMS(KP866025403, T1U, T1P);
					T3c = FNMS(KP618033988, T3b, T38);
					T3e = FMA(KP618033988, T38, T3b);
					T3q = FNMS(KP618033988, T3p, T3o);
					T3s = FMA(KP618033988, T3o, T3p);
					T1A = TU + T1z;
					T34 = TU - T1z;
				   }
			      }
			      {
				   E T2W, T33, T3m, T3h, T2Z, T3d, T35, T3l;
				   T3J = T2U + T2V;
				   T2W = T2U - T2V;
				   rio[0] = Tf + T1A;
				   T33 = FNMS(KP250000000, T1A, Tf);
				   T3m = T3f - T3g;
				   T3h = T3f + T3g;
				   T2Z = T2X - T2Y;
				   T3I = T2X + T2Y;
				   T3d = FMA(KP559016994, T34, T33);
				   T35 = FNMS(KP559016994, T34, T33);
				   iio[0] = T3h + T3k;
				   T3l = FNMS(KP250000000, T3h, T3k);
				   rio[WS(ios, 3)] = FMA(KP951056516, T3c, T35);
				   iio[-WS(ios, 12)] = FNMS(KP951056516, T3c, T35);
				   rio[WS(ios, 6)] = FMA(KP951056516, T3e, T3d);
				   iio[-WS(ios, 9)] = FNMS(KP951056516, T3e, T3d);
				   T3r = FMA(KP559016994, T3m, T3l);
				   T3n = FNMS(KP559016994, T3m, T3l);
				   T32 = FMA(KP618033988, T2W, T2Z);
				   T30 = FNMS(KP618033988, T2Z, T2W);
			      }
			      iio[-WS(ios, 3)] = FMA(KP951056516, T3q, T3n);
			      rio[WS(ios, 12)] = FMS(KP951056516, T3q, T3n);
			      iio[-WS(ios, 6)] = FMA(KP951056516, T3s, T3r);
			      rio[WS(ios, 9)] = FMS(KP951056516, T3s, T3r);
			      T2d = FMA(KP866025403, T1G, T1B);
			      T1H = FNMS(KP866025403, T1G, T1B);
			      T1W = T1O + T1V;
			      T3O = T1O - T1V;
			      T3P = T23 - T2a;
			      T2b = T23 + T2a;
			 }
		    }
		    {
			 E T3H, T3y, T2S, T3Q, T3S, T2R, T2c;
			 T3H = FNMS(KP866025403, T3x, T3w);
			 T3y = FMA(KP866025403, T3x, T3w);
			 T2c = T1W + T2b;
			 T2S = T1W - T2b;
			 T3Q = FMA(KP618033988, T3P, T3O);
			 T3S = FNMS(KP618033988, T3O, T3P);
			 rio[WS(ios, 5)] = T1H + T2c;
			 T2R = FNMS(KP250000000, T2c, T1H);
			 {
			      E T2g, T2j, T3G, T3E, T2A, T2N, T3v, T3A, T3M, T3L, T3z, T3F, T3B;
			      {
				   E T3C, T3D, T31, T2T, T3K;
				   T2g = T2e + T2f;
				   T3C = T2e - T2f;
				   T3D = T2h - T2i;
				   T2j = T2h + T2i;
				   T31 = FMA(KP559016994, T2S, T2R);
				   T2T = FNMS(KP559016994, T2S, T2R);
				   T3K = T3I + T3J;
				   T3M = T3J - T3I;
开发者ID:exic,项目名称:last.fm-dbus,代码行数:67,代码来源:hf_15.c


示例13: hc2cf2_16


//.........这里部分代码省略.........
				   E T3m, T3a, T3J, T3H;
				   Ip[0] = T3u + T3C;
				   Im[WS(rs, 7)] = T3u - T3C;
				   Rp[WS(rs, 4)] = T3q + T3t;
				   Rm[WS(rs, 3)] = T3q - T3t;
				   T3m = T36 - T39;
				   T3a = T36 + T39;
				   T3J = T3G - T3F;
				   T3H = T3F + T3G;
				   {
					E T2Q, T20, T3N, T3T, T2J, T2C, T3O, T2f, T34, T30, T2W, T2V, T3U, T2T, T2N;
					E T2v;
					{
					     E T2R, T27, T2e, T2S;
					     {
						  E T3n, T3f, T3o, T3k;
						  T2Q = T1U + T1Z;
						  T20 = T1U - T1Z;
						  T3n = T3e - T3b;
						  T3f = T3b + T3e;
						  T3o = T3g + T3j;
						  T3k = T3g - T3j;
						  T3N = T3L - T3M;
						  T3T = T3M + T3L;
						  {
						       E T3p, T3I, T3K, T3l;
						       T3p = T3n - T3o;
						       T3I = T3n + T3o;
						       T3K = T3k - T3f;
						       T3l = T3f + T3k;
						       Rp[WS(rs, 6)] = FMA(KP707106781, T3p, T3m);
						       Rm[WS(rs, 1)] = FNMS(KP707106781, T3p, T3m);
						       Ip[WS(rs, 2)] = FMA(KP707106781, T3I, T3H);
						       Im[WS(rs, 5)] = FMS(KP707106781, T3I, T3H);
						       Ip[WS(rs, 6)] = FMA(KP707106781, T3K, T3J);
						       Im[WS(rs, 1)] = FMS(KP707106781, T3K, T3J);
						       Rp[WS(rs, 2)] = FMA(KP707106781, T3l, T3a);
						       Rm[WS(rs, 5)] = FNMS(KP707106781, T3l, T3a);
						       T2R = T26 + T25;
						       T27 = T25 - T26;
						       T2e = T28 + T2d;
						       T2S = T28 - T2d;
						  }
					     }
					     {
						  E T2Y, T2Z, T2n, T2u;
						  T2J = T2D - T2I;
						  T2Y = T2D + T2I;
						  T2Z = T2B + T2w;
						  T2C = T2w - T2B;
						  T3O = T27 + T2e;
						  T2f = T27 - T2e;
						  T34 = FMA(KP414213562, T2Y, T2Z);
						  T30 = FNMS(KP414213562, T2Z, T2Y);
						  T2W = T2l - T2m;
						  T2n = T2l + T2m;
						  T2u = T2o - T2t;
						  T2V = T2o + T2t;
						  T3U = T2S - T2R;
						  T2T = T2R + T2S;
						  T2N = FNMS(KP414213562, T2n, T2u);
						  T2v = FMA(KP414213562, T2u, T2n);
					     }
					}
					{
					     E T33, T2X, T3X, T3Y;
开发者ID:BGCX261,项目名称:zpr-rozpoznawanie-dzwieku-git,代码行数:67,代码来源:hc2cf2_16.c


示例14: DK


//.........这里部分代码省略.........
			      TQ = iio[-WS(ios, 4)];
			      Td = FMA(Tb, Tc, T9);
			      T1q = T7 * Tc;
			      TL = Tl * TK;
			      T1b = T3 * TQ;
			      TM = iio[0];
			      T1r = FNMS(Tb, T8, T1q);
			      TO = rio[WS(ios, 3)];
			      {
				   E TN, T1c, T1a, TR, T19, TP;
				   Tx = rio[WS(ios, 1)];
				   TN = FMA(Tn, TM, TL);
				   T19 = Tl * TM;
				   T1c = FNMS(T6, TO, T1b);
				   TP = T3 * TO;
				   Ty = T2 * Tx;
				   T1a = FNMS(Tn, TK, T19);
				   TH = iio[-WS(ios, 2)];
				   TR = FMA(T6, TQ, TP);
				   Tz = iio[-WS(ios, 6)];
				   T1m = T1a + T1c;
				   T1d = T1a - T1c;
				   T14 = TC * TH;
				   TS = TN + TR;
				   T18 = TN - TR;
				   T12 = T2 * Tz;
				   TD = rio[WS(ios, 5)];
			      }
			      TA = FMA(T5, Tz, Ty);
			 }
		    }
		    {
			 E Te, T1p, T1t, T10, T1g, Tv, TJ, T1i, T1e, T1z, T1B, T17, T1h;
			 {
			      E T1x, T16, T1y, T11;
			      {
				   E TU, T13, T15, TE, TZ, TI;
				   Te = T1 + Td;
				   TU = T1 - Td;
				   T13 = FNMS(T5, Tx, T12);
				   T15 = FNMS(TG, TD, T14);
				   TE = TC * TD;
				   T1p = TW + TY;
				   TZ = TW - TY;
				   T1x = T1s - T1r;
				   T1t = T1r + T1s;
				   T1l = T13 + T15;
				   T16 = T13 - T15;
				   TI = FMA(TG, TH, TE);
				   T10 = TU + TZ;
				   T1g = TU - TZ;
				   Tv = Tk + Tu;
				   T1y = Tk - Tu;
				   T11 = TA - TI;
				   TJ = TA + TI;
			      }
			      T1i = T18 + T1d;
			      T1e = T18 - T1d;
			      T1z = T1x - T1y;
			      T1B = T1y + T1x;
			      T17 = T11 + T16;
			      T1h = T16 - T11;
			 }
			 {
			      E T1j, T1A, T1f, T1C;
			      T1j = T1h - T1i;
			      T1A = T1h + T1i;
			      T1f = T17 + T1e;
			      T1C = T1e - T17;
			      iio[-WS(ios, 1)] = FMA(KP707106781, T1A, T1z);
			      rio[WS(ios, 5)] = FMS(KP707106781, T1A, T1z);
			      rio[WS(ios, 3)] = FMA(KP707106781, T1j, T1g);
			      iio[-WS(ios, 7)] = FNMS(KP707106781, T1j, T1g);
			      iio[-WS(ios, 3)] = FMA(KP707106781, T1C, T1B);
			      rio[WS(ios, 7)] = FMS(KP707106781, T1C, T1B);
			      rio[WS(ios, 1)] = FMA(KP707106781, T1f, T10);
			      iio[-WS(ios, 5)] = FNMS(KP707106781, T1f, T10);
			      T1k = Te - Tv;
			      Tw = Te + Tv;
			 }
			 T1w = T1t - T1p;
			 T1u = T1p + T1t;
			 T1v = TS - TJ;
			 TT = TJ + TS;
		    }
	       }
	  }
	  iio[-WS(ios, 2)] = T1v + T1w;
	  rio[WS(ios, 6)] = T1v - T1w;
	  rio[0] = Tw + TT;
	  iio[-WS(ios, 4)] = Tw - TT;
	  T1n = T1l - T1m;
	  T1o = T1l + T1m;
	  iio[0] = T1o + T1u;
	  rio[WS(ios, 4)] = T1o - T1u;
	  rio[WS(ios, 2)] = T1k + T1n;
	  iio[-WS(ios, 6)] = T1k - T1n;
     }
     return W;
}
开发者ID:exic,项目名称:last.fm-dbus,代码行数:101,代码来源:hf2_8.c


示例15: r2cfII_25


//.........这里部分代码省略.........
			 T1t = FMA(KP481753674, T1b, KP1_752613360 * T1e);
			 T1u = FMA(KP684547105, T1i, KP1_457937254 * T1l);
			 T1B = T1t + T1u;
			 T1r = T1p - T1q;
			 T1C = T1A + T1B;
			 T1v = T1t - T1u;
		    }
		    {
			 E T10, T17, T1f, T1m;
			 T10 = FNMS(KP497379774, TZ, KP968583161 * TW);
			 T17 = FNMS(KP1_688655851, T16, KP535826794 * T13);
			 T18 = T10 + T17;
			 T1f = FNMS(KP963507348, T1e, KP876306680 * T1b);
			 T1m = FNMS(KP1_369094211, T1l, KP728968627 * T1i);
			 T1n = T1f + T1m;
			 T1o = T18 + T1n;
			 T1G = T10 - T17;
			 T1D = T1f - T1m;
		    }
		    {
			 E T1R, T1Q, T20, T1Z;
			 Cr[0] = TT + T1o;
			 Ci[0] = -(T1z + T1C);
			 T1R = KP559016994 * (T1P + T1M);
			 T1Q = FMA(KP250000000, T1M - T1P, TT);
			 Cr[WS(csr, 4)] = FMA(KP951056516, T1J, T1Q) + FMA(KP587785252, T1U, T1R);
			 Cr[WS(csr, 9)] = FMA(KP951056516, T1U, T1Q) + FNMA(KP587785252, T1J, T1R);
			 T20 = KP559016994 * (T1Y + T1X);
			 T1Z = FMA(KP250000000, T1X - T1Y, T1z);
			 Ci[WS(csi, 9)] = FMA(KP587785252, T1V, KP951056516 * T1W) + T1Z - T20;
			 Ci[WS(csi, 4)] = FMA(KP587785252, T1W, T1Z) + FNMS(KP951056516, T1V, T20);
			 {
			      E T1E, T1F, T1s, T1w;
			      T1E = FMS(KP250000000, T1C, T1z);
			      T1F = KP559016994 * (T1B - T1A);
			      Ci[WS(csi, 5)] = FMA(KP951056516, T1D, T1E) + FNMA(KP587785252, T1G, T1F);
			      Ci[WS(csi, 10)] = FMA(KP951056516, T1G, KP587785252 * T1D) + T1E + T1F;
			      T1s = FNMS(KP250000000, T1o, TT);
			      T1w = KP559016994 * (T18 - T1n);
			      Cr[WS(csr, 5)] = FMA(KP587785252, T1r, T1s) + FMS(KP951056516, T1v, T1w);
			      Cr[WS(csr, 10)] = T1w + FMA(KP587785252, T1v, T1s) - (KP951056516 * T1r);
			 }
		    }
	       }
	       {
		    E T21, T2z, T2L, T2K, T2M, T2F, T2P, T2C, T2Q, T2l, T2o, T2p, T2w, T2u, T28;
		    E T2f, T2g, T2s, T2h;
		    T21 = TS - TR;
		    {
			 E T2x, T2y, T2I, T2J;
			 T2x = FNMS(KP844327925, T29, KP1_071653589 * T2a);
			 T2y = FNMS(KP125581039, T2d, KP998026728 * T2c);
			 T2z = T2x + T2y;
			 T2L = T2y - T2x;
			 T2I = FNMS(KP481753674, T22, KP1_752613360 * T23);
			 T2J = FMA(KP904827052, T25, KP851558583 * T26);
			 T2K = T2I + T2J;
			 T2M = T2I - T2J;
		    }
		    {
			 E T2D, T2E, T2A, T2B;
			 T2D = FMA(KP535826794, T29, KP1_688655851 * T2a);
			 T2E = FMA(KP062790519, T2c, KP1_996053456 * T2d);
			 T2F = T2D + T2E;
			 T2P = T2E - T2D;
			 T2A = FMA(KP876306680, T22, KP963507348 * T23);
开发者ID:DougBurke,项目名称:sherpa,代码行数:67,代码来源:r2cfII_25.c


示例16: hf_10


//.........这里部分代码省略.........
					T1i = Tf * Tj;
					Th = Tf * Tg;
					Tm = W[14];
					Tp = W[15];
					T1j = FNMS(Ti, Tg, T1i);
					Tk = FMA(Ti, Tj, Th);
					T1f = Tm * Tq;
					To = Tm * Tn;
				   }
			      }
			      {
				   E T1m, T1g, Tr, T1h;
				   T1m = T1j - T1l;
				   T1I = T1l + T1j;
				   Tl = Te - Tk;
				   T13 = Te + Tk;
				   T1g = FNMS(Tp, Tn, T1f);
				   Tr = FMA(Tp, Tq, To);
				   T1J = T1g + T1e;
				   T1h = T1e - T1g;
				   Ty = Tr - Tx;
				   T14 = Tr + Tx;
				   T1n = T1h - T1m;
				   T1O = T1m + T1h;
			      }
			 }
		    }
	       }
	       T24 = T1I + T1J;
	       T1K = T1I - T1J;
	       {
		    E T2a, T15, Tz, T1Z;
		    T2a = T13 - T14;
		    T15 = T13 + T14;
		    Tz = Tl + Ty;
		    T1Z = Ty - Tl;
		    {
			 E T1L, T1N, T1E, T1M;
			 {
			      E T19, T1D, T1C, T11, T1b;
			      T19 = T15 + T18;
			      T1D = T15 - T18;
			      T11 = Tz + T10;
			      T1b = Tz - T10;
			      {
				   E T1B, T1z, T1a, T1A, T1c;
				   T1B = FNMS(KP618033988, T1n, T1y);
				   T1z = FMA(KP618033988, T1y, T1n);
				   ci[WS(rs, 4)] = T8 + T11;
				   T1a = FNMS(KP250000000, T11, T8);
				   T1A = FNMS(KP559016994, T1b, T1a);
				   T1c = FMA(KP559016994, T1b, T1a);
				   T1C = FNMS(KP250000000, T19, T12);
				   T1L = FNMS(KP618033988, T1K, T1H);
				   T1N = FMA(KP618033988, T1H, T1K);
				   cr[WS(rs, 1)] = FMA(KP951056516, T1z, T1c);
				   ci[0] = FNMS(KP951056516, T1z, T1c);
				   cr[WS(rs, 3)] = FMA(KP951056516, T1B, T1A);
				   ci[WS(rs, 2)] = FNMS(KP951056516, T1B, T1A);
			      }
			      cr[0] = T12 + T19;
			      T1E = FNMS(KP559016994, T1D, T1C);
			      T1M = FMA(KP559016994, T1D, T1C);
			 }
			 {
			      E T1X, T21, T20, T22, T1Q, T1W, T1V, T26, T28, T27;
			      T1Q = T1O + T1P;
			      T1W = T1P - T1O;
			      ci[WS(rs, 3)] = FMA(KP951056516, T1N, T1M);
			      cr[WS(rs, 4)] = FNMS(KP951056516, T1N, T1M);
			      ci[WS(rs, 1)] = FMA(KP951056516, T1L, T1E);
			      cr[WS(rs, 2)] = FNMS(KP951056516, T1L, T1E);
			      T1V = FMA(KP250000000, T1Q, T1U);
			      cr[WS(rs, 5)] = T1Q - T1U;
			      T1X = FNMS(KP559016994, T1W, T1V);
			      T21 = FMA(KP559016994, T1W, T1V);
			      T20 = FNMS(KP618033988, T1Z, T1Y);
			      T22 = FMA(KP618033988, T1Y, T1Z);
			      T26 = T24 + T25;
			      T28 = T24 - T25;
			      ci[WS(rs, 8)] = FMA(KP951056516, T22, T21);
			      cr[WS(rs, 9)] = FMS(KP951056516, T22, T21);
			      ci[WS(rs, 6)] = FMA(KP951056516, T20, T1X);
			      cr[WS(rs, 7)] = FMS(KP951056516, T20, T1X);
			      T27 = FNMS(KP250000000, T26, T23);
			      ci[WS(rs, 9)] = T26 + T23;
			      T29 = FMA(KP559016994, T28, T27);
			      T2d = FNMS(KP559016994, T28, T27);
			      T2c = FMA(KP618033988, T2b, T2a);
			      T2e = FNMS(KP618033988, T2a, T2b);
			 }
		    }
	       }
	  }
	  ci[WS(rs, 7)] = FMA(KP951056516, T2e, T2d);
	  cr[WS(rs, 8)] = FMS(KP951056516, T2e, T2d);
	  ci[WS(rs, 5)] = FMA(KP951056516, T2c, T29);
	  cr[WS(rs, 6)] = FMS(KP951056516, T2c, T29);
     }
}
开发者ID:BackupTheBerlios,项目名称:openvsipl,代码行数:101,代码来源:hf_10.c


示例17: r2cf_25


//.........这里部分代码省略.........
				   T17 = Tx - Ty;
				   Tz = Tx + Ty;
				   T3 = R0[WS(rs, 12)];
				   T5 = R0[WS(rs, 7)];
				   T1L = FMA(KP618033988, T16, T17);
				   T18 = FNMS(KP618033988, T17, T16);
				   T1a = Tz - Tw;
				   TA = Tw + Tz;
				   T4 = T2 + T3;
				   TU = T3 - T2;
				   T6 = R1[WS(rs, 9)];
			      }
			      TB = Tt + TA;
			      T19 = FNMS(KP250000000, TA, Tt);
			      {
				   E T7, TV, T1b, T1K, T8;
				   T7 = T5 + T6;
				   TV = T5 - T6;
				   T1b = FNMS(KP559016994, T1a, T19);
				   T1K = FMA(KP559016994, T1a, T19);
				   T1P = FMA(KP618033988, TU, TV);
				   TW = FNMS(KP618033988, TV, TU);
				   TS = T4 - T7;
				   T8 = T4 + T7;
				   T21 = FMA(KP869845200, T1K, T1L);
				   T1M = FNMS(KP786782374, T1L, T1K);
				   T2e = FMA(KP066152395, T1K, T1L);
				   T26 = FNMS(KP059835404, T1L, T1K);
				   T1B = FMA(KP120146378, T18, T1b);
				   T1r = FNMS(KP132830569, T1b, T18);
				   T1k = FMA(KP893101515, T18, T1b);
				   T1c = FNMS(KP987388751, T1b, T18);
				   T9 = T1 + T8;
				   TR = FMS(KP250000000, T8, T1);
			      }
			 }
			 {
			      E Ta, Te, TK, Td, Tf;
			      Ta = R1[0];
			      {
				   E Tb, Tc, T1Q, TT;
				   Tb = R0[WS(rs, 3)];
				   Tc = R1[WS(rs, 10)];
				   T1Q = FMA(KP559016994, TS, TR);
				   TT = FNMS(KP559016994, TS, TR);
				   Te = R1[WS(rs, 5)];
				   TK = Tb - Tc;
				   Td = Tb + Tc;
				   T1X = FNMS(KP120146378, T1P, T1Q);
				   T1R = FMA(KP132830569, T1Q, T1P);
				   T2k = FMA(KP494780565, T1Q, T1P);
				   T29 = FNMS(KP447533225, T1P, T1Q);
				   T1z = FMA(KP869845200, TT, TW);
				   T1v = FNMS(KP786782374, TW, TT);
				   T1h = FNMS(KP667278218, TT, TW);
				   TX = FMA(KP603558818, TW, TT);
				   Tf = R0[WS(rs, 8)];
			      }
			      {
				   E Tk, T1S, TM, TO, Tn, TZ, TN, T10, Tq, To, Th, Tp, TP, T1T, Tr;
				   Tk = R0[WS(rs, 1)];
				   {
					E Tl, TL, Tg, Tm;
					Tl = R1[WS(rs, 3)];
					TL = Tf - Te;
					Tg = Te + Tf;
开发者ID:phillipstanleymarbell,项目名称:sunflower-simulator,代码行数:67,代码来源:r2cf_25.c


示例18: hf_9


//.........这里部分代码省略.........
		    Tj = W[6];
		    Tl = W[7];
		    Tn = FMA(Tj, Tk, Tl * Tm);
		    TT = FNMS(Tl, Tk, Tj * Tm);
	       }
	       {
		    E Tp, Tr, To, Tq;
		    Tp = cr[WS(rs, 7)];
		    Tr = ci[WS(rs, 7)];
		    To = W[12];
		    Tq = W[13];
		    Ts = FMA(To, Tp, Tq * Tr);
		    TU = FNMS(Tq, Tp, To * Tr);
	       }
	       Tt = Tn + Ts;
	       TY = TT + TU;
	       Tu = Ti + Tt;
	       T1w = TX + TY;
	       {
		    E TS, TV, TZ, T10;
		    TS = FNMS(KP500000000, Tt, Ti);
		    TV = KP866025403 * (TT - TU);
		    TW = TS - TV;
		    T1k = TS + TV;
		    TZ = FNMS(KP500000000, TY, TX);
		    T10 = KP866025403 * (Ts - Tn);
		    T11 = TZ - T10;
		    T1l = T10 + TZ;
	       }
	  }
	  {
	       E T1y, Td, TM, T1v;
	       T1y = KP866025403 * (T1w - T1x);
	       Td = T1 + Tc;
	       TM = Tu + TL;
	       T1v = FNMS(KP500000000, TM, Td);
	       cr[0] = Td + TM;
	       cr[WS(rs, 3)] = T1v + T1y;
	       ci[WS(rs, 2)] = T1v - T1y;
	  }
	  {
	       E TR, T1I, T1e, T1K, T1i, T1H, T1f, T1J;
	       TR = TN - TQ;
	       T1I = T1D - T1A;
	       {
		    E T12, T1d, T1g, T1h;
		    T12 = FMA(KP173648177, TW, KP984807753 * T11);
		    T1d = FNMS(KP939692620, T1c, KP342020143 * T19);
		    T1e = T12 + T1d;
		    T1K = KP866025403 * (T1d - T12);
		    T1g = FNMS(KP984807753, TW, KP173648177 * T11);
		    T1h = FMA(KP342020143, T1c, KP939692620 * T19);
		    T1i = KP866025403 * (T1g + T1h);
		    T1H = T1g - T1h;
	       }
	       cr[WS(rs, 2)] = TR + T1e;
	       ci[WS(rs, 6)] = T1H + T1I;
	       T1f = FNMS(KP500000000, T1e, TR);
	       ci[0] = T1f - T1i;
	       ci[WS(rs, 3)] = T1f + T1i;
	       T1J = FMS(KP500000000, T1H, T1I);
	       cr[WS(rs, 5)] = T1J - T1K;
	       cr[WS(rs, 8)] = T1K + T1J;
	  }
	  {
	       E T1L, T1M, T1N, T1O;
	       T1L = KP866025403 * (TL - Tu);
	       T1M = T1C + T1B;
	       T1N = T1w + T1x;
	       T1O = FNMS(KP500000000, T1N, T1M);
	       cr[WS(rs, 6)] = T1L - T1O;
	       ci[WS(rs, 8)] = T1N + T1M;
	       ci[WS(rs, 5)] = T1L + T1O;
	  }
	  {
	       E T1j, T1E, T1q, T1z, T1u, T1F, T1r, T1G;
	       T1j = TN + TQ;
	       T1E = T1A + T1D;
	       {
		    E T1m, T1p, T1s, T1t;
		    T1m = FMA(KP766044443, T1k, KP642787609 * T1l);
		    T1p = FMA(KP173648177, T1n, KP984807753 * T1o);
		    T1q = T1m + T1p;
		    T1z = KP866025403 * (T1p - T1m);
		    T1s = FNMS(KP642787609, T1k, KP766044443 * T1l);
		    T1t = FNMS(KP984807753, T1n, KP173648177 * T1o);
		    T1u = KP866025403 * (T1s - T1t);
		    T1F = T1s + T1t;
	       }
	       cr[WS(rs, 1)] = T1j + T1q;
	       T1r = FNMS(KP500000000, T1q, T1j);
	       ci[WS(rs, 1)] = T1r - T1u;
	       cr[WS(rs, 4)] = T1r + T1u;
	       ci[WS(rs, 7)] = T1F + T1E;
	       T1G = FNMS(KP500000000, T1F, T1E);
	       cr[WS(rs, 7)] = T1z - T1G;
	       ci[WS(rs, 4)] = T1z + T1G;
	  }
     }
}
开发者ID:bambang,项目名称:vsipl,代码行数:101,代码来源:hf_9.c


示例19: r2cbIII_25

static void r2cbIII_25(R *R0, R *R1, R *Cr, R *Ci, stride rs, stride csr, stride csi, INT v, INT ivs, INT ovs)
{
     DK(KP968583161, +0.968583161128631119490168375464735813836012403);
     DK(KP248689887, +0.248689887164854788242283746006447968417567406);
     DK(KP684547105, +0.684547105928688673732283357621209269889519233);
     DK(KP728968627, +0.728968627421411523146730319055259111372571664);
     DK(KP062790519, +0.062790519529313376076178224565631133122484832);
     DK(KP998026728, +0.998026728428271561952336806863450553336905220);
     DK(KP876306680, +0.876306680043863587308115903922062583399064238);
     DK(KP481753674, +0.481753674101715274987191502872129653528542010);
     DK(KP535826794, +0.535826794978996618271308767867639978063575346);
     DK(KP844327925, +0.844327925502015078548558063966681505381659241);
     DK(KP904827052, +0.904827052466019527713668647932697593970413911);
     DK(KP425779291, +0.425779291565072648862502445744251703979973042);
     DK(KP250000000, +0.250000000000000000000000000000000000000000000);
     DK(KP951056516, +0.951056516295153572116439333379382143405698634);
     DK(KP587785252, +0.587785252292473129168705954639072768597652438);
     DK(KP559016994, +0.559016994374947424102293417182819058860154590);
     DK(KP500000000, +0.500000000000000000000000000000000000000000000);
     DK(KP2_000000000, +2.000000000000000000000000000000000000000000000);
     DK(KP1_118033988, +1.118033988749894848204586834365638117720309180);
     DK(KP1_175570504, +1.175570504584946258337411909278145537195304875);
     DK(KP1_902113032, +1.902113032590307144232878666758764286811397268);
     {
	  INT i;
	  for (i = v; i > 0; i = i - 1, R0 = R0 + ovs, R1 = R1 + ovs, Cr = Cr + ivs, Ci = Ci + ivs, MAKE_VOLATILE_STRIDE(100, rs), MAKE_VOLATILE_STRIDE(100, csr), MAKE_VOLATILE_STRIDE(100, csi)) {
	       E TS, T1O, T5, TP, T1N, TI, TH, Te, T17, T2h, T1y, T1V, T10, T2g, T1x;
	       E T1S, Tz, Ty, Tn, T1m, T2e, T1B, T22, T1f, T2d, T1A, T1Z, TQ, TR;
	       TQ = Ci[WS(csi, 2)];
	       TR = Ci[WS(csi, 7)];
	       TS = FNMS(KP1_175570504, TR, KP1_902113032 * TQ);
	       T1O = FMA(KP1_902113032, TR, KP1_175570504 * TQ);
	       {
		    E T1, T4, TN, T2, T3, TO;
		    T1 = Cr[WS(csr, 12)];
		    T2 = Cr[WS(csr, 7)];
		    T3 = Cr[WS(csr, 2)];
		    T4 = T2 + T3;
		    TN = KP1_118033988 * (T3 - T2);
		    T5 = FMA(KP2_000000000, T4, T1);
		    TO = FMS(KP500000000, T4, T1);
		    TP = TN - TO;
		    T1N = TO + TN;
	       }
	       {
		    E T6, Td, T15, TU, T14, T11, TX, TY;
		    T6 = Cr[WS(csr, 11)];
		    TI = Ci[WS(csi, 11)];
		    {
			 E T7, T8, T9, Ta, Tb, Tc;
			 T7 = Cr[WS(csr, 6)];
			 T8 = Cr[WS(csr, 8)];
			 T9 = T7 + T8;
			 Ta = Cr[WS(csr, 1)];
			 Tb = Cr[WS(csr, 3)];
			 Tc = Ta + Tb;
			 Td = T9 + Tc;
			 T15 = Ta - Tb;
			 TU = KP559016994 * (Tc - T9);
			 T14 = T8 - T7;
		    }
		    {
			 E TB, TC, TD, TE, TF, TG;
			 TB = Ci[WS(csi, 6)];
			 TC = Ci[WS(csi, 8)];
			 TD = TB - TC;
			 TE = Ci[WS(csi, 1)];
			 TF = Ci[WS(csi, 3)];
			 TG = TE - TF;
			 TH = TD + TG;
			 T11 = KP559016994 * (TD - TG);
			 TX = TB + TC;
			 TY = TE + TF;
		    }
		    Te = T6 + Td;
		    {
			 E T16, T1T, T13, T1U, T12;
			 T16 = FMA(KP587785252, T14, KP951056516 * T15);
			 T1T = FNMS(KP587785252, T15, KP951056516 * T14);
			 T12 = FNMS(KP250000000, TH, TI);
			 T13 = T11 - T12;
			 T1U = T11 + T12;
			 T17 = T13 - T16;
			 T2h = T1T - T1U;
			 T1y = T16 + T13;
			 T1V = T1T + T1U;
		    }
		    {
			 E TZ, T1R, TW, T1Q, TV;
			 TZ = FNMS(KP951056516, TY, KP587785252 * TX);
			 T1R = FMA(KP951056516, TX, KP587785252 * TY);
			 TV = FMS(KP250000000, Td, T6);
			 TW = TU - TV;
			 T1Q = TV + TU;
			 T10 = TW + TZ;
			 T2g = T1Q + T1R;
			 T1x = TZ - TW;
			 T1S = T1Q - T1R;
		    }
	       }
//.........这里部分代码省略.........
开发者ID:DougBurke,项目名称:sherpa,代码行数:101,代码来源:r2cbIII_25.c


示例20: hc2cf_8


//.........这里部分代码省略.........
			      TQ = FNMS(Tc, Ta, TP);
			      Te = FMA(Tc, Td, Tb);
			      T13 = TA * TE;
			      TC = TA * TB;
			      TK = Im[WS(rs, 1)];
			      TG = W[4];
			      TD = W[13];
			      TJ = W[5];
			      {
				   E T14, TF, T16, TL, T15, TI;
				   To = Ip[0];
				   T15 = TG * TK;
				   TI = TG * TH;
				   T14 = FNMS(TD, TB, T13);
				   TF = FMA(TD, TE, TC);
				   T16 = FNMS(TJ, TH, T15);
				   TL = FMA(TJ, TK, TI);
				   Tr = Im[0];
				   Tn = W[0];
				   T17 = T14 - T16;
				   T1g = T14 + T16;
				   TM = TF + TL;
				   T12 = TF - TL;
			      }
			      Tu = Ip[WS(rs, 2)];
			      TW = Tn * Tr;
			      Tp = Tn * To;
			      Tx = Im[WS(rs, 2)];
			      Tt = W[8];
			      Tq = W[1];
			      Tw = W[9];
			 }
		    }
	       }
	       {
		    E T8, T1j, T1n, Tz, T1a, TU, Tl, T1b, T1c, T1v, T1t, T1w, T19, T1u, T1d;
		    {
			 E T1r, T10, TV, T1s, T11, T18;
			 {
			      E TO, TX, Ts, TZ, Ty, TT, TY, Tv;
			      T8 = T1 + T7;
			      TO = T1 - T7;
			      TY = Tt * Tx;
			      Tv = Tt * Tu;
			      TX = FNMS(Tq, To, TW);
			      Ts = FMA(Tq, Tr, Tp);
			      TZ = FNMS(Tw, Tu, TY);
			      Ty = FMA(Tw, Tx, Tv);
			      TT = TQ - TS;
			      T1j = TQ + TS;
			      T1n = T1l + T1m;
			      T1r = T1m - T1l;
			      T10 = TX - TZ;
			      T1f = TX + TZ;
			      Tz = Ts + Ty;
			      TV = Ts - Ty;
			      T1a = TO - TT;
			      TU = TO + TT;
			      T1s = Te - Tk;
			      Tl = Te + Tk;
			 }
			 T1b = T10 - TV;
			 T11 = TV + T10;
			 T18 = T12 - T17;
			 T1c = T12 + T17;
			 T1v = T1s + T1r;
			 T1t = T1r - T1s;
			 T1w = T18 - T11;
			 T19 = T11 + T18;
		    }
		    Ip[WS(rs, 3)] = FMA(KP707106781, T1w, T1v);
		    Im[0] = FMS(KP707106781, T1w, T1v);
		    Rp[WS(rs, 1)] = FMA(KP707106781, T19, TU);
		    Rm[WS(rs, 2)] = FNMS(KP707106781, T19, TU);
		    T1u = T1b + T1c;
		    T1d = T1b - T1c;
		    Ip[WS(rs, 1)] = FMA(KP707106781, T1u, T1t);
		    Im[WS(rs, 2)] = FMS(KP707106781, T1u, T1t);
		    Rp[WS(rs, 3)] = FMA(KP707106781, T1d, T1a);
		    Rm[0 

鲜花

握手

雷人

路过

鸡蛋
该文章已有0人参与评论

请发表评论

全部评论

专题导读
上一篇:
C++ FMargin函数代码示例发布时间:2022-05-30
下一篇:
C++ FMR_ASSERT函数代码示例发布时间:2022-05-30
热门推荐
阅读排行榜

扫描微信二维码

查看手机版网站

随时了解更新最新资讯

139-2527-9053

在线客服(服务时间 9:00~18:00)

在线QQ客服
地址:深圳市南山区西丽大学城创智工业园
电邮:jeky_zhao#qq.com
移动电话:139-2527-9053

Powered by 互联科技 X3.4© 2001-2213 极客世界.|Sitemap