Lines Matching full:1
21 qref=(/2.d-1,6.d-1,1.d0/3.d0,2.d-1,2.d-1,2.d-1,1.d0/3.d0,6.d-1/)
24 do i=0,q-1
25 x1 = qref(0*q+i+1)
26 x2 = qref(1*q+i+1);
28 interp(i*P+1)=2.*(x1+x2-1.)*(x1+x2-1./2.);
29 interp(i*P+2)=-4.*x1*(x1+x2-1.);
30 interp(i*P+3)=2.*x1*(x1-1./2.);
31 interp(i*P+4)=-4.*x2*(x1+x2-1.);
33 interp(i*P+6)=2.*x2*(x2-1./2.);
35 grad((i+0)*P+1)=2.*(1.*(x1+x2-1./2.)+(x1+x2-1.)*1.);
36 grad((i+Q)*P+1)=2.*(1.*(x1+x2-1./2.)+(x1+x2-1.)*1.);
37 grad((i+0)*P+2)=-4.*(1.*(x1+x2-1.)+x1*1.);
38 grad((i+Q)*P+2)=-4.*(x1*1.);
39 grad((i+0)*P+3)=2.*(1.*(x1-1./2.)+x1*1.);
41 grad((i+0)*P+4)=-4.*(x2*1.);
42 grad((i+Q)*P+4)=-4.*(1.*(x1+x2-1.)+x2*1.);
43 grad((i+0)*P+5)=4.*(1.*x2);
44 grad((i+Q)*P+5)=4.*(x1*1.);
46 grad((i+Q)*P+6)=2.*(1.*(x2-1./2.)+x2*1.);