computes floating point errors

floating_point_01.R

library(ggplot2)
library(dplyr)
library(paletteer)

# specify the range
lo <- 0
hi <- 4
by <- .01

# create data structure
val <- seq(lo, hi, by)
dat <- expand.grid(
  x = val, 
  y = val, 
  z = seq(lo, 2*hi, by)
)
dat <- dat %>% 
  mutate(s = x + y - z) %>%   # compute all sums
  filter(abs(s) < by/10) %>%  # retain ostensibly zero rows
  arrange(x, y)               # be pretty


# create plot object
pic <- ggplot(dat, aes(x, y, fill = s)) + 
  geom_raster(show.legend = FALSE) + 
  theme_void() + 
  coord_fixed() + 
  scale_fill_paletteer_c("scico::broc") + 
  scale_x_continuous(expand = c(0, 0)) +
  scale_y_continuous(expand = c(0, 0)) +
  NULL

# draw plot
ggsave(
  filename = here::here("image", "floating_point_01.png"),
  plot = pic,
  height = 5000/300,
  width = 5000/300,
  dpi = 300
)

This is so beautiful 😍 thanks for sharing it Danielle!

❤️