Preliminaries

Go back to the Contents page.

---

Press Show to reveal the code chunks.

---

# Create a clipboard button on the rendered HTML page
source(here::here("clipboard.R")); clipboard

# Set seed for reproducibility
set.seed(1982) 
# Set global options for all code chunks
knitr::opts_chunk$set(
  # Disable messages printed by R code chunks
  message = FALSE,    
  # Disable warnings printed by R code chunks
  warning = FALSE,    
  # Show R code within code chunks in output
  echo = TRUE,        
  # Include both R code and its results in output
  include = TRUE,     
  # Evaluate R code chunks
  eval = TRUE,       
  # Enable caching of R code chunks for faster rendering
  cache = FALSE,      
  # Align figures in the center of the output
  fig.align = "center",
  # Enable retina display for high-resolution figures
  retina = 2,
  # Show errors in the output instead of stopping rendering
  error = TRUE,
  # Do not collapse code and output into a single block
  collapse = FALSE
)
# Start the figure counter
fig_count <- 0
# Define the captioner function
captioner <- function(caption) {
  fig_count <<- fig_count + 1
  paste0("Figure ", fig_count, ": ", caption)
}

library(MetricGraph)
library(ggplot2)
library(reshape2)
library(plotly)

library(fmesher)

capture.output(
  knitr::purl(here::here("functionality.Rmd"), output = here::here("functionality.R")),
  file = here::here("purl_log.txt")
)
source(here::here("functionality.R"))

load(here::here("data_files/fem_basis_tadpole.RData"))
p_b1

Figure 1: Time evolution of the absolute difference between the exact and approximate solution.

#plot_3d_square_mesh(gets.mesh.and.FEM.on.rectangle(1, 1, 8)$mesh)
plot_3d_square_mesh_with_hat(gets.mesh.and.FEM.on.rectangle(1, 1, 6)$mesh, hat_nodes = c(1, 20, 43), hat_height = 1, hat_color = c("red", "red", "blue"), hat_alpha = 0.2)

Figure 2: Time evolution of the absolute difference between the exact and approximate solution.

#plot_3d_mesh_edges_faces(fm_rcdt_2d(globe = 3))
plot_3d_mesh_edges_faces_with_hat(
  fm_rcdt_2d(globe = 3),
  hat_nodes = c(1, 17, 37),
  hat_color = c("blue", "blue", "blue"),
  hat_height = 1,
  hat_alpha = 0.2
)

Figure 3: Time evolution of the absolute difference between the exact and approximate solution.

1 For the error

load(here::here("data_files/exp1_error_tadpole.RData"))
error_tadpole

Figure 4: Time evolution of the absolute difference between the exact and approximate solution.

load(here::here("data_files/exp_1_error_rectangle.RData"))
error_rectangle

Figure 5: Time evolution of the absolute difference between the exact and approximate solution.

load(here::here("data_files/exp_1_error_sphere.RData"))
error_sphere

Figure 6: Time evolution of the absolute difference between the exact and approximate solution.

2 Functions values

3 For the error

load(here::here("data_files/exp1_f_tadpole.RData"))
f_tadpole

Figure 7: Time evolution of the absolute difference between the exact and approximate solution.

load(here::here("data_files/exp_1_f_rect.RData"))
f_rect

Figure 8: Time evolution of the absolute difference between the exact and approximate solution.

load(here::here("data_files/exp_1_f_sphere.RData"))
f_sphere

Figure 9: Time evolution of the absolute difference between the exact and approximate solution.