Research Article
Efficient Alternative for Construction of the Linear System Stemming from Numerical Solution of Heat Transfer Problems via FEM
Code 3
Code used for the construction of KNN for CASE 3.
Mirow = 0; Mjcol = 0; KNN = 0; k = 1; ka = 1 | do i = 1,NNost | ki = 1; iel = 1 | do j = 1,Nelem | do ii = 1,NNos | if(KCONEC(j,ii).eq.i) then | iAUX = j | KNN(i,ki) = j | ki = ki + 1 | iel = iel + 1 | go to (1) | end if | if(ii.eq.NNos) then | go to (3) | end if | end do | (1) continue | if(iel.eq.2) then | do ii = 1,NNos | Mirow(k) = i | Mjcol(k) = KCONEC(iAUX,ii) | k = k + 1 | end do | else | do ii = 1,NNos | do kk = ka,k | if(Mirow(kk).eq.i.and.Mjcol(kk).eq.KCONEC(iAUX,ii)) then | go to (2) | end if | end do | Mirow(k) = i | Mjcol(k) = KCONEC(iAUX,ii) | k = k + 1 | (2) continue | end do | end if | (3) continue | end do | ka = k - 1 | end do | k = k - 1 | continue |
|