camm_strat1.h

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

      pm(2,3) \
      pa(3,0) \
      pm(2,1) \
      pa(1,6)
#undef p2_dot_1c
#define p2_dot_1c(a_) \
      f(nta,SS(a_,MM((SS(2,CL)),RS4)),ax) \
      plq(SS(a_,MM(1,RS4)),ax,1) \
      pl(SS(a_,MM(1,RS4)),cx,2) \
      pc(3,5) \
      ps(HSHUF,3,3) \
      ps(LSHUF,5,5) \
      pm(7,3) \
      pm(4,5) \
      pa(5,0) \
      pm(4,3) \
      pa(3,6) \
      f(nta,SS(a_,MM((SS(2,CL)),RS4)),cx) \
      pl(SS(a_,MM(2,RS4)),cx,4) \
      plq(SS(a_,MM(2,RS4)),ax,3) \
      pc(1,5) \
      ps(HSHUF,1,1) \
      ps(LSHUF,5,5) \
      pm(7,1) \
      pm(2,5) \
      pa(5,0) \
      pm(2,1) \
      pa(1,6)
#undef lpdot_1c
#define lpdot_1c(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_1c
#define dpdot_1c(a_) \
      plq(SS(a_,MM(1,RS4)),ax,1) \
      pl(SS(a_,MM(1,RS4)),cx,2) \
      pc(3,5) \
      ps(HSHUF,3,3) \
      ps(LSHUF,5,5) \
      pm(7,3) \
      pm(4,5) \
      pa(5,0) \
      pm(4,3) \
      pa(3,6) \
      pc(1,5) \
      ps(HSHUF,1,1) \
      ps(LSHUF,5,5) \
      pm(7,1) \
      pm(2,5) \
      pa(5,0) \
      pm(2,1) \
      pa(1,6)
#undef pldot_1c
#define pldot_1c 8

#undef p1_4_dot_2c
#define p1_4_dot_2c(a_)
#undef p1_2_dot_2c
#define p1_2_dot_2c(a_) \
      px(1) \
      pld(a_,ax,1) \
      px(2) \
      pld(a_,cx,2) \
      pc(1,3) \
      ps(CSHUF,1,1) \
      pm(2,3) \
      pa(3,0) \
      pm(2,1) \
      pa(1,6)
#undef p1_dot_2c
#define p1_dot_2c(a_) \
      plq(a_,ax,1) \
      pl(a_,cx,2) \
      pc(1,3) \
      ps(CSHUF,1,1) \
      pm(2,3) \
      pa(3,0) \
      pm(2,1) \
      pa(1,6)
#undef p2_dot_2c
#define p2_dot_2c(a_) \
      f(nta,SS(a_,MM((SS(2,CL)),RS4)),ax) \
      plq(SS(a_,MM(1,RS4)),ax,1) \
      pl(SS(a_,MM(1,RS4)),cx,2) \
      pc(3,5) \
      ps(CSHUF,3,3) \
      pm(4,5) \
      pa(5,0) \
      pm(4,3) \
      pa(3,6) \
      f(nta,SS(a_,MM((SS(2,CL)),RS4)),cx) \
      pl(SS(a_,MM(2,RS4)),cx,4) \
      plq(SS(a_,MM(2,RS4)),ax,3) \
      pc(1,5) \
      ps(CSHUF,1,1) \
      pm(2,5) \
      pa(5,0) \
      pm(2,1) \
      pa(1,6)
#undef lpdot_2c
#define lpdot_2c(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_2c
#define dpdot_2c(a_) \
      plq(SS(a_,MM(1,RS4)),ax,1) \
      pl(SS(a_,MM(1,RS4)),cx,2) \
      pc(3,5) \
      ps(CSHUF,3,3) \
      pm(4,5) \
      pa(5,0) \
      pm(4,3) \
      pa(3,6) \
      pc(1,5) \
      ps(CSHUF,1,1) \
      pm(2,5) \
      pa(5,0) \
      pm(2,1) \
      pa(1,6)
#undef pldot_2c
#define pldot_2c 8

#undef p1_4_axpby_3
#define p1_4_axpby_3(a_) \
      pls(a_,ax,0) \
      pls(a_,cx,3) \
      pmsr(5,0) \
      pmsr(6,3) \
      pasr(3,0) \
      pus(0,a_,ax)
#undef p1_2_axpby_3
#define p1_2_axpby_3(a_) \
      pld(a_,ax,0) \
      pld(a_,cx,3) \
      pm(5,0) \
      pm(6,3) \
      pa(3,0) \
      pud(0,a_,ax)
#undef p1_axpby_3
#define p1_axpby_3(a_) \
      plq(a_,ax,0) \
      pl(a_,cx,3) \
      pm(5,0) \
      pm(6,3) \
      pa(3,0) \
      punt(0,a_,ax)
#undef p2_axpby_3
#define p2_axpby_3(a_) \
      plq(a_,ax,0) \
      pl(a_,cx,3) \
      plq(SS(a_,RS4),ax,1) \
      pm(5,0) \
      pm(6,3) \
      pa(3,0) \
      pl(SS(a_,RS4),cx,3) \
      punt(0,a_,ax) \
      pm(5,1) \
      pm(6,3) \
      pa(3,1) \
      punt(1,SS(a_,RS4),ax)
#undef p4_axpby_3
#define p4_axpby_3(a_) \
      plq(SS(a_,MM(3,RS4)),ax,3) \
      pm(5,2) \
      pl(SS(a_,MM(3,RS4)),cx,7) \
      pm(6,4) \
      pa(4,2) \
      punt(0,a_,ax) \
      f(nta,SS(a_,MM((SS(2,CL)),RS4)),cx) \
      pl(SS(a_,MM(4,RS4)),cx,4) \
      pm(5,3) \
      plq(SS(a_,MM(4,RS4)),ax,0) \
      pm(6,7) \
      pa(7,3) \
      punt(1,SS(a_,RS4),ax) \
      f(nta,SS(a_,MM((SS(2,CL)),RS4)),ax) \
      plq(SS(a_,MM(5,RS4)),ax,1) \
      pm(5,0) \
      pl(SS(a_,MM(5,RS4)),cx,7) \
      pm(6,4) \
      pa(4,0) \
      punt(2,SS(a_,MM(2,RS4)),ax) \
      f(nta,SS(a_,MM((SS(4,CL)),RS4)),cx) \
      pl(SS(a_,MM(6,RS4)),cx,4) \
      pm(5,1) \
      plq(SS(a_,MM(6,RS4)),ax,2) \
      pm(6,7) \
      pa(7,1) \
      punt(3,SS(a_,MM(3,RS4)),ax) \
      f(nta,SS(a_,MM((SS(4,CL)),RS4)),ax) 
#undef lpaxpby_3
#define lpaxpby_3(a_) \
      f(nta,SS(a_,MM((SS(0,CL)),RS4)),cx) \
      pl(SS(a_,MM(0,RS4)),cx,4) \
      plq(SS(a_,MM(0,RS4)),ax,0) \
      pl(SS(a_,MM(1,RS4)),cx,7) \
      pm(5,0) \
      plq(SS(a_,MM(1,RS4)),ax,1) \
      pm(6,4) \
      pa(4,0) \
      f(nta,SS(a_,MM((SS(0,CL)),RS4)),ax) \
      plq(SS(a_,MM(2,RS4)),ax,2) \
      pm(5,1) \
      pl(SS(a_,MM(2,RS4)),cx,4) \
      pm(6,7) \
      pa(7,1)
#undef dpaxpby_3
#define dpaxpby_3(a_) \
      pl(SS(a_,MM(3,RS4)),cx,7) \
      pm(5,2) \
      plq(SS(a_,MM(3,RS4)),ax,3) \
      pm(6,4) \
      pa(4,2) \
      pm(5,3) \
      punt(0,a_,ax) \
      pm(6,7) \
      pa(7,3) \
      punt(1,SS(a_,RS4),ax) \
      punt(2,SS(a_,MM(2,RS4)),ax) \
      punt(3,SS(a_,MM(3,RS4)),ax)
#undef plaxpby_3
#define plaxpby_3 16

#undef p1_4_axpby_3c
#define p1_4_axpby_3c(a_) 
#undef p1_2_axpby_3c
#define p1_2_axpby_3c(a_) \
      pld(a_,ax,0) \
      pld(a_,cx,2) \
      pc(0,3) \
      pm(5,0) \
      ps(CSHUF,3,3) \
      pm(4,3) \
      pa(3,0) \
      pc(2,3) \
      pm(6,2) \
      pa(2,0) \
      ps(CSHUF,3,3) \
      pm(7,3) \
      pa(3,0) \
      pud(0,a_,ax)
#undef p1_axpby_3c
#define p1_axpby_3c(a_) \
      plq(a_,ax,0) \
      pl(a_,cx,2) \
      pc(0,3) \
      pm(5,0) \
      ps(CSHUF,3,3) \
      pm(4,3) \
      pa(3,0) \
      pc(2,3) \
      pm(6,2) \
      pa(2,0) \
      ps(CSHUF,3,3) \
      pm(7,3) \
      pa(3,0) \
      puq(0,a_,ax)
#undef p2_axpby_3c
#define p2_axpby_3c(a_) \
      f(nta,SS(a_,MM((SS(2,CL)),RS4)),ax) \
      plq(SS(a_,MM(1,RS4)),ax,1) \
      pl(SS(a_,MM(1,RS4)),cx,3) \
      pc(1,2) \
      pm(5,1) \
      ps(CSHUF,2,2) \
      pm(4,2) \
      pa(2,1) \
      pc(3,2) \
      pm(6,3) \
      pa(3,1) \
      ps(CSHUF,2,2) \
      pm(7,2) \
      pa(2,1) \
      puq(0,a_,ax) \
      plq(SS(a_,MM(2,RS4)),ax,0) \
      f(nta,SS(a_,MM((SS(2,CL)),RS4)),cx) \
      pl(SS(a_,MM(2,RS4)),cx,2) \
      pc(0,3) \
      pm(5,0) \
      ps(CSHUF,3,3) \
      pm(4,3) \
      pa(3,0) \
      pc(2,3) \
      pm(6,2) \
      pa(2,0) \
      ps(CSHUF,3,3) \
      pm(7,3) \
      pa(3,0) \
      puq(1,SS(a_,RS4),ax) 
#undef lpaxpby_3c
#define lpaxpby_3c(a_) \
      f(nta,SS(a_,MM((SS(0,CL)),RS4)),ax) \
      plq(SS(a_,MM(0,RS4)),ax,0) \
      f(nta,SS(a_,MM((SS(0,CL)),RS4)),cx) \
      pl(SS(a_,MM(0,RS4)),cx,2) \
      pc(0,3) \
      pm(5,0) \
      ps(CSHUF,3,3) \
      pm(4,3) \
      pa(3,0) \
      pc(2,3) \
      pm(6,2) \
      pa(2,0) \
      ps(CSHUF,3,3) \
      pm(7,3) \
      pa(3,0) 
#undef dpaxpby_3c
#define dpaxpby_3c(a_) \
      plq(SS(a_,MM(1,RS4)),ax,1) \
      pl(SS(a_,MM(1,RS4)),cx,3) \
      pc(1,2) \
      pm(5,1) \
      ps(CSHUF,2,2) \
      pm(4,2) \
      pa(2,1) \
      pc(3,2) \
      pm(6,3) \
      pa(3,1) \
      ps(CSHUF,2,2) \
      pm(7,2) \
      pa(2,1) \
      puq(0,a_,ax) \
      puq(1,SS(a_,RS4),ax) 
#undef plaxpby_3c
#define plaxpby_3c 8

#undef p1_4_axpby_2
#define p1_4_axpby_2(a_) \
      pls(a_,cx,5) \
      pls(a_,ax,0) \
      pmsr(6,5) \
      pasr(5,0) \
      pus(0,a_,ax)
#undef p1_2_axpby_2
#define p1_2_axpby_2(a_) \
      pld(a_,cx,5) \
      pld(a_,ax,0) \
      pm(6,5) \
      pa(5,0) \
      pud(0,a_,ax)
#undef p1_axpby_2
#define p1_axpby_2(a_) \
      pl(a_,cx,5) \
      plq(a_,ax,0) \
      pm(6,5) \
      pa(5,0) \
      puq(0,a_,ax)
#undef p2_axpby_2
#define p2_axpby_2(a_) \
      pl(a_,cx,5) \
      plq(a_,ax,0) \
      pl(SS(a_,RS4),cx,4) \
      pm(6,5) \
      pa(5,0) \
      plq(SS(a_,RS4),ax,1) \
      puq(0,a_,ax) \
      pm(6,4) \
      pa(4,1) \
      puq(1,SS(a_,RS4),ax) 
#undef p4_axpby_2
#define p4_axpby_2(a_) \
      plq(SS(a_,MM(3,RS4)),ax,3) \
      pl(SS(a_,MM(3,RS4)),cx,5) \
      pm(6,4) \
      pa(4,2) \
      puq(0,a_,ax) \
      f(nta,SS(a_,MM((SS(2,CL)),RS4)),cx) \
      pl(SS(a_,MM(4,RS4)),cx,4) \
      plq(SS(a_,MM(4,RS4)),ax,0) \
      pm(6,5) \
      pa(5,3) \
      puq(1,SS(a_,RS4),ax) \
      f(nta,SS(a_,MM((SS(2,CL)),RS4)),ax) \
      plq(SS(a_,MM(5,RS4)),ax,1) \
      pl(SS(a_,MM(5,RS4)),cx,5) \
      pm(6,4) \
      pa(4,0) \
      puq(2,SS(a_,MM(2,RS4)),ax) \
      f(nta,SS(a_,MM((SS(4,CL)),RS4)),cx) \
      pl(SS(a_,MM(6,RS4)),cx,4) \
      plq(SS(a_,MM(6,RS4)),ax,2) \
      pm(6,5) \
      pa(5,1) \
      puq(3,SS(a_,MM(3,RS4)),ax) \
      f(nta,SS(a_,MM((SS(4,CL)),RS4)),ax) 
#undef lpaxpby_2
#define lpaxpby_2(a_) \
      f(nta,SS(a_,MM((SS(0,CL)),RS4)),cx) \
      pl(SS(a_,MM(0,RS4)),cx,4) \
      plq(SS(a_,MM(0,RS4)),ax,0) \
      pl(SS(a_,MM(1,RS4)),cx,5) \
      plq(SS(a_,MM(1,RS4)),ax,1) \
      pm(6,4) \
      pa(4,0) \
      f(nta,SS(a_,MM((SS(0,CL)),RS4)),ax) \
      plq(SS(a_,MM(2,RS4)),ax,2) \
      pl(SS(a_,MM(2,RS4)),cx,4) \
      pm(6,5) \
      pa(5,1)
#undef dpaxpby_2
#define dpaxpby_2(a_) \
      pl(SS(a_,MM(3,RS4)),cx,5) \
      plq(SS(a_,MM(3,RS4)),ax,3) \
      pm(6,4) \
      pa(4,2) \
      puq(0,a_,ax) \
      pm(6,5) \
      pa(5,3) \
      puq(1,SS(a_,RS4),ax) \
      puq(2,SS(a_,MM(2,RS4)),ax) \
      puq(3,SS(a_,MM(3,RS4)),ax)
#undef plaxpby_2
#define plaxpby_2 16

#undef p1_4_axpby_2c
#define p1_4_axpby_2c(a_) 
#undef p1_2_axpby_2c
#define p1_2_axpby_2c(a_) \
      pld(a_,cx,5) \
      pld(a_,ax,0) \
      pc(5,1) \
      pm(6,5) \
      pa(5,0) \
      ps(CSHUF,1,1) \
      pm(7,1) \
      pa(1,0) \
      pud(0,a_,ax)
#undef p1_axpby_2c
#define p1_axpby_2c(a_) \
      pl(a_,cx,5) \
      plq(a_,ax,0) \
      pc(5,1) \
      pm(6,5) \
      pa(5,0) \
      ps(CSHUF,1,1) \
      pm(7,1) \
      pa(1,0) \
      puq(0,a_,ax)
#undef p2_axpby_2c
#define p2_axpby_2c(a_) \
      pl(a_,cx,5) \
      plq(a_,ax,0) \
      pl(SS(a_,RS4),cx,4) \
      pc(5,1) \
      pm(6,5) \
      pa(5,0) \
      ps(CSHUF,2,2) \
      pm(7,2) \
      pa(2,0) \
      plq(SS(a_,RS4),ax,1) \
      puq(0,a_,ax) \
      pc(4,3) \
      pm(6,4) \
      pa(4,1) \
      ps(CSHUF,3,3) \
      pm(7,3) \
      pa(3,1) \
      puq(1,SS(a_,RS4),ax) 
#undef p4_axpby_2c
#define p4_axpby_2c(a_) \
      plq(SS(a_,MM(3,RS4)),ax,3) \
      puq(0,a_,ax) \
      pc(4,0) \
      pm(6,4) \
      pa(4,2) \
      ps(CSHUF,0,0) \
      f(nta,SS(a_,MM((SS(2,CL)),RS4)),cx) \
      pl(SS(a_,MM(4,RS4)),cx,4) \
      pm(7,0) \
      pa(0,2) \
      plq(SS(a_,MM(4,RS4)),ax,0) \
      puq(1,SS(a_,RS4),ax) \
      pc(5,1) \
      pm(6,5) \
      pa(5,3) \
      ps(CSHUF,1,1) \
      pl(SS(a_,MM(5,RS4)),cx,5) \
      pm(7,1) \
      pa(1,3) \
      f(nta,SS(a_,MM((SS(2,CL)),RS4)),ax) \
      plq(SS(a_,MM(5,RS4)),ax,1) \
      puq(2,SS(a_,MM(2,RS4)),ax) \
      pc(4,2) \
      pm(6,4) \
      pa(4,0) \
      ps(CSHUF,2,2) \
      f(nta,SS(a_,MM((SS(4,CL)),RS4)),cx) \
      pl(SS(a_,MM(6,RS4)),cx,4) \
      pm(7,2) \
      pa(2,0) \
      plq(SS(a_,MM(6,RS4)),ax,2) \
      puq(3,SS(a_,MM(3,RS4)),ax) \
      pc(5,3) \
      pm(6,5) \
      pa(5,1) \
      ps(CSHUF,3,3) \
      pl(SS(a_,MM(7,RS4)),cx,5) \
      pm(7,3) \
      pa(3,1) \
      f(nta,SS(a_,MM((SS(4,CL)),RS4)),ax) 
#undef lpaxpby_2c
#define lpaxpby_2c(a_) \
      f(nta,SS(a_,MM((SS(0,CL)),RS4)),cx) \
      pl(SS(a_,MM(0,RS4)),cx,4) \
      plq(SS(a_,MM(0,RS4)),ax,0) \
      pl(SS(a_,MM(1,RS4)),cx,5) \
      plq(SS(a_,MM(1,RS4)),ax,1) \
      pc(4,2) \
      pm(6,4) \
      pa(4,0) \
      ps(CSHUF,2,2) \
      pl(SS(a_,MM(2,RS4)),cx,4) \
      pm(7,2) \
      pa(2,0) \
      f(nta,SS(a_,MM((SS(0,CL)),RS4)),ax) \
      plq(SS(a_,MM(2,RS4)),ax,2) \
      pc(5,3) \
      pm(6,5) \
      pa(5,1) \
      ps(CSHUF,3,3) \
      pl(SS(a_,MM(3,RS4)),cx,5) \
      pm(7,3) \
      pa(3,1)
#undef dpaxpby_2c
#define dpaxpby_2c(a_) \
      plq(SS(a_,MM(3,RS4)),ax,3) \
      puq(0,a_,ax) \
      pc(4,0) \
      pm(6,4) \
      pa(4,2) \
      ps(CSHUF,0,0) \
      puq(1,SS(a_,RS4),ax) \
      pm(7,0) \
      pa(0,2) \
      pc(5,1) \
      pm(6,5) \
      pa(5,3) \
      ps(CSHUF,1,1) \
      puq(2,SS(a_,MM(2,RS4)),ax) \
      pm(7,1) \
      pa(1,3) \
      puq(3,SS(a_,MM(3,RS4)),ax)
#undef plaxpby_2c
#define plaxpby_2c 16

#undef p1_4_axpby_1
#define p1_4_axpby_1(a_) \
      pls(a_,ax,1) \
      pls(a_,cx,2) \
      pmsr(5,1) \
      pmsr(6,2) \
      pasr(2,1) \
      pus(1,a_,ax)
#undef p1_2_axpby_1
#define p1_2_axpby_1(a_) \
      pld(a_,ax,1) \
      pld(a_,cx,2) \
      pm(5,1) \
      pm(6,2) \
      pa(2,1) \
      pud(1,a_,ax)
#undef p1_axpby_1
#define p1_axpby_1(a_) \
      plq(a_,ax,1) \
      pl(a_,cx,2) \
      pm(5,1) \
      pm(6,2) \
      pa(2,1) \
      puq(1,a_,ax)