Mathematical Problems in Engineering

Research Article

Efficient Alternative for Construction of the Linear System Stemming from Numerical Solution of Heat Transfer Problems via FEM

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