MovementMatrix_Fn: Calculate movement matrices on triangulated mesh

In James-Thorson/movement_tools: Tools for movement operators in GMRF models

Usage

MovementMatrix_Fn(mesh, TriList)

Arguments

mesh
TriList

Examples

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (mesh, TriList) 
{
    M1 = M2 = M3 = M4 = matrix(0, ncol = n_r, nrow = n_r)
    for (i in 1:n_r) {
        for (j in 1:3) {
            if (!is.na(mesh$graph$tt[i, j])) {
                v2 <- list(TriList$E0, TriList$E1, TriList$E2)[[j]][i, 
                  ]
                M1[mesh$graph$tt[i, j], i] = M1[mesh$graph$tt[i, 
                  j], i] + posflux(c(1, 0), v2)/TriList$Tri_Area[i]
                M2[mesh$graph$tt[i, j], i] = M2[mesh$graph$tt[i, 
                  j], i] + posflux(c(0, 1), v2)/TriList$Tri_Area[i]
                M3[mesh$graph$tt[i, j], i] = M3[mesh$graph$tt[i, 
                  j], i] + posflux(c(-1, 0), v2)/TriList$Tri_Area[i]
                M4[mesh$graph$tt[i, j], i] = M4[mesh$graph$tt[i, 
                  j], i] + posflux(c(0, -1), v2)/TriList$Tri_Area[i]
            }
        }
        M1[i, i] = -1 * sum(M1[, i])
        M2[i, i] = -1 * sum(M2[, i])
        M3[i, i] = -1 * sum(M3[, i])
        M4[i, i] = -1 * sum(M4[, i])
    }
    maxout = max(abs(diag(M1 + M2 + M3 + M4)))
    Return = list(maxout = maxout, M1 = M1/maxout, M2 = M2/maxout, 
        M3 = M3/maxout, M4 = M4/maxout)
  }

James-Thorson/movement_tools documentation built on May 7, 2019, 10:19 a.m.