C++ BYTW函数代码示例

OStack程序员社区-中国程序员成长平台 › 门户 › 编程› C++›C++教程

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

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

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

示例1: t1bv_4

static void t1bv_4(R *ri, R *ii, const R *W, stride rs, INT mb, INT me, INT ms)
{
     {
	  INT m;
	  R *x;
	  x = ii;
	  for (m = mb, W = W + (mb * ((TWVL / VL) * 6)); m < me; m = m + VL, x = x + (VL * ms), W = W + (TWVL * 6), MAKE_VOLATILE_STRIDE(4, rs)) {
	       V T1, T8, T3, T6, T7, T2, T5;
	       T1 = LD(&(x[0]), ms, &(x[0]));
	       T7 = LD(&(x[WS(rs, 3)]), ms, &(x[WS(rs, 1)]));
	       T8 = BYTW(&(W[TWVL * 4]), T7);
	       T2 = LD(&(x[WS(rs, 2)]), ms, &(x[0]));
	       T3 = BYTW(&(W[TWVL * 2]), T2);
	       T5 = LD(&(x[WS(rs, 1)]), ms, &(x[WS(rs, 1)]));
	       T6 = BYTW(&(W[0]), T5);
	       {
		    V T4, T9, Ta, Tb;
		    T4 = VSUB(T1, T3);
		    T9 = VBYI(VSUB(T6, T8));
		    ST(&(x[WS(rs, 3)]), VSUB(T4, T9), ms, &(x[WS(rs, 1)]));
		    ST(&(x[WS(rs, 1)]), VADD(T4, T9), ms, &(x[WS(rs, 1)]));
		    Ta = VADD(T1, T3);
		    Tb = VADD(T6, T8);
		    ST(&(x[WS(rs, 2)]), VSUB(Ta, Tb), ms, &(x[0]));
		    ST(&(x[0]), VADD(Ta, Tb), ms, &(x[0]));
	       }
	  }
     }
     VLEAVE();
}

开发者ID:barnex，项目名称:fftw，代码行数:30，代码来源:dft_simd_sse2_t1bv_4.c

示例2: t1buv_3

static void t1buv_3(R *ri, R *ii, const R *W, stride rs, INT mb, INT me, INT ms)
{
    DVK(KP500000000, +0.500000000000000000000000000000000000000000000);
    DVK(KP866025403, +0.866025403784438646763723170752936183471402627);
    INT m;
    R *x;
    x = ii;
    for (m = mb, W = W + (mb * ((TWVL / VL) * 4)); m < me; m = m + VL, x = x + (VL * ms), W = W + (TWVL * 4), MAKE_VOLATILE_STRIDE(rs)) {
        V T1, T2, T4;
        T1 = LD(&(x[0]), ms, &(x[0]));
        T2 = LD(&(x[WS(rs, 1)]), ms, &(x[WS(rs, 1)]));
        T4 = LD(&(x[WS(rs, 2)]), ms, &(x[0]));
        {
            V T3, T5, T8, T6, T7;
            T3 = BYTW(&(W[0]), T2);
            T5 = BYTW(&(W[TWVL * 2]), T4);
            T8 = VMUL(LDK(KP866025403), VSUB(T3, T5));
            T6 = VADD(T3, T5);
            T7 = VFNMS(LDK(KP500000000), T6, T1);
            ST(&(x[0]), VADD(T1, T6), ms, &(x[0]));
            ST(&(x[WS(rs, 2)]), VFNMSI(T8, T7), ms, &(x[0]));
            ST(&(x[WS(rs, 1)]), VFMAI(T8, T7), ms, &(x[WS(rs, 1)]));
        }
    }
}

开发者ID:simonmr，项目名称:SuperDARN_MSI_ROS，代码行数:25，代码来源:t1buv_3.c

示例3: t1bv_6

static void t1bv_6(R *ri, R *ii, const R *W, stride rs, INT mb, INT me, INT ms)
{
     DVK(KP500000000, +0.500000000000000000000000000000000000000000000);
     DVK(KP866025403, +0.866025403784438646763723170752936183471402627);
     {
	  INT m;
	  R *x;
	  x = ii;
	  for (m = mb, W = W + (mb * ((TWVL / VL) * 10)); m < me; m = m + VL, x = x + (VL * ms), W = W + (TWVL * 10), MAKE_VOLATILE_STRIDE(6, rs)) {
	       V T1, T2, Ta, Tc, T5, T7;
	       T1 = LD(&(x[0]), ms, &(x[0]));
	       T2 = LD(&(x[WS(rs, 3)]), ms, &(x[WS(rs, 1)]));
	       Ta = LD(&(x[WS(rs, 4)]), ms, &(x[0]));
	       Tc = LD(&(x[WS(rs, 1)]), ms, &(x[WS(rs, 1)]));
	       T5 = LD(&(x[WS(rs, 2)]), ms, &(x[0]));
	       T7 = LD(&(x[WS(rs, 5)]), ms, &(x[WS(rs, 1)]));
	       {
		    V T3, Tb, Td, T6, T8;
		    T3 = BYTW(&(W[TWVL * 4]), T2);
		    Tb = BYTW(&(W[TWVL * 6]), Ta);
		    Td = BYTW(&(W[0]), Tc);
		    T6 = BYTW(&(W[TWVL * 2]), T5);
		    T8 = BYTW(&(W[TWVL * 8]), T7);
		    {
			 V Ti, T4, Tk, Te, Tj, T9;
			 Ti = VADD(T1, T3);
			 T4 = VSUB(T1, T3);
			 Tk = VADD(Tb, Td);
			 Te = VSUB(Tb, Td);
			 Tj = VADD(T6, T8);
			 T9 = VSUB(T6, T8);
			 {
			      V Tl, Tn, Tf, Th, Tm, Tg;
			      Tl = VADD(Tj, Tk);
			      Tn = VMUL(LDK(KP866025403), VSUB(Tj, Tk));
			      Tf = VADD(T9, Te);
			      Th = VMUL(LDK(KP866025403), VSUB(T9, Te));
			      ST(&(x[0]), VADD(Ti, Tl), ms, &(x[0]));
			      Tm = VFNMS(LDK(KP500000000), Tl, Ti);
			      ST(&(x[WS(rs, 3)]), VADD(T4, Tf), ms, &(x[WS(rs, 1)]));
			      Tg = VFNMS(LDK(KP500000000), Tf, T4);
			      ST(&(x[WS(rs, 4)]), VFMAI(Tn, Tm), ms, &(x[0]));
			      ST(&(x[WS(rs, 2)]), VFNMSI(Tn, Tm), ms, &(x[0]));
			      ST(&(x[WS(rs, 5)]), VFNMSI(Th, Tg), ms, &(x[WS(rs, 1)]));
			      ST(&(x[WS(rs, 1)]), VFMAI(Th, Tg), ms, &(x[WS(rs, 1)]));
			 }
		    }
	       }
	  }
     }
     VLEAVE();
}

开发者ID:Aegisub，项目名称:fftw3，代码行数:52，代码来源:t1bv_6.c

示例4: DVK

static const R *t1bv_6(R *ri, R *ii, const R *W, stride ios, int m, int dist)
{
     DVK(KP500000000, +0.500000000000000000000000000000000000000000000);
     DVK(KP866025403, +0.866025403784438646763723170752936183471402627);
     int i;
     R *x;
     x = ii;
     BEGIN_SIMD();
     for (i = m; i > 0; i = i - VL, x = x + (VL * dist), W = W + (TWVL * 10)) {
	  V Tf, Ti, Ta, Tk, T5, Tj, Tc, Te, Td;
	  Tc = LD(&(x[0]), dist, &(x[0]));
	  Td = LD(&(x[WS(ios, 3)]), dist, &(x[WS(ios, 1)]));
	  Te = BYTW(&(W[TWVL * 4]), Td);
	  Tf = VSUB(Tc, Te);
	  Ti = VADD(Tc, Te);
	  {
	       V T7, T9, T6, T8;
	       T6 = LD(&(x[WS(ios, 4)]), dist, &(x[0]));
	       T7 = BYTW(&(W[TWVL * 6]), T6);
	       T8 = LD(&(x[WS(ios, 1)]), dist, &(x[WS(ios, 1)]));
	       T9 = BYTW(&(W[0]), T8);
	       Ta = VSUB(T7, T9);
	       Tk = VADD(T7, T9);
	  }
	  {
	       V T2, T4, T1, T3;
	       T1 = LD(&(x[WS(ios, 2)]), dist, &(x[0]));
	       T2 = BYTW(&(W[TWVL * 2]), T1);
	       T3 = LD(&(x[WS(ios, 5)]), dist, &(x[WS(ios, 1)]));
	       T4 = BYTW(&(W[TWVL * 8]), T3);
	       T5 = VSUB(T2, T4);
	       Tj = VADD(T2, T4);
	  }
	  {
	       V Tb, Tg, Th, Tn, Tl, Tm;
	       Tb = VBYI(VMUL(LDK(KP866025403), VSUB(T5, Ta)));
	       Tg = VADD(T5, Ta);
	       Th = VFNMS(LDK(KP500000000), Tg, Tf);
	       ST(&(x[WS(ios, 1)]), VADD(Tb, Th), dist, &(x[WS(ios, 1)]));
	       ST(&(x[WS(ios, 3)]), VADD(Tf, Tg), dist, &(x[WS(ios, 1)]));
	       ST(&(x[WS(ios, 5)]), VSUB(Th, Tb), dist, &(x[WS(ios, 1)]));
	       Tn = VBYI(VMUL(LDK(KP866025403), VSUB(Tj, Tk)));
	       Tl = VADD(Tj, Tk);
	       Tm = VFNMS(LDK(KP500000000), Tl, Ti);
	       ST(&(x[WS(ios, 2)]), VSUB(Tm, Tn), dist, &(x[0]));
	       ST(&(x[0]), VADD(Ti, Tl), dist, &(x[0]));
	       ST(&(x[WS(ios, 4)]), VADD(Tn, Tm), dist, &(x[0]));
	  }
     }
     END_SIMD();
     return W;
}

开发者ID:abrahamneben，项目名称:orbcomm_beam_mapping，代码行数:52，代码来源:t1bv_6.c

示例5: t1buv_5

static void t1buv_5(R *ri, R *ii, const R *W, stride rs, INT mb, INT me, INT ms)
{
     DVK(KP559016994, +0.559016994374947424102293417182819058860154590);
     DVK(KP250000000, +0.250000000000000000000000000000000000000000000);
     DVK(KP618033988, +0.618033988749894848204586834365638117720309180);
     DVK(KP951056516, +0.951056516295153572116439333379382143405698634);
     {
	  INT m;
	  R *x;
	  x = ii;
	  for (m = mb, W = W + (mb * ((TWVL / VL) * 8)); m < me; m = m + VL, x = x + (VL * ms), W = W + (TWVL * 8), MAKE_VOLATILE_STRIDE(5, rs)) {
	       V T1, T2, T9, T4, T7;
	       T1 = LD(&(x[0]), ms, &(x[0]));
	       T2 = LD(&(x[WS(rs, 1)]), ms, &(x[WS(rs, 1)]));
	       T9 = LD(&(x[WS(rs, 3)]), ms, &(x[WS(rs, 1)]));
	       T4 = LD(&(x[WS(rs, 4)]), ms, &(x[0]));
	       T7 = LD(&(x[WS(rs, 2)]), ms, &(x[0]));
	       {
		    V T3, Ta, T5, T8;
		    T3 = BYTW(&(W[0]), T2);
		    Ta = BYTW(&(W[TWVL * 4]), T9);
		    T5 = BYTW(&(W[TWVL * 6]), T4);
		    T8 = BYTW(&(W[TWVL * 2]), T7);
		    {
			 V T6, Tg, Tb, Th;
			 T6 = VADD(T3, T5);
			 Tg = VSUB(T3, T5);
			 Tb = VADD(T8, Ta);
			 Th = VSUB(T8, Ta);
			 {
			      V Te, Tc, Tk, Ti, Td, Tj, Tf;
			      Te = VSUB(T6, Tb);
			      Tc = VADD(T6, Tb);
			      Tk = VMUL(LDK(KP951056516), VFNMS(LDK(KP618033988), Tg, Th));
			      Ti = VMUL(LDK(KP951056516), VFMA(LDK(KP618033988), Th, Tg));
			      Td = VFNMS(LDK(KP250000000), Tc, T1);
			      ST(&(x[0]), VADD(T1, Tc), ms, &(x[0]));
			      Tj = VFNMS(LDK(KP559016994), Te, Td);
			      Tf = VFMA(LDK(KP559016994), Te, Td);
			      ST(&(x[WS(rs, 2)]), VFNMSI(Tk, Tj), ms, &(x[0]));
			      ST(&(x[WS(rs, 3)]), VFMAI(Tk, Tj), ms, &(x[WS(rs, 1)]));
			      ST(&(x[WS(rs, 4)]), VFNMSI(Ti, Tf), ms, &(x[0]));
			      ST(&(x[WS(rs, 1)]), VFMAI(Ti, Tf), ms, &(x[WS(rs, 1)]));
			 }
		    }
	       }
	  }
     }
     VLEAVE();
}

开发者ID:Enny1991，项目名称:ExApp，代码行数:50，代码来源:t1buv_5.c

示例6: t1buv_6

static void t1buv_6(R *ri, R *ii, const R *W, stride rs, INT mb, INT me, INT ms)
{
     DVK(KP500000000, +0.500000000000000000000000000000000000000000000);
     DVK(KP866025403, +0.866025403784438646763723170752936183471402627);
     INT m;
     R *x;
     x = ii;
     for (m = mb, W = W + (mb * ((TWVL / VL) * 10)); m < me; m = m + VL, x = x + (VL * ms), W = W + (TWVL * 10), MAKE_VOLATILE_STRIDE(rs)) {
	  V Tf, Ti, Ta, Tk, T5, Tj, Tc, Te, Td;
	  Tc = LD(&(x[0]), ms, &(x[0]));
	  Td = LD(&(x[WS(rs, 3)]), ms, &(x[WS(rs, 1)]));
	  Te = BYTW(&(W[TWVL * 4]), Td);
	  Tf = VSUB(Tc, Te);
	  Ti = VADD(Tc, Te);
	  {
	       V T7, T9, T6, T8;
	       T6 = LD(&(x[WS(rs, 4)]), ms, &(x[0]));
	       T7 = BYTW(&(W[TWVL * 6]), T6);
	       T8 = LD(&(x[WS(rs, 1)]), ms, &(x[WS(rs, 1)]));
	       T9 = BYTW(&(W[0]), T8);
	       Ta = VSUB(T7, T9);
	       Tk = VADD(T7, T9);
	  }
	  {
	       V T2, T4, T1, T3;
	       T1 = LD(&(x[WS(rs, 2)]), ms, &(x[0]));
	       T2 = BYTW(&(W[TWVL * 2]), T1);
	       T3 = LD(&(x[WS(rs, 5)]), ms, &(x[WS(rs, 1)]));
	       T4 = BYTW(&(W[TWVL * 8]), T3);
	       T5 = VSUB(T2, T4);
	       Tj = VADD(T2, T4);
	  }
	  {
	       V Tb, Tg, Th, Tn, Tl, Tm;
	       Tb = VBYI(VMUL(LDK(KP866025403), VSUB(T5, Ta)));
	       Tg = VADD(T5, Ta);
	       Th = VFNMS(LDK(KP500000000), Tg, Tf);
	       ST(&(x[WS(rs, 1)]), VADD(Tb, Th), ms, &(x[WS(rs, 1)]));
	       ST(&(x[WS(rs, 3)]), VADD(Tf, Tg), ms, &(x[WS(rs, 1)]));
	       ST(&(x[WS(rs, 5)]), VSUB(Th, Tb), ms, &(x[WS(rs, 1)]));
	       Tn = VBYI(VMUL(LDK(KP866025403), VSUB(Tj, Tk)));
	       Tl = VADD(Tj, Tk);
	       Tm = VFNMS(LDK(KP500000000), Tl, Ti);
	       ST(&(x[WS(rs, 2)]), VSUB(Tm, Tn), ms, &(x[0]));
	       ST(&(x[0]), VADD(Ti, Tl), ms, &(x[0]));
	       ST(&(x[WS(rs, 4)]), VADD(Tn, Tm), ms, &(x[0]));
	  }
     }
}

开发者ID:8cH9azbsFifZ，项目名称:wspr，代码行数:49，代码来源:t1buv_6.c

示例7: t1bv_7

static void t1bv_7(R *ri, R *ii, const R *W, stride rs, INT mb, INT me, INT ms)
{
     DVK(KP222520933, +0.222520933956314404288902564496794759466355569);
     DVK(KP900968867, +0.900968867902419126236102319507445051165919162);
     DVK(KP623489801, +0.623489801858733530525004884004239810632274731);
     DVK(KP433883739, +0.433883739117558120475768332848358754609990728);
     DVK(KP781831482, +0.781831482468029808708444526674057750232334519);
     DVK(KP974927912, +0.974927912181823607018131682993931217232785801);
     INT m;
     R *x;
     x = ii;
     for (m = mb, W = W + (mb * ((TWVL / VL) * 12)); m < me; m = m + VL, x = x + (VL * ms), W = W + (TWVL * 12), MAKE_VOLATILE_STRIDE(rs)) {
	  V Th, Tf, Ti, T5, Tk, Ta, Tj, To, Tp;
	  Th = LD(&(x[0]), ms, &(x[0]));
	  {
	       V Tc, Te, Tb, Td;
	       Tb = LD(&(x[WS(rs, 2)]), ms, &(x[0]));
	       Tc = BYTW(&(W[TWVL * 2]), Tb);
	       Td = LD(&(x[WS(rs, 5)]), ms, &(x[WS(rs, 1)]));
	       Te = BYTW(&(W[TWVL * 8]), Td);
	       Tf = VSUB(Tc, Te);
	       Ti = VADD(Tc, Te);
	  }
	  {
	       V T2, T4, T1, T3;
	       T1 = LD(&(x[WS(rs, 1)]), ms, &(x[WS(rs, 1)]));
	       T2 = BYTW(&(W[0]), T1);
	       T3 = LD(&(x[WS(rs, 6)]), ms, &(x[0]));
	       T4 = BYTW(&(W[TWVL * 10]), T3);
	       T5 = VSUB(T2, T4);
	       Tk = VADD(T2, T4);
	  }
	  {
	       V T7, T9, T6, T8;
	       T6 = LD(&(x[WS(rs, 3)]), ms, &(x[WS(rs, 1)]));
	       T7 = BYTW(&(W[TWVL * 4]), T6);
	       T8 = LD(&(x[WS(rs, 4)]), ms, &(x[0]));
	       T9 = BYTW(&(W[TWVL * 6]), T8);
	       Ta = VSUB(T7, T9);
	       Tj = VADD(T7, T9);
	  }
	  ST(&(x[0]), VADD(Th, VADD(Tk, VADD(Ti, Tj))), ms, &(x[0]));
	  To = VBYI(VFNMS(LDK(KP781831482), Ta, VFNMS(LDK(KP433883739), Tf, VMUL(LDK(KP974927912), T5))));
	  Tp = VFMA(LDK(KP623489801), Tj, VFNMS(LDK(KP900968867), Ti, VFNMS(LDK(KP222520933), Tk, Th)));
	  ST(&(x[WS(rs, 2)]), VADD(To, Tp), ms, &(x[0]));
	  ST(&(x[WS(rs, 5)]), VSUB(Tp, To), ms, &(x[WS(rs, 1)]));
	  {
	       V Tg, Tl, Tm, Tn;
	       Tg = VBYI(VFMA(LDK(KP433883739), T5, VFNMS(LDK(KP781831482), Tf, VMUL(LDK(KP974927912), Ta))));
	       Tl = VFMA(LDK(KP623489801), Ti, VFNMS(LDK(KP222520933), Tj, VFNMS(LDK(KP900968867), Tk, Th)));
	       ST(&(x[WS(rs, 3)]), VADD(Tg, Tl), ms, &(x[WS(rs, 1)]));
	       ST(&(x[WS(rs, 4)]), VSUB(Tl, Tg), ms, &(x[0]));
	       Tm = VBYI(VFMA(LDK(KP781831482), T5, VFMA(LDK(KP974927912), Tf, VMUL(LDK(KP433883739), Ta))));
	       Tn = VFMA(LDK(KP623489801), Tk, VFNMS(LDK(KP900968867), Tj, VFNMS(LDK(KP222520933), Ti, Th)));
	       ST(&(x[WS(rs, 1)]), VADD(Tm, Tn), ms, &(x[WS(rs, 1)]));
	       ST(&(x[WS(rs, 6)]), VSUB(Tn, Tm), ms, &(x[0]));
	  }
     }
}

开发者ID:8cH9azbsFifZ，项目名称:wspr，代码行数:59，代码来源:t1bv_7.c

示例8: q1bv_2

static void q1bv_2(R *ri, R *ii, const R *W, stride rs, stride vs, INT mb, INT me, INT ms)
{
     {
	  INT m;
	  R *x;
	  x = ii;
	  for (m = mb, W = W + (mb * ((TWVL / VL) * 2)); m < me; m = m + VL, x = x + (VL * ms), W = W + (TWVL * 2), MAKE_VOLATILE_STRIDE(4, rs), MAKE_VOLATILE_STRIDE(4, vs)) {
	       V T1, T2, T3, T4, T5, T6;
	       T1 = LD(&(x[0]), ms, &(x[0]));
	       T2 = LD(&(x[WS(rs, 1)]), ms, &(x[WS(rs, 1)]));
	       T3 = BYTW(&(W[0]), VSUB(T1, T2));
	       T4 = LD(&(x[WS(vs, 1)]), ms, &(x[WS(vs, 1)]));
	       T5 = LD(&(x[WS(vs, 1) + WS(rs, 1)]), ms, &(x[WS(vs, 1) + WS(rs, 1)]));
	       T6 = BYTW(&(W[0]), VSUB(T4, T5));
	       ST(&(x[WS(vs, 1)]), T3, ms, &(x[WS(vs, 1)]));
	       ST(&(x[WS(vs, 1) + WS(rs, 1)]), T6, ms, &(x[WS(vs, 1) + WS(rs, 1)]));
	       ST(&(x[0]), VADD(T1, T2), ms, &(x[0]));
	       ST(&(x[WS(rs, 1)]), VADD(T4, T5), ms, &(x[WS(rs, 1)]));
	  }
     }
     VLEAVE();
}

开发者ID:Aegisub，项目名称:fftw3，代码行数:22，代码来源:q1bv_2.c

示例9: t1bv_2

static void t1bv_2(R *ri, R *ii, const R *W, stride rs, INT mb, INT me, INT ms)
{
     INT m;
     R *x;
     x = ii;
     for (m = mb, W = W + (mb * ((TWVL / VL) * 2)); m < me; m = m + VL, x = x + (VL * ms), W = W + (TWVL * 2), MAKE_VOLATILE_STRIDE(rs)) {
	  V T1, T3, T2;
	  T1 = LD(&(x[0]), ms, &(x[0]));
	  T2 = LD(&(x[WS(rs, 1)]), ms, &(x[WS(rs, 1)]));
	  T3 = BYTW(&(W[0]), T2);
	  ST(&(x[WS(rs, 1)]), VSUB(T1, T3), ms, &(x[WS(rs, 1)]));
	  ST(&(x[0]), VADD(T1, T3), ms, &(x[0]));
     }
}

开发者ID:phillipstanleymarbell，项目名称:sunflower-simulator，代码行数:14，代码来源:t1bv_2.c

示例10: LD

static const R *t2bv_2(R *ri, R *ii, const R *W, stride ios, INT m, INT dist)
{
     INT i;
     R *x;
     x = ii;
     for (i = m; i > 0; i = i - VL, x = x + (VL * dist), W = W + (TWVL * 2), MAKE_VOLATILE_STRIDE(ios)) {
	  V T1, T3, T2;
	  T1 = LD(&(x[0]), dist, &(x[0]));
	  T2 = LD(&(x[WS(ios, 1)]), dist, &(x[WS(ios, 1)]));
	  T3 = BYTW(&(W[0]), T2);
	  ST(&(x[WS(ios, 1)]), VSUB(T1, T3), dist, &(x[WS(ios, 1)]));
	  ST(&(x[0]), VADD(T1, T3), dist, &(x[0]));
     }
     return W;
}

开发者ID:exic，项目名称:last.fm-dbus，代码行数:15，代码来源:t2bv_2.c

示例11: t1buv_8

static void t1buv_8(R *ri, R *ii, const R *W, stride rs, INT mb, INT me, INT ms)
{
     DVK(KP707106781, +0.707106781186547524400844362104849039284835938);
     INT m;
     R *x;
     x = ii;
     for (m = mb, W = W + (mb * ((TWVL / VL) * 14)); m < me; m = m + VL, x = x + (VL * ms), W = W + (TWVL * 14), MAKE_VOLATILE_STRIDE(rs)) {
	  V Tl, Tq, Tg, Tr, T5, Tt, Ta, Tu, Ti, Tk, Tj;
	  Ti = LD(&(x[0]), ms, &(x[0]));
	  Tj = LD(&(x[WS(rs, 4)]), ms, &(x[0]));
	  Tk = BYTW(&(W[TWVL * 6]), Tj);
	  Tl = VSUB(Ti, Tk);
	  Tq = VADD(Ti, Tk);
	  {
	       V Td, Tf, Tc, Te;
	       Tc = LD(&(x[WS(rs, 2)]), ms, &(x[0]));
	       Td = BYTW(&(W[TWVL * 2]), Tc);
	       Te = LD(&(x[WS(rs, 6)]), ms, &(x[0]));
	       Tf = BYTW(&(W[TWVL * 10]), Te);
	       Tg = VSUB(Td, Tf);
	       Tr = VADD(Td, Tf);
	  }
	  {
	       V T2, T4, T1, T3;
	       T1 = LD(&(x[WS(rs, 1)]), ms, &(x[WS(rs, 1)]));
	       T2 = BYTW(&(W[0]), T1);
	       T3 = LD(&(x[WS(rs, 5)]), ms, &(x[WS(rs, 1)]));
	       T4 = BYTW(&(W[TWVL * 8]), T3);
	       T5 = VSUB(T2, T4);
	       Tt = VADD(T2, T4);
	  }
	  {
	       V T7, T9, T6, T8;
	       T6 = LD(&(x[WS(rs, 7)]), ms, &(x[WS(rs, 1)]));
	       T7 = BYTW(&(W[TWVL * 12]), T6);
	       T8 = LD(&(x[WS(rs, 3)]), ms, &(x[WS(rs, 1)]));
	       T9 = BYTW(&(W[TWVL * 4]), T8);
	       Ta = VSUB(T7, T9);
	       Tu = VADD(T7, T9);
	  }
	  {
	       V Ts, Tv, Tw, Tx;
	       Ts = VSUB(Tq, Tr);
	       Tv = VBYI(VSUB(Tt, Tu));
	       ST(&(x[WS(rs, 6)]), VSUB(Ts, Tv), ms, &(x[0]));
	       ST(&(x[WS(rs, 2)]), VADD(Ts, Tv), ms, &(x[0]));
	       Tw = VADD(Tq, Tr);
	       Tx = VADD(Tt, Tu);
	       ST(&(x[WS(rs, 4)]), VSUB(Tw, Tx), ms, &(x[0]));
	       ST(&(x[0]), VADD(Tw, Tx), ms, &(x[0]));
	       {
		    V Th, To, Tn, Tp, Tb, Tm;
		    Tb = VMUL(LDK(KP707106781), VSUB(T5, Ta));
		    Th = VBYI(VSUB(Tb, Tg));
		    To = VBYI(VADD(Tg, Tb));
		    Tm = VMUL(LDK(KP707106781), VADD(T5, Ta));
		    Tn = VSUB(Tl, Tm);
		    Tp = VADD(Tl, Tm);
		    ST(&(x[WS(rs, 3)]), VADD(Th, Tn), ms, &(x[WS(rs, 1)]));
		    ST(&(x[WS(rs, 7)]), VSUB(Tp, To), ms, &(x[WS(rs, 1)]));
		    ST(&(x[WS(rs, 5)]), VSUB(Tn, Th), ms, &(x[WS(rs, 1)]));
		    ST(&(x[WS(rs, 1)]), VADD(To, Tp), ms, &(x[WS(rs, 1)]));
	       }
	  }
     }
}

开发者ID:BackupTheBerlios，项目名称:openvsipl，代码行数:66，代码来源:t1buv_8.c

示例12: DVK

static const R *t1bv_9(R *ri, R *ii, const R *W, stride ios, INT m, INT dist)
{
     DVK(KP939692620, +0.939692620785908384054109277324731469936208134);
     DVK(KP907603734, +0.907603734547952313649323976213898122064543220);
     DVK(KP666666666, +0.666666666666666666666666666666666666666666667);
     DVK(KP852868531, +0.852868531952443209628250963940074071936020296);
     DVK(KP879385241, +0.879385241571816768108218554649462939872416269);
     DVK(KP984807753, +0.984807753012208059366743024589523013670643252);
     DVK(KP826351822, +0.826351822333069651148283373230685203999624323);
     DVK(KP347296355, +0.347296355333860697703433253538629592000751354);
     DVK(KP898197570, +0.898197570222573798468955502359086394667167570);
     DVK(KP673648177, +0.673648177666930348851716626769314796000375677);
     DVK(KP420276625, +0.420276625461206169731530603237061658838781920);
     DVK(KP866025403, +0.866025403784438646763723170752936183471402627);
     DVK(KP586256827, +0.586256827714544512072145703099641959914944179);
     DVK(KP968908795, +0.968908795874236621082202410917456709164223497);
     DVK(KP726681596, +0.726681596905677465811651808188092531873167623);
     DVK(KP439692620, +0.439692620785908384054109277324731469936208134);
     DVK(KP203604859, +0.203604859554852403062088995281827210665664861);
     DVK(KP152703644, +0.152703644666139302296566746461370407999248646);
     DVK(KP500000000, +0.500000000000000000000000000000000000000000000);
     INT i;
     R *x;
     x = ii;
     for (i = m; i > 0; i = i - VL, x = x + (VL * dist), W = W + (TWVL * 16), MAKE_VOLATILE_STRIDE(ios)) {
	  V T1, T3, T5, T9, Tn, Tb, Td, Th, Tj, Tx, T6;
	  T1 = LD(&(x[0]), dist, &(x[0]));
	  {
	       V T2, T4, T8, Tm;
	       T2 = LD(&(x[WS(ios, 3)]), dist, &(x[WS(ios, 1)]));
	       T4 = LD(&(x[WS(ios, 6)]), dist, &(x[0]));
	       T8 = LD(&(x[WS(ios, 2)]), dist, &(x[0]));
	       Tm = LD(&(x[WS(ios, 1)]), dist, &(x[WS(ios, 1)]));
	       {
		    V Ta, Tc, Tg, Ti;
		    Ta = LD(&(x[WS(ios, 5)]), dist, &(x[WS(ios, 1)]));
		    Tc = LD(&(x[WS(ios, 8)]), dist, &(x[0]));
		    Tg = LD(&(x[WS(ios, 4)]), dist, &(x[0]));
		    Ti = LD(&(x[WS(ios, 7)]), dist, &(x[WS(ios, 1)]));
		    T3 = BYTW(&(W[TWVL * 4]), T2);
		    T5 = BYTW(&(W[TWVL * 10]), T4);
		    T9 = BYTW(&(W[TWVL * 2]), T8);
		    Tn = BYTW(&(W[0]), Tm);
		    Tb = BYTW(&(W[TWVL * 8]), Ta);
		    Td = BYTW(&(W[TWVL * 14]), Tc);
		    Th = BYTW(&(W[TWVL * 6]), Tg);
		    Tj = BYTW(&(W[TWVL * 12]), Ti);
	       }
	  }
	  Tx = VSUB(T3, T5);
	  T6 = VADD(T3, T5);
	  {
	       V Tl, Te, Tk, To, T7, TN;
	       Tl = VSUB(Td, Tb);
	       Te = VADD(Tb, Td);
	       Tk = VSUB(Th, Tj);
	       To = VADD(Th, Tj);
	       T7 = VFNMS(LDK(KP500000000), T6, T1);
	       TN = VADD(T1, T6);
	       {
		    V Tf, TP, Tp, TO;
		    Tf = VFNMS(LDK(KP500000000), Te, T9);
		    TP = VADD(T9, Te);
		    Tp = VFNMS(LDK(KP500000000), To, Tn);
		    TO = VADD(Tn, To);
		    {
			 V Tz, TC, Tu, TD, TA, Tq, TQ, TS;
			 Tz = VFNMS(LDK(KP152703644), Tl, Tf);
			 TC = VFMA(LDK(KP203604859), Tf, Tl);
			 Tu = VFNMS(LDK(KP439692620), Tk, Tf);
			 TD = VFNMS(LDK(KP726681596), Tk, Tp);
			 TA = VFMA(LDK(KP968908795), Tp, Tk);
			 Tq = VFNMS(LDK(KP586256827), Tp, Tl);
			 TQ = VADD(TO, TP);
			 TS = VMUL(LDK(KP866025403), VSUB(TO, TP));
			 {
			      V TI, TB, TH, TE, Tr, TR, Tw, Tv;
			      Tv = VFNMS(LDK(KP420276625), Tu, Tl);
			      TI = VFMA(LDK(KP673648177), TA, Tz);
			      TB = VFNMS(LDK(KP673648177), TA, Tz);
			      TH = VFNMS(LDK(KP898197570), TD, TC);
			      TE = VFMA(LDK(KP898197570), TD, TC);
			      Tr = VFNMS(LDK(KP347296355), Tq, Tk);
			      ST(&(x[0]), VADD(TQ, TN), dist, &(x[0]));
			      TR = VFNMS(LDK(KP500000000), TQ, TN);
			      Tw = VFNMS(LDK(KP826351822), Tv, Tp);
			      {
				   V TM, TL, TF, TJ, Ts, Ty, TG, TK, Tt;
				   TM = VMUL(LDK(KP984807753), VFMA(LDK(KP879385241), Tx, TI));
				   TL = VFMA(LDK(KP852868531), TE, T7);
				   TF = VFNMS(LDK(KP500000000), TE, TB);
				   TJ = VFMA(LDK(KP666666666), TI, TH);
				   Ts = VFNMS(LDK(KP907603734), Tr, Tf);
				   ST(&(x[WS(ios, 6)]), VFNMSI(TS, TR), dist, &(x[0]));
				   ST(&(x[WS(ios, 3)]), VFMAI(TS, TR), dist, &(x[WS(ios, 1)]));
				   Ty = VMUL(LDK(KP984807753), VFNMS(LDK(KP879385241), Tx, Tw));
				   ST(&(x[WS(ios, 8)]), VFNMSI(TM, TL), dist, &(x[0]));
				   ST(&(x[WS(ios, 1)]), VFMAI(TM, TL), dist, &(x[WS(ios, 1)]));
				   TG = VFMA(LDK(KP852868531), TF, T7);
				   TK = VMUL(LDK(KP866025403), VFNMS(LDK(KP852868531), TJ, Tx));
//.........这里部分代码省略.........

开发者ID:exic，项目名称:last.fm-dbus，代码行数:101，代码来源:t1bv_9.c

示例13: t1bv_15

static void t1bv_15(R *ri, R *ii, const R *W, stride rs, INT mb, INT me, INT ms)
{
     DVK(KP216506350, +0.216506350946109661690930792688234045867850657);
     DVK(KP484122918, +0.484122918275927110647408174972799951354115213);
     DVK(KP866025403, +0.866025403784438646763723170752936183471402627);
     DVK(KP509036960, +0.509036960455127183450980863393907648510733164);
     DVK(KP823639103, +0.823639103546331925877420039278190003029660514);
     DVK(KP951056516, +0.951056516295153572116439333379382143405698634);
     DVK(KP587785252, +0.587785252292473129168705954639072768597652438);
     DVK(KP250000000, +0.250000000000000000000000000000000000000000000);
     DVK(KP559016994, +0.559016994374947424102293417182819058860154590);
     DVK(KP500000000, +0.500000000000000000000000000000000000000000000);
     {
	  INT m;
	  R *x;
	  x = ii;
	  for (m = mb, W = W + (mb * ((TWVL / VL) * 28)); m < me; m = m + VL, x = x + (VL * ms), W = W + (TWVL * 28), MAKE_VOLATILE_STRIDE(rs)) {
	       V Ts, TV, T1f, TZ, T10, Tb, Tm, Tt, T1j, T1k, T1l, TI, TM, TR, Tz;
	       V TD, TQ, T1g, T1h, T1i;
	       {
		    V TT, Tr, Tp, Tq, To, TU;
		    TT = LD(&(x[0]), ms, &(x[0]));
		    Tq = LD(&(x[WS(rs, 10)]), ms, &(x[0]));
		    Tr = BYTW(&(W[TWVL * 18]), Tq);
		    To = LD(&(x[WS(rs, 5)]), ms, &(x[WS(rs, 1)]));
		    Tp = BYTW(&(W[TWVL * 8]), To);
		    Ts = VSUB(Tp, Tr);
		    TU = VADD(Tp, Tr);
		    TV = VFNMS(LDK(KP500000000), TU, TT);
		    T1f = VADD(TT, TU);
	       }
	       {
		    V Tx, TG, TK, TB, T5, Ty, Tg, TH, Tl, TL, Ta, TC;
		    {
			 V Tw, TF, TJ, TA;
			 Tw = LD(&(x[WS(rs, 3)]), ms, &(x[WS(rs, 1)]));
			 Tx = BYTW(&(W[TWVL * 4]), Tw);
			 TF = LD(&(x[WS(rs, 6)]), ms, &(x[0]));
			 TG = BYTW(&(W[TWVL * 10]), TF);
			 TJ = LD(&(x[WS(rs, 9)]), ms, &(x[WS(rs, 1)]));
			 TK = BYTW(&(W[TWVL * 16]), TJ);
			 TA = LD(&(x[WS(rs, 12)]), ms, &(x[0]));
			 TB = BYTW(&(W[TWVL * 22]), TA);
		    }
		    {
			 V T2, T4, T1, T3;
			 T1 = LD(&(x[WS(rs, 8)]), ms, &(x[0]));
			 T2 = BYTW(&(W[TWVL * 14]), T1);
			 T3 = LD(&(x[WS(rs, 13)]), ms, &(x[WS(rs, 1)]));
			 T4 = BYTW(&(W[TWVL * 24]), T3);
			 T5 = VSUB(T2, T4);
			 Ty = VADD(T2, T4);
		    }
		    {
			 V Td, Tf, Tc, Te;
			 Tc = LD(&(x[WS(rs, 11)]), ms, &(x[WS(rs, 1)]));
			 Td = BYTW(&(W[TWVL * 20]), Tc);
			 Te = LD(&(x[WS(rs, 1)]), ms, &(x[WS(rs, 1)]));
			 Tf = BYTW(&(W[0]), Te);
			 Tg = VSUB(Td, Tf);
			 TH = VADD(Td, Tf);
		    }
		    {
			 V Ti, Tk, Th, Tj;
			 Th = LD(&(x[WS(rs, 14)]), ms, &(x[0]));
			 Ti = BYTW(&(W[TWVL * 26]), Th);
			 Tj = LD(&(x[WS(rs, 4)]), ms, &(x[0]));
			 Tk = BYTW(&(W[TWVL * 6]), Tj);
			 Tl = VSUB(Ti, Tk);
			 TL = VADD(Ti, Tk);
		    }
		    {
			 V T7, T9, T6, T8;
			 T6 = LD(&(x[WS(rs, 2)]), ms, &(x[0]));
			 T7 = BYTW(&(W[TWVL * 2]), T6);
			 T8 = LD(&(x[WS(rs, 7)]), ms, &(x[WS(rs, 1)]));
			 T9 = BYTW(&(W[TWVL * 12]), T8);
			 Ta = VSUB(T7, T9);
			 TC = VADD(T7, T9);
		    }
		    TZ = VSUB(T5, Ta);
		    T10 = VSUB(Tg, Tl);
		    Tb = VADD(T5, Ta);
		    Tm = VADD(Tg, Tl);
		    Tt = VADD(Tb, Tm);
		    T1j = VADD(TG, TH);
		    T1k = VADD(TK, TL);
		    T1l = VADD(T1j, T1k);
		    TI = VFNMS(LDK(KP500000000), TH, TG);
		    TM = VFNMS(LDK(KP500000000), TL, TK);
		    TR = VADD(TI, TM);
		    Tz = VFNMS(LDK(KP500000000), Ty, Tx);
		    TD = VFNMS(LDK(KP500000000), TC, TB);
		    TQ = VADD(Tz, TD);
		    T1g = VADD(Tx, Ty);
		    T1h = VADD(TB, TC);
		    T1i = VADD(T1g, T1h);
	       }
	       {
		    V T1o, T1m, T1n, T1s, T1t, T1q, T1r, T1u, T1p;
//.........这里部分代码省略.........

开发者ID:AlexanderB1290，项目名称:scafacos，代码行数:101，代码来源:t1bv_15.c

示例14: t1bv_16

static void t1bv_16(R *ri, R *ii, const R *W, stride rs, INT mb, INT me, INT ms)
{
     DVK(KP382683432, +0.382683432365089771728459984030398866761344562);
     DVK(KP923879532, +0.923879532511286756128183189396788286822416626);
     DVK(KP707106781, +0.707106781186547524400844362104849039284835938);
     {
	  INT m;
	  R *x;
	  x = ii;
	  for (m = mb, W = W + (mb * ((TWVL / VL) * 30)); m < me; m = m + VL, x = x + (VL * ms), W = W + (TWVL * 30), MAKE_VOLATILE_STRIDE(rs)) {
	       V TJ, T1b, TD, T1c, T17, T18, Ty, TK, T10, T11, T12, Tb, TM, T13, T14;
	       V T15, Tm, TN, TG, TI, TH;
	       TG = LD(&(x[0]), ms, &(x[0]));
	       TH = LD(&(x[WS(rs, 8)]), ms, &(x[0]));
	       TI = BYTW(&(W[TWVL * 14]), TH);
	       TJ = VSUB(TG, TI);
	       T1b = VADD(TG, TI);
	       {
		    V TA, TC, Tz, TB;
		    Tz = LD(&(x[WS(rs, 4)]), ms, &(x[0]));
		    TA = BYTW(&(W[TWVL * 6]), Tz);
		    TB = LD(&(x[WS(rs, 12)]), ms, &(x[0]));
		    TC = BYTW(&(W[TWVL * 22]), TB);
		    TD = VSUB(TA, TC);
		    T1c = VADD(TA, TC);
	       }
	       {
		    V Tp, Tw, Tr, Tu, Ts, Tx;
		    {
			 V To, Tv, Tq, Tt;
			 To = LD(&(x[WS(rs, 2)]), ms, &(x[0]));
			 Tp = BYTW(&(W[TWVL * 2]), To);
			 Tv = LD(&(x[WS(rs, 6)]), ms, &(x[0]));
			 Tw = BYTW(&(W[TWVL * 10]), Tv);
			 Tq = LD(&(x[WS(rs, 10)]), ms, &(x[0]));
			 Tr = BYTW(&(W[TWVL * 18]), Tq);
			 Tt = LD(&(x[WS(rs, 14)]), ms, &(x[0]));
			 Tu = BYTW(&(W[TWVL * 26]), Tt);
		    }
		    T17 = VADD(Tp, Tr);
		    T18 = VADD(Tu, Tw);
		    Ts = VSUB(Tp, Tr);
		    Tx = VSUB(Tu, Tw);
		    Ty = VMUL(LDK(KP707106781), VSUB(Ts, Tx));
		    TK = VMUL(LDK(KP707106781), VADD(Ts, Tx));
	       }
	       {
		    V T2, T9, T4, T7, T5, Ta;
		    {
			 V T1, T8, T3, T6;
			 T1 = LD(&(x[WS(rs, 1)]), ms, &(x[WS(rs, 1)]));
			 T2 = BYTW(&(W[0]), T1);
			 T8 = LD(&(x[WS(rs, 13)]), ms, &(x[WS(rs, 1)]));
			 T9 = BYTW(&(W[TWVL * 24]), T8);
			 T3 = LD(&(x[WS(rs, 9)]), ms, &(x[WS(rs, 1)]));
			 T4 = BYTW(&(W[TWVL * 16]), T3);
			 T6 = LD(&(x[WS(rs, 5)]), ms, &(x[WS(rs, 1)]));
			 T7 = BYTW(&(W[TWVL * 8]), T6);
		    }
		    T10 = VADD(T2, T4);
		    T11 = VADD(T7, T9);
		    T12 = VSUB(T10, T11);
		    T5 = VSUB(T2, T4);
		    Ta = VSUB(T7, T9);
		    Tb = VFNMS(LDK(KP382683432), Ta, VMUL(LDK(KP923879532), T5));
		    TM = VFMA(LDK(KP382683432), T5, VMUL(LDK(KP923879532), Ta));
	       }
	       {
		    V Td, Tk, Tf, Ti, Tg, Tl;
		    {
			 V Tc, Tj, Te, Th;
			 Tc = LD(&(x[WS(rs, 15)]), ms, &(x[WS(rs, 1)]));
			 Td = BYTW(&(W[TWVL * 28]), Tc);
			 Tj = LD(&(x[WS(rs, 11)]), ms, &(x[WS(rs, 1)]));
			 Tk = BYTW(&(W[TWVL * 20]), Tj);
			 Te = LD(&(x[WS(rs, 7)]), ms, &(x[WS(rs, 1)]));
			 Tf = BYTW(&(W[TWVL * 12]), Te);
			 Th = LD(&(x[WS(rs, 3)]), ms, &(x[WS(rs, 1)]));
			 Ti = BYTW(&(W[TWVL * 4]), Th);
		    }
		    T13 = VADD(Td, Tf);
		    T14 = VADD(Ti, Tk);
		    T15 = VSUB(T13, T14);
		    Tg = VSUB(Td, Tf);
		    Tl = VSUB(Ti, Tk);
		    Tm = VFMA(LDK(KP923879532), Tg, VMUL(LDK(KP382683432), Tl));
		    TN = VFNMS(LDK(KP382683432), Tg, VMUL(LDK(KP923879532), Tl));
	       }
	       {
		    V T1a, T1g, T1f, T1h;
		    {
			 V T16, T19, T1d, T1e;
			 T16 = VMUL(LDK(KP707106781), VSUB(T12, T15));
			 T19 = VSUB(T17, T18);
			 T1a = VBYI(VSUB(T16, T19));
			 T1g = VBYI(VADD(T19, T16));
			 T1d = VSUB(T1b, T1c);
			 T1e = VMUL(LDK(KP707106781), VADD(T12, T15));
			 T1f = VSUB(T1d, T1e);
			 T1h = VADD(T1d, T1e);
//.........这里部分代码省略.........

开发者ID:AlexanderB1290，项目名称:scafacos，代码行数:101，代码来源:t1bv_16.c

示例15: t1bv_32

static void t1bv_32(R *ri, R *ii, const R *W, stride rs, INT mb, INT me, INT ms)
{
     DVK(KP195090322, +0.195090322016128267848284868477022240927691618);
     DVK(KP980785280, +0.980785280403230449126182236134239036973933731);
     DVK(KP555570233, +0.555570233019602224742830813948532874374937191);
     DVK(KP831469612, +0.831469612302545237078788377617905756738560812);
     DVK(KP382683432, +0.382683432365089771728459984030398866761344562);
     DVK(KP923879532, +0.923879532511286756128183189396788286822416626);
     DVK(KP707106781, +0.707106781186547524400844362104849039284835938);
     INT m;
     R *x;
     x = ii;
     for (m = mb, W = W + (mb * ((TWVL / VL) * 62)); m < me; m = m + VL, x = x + (VL * ms), W = W + (TWVL * 62), MAKE_VOLATILE_STRIDE(rs)) {
	  V T4, T1D, T2P, T3h, Tf, T1y, T2K, T3i, TC, T1w, T2G, T3e, Tr, T1v, T2D;
	  V T3d, T1k, T20, T2y, T3a, T1r, T21, T2v, T39, TV, T1X, T2r, T37, T12, T1Y;
	  V T2o, T36;
	  {
	       V T1, T1C, T3, T1A, T1B, T2, T1z, T2N, T2O;
	       T1 = LD(&(x[0]), ms, &(x[0]));
	       T1B = LD(&(x[WS(rs, 24)]), ms, &(x[0]));
	       T1C = BYTW(&(W[TWVL * 46]), T1B);
	       T2 = LD(&(x[WS(rs, 16)]), ms, &(x[0]));
	       T3 = BYTW(&(W[TWVL * 30]), T2);
	       T1z = LD(&(x[WS(rs, 8)]), ms, &(x[0]));
	       T1A = BYTW(&(W[TWVL * 14]), T1z);
	       T4 = VSUB(T1, T3);
	       T1D = VSUB(T1A, T1C);
	       T2N = VADD(T1, T3);
	       T2O = VADD(T1A, T1C);
	       T2P = VSUB(T2N, T2O);
	       T3h = VADD(T2N, T2O);
	  }
	  {
	       V T6, Td, T8, Tb;
	       {
		    V T5, Tc, T7, Ta;
		    T5 = LD(&(x[WS(rs, 4)]), ms, &(x[0]));
		    T6 = BYTW(&(W[TWVL * 6]), T5);
		    Tc = LD(&(x[WS(rs, 12)]), ms, &(x[0]));
		    Td = BYTW(&(W[TWVL * 22]), Tc);
		    T7 = LD(&(x[WS(rs, 20)]), ms, &(x[0]));
		    T8 = BYTW(&(W[TWVL * 38]), T7);
		    Ta = LD(&(x[WS(rs, 28)]), ms, &(x[0]));
		    Tb = BYTW(&(W[TWVL * 54]), Ta);
	       }
	       {
		    V T9, Te, T2I, T2J;
		    T9 = VSUB(T6, T8);
		    Te = VSUB(Tb, Td);
		    Tf = VMUL(LDK(KP707106781), VADD(T9, Te));
		    T1y = VMUL(LDK(KP707106781), VSUB(T9, Te));
		    T2I = VADD(T6, T8);
		    T2J = VADD(Tb, Td);
		    T2K = VSUB(T2I, T2J);
		    T3i = VADD(T2I, T2J);
	       }
	  }
	  {
	       V Tt, TA, Tv, Ty;
	       {
		    V Ts, Tz, Tu, Tx;
		    Ts = LD(&(x[WS(rs, 6)]), ms, &(x[0]));
		    Tt = BYTW(&(W[TWVL * 10]), Ts);
		    Tz = LD(&(x[WS(rs, 14)]), ms, &(x[0]));
		    TA = BYTW(&(W[TWVL * 26]), Tz);
		    Tu = LD(&(x[WS(rs, 22)]), ms, &(x[0]));
		    Tv = BYTW(&(W[TWVL * 42]), Tu);
		    Tx = LD(&(x[WS(rs, 30)]), ms, &(x[0]));
		    Ty = BYTW(&(W[TWVL * 58]), Tx);
	       }
	       {
		    V Tw, TB, T2E, T2F;
		    Tw = VSUB(Tt, Tv);
		    TB = VSUB(Ty, TA);
		    TC = VFNMS(LDK(KP382683432), TB, VMUL(LDK(KP923879532), Tw));
		    T1w = VFMA(LDK(KP923879532), TB, VMUL(LDK(KP382683432), Tw));
		    T2E = VADD(Ty, TA);
		    T2F = VADD(Tt, Tv);
		    T2G = VSUB(T2E, T2F);
		    T3e = VADD(T2E, T2F);
	       }
	  }
	  {
	       V Ti, Tp, Tk, Tn;
	       {
		    V Th, To, Tj, Tm;
		    Th = LD(&(x[WS(rs, 2)]), ms, &(x[0]));
		    Ti = BYTW(&(W[TWVL * 2]), Th);
		    To = LD(&(x[WS(rs, 26)]), ms, &(x[0]));
		    Tp = BYTW(&(W[TWVL * 50]), To);
		    Tj = LD(&(x[WS(rs, 18)]), ms, &(x[0]));
		    Tk = BYTW(&(W[TWVL * 34]), Tj);
		    Tm = LD(&(x[WS(rs, 10)]), ms, &(x[0]));
		    Tn = BYTW(&(W[TWVL * 18]), Tm);
	       }
	       {
		    V Tl, Tq, T2B, T2C;
		    Tl = VSUB(Ti, Tk);
		    Tq = VSUB(Tn, Tp);
		    Tr = VFMA(LDK(KP382683432), Tl, VMUL(LDK(KP923879532), Tq));
//.........这里部分代码省略.........

开发者ID:phillipstanleymarbell，项目名称:sunflower-simulator，代码行数:101，代码来源:t1bv_32.c

示例16: t2bv_10

static void t2bv_10(R *ri, R *ii, const R *W, stride rs, INT mb, INT me, INT ms)
{
     DVK(KP587785252, +0.587785252292473129168705954639072768597652438);
     DVK(KP951056516, +0.951056516295153572116439333379382143405698634);
     DVK(KP250000000, +0.250000000000000000000000000000000000000000000);
     DVK(KP559016994, +0.559016994374947424102293417182819058860154590);
     {
	  INT m;
	  R *x;
	  x = ii;
	  for (m = mb, W = W + (mb * ((TWVL / VL) * 18)); m < me; m = m + VL, x = x + (VL * ms), W = W + (TWVL * 18), MAKE_VOLATILE_STRIDE(10, rs)) {
	       V Tu, TH, Tg, Tl, Tp, TD, TE, TJ, T5, Ta, To, TA, TB, TI, Tr;
	       V Tt, Ts;
	       Tr = LD(&(x[0]), ms, &(x[0]));
	       Ts = LD(&(x[WS(rs, 5)]), ms, &(x[WS(rs, 1)]));
	       Tt = BYTW(&(W[TWVL * 8]), Ts);
	       Tu = VSUB(Tr, Tt);
	       TH = VADD(Tr, Tt);
	       {
		    V Td, Tk, Tf, Ti;
		    {
			 V Tc, Tj, Te, Th;
			 Tc = LD(&(x[WS(rs, 4)]), ms, &(x[0]));
			 Td = BYTW(&(W[TWVL * 6]), Tc);
			 Tj = LD(&(x[WS(rs, 1)]), ms, &(x[WS(rs, 1)]));
			 Tk = BYTW(&(W[0]), Tj);
			 Te = LD(&(x[WS(rs, 9)]), ms, &(x[WS(rs, 1)]));
			 Tf = BYTW(&(W[TWVL * 16]), Te);
			 Th = LD(&(x[WS(rs, 6)]), ms, &(x[0]));
			 Ti = BYTW(&(W[TWVL * 10]), Th);
		    }
		    Tg = VSUB(Td, Tf);
		    Tl = VSUB(Ti, Tk);
		    Tp = VADD(Tg, Tl);
		    TD = VADD(Td, Tf);
		    TE = VADD(Ti, Tk);
		    TJ = VADD(TD, TE);
	       }
	       {
		    V T2, T9, T4, T7;
		    {
			 V T1, T8, T3, T6;
			 T1 = LD(&(x[WS(rs, 2)]), ms, &(x[0]));
			 T2 = BYTW(&(W[TWVL * 2]), T1);
			 T8 = LD(&(x[WS(rs, 3)]), ms, &(x[WS(rs, 1)]));
			 T9 = BYTW(&(W[TWVL * 4]), T8);
			 T3 = LD(&(x[WS(rs, 7)]), ms, &(x[WS(rs, 1)]));
			 T4 = BYTW(&(W[TWVL * 12]), T3);
			 T6 = LD(&(x[WS(rs, 8)]), ms, &(x[0]));
			 T7 = BYTW(&(W[TWVL * 14]), T6);
		    }
		    T5 = VSUB(T2, T4);
		    Ta = VSUB(T7, T9);
		    To = VADD(T5, Ta);
		    TA = VADD(T2, T4);
		    TB = VADD(T7, T9);
		    TI = VADD(TA, TB);
	       }
	       {
		    V Tq, Tv, Tw, Tn, Tz, Tb, Tm, Ty, Tx;
		    Tq = VMUL(LDK(KP559016994), VSUB(To, Tp));
		    Tv = VADD(To, Tp);
		    Tw = VFNMS(LDK(KP250000000), Tv, Tu);
		    Tb = VSUB(T5, Ta);
		    Tm = VSUB(Tg, Tl);
		    Tn = VBYI(VFMA(LDK(KP951056516), Tb, VMUL(LDK(KP587785252), Tm)));
		    Tz = VBYI(VFNMS(LDK(KP951056516), Tm, VMUL(LDK(KP587785252), Tb)));
		    ST(&(x[WS(rs, 5)]), VADD(Tu, Tv), ms, &(x[WS(rs, 1)]));
		    Ty = VSUB(Tw, Tq);
		    ST(&(x[WS(rs, 3)]), VSUB(Ty, Tz), ms, &(x[WS(rs, 1)]));
		    ST(&(x[WS(rs, 7)]), VADD(Tz, Ty), ms, &(x[WS(rs, 1)]));
		    Tx = VADD(Tq, Tw);
		    ST(&(x[WS(rs, 1)]), VADD(Tn, Tx), ms, &(x[WS(rs, 1)]));
		    ST(&(x[WS(rs, 9)]), VSUB(Tx, Tn), ms, &(x[WS(rs, 1)]));
	       }
	       {
		    V TM, TK, TL, TG, TP, TC, TF, TO, TN;
		    TM = VMUL(LDK(KP559016994), VSUB(TI, TJ));
		    TK = VADD(TI, TJ);
		    TL = VFNMS(LDK(KP250000000), TK, TH);
		    TC = VSUB(TA, TB);
		    TF = VSUB(TD, TE);
		    TG = VBYI(VFNMS(LDK(KP951056516), TF, VMUL(LDK(KP587785252), TC)));
		    TP = VBYI(VFMA(LDK(KP951056516), TC, VMUL(LDK(KP587785252), TF)));
		    ST(&(x[0]), VADD(TH, TK), ms, &(x[0]));
		    TO = VADD(TM, TL);
		    ST(&(x[WS(rs, 4)]), VSUB(TO, TP), ms, &(x[0]));
		    ST(&(x[WS(rs, 6)]), VADD(TP, TO), ms, &(x[0]));
		    TN = VSUB(TL, TM);
		    ST(&(x[WS(rs, 2)]), VADD(TG, TN), ms, &(x[0]));
		    ST(&(x[WS(rs, 8)]), VSUB(TN, TG), ms, &(x[0]));
	       }
	  }
     }
     VLEAVE();
}

开发者ID:SKA-ScienceDataProcessor，项目名称:FastImaging，代码行数:96，代码来源:t2bv_10.c

示例17: t1bv_12

static void t1bv_12(R *ri, R *ii, const R *W, stride rs, INT mb, INT me, INT ms)
{
     DVK(KP866025403, +0.866025403784438646763723170752936183471402627);
     DVK(KP500000000, +0.500000000000000000000000000000000000000000000);
     INT m;
     R *x;
     x = ii;
     for (m = mb, W = W + (mb * ((TWVL / VL) * 22)); m < me; m = m + VL, x = x + (VL * ms), W = W + (TWVL * 22), MAKE_VOLATILE_STRIDE(rs)) {
	  V TI, Ti, TA, T7, Tm, TE, Tw, Tk, Tf, TB, TU, TM;
	  {
	       V T9, TK, Tj, TL, Te;
	       {
		    V T1, T4, T2, Tp, Tt, Tr;
		    T1 = LD(&(x[0]), ms, &(x[0]));
		    T4 = LD(&(x[WS(rs, 8)]), ms, &(x[0]));
		    T2 = LD(&(x[WS(rs, 4)]), ms, &(x[0]));
		    Tp = LD(&(x[WS(rs, 9)]), ms, &(x[WS(rs, 1)]));
		    Tt = LD(&(x[WS(rs, 5)]), ms, &(x[WS(rs, 1)]));
		    Tr = LD(&(x[WS(rs, 1)]), ms, &(x[WS(rs, 1)]));
		    {
			 V T5, T3, Tq, Tu, Ts, Td, Tb, T8, Tc, Ta;
			 T8 = LD(&(x[WS(rs, 6)]), ms, &(x[0]));
			 Tc = LD(&(x[WS(rs, 2)]), ms, &(x[0]));
			 Ta = LD(&(x[WS(rs, 10)]), ms, &(x[0]));
			 T5 = BYTW(&(W[TWVL * 14]), T4);
			 T3 = BYTW(&(W[TWVL * 6]), T2);
			 Tq = BYTW(&(W[TWVL * 16]), Tp);
			 Tu = BYTW(&(W[TWVL * 8]), Tt);
			 Ts = BYTW(&(W[0]), Tr);
			 T9 = BYTW(&(W[TWVL * 10]), T8);
			 Td = BYTW(&(W[TWVL * 2]), Tc);
			 Tb = BYTW(&(W[TWVL * 18]), Ta);
			 {
			      V Th, T6, Tl, Tv;
			      Th = LD(&(x[WS(rs, 3)]), ms, &(x[WS(rs, 1)]));
			      TK = VSUB(T3, T5);
			      T6 = VADD(T3, T5);
			      Tl = LD(&(x[WS(rs, 11)]), ms, &(x[WS(rs, 1)]));
			      Tv = VADD(Ts, Tu);
			      TI = VSUB(Tu, Ts);
			      Tj = LD(&(x[WS(rs, 7)]), ms, &(x[WS(rs, 1)]));
			      TL = VSUB(Tb, Td);
			      Te = VADD(Tb, Td);
			      Ti = BYTW(&(W[TWVL * 4]), Th);
			      TA = VFNMS(LDK(KP500000000), T6, T1);
			      T7 = VADD(T1, T6);
			      Tm = BYTW(&(W[TWVL * 20]), Tl);
			      TE = VFNMS(LDK(KP500000000), Tv, Tq);
			      Tw = VADD(Tq, Tv);
			 }
		    }
	       }
	       Tk = BYTW(&(W[TWVL * 12]), Tj);
	       Tf = VADD(T9, Te);
	       TB = VFNMS(LDK(KP500000000), Te, T9);
	       TU = VSUB(TK, TL);
	       TM = VADD(TK, TL);
	  }
	  {
	       V Tn, TH, TC, TQ, Ty, Tg;
	       Tn = VADD(Tk, Tm);
	       TH = VSUB(Tk, Tm);
	       TC = VADD(TA, TB);
	       TQ = VSUB(TA, TB);
	       Ty = VADD(T7, Tf);
	       Tg = VSUB(T7, Tf);
	       {
		    V To, TD, TJ, TR;
		    To = VADD(Ti, Tn);
		    TD = VFNMS(LDK(KP500000000), Tn, Ti);
		    TJ = VSUB(TH, TI);
		    TR = VADD(TH, TI);
		    {
			 V TP, TN, TW, TS, TO, TG, TX, TV;
			 {
			      V Tz, Tx, TF, TT;
			      Tz = VADD(To, Tw);
			      Tx = VSUB(To, Tw);
			      TF = VADD(TD, TE);
			      TT = VSUB(TD, TE);
			      TP = VMUL(LDK(KP866025403), VADD(TM, TJ));
			      TN = VMUL(LDK(KP866025403), VSUB(TJ, TM));
			      TW = VFMA(LDK(KP866025403), TR, TQ);
			      TS = VFNMS(LDK(KP866025403), TR, TQ);
			      ST(&(x[WS(rs, 6)]), VSUB(Ty, Tz), ms, &(x[0]));
			      ST(&(x[0]), VADD(Ty, Tz), ms, &(x[0]));
			      ST(&(x[WS(rs, 9)]), VFMAI(Tx, Tg), ms, &(x[WS(rs, 1)]));
			      ST(&(x[WS(rs, 3)]), VFNMSI(Tx, Tg), ms, &(x[WS(rs, 1)]));
			      TO = VADD(TC, TF);
			      TG = VSUB(TC, TF);
			      TX = VFNMS(LDK(KP866025403), TU, TT);
			      TV = VFMA(LDK(KP866025403), TU, TT);
			 }
			 ST(&(x[WS(rs, 8)]), VFNMSI(TP, TO), ms, &(x[0]));
			 ST(&(x[WS(rs, 4)]), VFMAI(TP, TO), ms, &(x[0]));
			 ST(&(x[WS(rs, 2)]), VFMAI(TN, TG), ms, &(x[0]));
			 ST(&(x[WS(rs, 10)]), VFNMSI(TN, TG), ms, &(x[0]));
			 ST(&(x[WS(rs, 5)]), VFMAI(TX, TW), ms, &(x[WS(rs, 1)]));
			 ST(&(x[WS(rs, 7)]), VFNMSI(TX, TW), ms, &(x[WS(rs, 1)]));
			 ST(&(x[WS(rs, 11)]), VFNMSI(TV, TS), ms, &(x[WS(rs, 1)]));
//.........这里部分代码省略.........

开发者ID:phillipstanleymarbell，项目名称:sunflower-simulator，代码行数:101，代码来源:t1bv_12.c

示例18: t2bv_20

该文章已有0人参与评论

请发表评论

全部评论

专题导读

More+

10-27 六六分期app的软件客服如何联系？(六六分期

11-06 可心卡盟:win10系统火狐flash插件崩溃怎么

11-06 亲亲特价:怎么删除回收站图标

11-06 济南大学虚拟社区:鲁大师节能降温的具体办

11-06 xlueops.exe:无线网络安装向导

11-06 女斗合众国:win7系统cf与主机连接不稳定怎

11-06 0xc000022-[cf烟雾头]cf怎么调烟雾头

11-06 qizideyouhuo:应用程序无法正常启动0xc0000

11-06 ipz-185:win7系统vcf文件怎么打开

11-06 傻哥蹦迪:win10系统s4怎么打开usb调试

11-06 八神浩树gtaste:回收站清空了怎么恢复

11-06 妖尾之黑色守护:win10系统电脑没有1440x900

11-06 校园至尊魔王小说:win7系统浏览网页时字体

11-06 女斗合众国:win10系统访问共享文件夹提示请

11-06 tokyo hot n0654:恢复win7系统默认字体一招

11-06 雨酷仙境:设置win7系统转移临时文件夹腾出

11-06 阿穆纳伊之杖:win7系统开始菜单在右边还原

11-06 tunespotting:win10系统火狐flash插件总是

11-06 甘尔葛分析师：计谋网站seo关键词暴涨有什

11-06 蔡贵霖: 计谋网站seo关键词暴涨有什么秘密

11-06 博益网首页:ao3网页版进入不了解决方法

11-06 漏斗子专栏: 网站数据分析小白易懂精华篇

11-06 见证双虹怎么做:win7系统开启telnet命令的

11-06 颾狐蝶蜋:系统资源不足无法完成请求的服务

11-06 国光中学校歌:提交网站到alexa查询详细步骤

11-06 西安有情天:静态网页和动态网页的区别

11-06 红木雅尚斋:外部链接构造对网站的好处

11-06 前官礼遇：防止域名劫持–增强域安全性的10

11-06 密传二转答案: 中文分词算法有哪些

11-06 金泉家园邮编:百度快照劫持的表现及应对方

C++ BYTWJ函数代码示例发布时间：2022-05-30

C++ BYTE_TYPE函数代码示例发布时间：2022-05-30

剪的笔顺,诠释剪的笔画,认识剪的部首

六六分期app的软件客服如何联系？(六六分期

2023-10-27

florent37/ViewAnimator: A fluent Android

2022-08-15

florent37/Shrine-MaterialDesign2: implem

2022-08-17

CVE-2020-36276

2022-09-23

SimpleSoftwareIO/simple-sms: Send and re

2022-08-13

阅读排行榜

1 六六分期app的软件客服如何联系？(六六分期

六六分期app的软件客服如何联系？不知道吗？加qq群【895510560】即可！标题：六六分期

阅读：18860|2023-10-27

2 可心卡盟:win10系统火狐flash插件崩溃怎么

今天小编告诉大家如何处理win10系统火狐flash插件总是崩溃的问题，可能很多用户都不知

阅读：9882|2022-11-06

3 亲亲特价:怎么删除回收站图标

今天小编告诉大家如何对win10系统删除桌面回收站图标进行设置，可能很多用户都不知道

阅读：8285|2022-11-06

4 济南大学虚拟社区:鲁大师节能降温的具体办

今天小编告诉大家如何对win10系统电脑设置节能降温的设置方法，想必大家都遇到过需要

阅读：8643|2022-11-06

5 xlueops.exe:无线网络安装向导

我们在使用xp系统的过程中,经常需要对xp系统无线网络安装向导设置进行设置，可能很多

阅读：8574|2022-11-06

6 女斗合众国:win7系统cf与主机连接不稳定怎

今天小编告诉大家如何处理win7系统玩cf老是与主机连接不稳定的问题，可能很多用户都不

阅读：9578|2022-11-06

7 0xc000022-[cf烟雾头]cf怎么调烟雾头

电脑对日常生活的重要性小编就不多说了，可是一旦碰到win7系统设置cf烟雾头的问题，很

阅读：8570|2022-11-06

8 qizideyouhuo:应用程序无法正常启动0xc0000

我们在日常使用电脑的时候，有的小伙伴们可能在打开应用的时候会遇见提示应用程序无法

阅读：7960|2022-11-06

9 ipz-185:win7系统vcf文件怎么打开

今天小编告诉大家如何对win7系统打开vcf文件进行设置，可能很多用户都不知道怎么对win

阅读：8575|2022-11-06

10 傻哥蹦迪:win10系统s4怎么打开usb调试

今天小编告诉大家如何对win10系统s4开启USB调试模式进行设置，可能很多用户都不知道怎

阅读：7499|2022-11-06

扫描微信二维码

查看手机版网站

随时了解更新最新资讯

139-2527-9053

在线客服（服务时间 9:00～18:00）

在线QQ客服

地址：深圳市南山区西丽大学城创智工业园

电邮：jeky_zhao#qq.com

移动电话：139-2527-9053

客服电话

电子邮件