From 04dbe1e0acc587c6026c9f31f15942751b579077 Mon Sep 17 00:00:00 2001
From: Andrea Manica <am315@cam.ac.uk>
Date: Tue, 27 Feb 2024 13:08:57 +0000
Subject: [PATCH] speed test file

---
 data-raw/temp/drop_nodes_speed.R | 68 ++++++++++++++++++++++++++++++++
 1 file changed, 68 insertions(+)
 create mode 100644 data-raw/temp/drop_nodes_speed.R

diff --git a/data-raw/temp/drop_nodes_speed.R b/data-raw/temp/drop_nodes_speed.R
new file mode 100644
index 0000000..c8e0e7f
--- /dev/null
+++ b/data-raw/temp/drop_nodes_speed.R
@@ -0,0 +1,68 @@
+library(geoGraph)
+library(tictoc)
+devtools::load_all()
+
+x<-worldgraph.40k
+myGraph <- getGraph(x)
+connected_sets <- RBGL::connectedComp(myGraph)
+# find the largest set
+max_set <- connected_sets[[which.max(lapply(connected_sets, length))]]
+max_set <- as.numeric(max_set)
+# all cells NOT in the largest set need to be removed
+
+# get the edges from the graph
+edgeW <- edgeWeights(myGraph)
+edgeL <- edgeL(myGraph)
+
+tic()
+# We create a new list and then copy over only the edges for which we have a node
+# from the largest set
+newEdgeL <- list()
+for (i in 1:length(edgeL)) {
+  newEdgeL[[i]] <- list()
+  # if the source is in the set, we keep its edges but remove any destination not in the set
+  if (i %in% max_set){
+    newEdgeL[[i]]$edges <- edgeL[[i]]$edges[edgeL[[i]]$edges %in% max_set]
+    newEdgeL[[i]]$weights <- edgeW[[i]][edgeL[[i]]$edges %in% max_set]
+  } else { #we remove this edge
+    newEdgeL[[i]]$edges <- numeric(0)
+    newEdgeL[[i]]$weights <- numeric(0)
+  }
+}
+toc()
+# 6.974 secs
+
+## use an apply function to loop over all elements
+max_set
+keep_selected_nodes <-function(i, edgeL, max_set){
+  this_edge <- list()
+  # if the source is in the set, we keep its edges but remove any destination not in the set
+  if (i %in% max_set){
+    this_edge$edges <- edgeL[[i]]$edges[edgeL[[i]]$edges %in% max_set]
+    this_edge$weights <- edgeW[[i]][edgeL[[i]]$edges %in% max_set]
+  } else { #we remove this edge
+    this_edge$edges <- numeric(0)
+    this_edge$weights <- numeric(0)
+  }
+  this_edge
+}
+tic()
+newEdgeL2<-lapply(1:length(edgeL),FUN=keep_selected_nodes,edgeL=edgeL,max_set=max_set)
+toc()
+identical(newEdgeL,newEdgeL2)
+
+
+## break down in two steps # not working yet????
+tic()
+newEdgeL3 <-lapply(1:length(edgeL),FUN=function(i) {list(edges=numeric(0), weights=numeric(0))})
+for (i in max_set) {
+  newEdgeL3[[i]] <- list()
+  # if the source is in the set, we keep its edges but remove any destination not in the set
+    newEdgeL3[[i]]$edges <- edgeL[[i]]$edges[edgeL[[i]]$edges %in% max_set]
+    newEdgeL3[[i]]$weights <- edgeW[[i]][edgeL[[i]]$edges %in% max_set]
+}
+toc()
+identical(newEdgeL,newEdgeL3)
+
+
+##