camm_strat1.h

基于Blas CLapck的.用过的人知道是干啥的

      pm(1,1) \
      pm(2,2) \
      pa(1,0) \
      pm(2,0) 
#undef p4_0x1_nrm2_1
#define p4_0x1_nrm2_1(a_) \
      f(nta,SS(a_,MM((SS(2,CL)),RS4)),ax) \
      plq(SS(a_,MM(2,RS4)),ax,1) \
      pm(3,3) \
      pa(7,0) \
      plq(SS(a_,MM(3,RS4)),ax,2) \
      pm(1,1) \
      pa(3,0) \
      f(nta,SS(a_,MM((SS(4,CL)),RS4)),ax) \
      plq(SS(a_,MM(4,RS4)),ax,7) \
      pm(2,2) \
      pa(1,0) \
      plq(SS(a_,MM(5,RS4)),ax,3) \
      pm(7,7) \
      pa(2,0) 
#undef lp0x1_nrm2_1
#define lp0x1_nrm2_1(a_) \
      plq(a_,ax,7) \
      plq(SS(a_,MM(1,RS4)),ax,3) \
      pm(7,7) 
#undef dp0x1_nrm2_1
#define dp0x1_nrm2_1(a_) \
      plq(SS(a_,MM(2,RS4)),ax,1) \
      pm(3,3) \
      pa(7,0) \
      plq(SS(a_,MM(3,RS4)),ax,2) \
      pm(1,1) \
      pa(3,0) \
      pm(2,2) \
      pa(1,0) \
      pa(2,0) 
#undef pl0x1_nrm2_1
#define pl0x1_nrm2_1 RS4


#undef p1_4_nrm2_2
#define p1_4_nrm2_2(a_) \
      pls(a_,ax,1) dbg(1) \
      pan(4,1) dbg(1) \
      pcs(5,6) dbg(6) \
      pcs(5,7) dbg(7) \
      paxs(1,5) dbg(5) \
      prps(5,2) dbg(2) \
      px(3) \
      pcms(0,2,3) dbg(3) \
      pan(3,7) dbg(7) \
      pann(5,3) dbg(3) \
      pasr(3,7) dbg(7) \
      pcs(7,5) dbg(5) \
      pdsr(5,6) dbg(6) \
      pdsr(5,1) dbg(1) \
      pmsr(6,6) dbg(6) \
      pmsr(1,1) dbg(1) \
      pm(6,0) dbg(0) \
      pasr(1,0) dbg(0)
#undef p1_2_nrm2_2
#define p1_2_nrm2_2(a_) \
      px(1) pld(a_,ax,1) dbg(1) \
      pan(4,1) dbg(1) \
      pc(5,6) dbg(6) \
      pc(5,7) dbg(7) \
      pax(1,5) dbg(5) \
      prp(5,2) dbg(2) \
      px(3) \
      pcm(0,2,3)dbg(3) \
      pan(3,7) dbg(7) \
      pann(5,3) dbg(3) \
      pa(3,7) dbg(7) \
      pc(7,5) dbg(5) \
      pd(5,6) dbg(6) \
      pd(5,1) dbg(1) \
      pm(6,6) dbg(6) \
      pm(1,1) dbg(1) \
      pm(6,0) dbg(0) \
      pa(1,0) dbg(0)
#undef p1_nrm2_2
#define p1_nrm2_2(a_) \
      plq(a_,ax,1) dbg(1) \
      pan(4,1) dbg(1) \
      pc(5,6) dbg(6) \
      pc(5,7) dbg(7) \
      pax(1,5) dbg(5) \
      prp(5,2) dbg(2) \
      px(3) \
      pcm(0,2,3)dbg(3) \
      pan(3,7) dbg(7) \
      pann(5,3) dbg(3) \
      pa(3,7) dbg(7) \
      pc(7,5) dbg(5) \
      pd(5,6) dbg(6) \
      pd(5,1) dbg(1) \
      pm(6,6) dbg(6) \
      pm(1,1) dbg(1) \
      pm(6,0) dbg(0) \
      pa(1,0) dbg(0)
#define p2_nrm2_2(a_) \
      plq(SS(a_,RS4),ax,1) dbg(1) \
      pan(4,1) dbg(1) \
      pc(5,6) dbg(6) \
      pc(5,7) dbg(7) \
      pax(1,5) dbg(5) \
      prp(5,2) dbg(2) \
      px(3) \
      pcm(0,2,3)dbg(3) \
      pan(3,7) dbg(7) \
      pann(5,3) dbg(3) \
      pa(3,7) dbg(7) \
      pc(7,5) dbg(5) \
      pd(5,6) dbg(6) \
      pd(5,1) dbg(1) \
      pm(6,6) dbg(6) \
      pm(1,1) dbg(1) \
      pm(6,0) dbg(0) \
      pa(1,0) dbg(0) \
      f(nta,SS(a_,MM((SS(2,CL)),RS4)),ax) \
      plq(SS(a_,MM(2,RS4)),ax,1) dbg(1) \
      pan(4,1) dbg(1) \
      pc(5,6) dbg(6) \
      pc(5,7) dbg(7) \
      pax(1,5) dbg(5) \
      prp(5,2) dbg(2) \
      px(3) \
      pcm(0,2,3)dbg(3) \
      pan(3,7) dbg(7) \
      pann(5,3) dbg(3) \
      pa(3,7) dbg(7) \
      pc(7,5) dbg(5) \
      pd(5,6) dbg(6) \
      pd(5,1) dbg(1) \
      pm(6,6) dbg(6) \
      pm(1,1) dbg(1) \
      pm(6,0) dbg(0) \
      pa(1,0) dbg(0)
#undef lpnrm2_2
#define lpnrm2_2(a_) \
      f(nta,SS(a_,MM((SS(0,CL)),RS4)),ax) \
      plq(SS(a_,MM(0,RS4)),ax,1) dbg(1) \
      pan(4,1) dbg(1) \
      pc(5,6) dbg(6) \
      pc(5,7) dbg(7) \
      pax(1,5) dbg(5) \
      prp(5,2) dbg(2) \
      px(3) \
      pcm(0,2,3)dbg(3) \
      pan(3,7) dbg(7) \
      pann(5,3) dbg(3) \
      pa(3,7) dbg(7) \
      pc(7,5) dbg(5) \
      pd(5,6) dbg(6) \
      pd(5,1) dbg(1) \
      pm(6,6) dbg(6) \
      pm(1,1) dbg(1) \
      pm(6,0) dbg(0) \
      pa(1,0) dbg(0)
#undef dpnrm2_2
#define dpnrm2_2(a_) \
      plq(SS(a_,RS4),ax,1) dbg(1) \
      pan(4,1) dbg(1) \
      pc(5,6) dbg(6) \
      pc(5,7) dbg(7) \
      pax(1,5) dbg(5) \
      prp(5,2) dbg(2) \
      px(3) \
      pcm(0,2,3)dbg(3) \
      pan(3,7) dbg(7) \
      pann(5,3) dbg(3) \
      pa(3,7) dbg(7) \
      pc(7,5) dbg(5) \
      pd(5,6) dbg(6) \
      pd(5,1) dbg(1) \
      pm(6,6) dbg(6) \
      pm(1,1) dbg(1) \
      pm(6,0) dbg(0) \
      pa(1,0) dbg(0) 
#undef plnrm2_2
#define plnrm2_2 8


#undef p1_4_nrm2_3
#define p1_4_nrm2_3(a_) \
      pls(a_,ax,1) dbg(1) \
      pcs(5,6) dbg(6) \
      pan(4,1) dbg(1) \
      paxs(1,5) dbg(5) \
      pdsr(5,6) dbg(6) \
      pdsr(5,1) dbg(1) \
      pmsr(6,6) dbg(6) \
      pmsr(1,1) dbg(1) \
      pm(6,0) dbg(0) \
      pasr(1,0) dbg(0)
#undef p1_2_nrm2_3
#define p1_2_nrm2_3(a_) \
      px(1) pld(a_,ax,1) dbg(1) \
      pc(5,6) dbg(6) \
      pan(4,1) dbg(1) \
      pax(1,5) dbg(5) \
      pd(5,6) dbg(6) \
      pd(5,1) dbg(1) \
      pm(6,6) dbg(6) \
      pm(1,1) dbg(1) \
      pm(6,0) dbg(0) \
      pa(1,0) dbg(0)
#undef p1_nrm2_3
#define p1_nrm2_3(a_) \
      plq(a_,ax,1) dbg(1) \
      pc(5,6) dbg(6) \
      pan(4,1) dbg(1) \
      pax(1,5) dbg(5) \
      pd(5,6) dbg(6) \
      pd(5,1) dbg(1) \
      pm(6,6) dbg(6) \
      pm(1,1) dbg(1) \
      pm(6,0) dbg(0) \
      pa(1,0) dbg(0)
#define p2_nrm2_3(a_) \
      plq(SS(a_,RS4),ax,1) dbg(1) \
      pc(5,6) dbg(6) \
      pan(4,1) dbg(1) \
      pax(1,5) dbg(5) \
      pd(5,6) dbg(6) \
      pd(5,1) dbg(1) \
      pm(6,6) dbg(6) \
      pm(1,1) dbg(1) \
      pm(6,0) dbg(0) \
      pa(1,0) dbg(0) \
      f(nta,SS(a_,MM((SS(2,CL)),RS4)),ax) \
      plq(SS(a_,MM(2,RS4)),ax,1) dbg(1) \
      pc(5,6) dbg(6) \
      pan(4,1) dbg(1) \
      pax(1,5) dbg(5) \
      pd(5,6) dbg(6) \
      pd(5,1) dbg(1) \
      pm(6,6) dbg(6) \
      pm(1,1) dbg(1) \
      pm(6,0) dbg(0) \
      pa(1,0) dbg(0)
#undef lpnrm2_3
#define lpnrm2_3(a_) \
      f(nta,SS(a_,MM((SS(0,CL)),RS4)),ax) \
      plq(SS(a_,MM(0,RS4)),ax,1) dbg(1) \
      pc(5,6) dbg(6) \
      pan(4,1) dbg(1) \
      pax(1,5) dbg(5) \
      pd(5,6) dbg(6) \
      pd(5,1) dbg(1) \
      pm(6,6) dbg(6) \
      pm(1,1) dbg(1) \
      pm(6,0) dbg(0) \
      pa(1,0) dbg(0)
#undef dpnrm2_3
#define dpnrm2_3(a_) \
      plq(SS(a_,RS4),ax,1) dbg(1) \
      pc(5,6) dbg(6) \
      pan(4,1) dbg(1) \
      pax(1,5) dbg(5) \
      pd(5,6) dbg(6) \
      pd(5,1) dbg(1) \
      pm(6,6) dbg(6) \
      pm(1,1) dbg(1) \
      pm(6,0) dbg(0) \
      pa(1,0) dbg(0) 
#undef plnrm2_3
#define plnrm2_3 8

#define block_nrm2_4(a_,b_) \
      Mjoin(pc,a_)(5,6) dbg(6) \
      pan(4,1) dbg(1) \
      Mjoin(pax,a_)(1,5) dbg(5) \
      Mjoin(pc,a_)(2,7) dbg(7) \
      Mjoin(pd,b_)(5,7) dbg(7) \
      Mjoin(pm,b_)(7,6) dbg(6) \
      Mjoin(pm,b_)(7,1) dbg(1) \
      Mjoin(pm,b_)(6,6) dbg(6) \
      Mjoin(pm,b_)(6,0) dbg(0) \
      Mjoin(pm,b_)(1,1) dbg(1) \
      Mjoin(pa,b_)(1,0) dbg(0)


/*  #undef p1_4_nrm2_4 */
/*  #define p1_4_nrm2_4(a_) \ */
/*        pls(a_,ax,1) dbg(1) \ */
/*        pcs(5,6) dbg(6) \ */
/*        pan(4,1) dbg(1) \ */
/*        paxs(1,5) dbg(5) \ */
/*        pcs(2,7) dbg(7) \ */
/*        pdsr(5,7) dbg(7) \ */
/*        pmsr(7,6) dbg(6) \ */
/*        pmsr(7,1) dbg(1) \ */
/*        pmsr(6,6) dbg(6) \ */
/*        pmsr(6,0) dbg(0) \ */
/*        pmsr(1,1) dbg(1) \ */
/*        pasr(1,0) dbg(0) */
#undef p1_4_nrm2_4
#define p1_4_nrm2_4(a_) \
      pls(a_,ax,1) dbg(1) \
      block_nrm2_4(s,sr)
#undef p1_2_nrm2_4
#define p1_2_nrm2_4(a_) \
      px(1) pld(a_,ax,1) dbg(1) \
      block_nrm2_4(,)
#undef p1_nrm2_4
#define p1_nrm2_4(a_) \
      plq(a_,ax,1) dbg(1) \
      block_nrm2_4(,)
#define p2_nrm2_4(a_) \
      plq(SS(a_,RS4),ax,1) dbg(1) \
      block_nrm2_4(,) \
      plq(SS(a_,MM(2,RS4)),ax,1) dbg(1) \
      f(nta,SS(a_,MM((SS(2,CL)),RS4)),ax) \
      block_nrm2_4(,)
#undef lpnrm2_4
#define lpnrm2_4(a_) \
      plq(SS(a_,MM(0,RS4)),ax,1) dbg(1) \
      f(nta,SS(a_,MM((SS(0,CL)),RS4)),ax) \
      block_nrm2_4(,)
#undef dpnrm2_4
#define dpnrm2_4(a_) \
      plq(SS(a_,RS4),ax,1) dbg(1) \
      block_nrm2_4(,)
#undef plnrm2_4
#define plnrm2_4 8


#undef p1_4_1x1_1
#define p1_4_1x1_1(a_) \
      pls(a_,ax,1) \
      pls(a_,bx,0) \
      pm(0,1) \
      pa(1,6) 
#undef p1_2_1x1_1
#define p1_2_1x1_1(a_) \
      pld(a_,ax,1) \
      pld(a_,bx,0) \
      pm(0,1) \
      pa(1,6) 
#undef p1_1x1_1
#define p1_1x1_1(a_) \
      plq(a_,ax,1) \
      plq(a_,bx,0) \
      pm(0,1) \
      pa(0,6) 
#undef p2_1x1_1
#define p2_1x1_1(a_) \
      plq(a_,ax,1) \
      plq(a_,bx,0) \
      plq(SS(a_,RS4),ax,2) \
      plq(SS(a_,RS4),bx,3) \
      pm(0,1) \
      pm(2,3) \
      pa(1,6) \
      pa(3,6) 
#undef p4_1x1_1
#define p4_1x1_1(a_) \
      f(nta,SS(a_,MM(4,RS4)),ax) \
      plq(SS(a_,MM(2,RS4)),ax,1) \
      plq(SS(a_,MM(2,RS4)),ax,1) \
      pm(0,3) \
      puq(7,a_,ax) \
      plq(SS(a_,MM(3,RS4)),ax,2) \
      pm(0,1) \
      puq(3,SS(a_,RS4),ax) \
      f(nta,SS(a_,MM(6,RS4)),ax) \
      plq(SS(a_,MM(4,RS4)),ax,7) \
      pm(0,2) \
      puq(1,SS(a_,MM(2,RS4)),ax) \
      plq(SS(a_,MM(5,RS4)),ax,3) \
      pm(0,7) \
      puq(2,SS(a_,MM(3,RS4)),ax) 
#undef lp1x1_1
#define lp1x1_1(a_) \
      plq(a_,ax,7) \
      plq(SS(a_,RS4),ax,3) \
      pm(0,7) 
#undef dp1x1_1
#define dp1x1_1(a_) \
      plq(SS(,a_,MM(2,RS4)),ax,1) \
      pm(0,3) \
      puq(7,a_,ax) \
      plq(SS(a_,MM(3,RS4)),ax,2) \
      pm(0,1) \
      puq(3,SS(a_,RS4),ax) \
      pm(0,2) \
      puq(1,SS(a_,MM(2,RS4)),ax) \
      puq(2,SS(a_,MM(3,RS4)),ax) 
#undef pl1x1_1
#define pl1x1_1 RS4


#undef p1_4_0x1_asum_1
#define p1_4_0x1_asum_1(a_) \
      pls(a_,ax,1) \
      pan(4,1) \
      pasr(1,0) 
#undef p1_2_0x1_asum_1
#define p1_2_0x1_asum_1(a_) \
      px(1) \
      pld(a_,ax,1) \
      pan(4,1) \
      pa(1,0) 
#undef p1_0x1_asum_1
#define p1_0x1_asum_1(a_) \
      plq(a_,ax,1) \
      pan(4,1) \
      pa(1,0) 
#undef p2_0x1_asum_1
#define p2_0x1_asum_1(a_) \
      plq(a_,ax,1) \
      plq(SS(a_,RS4),ax,2) \
      pan(4,1) \
      pan(4,2) \
      pa(1,0) \
      pa(2,0) 
#undef p4_0x1_asum_1
#define p4_0x1_asum_1(a_) \
      f(nta,SS(a_,MM((SS(2,CL)),RS4)),ax) \
      plq(SS(a_,MM(2,RS4)),ax,1) \
      pan(4,3) \
      pa(7,0) \
      plq(SS(a_,MM(3,RS4)),ax,2) \
      pan(4,1) \
      pa(3,0) \
      f(nta,SS(a_,MM((SS(4,CL)),RS4)),ax) \
      plq(SS(a_,MM(4,RS4)),ax,7) \
      pan(4,2) \
      pa(1,0) \
      plq(SS(a_,MM(5,RS4)),ax,3) \
      pan(4,7) \
      pa(2,0) 
#undef lp0x1_asum_1
#define lp0x1_asum_1(a_) \
      plq(a_,ax,7) \
      plq(SS(a_,MM(1,RS4)),ax,3) \
      pan(4,7) 
#undef dp0x1_asum_1
#define dp0x1_asum_1(a_) \
      plq(SS(a_,MM(2,RS4)),ax,1) \
      pan(4,3) \
      pa(7,0) \
      plq(SS(a_,MM(3,RS4)),ax,2) \
      pan(4,1) \
      pa(3,0) \
      pan(4,2) \
      pa(1,0) \
      pa(2,0) 
#undef pl0x1_asum_1
#define pl0x1_asum_1 RS4


#undef p1_4_sum_1
#define p1_4_sum_1(a_) \
      pls(a_,ax,1) \
      pasr(1,0) 
#undef p1_2_sum_1
#define p1_2_sum_1(a_) \
      px(1) \
      pld(a_,ax,1) \
      pa(1,0) 
#undef p1_sum_1
#define p1_sum_1(a_) \
      plq(a_,ax,1) \
      pa(1,0) 
#undef p2_sum_1
#define p2_sum_1(a_) \
      plq(a_,ax,1) \
      plq(SS(a_,RS4),ax,2) \
      pa(1,0) \
      pa(2,0) 
#undef p4_sum_1
#define p4_sum_1(a_) \
      f(nta,SS(a_,MM((SS(2,CL)),RS4)),ax) \
      plq(SS(a_,MM(2,RS4)),ax,1) \
      pa(7,0) \
      plq(SS(a_,MM(3,RS4)),ax,2) \
      pa(3,0) \
      f(nta,SS(a_,MM((SS(4,CL)),RS4)),ax) \
      plq(SS(a_,MM(4,RS4)),ax,7) \
      pa(1,0) \
      plq(SS(a_,MM(5,RS4)),ax,3) \
      pa(2,0) 
#undef lpsum_1
#define lpsum_1(a_) \
      plq(a_,ax,7) \
      plq(SS(a_,MM(1,RS4)),ax,3) 
#undef dpsum_1
#define dpsum_1(a_) \
      plq(SS(a_,MM(2,RS4)),ax,1) \
      pa(7,0) \
      plq(SS(a_,MM(3,RS4)),ax,2) \
      pa(3,0) \
      pa(1,0) \
      pa(2,0) 
#undef plsum_1
#define plsum_1 RS4


#undef p1_4_dot_1
#define p1_4_dot_1(a_) \
      pls(a_,ax,1) \
      pls(a_,cx,2) \
      pmsr(2,1) \
      pasr(1,0) 
#undef p1_2_dot_1
#define p1_2_dot_1(a_) \
      px(1) \
      pld(a_,ax,1) \
      px(2) \
      pld(a_,cx,2) \
      pm(2,1) \
      pa(1,0) 
#undef p1_dot_1
#define p1_dot_1(a_) \
      plq(a_,ax,1) \
      pl(a_,cx,2) \
      pm(2,1) \
      pa(1,0) 
#undef p2_dot_1
#define p2_dot_1(a_) \
      plq(SS(a_,MM(1,RS4)),ax,1) \
      pl(SS(a_,MM(1,RS4)),cx,2) \
      pm(4,3) \
      pa(3,0) \
      f(nta,SS(a_,MM((SS(2,CL)),RS4)),ax) \
      plq(SS(a_,MM(2,RS4)),ax,3) \
      f(nta,SS(a_,MM((SS(2,CL)),RS4)),cx) \
      pl(SS(a_,MM(2,RS4)),cx,4) \
      pm(2,1) \
      pa(1,0) 
#undef lpdot_1
#define lpdot_1(a_) \
      f(nta,SS(a_,MM((SS(0,CL)),RS4)),ax) \
      plq(a_,ax,3) \
      f(nta,SS(a_,MM((SS(0,CL)),RS4)),cx) \
      pl(a_,cx,4) 
#undef dpdot_1
#define dpdot_1(a_) \
      plq(SS(a_,MM(1,RS4)),ax,1) \
      pl(SS(a_,MM(1,RS4)),cx,2) \
      pm(4,3) \
      pa(3,0) \
      pm(2,1) \
      pa(1,0)
#undef pldot_1
#define pldot_1 8

#undef p1_4_dot_1c
#define p1_4_dot_1c(a_)
#undef p1_2_dot_1c
#define p1_2_dot_1c(a_) \
      px(1) \
      pld(a_,ax,1) \
      px(2) \
      pld(a_,cx,2) \
      pc(1,3) \
      ps(HSHUF,1,1) \
      ps(LSHUF,3,3) \
      pm(7,1) \
      pm(2,3) \
      pa(3,0) \
      pm(2,1) \
      pa(1,6)
#undef p1_dot_1c
#define p1_dot_1c(a_) \
      plq(a_,ax,1) \
      pl(a_,cx,2) \
      pc(1,3) \
      ps(HSHUF,1,1) \
      ps(LSHUF,3,3) \
      pm(7,1) \