| /phasta/svLS/ ! |
| H A D | CGRAD.f | 42 2 ls, rowPtr, colPtr, D, G, L, R)
53 REAL*8, INTENT(INOUT) :: R(nNo) local
66 errO = NORMS(mynNo, commu, R)
74 P = R
99 R = R - alpha*SP
100 err = NORMS(mynNo, commu, R)
102 P = R + err/errO*P
104 R = X
118 2 rowPtr, colPtr, K, R)
128 REAL*8, INTENT(INOUT) :: R(dof,nNo) local
[all …]
|
| H A D | SOLVE.f | 55 REAL*8, ALLOCATABLE :: R(:,:), W(:,:) local
84 ALLOCATE(R(dof,nNo), W(dof,nNo))
86 R(:,lhs%map(a)) = Ri(:,a)
89 CALL COMMUV(dof, nNo, lhs%commu, lhs%cS, R)
91 2 lhs%face, lhs%rowPtr, lhs%colPtr, lhs%diagPtr, Val, R, W)
97 3 Val, R)
100 2 lhs%cS, lhs%face, ls%RI, lhs%rowPtr, lhs%colPtr, Val, R)
104 2 ls%RI, lhs%rowPtr, lhs%colPtr, Val, R)
107 2 ls%RI, lhs%rowPtr, lhs%colPtr, Val, R)
113 R = R*W
[all …]
|
| H A D | INCOMMU.f | 41 SUBROUTINE COMMUV(dof, nNo, commu, cS, R) argument
48 REAL*8, INTENT(INOUT) :: R(dof,nNo) local
74 CALL MPI_ISEND(R(:,s:e), cS(i)%n*dof, mpreal, i-1,
87 R(:,s:e) = R(:,s:e) + rTmp(:,k:k+e-s)
98 rTmp(:,k:k+e-s) = R(:,s:e)
109 CALL MPI_IRECV(R(:,s:e), cS(i)%n*dof, mpreal, i-1,
129 SUBROUTINE COMMUS(nNo, commu, cS, R) argument
136 REAL*8, INTENT(INOUT) :: R(nNo) local
162 CALL MPI_ISEND(R(s:e), cS(i)%n, mpreal, i-1,
175 R(s:e) = R(s:e) + rTmp(k:k+e-s)
[all …]
|
| H A D | GMRES.f | 42 2 ls, rowPtr, colPtr, Val, R, X)
52 REAL*8, INTENT(IN) :: Val(dof*dof,nnz), R(dof,nNo) local
70 u(:,:,1) = R
79 u(:,:,1) = R - u(:,:,1)
164 2 ls, rowPtr, colPtr, Val, R)
175 REAL*8, INTENT(INOUT) :: R(dof,nNo) local
188 eps = NORMV(dof, mynNo, commu, R)
208 u(:,:,1) = R - u(:,:,1)
257 R = X
|
| H A D | PRECOND.f | 42 2 rowPtr, colPtr, diagPtr, Val, R, W)
51 REAL*8, INTENT(INOUT) :: Val(dof*dof,nnz), R(dof,nNo) local
156 R = W*R
|
| /phasta/phSolver/common/ ! |
| H A D | symtri.c | 43 #define R 0.059715871789770 macro
75 {Q,Q,R},{Q,R,Q},{R,Q,Q}};
|