+ p->c2[0]=p->c[0]*p->c[0];
+ p->c2[1]=p->c[1]*p->c[1];
+ p->c2_mixed=p->c_mixed*p->c_mixed;
+ p->d2[0]=p->d[0]*p->d[0];
+ p->d2[1]=p->d[1]*p->d[1];
+ p->d2_mixed=p->d_mixed*p->d_mixed;
+ p->c2d2[0]=p->c2[0]/p->d2[0];
+ p->c2d2[1]=p->c2[1]/p->d2[1];
+ p->c2d2_m=p->c2_mixed/p->d2_mixed;