context("Sampling points from a sphere") test_that("Sampling sphere surface works", { library(igraph) library(digest) set.seed(42) s1 <- sample_sphere_surface(4, 100, positive=FALSE) expect_that(colSums(s1^2), equals(rep(1, 100))) s2 <- sample_sphere_surface(3, 100, radius=2, positive=FALSE) expect_that(sqrt(colSums(s2^2)), equals(rep(2, 100))) s3 <- sample_sphere_surface(2, 100, radius=1/2, positive=TRUE) expect_that(sqrt(colSums(s3^2)), equals(rep(1/2, 100))) expect_true(all(s3 >= 0)) }) test_that("Sampling sphere volume works", { library(igraph) library(digest) set.seed(42) s1 <- sample_sphere_volume(4, 10000, positive=FALSE) expect_true(all(colSums(s1^2) < 1)) s2 <- sample_sphere_volume(3, 100, radius=2, positive=FALSE) expect_true(all(sqrt(colSums(s2^2)) < 2)) s3 <- sample_sphere_volume(2, 100, radius=1/2, positive=TRUE) expect_true(all(sqrt(colSums(s3^2)) < 1/2)) expect_true(all(s3 >= 0)) })
Generated by dwww version 1.15 on Sun Jun 16 16:11:18 CEST 2024.