Lines Matching refs:u1
14 subroutine masssetupf(ctx,q,u1,u2,u3,u4,u5,u6,u7,
18 C Input: u1,u2,u3,q Output: v1,v2,ierr
21 real*8 u1(3*q)
63 $ *dsqrt(u1(i+q*0)*u1(i+q*0)
64 $ +u1(i+q*1)*u1(i+q*1)
65 $ +u1(i+q*2)*u1(i+q*2))
71 subroutine massf(ctx,q,u1,u2,u3,u4,u5,u6,u7,
75 C Input: u1,u2,q Output: v1,ierr
78 real*8 u1(q)
84 v1(i)=u2(i)*u1(i)
90 subroutine diffsetupf(ctx,q,u1,u2,u3,u4,u5,u6,u7,
94 C Input: u1,u2,u3,q Output: v1,v2,ierr
97 real*8 u1(3*q)
160 $ *dsin(pi*(c(1)+k(1)*u1(i+0*q)))
161 $ *dsin(pi*(c(2)+k(2)*u1(i+1*q)))
162 $ *dsin(pi*(c(3)+k(3)*u1(i+2*q)))
170 subroutine diffusionf(ctx,q,u1,u2,u3,u4,u5,u6,u7,
174 C Input: u1,u2,q Output: v1,ierr
177 real*8 u1(3*q)
184 $ u2(i+0*q)*u1(i)+u2(i+5*q)*u1(i+q)+u2(i+4*q)*u1(i+2*q)
186 $ u2(i+5*q)*u1(i)+u2(i+1*q)*u1(i+q)+u2(i+3*q)*u1(i+2*q)
188 $ u2(i+4*q)*u1(i)+u2(i+3*q)*u1(i+q)+u2(i+2*q)*u1(i+2*q)
710 common /vcrns/ u1(lt),r1(lt),r2(lt),r3(lt)
881 call cggos(u1,r1,h1,h2,vmult,binvm1,tol,ceed,op_mass,
888 er1 = glrdif(u1,e1,n)
916 call cggos(u1,r2,h1,h2,vmult,binvm1,tol,ceed,op_mass,
923 er1 = glrdif(u1,e1,n)
951 call cggos(u1,r3,h1,h2,vmult,binvm1,tol,ceed,op_mass,
958 er1 = glrdif(u1,e1,n)
1019 common /vcrns/ u1(lt),r1(lt),r2(lt),r3(lt)
1203 call cggos(u1,r1,h1,h2,vmult,binvm1,tol,ceed,op_diffusion,
1210 er1 = glrdif(u1,e1,n)
1238 call cggos(u1,r2,h1,h2,vmult,binvm1,tol,ceed,op_diffusion,
1245 er1 = glrdif(u1,e1,n)
1273 call cggos(u1,r3,h1,h2,vmult,binvm1,tol,ceed,op_diffusion,
1280 er1 = glrdif(u1,e1,n)
1327 subroutine cggos(u1,r1,h1,h2,rmult,binv,tin,ceed,ceed_op,vec_p1,
1332 C Output: u1,maxit
1344 real*8 u1(lt),r1(lt),h1(lt),h2(lt)
1367 call rzero (u1,n) ! Initialize solution
1388 u1(i)=u1(i)+alph(1)* p1(i)