poissons_eq_quadrilateral: Get sparse Poisson's matrix elements for quadrilateral domain

In GJMeijer/soilmech: R functions and plots for Soil Mechanics teaching

Get sparse Poisson's matrix elements for quadrilateral domain

Description

Function generates the sparse matrix elements for the Poisson's equation in all real nodes, for a single quadrilateral domain

Function generates the sparse matrix elements for the Poisson's equation in all real nodes, for a single quadrilateral domain

Usage

poissons_eq_quadrilateral(
  nx,
  ny,
  kx,
  ky,
  x1,
  x2,
  x3,
  x4,
  y1,
  y2,
  y3,
  y4,
  i0 = 0,
  ...
)

poissons_eq_quadrilateral(
  nx,
  ny,
  kx,
  ky,
  x1,
  x2,
  x3,
  x4,
  y1,
  y2,
  y3,
  y4,
  i0 = 0,
  ...
)

Arguments

nx, ny	number of real nodes in a and b-directions
kx, ky	horizontal and vertical permeability in domain
x1, x2, x3, x4, y1, y2, y3, y4	x and y positions of the corner points of the quadrilateral. Corners are numbered in clockwise order.
i0	index offset for nodes
...	extra arguments to pass

Value

a tibble with rows ('row'), columns ('col') and values ('val') of non-zero elements in the sparse matrix

a tibble with rows ('row'), columns ('col') and values ('val') of non-zero elements in the sparse matrix

Examples

nx <- 20
ny <- 10
kx <- 1
ky <- 1
x1 <- 0
x2 <- 0
x4 <- 1
x3 <- 1
y1 <- 0
y2 <- 1
y3 <- 2
y4 <- 0
poissons_eq_quadrilateral(nx, ny, kx, ky, x1, x2, x3, x4, y1, y2, y3, y4)
nx <- 20
ny <- 10
kx <- 1
ky <- 1
x1 <- 0
x2 <- 0
x4 <- 1
x3 <- 1
y1 <- 0
y2 <- 1
y3 <- 2
y4 <- 0
poissons_eq_quadrilateral(nx, ny, kx, ky, x1, x2, x3, x4, y1, y2, y3, y4)

GJMeijer/soilmech documentation built on May 22, 2022, 10:39 a.m.