first commit

Author: stla

.Rbuildignore

@@ -0,0 +1,3 @@
+^PolyhedralCubature\.Rproj$
+^\.Rproj\.user$
+^inst/essais$

DESCRIPTION

@@ -0,0 +1,16 @@
+Package: PolyhedralCubature
+Title: Integration Over Convex Polyhedra
+Version: 0.0.0.9000
+Authors@R: 
+    person("Stéphane", "Laurent", , "laurent_step@outlook.fr", role = c("aut", "cre"))
+Description: What the package does (one paragraph).
+License: GPL-3
+URL: https://github.com/stla/PolyhedralCubature
+BugReports: https://github.com/stla/PolyhedralCubature/issues
+Imports: 
+    gmp,
+    rcdd,
+    SimplicialCubature,
+    tessellation
+Encoding: UTF-8
+RoxygenNote: 7.2.3

NAMESPACE

@@ -0,0 +1,2 @@
+# Generated by roxygen2: do not edit by hand
+

PolyhedralCubature.Rproj

@@ -0,0 +1,22 @@
+Version: 1.0
+
+RestoreWorkspace: No
+SaveWorkspace: No
+AlwaysSaveHistory: Default
+
+EnableCodeIndexing: Yes
+UseSpacesForTab: Yes
+NumSpacesForTab: 2
+Encoding: UTF-8
+
+RnwWeave: Sweave
+LaTeX: pdfLaTeX
+
+AutoAppendNewline: Yes
+StripTrailingWhitespace: Yes
+LineEndingConversion: Posix
+
+BuildType: Package
+PackageUseDevtools: Yes
+PackageInstallArgs: --no-multiarch --with-keep.source
+PackageRoxygenize: rd,collate,namespace

inst/essais/essai01.R

@@ -0,0 +1,52 @@
+library(rcdd)
+library(tessellation)
+library(SimplicialCubature)
+
+#' @importFrom rcdd makeH validcdd scdd
+#' @importFrom tessellation delaunay getDelaunaySimplicies
+#' @importFrom SimplicialCubature adaptIntegrateSimplex
+
+A <- rbind(
+  c(-1, 0, 0), # -x
+  c( 1, 0, 0), # x
+  c( 0,-1, 0), # -y
+  c( 1, 1, 0), # x+y
+  c( 0, 0,-1), # -z
+  c( 1, 1, 1)  # x+y+z
+)
+b <- c(5, 4, 5, 3, 10, 6)
+
+H <- makeH(A, b)
+. <- validcdd(H)
+
+V <- scdd(H)[["output"]]
+if(any(V[, 1L] != 0) || any(V[, 2L] != 1)) {
+  stop("Invalid input.")
+}
+vertices <- V[, -c(1L, 2L)]
+
+dlny <- delaunay(vertices)
+simplices <- getDelaunaySimplicies(dlny)
+nsimplices <- length(simplices)
+
+# union of simplicies
+d <- ncol(A)
+S <- array(NA_real_, dim=c(d, d+1L, nsimplices))
+for(i in seq_len(nsimplices)) {
+  S[, , i] <- t(simplices[[i]])
+}
+
+# integration
+f <- function(x, y, z) {
+  x + y*z
+}
+stopifnot(length(formals(f)) == d)
+g <- function(v) NULL
+bdy <- sprintf(
+  "f(%s)", paste0(sprintf("v[%d]", 1:d), collapse = ",")
+)
+body(g) <- parse(text = bdy)
+
+adaptIntegrateSimplex(g, S)
+
+