Commit: b49d1bf49a379bfc69962211774c1f1918cca8a3 - myprojects-hg-reborn (hg)

Commit MetaInfo

• Revisión: b49d1bf49a379bfc69962211774c1f1918cca8a3 (tree)
• Tiempo: 2023-03-05 04:43:34
• Autor: Lorenzo Isella <lorenzo.isella@gmai...>
• Commiter: Lorenzo Isella

Log Message

I merged some heads.

Cambiar Resumen

• add: R-codes/attractor1.R

Diferencia incremental

diff -r 3f40b4513eb5 -r b49d1bf49a37 R-codes/attractor1.R
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/R-codes/attractor1.R	Sat Mar 04 20:43:34 2023 +0100
@@ -0,0 +1,87 @@
+rm(list=ls())
+
+library(Rcpp)
+library(ggplot2)
+library(dplyr)
+library(purrr)
+library(scattermore)
+library(tictoc())
+library(viridis)
+library(scales)
+
+tic()
+
+opt = theme(legend.position  = "none",
+            panel.background = element_rect(fill="black"),
+            axis.ticks       = element_blank(),
+            panel.grid       = element_blank(),
+            axis.title       = element_blank(),
+            axis.text        = element_blank())
+
+
+
+## #bedhead
+cppFunction('DataFrame createTrajectory(int n, double x0, double y0, 
+            double a, double b) {
+            // create the columns
+            NumericVector x(n);
+            NumericVector y(n);
+            x[0]=x0;
+            y[0]=y0;
+            for(int i = 1; i < n; ++i) {
+            x[i] = sin(x[i-1]*y[i-1]/b)*y[i-1]+cos(a*x[i-1]-y[i-1]);
+            y[i] = x[i-1]+sin(y[i-1])/b;
+            }
+            // return a new data frame
+            return DataFrame::create(_["x"]= x, _["y"]= y);
+            }
+            ')
+
+a=1.2
+b=0.15
+
+df3=createTrajectory(4000000, 1, 1, a, b)
+
+#something new
+#color by dist from origin
+eu_dist <- function(x1, y1, x2, y2) {
+  sqrt((x1-x2)^2 + (y1-y2)^2)
+}
+
+
+df3$dist <- map2_dbl(df3$x, df3$y, ~eu_dist(.x, .y, 0, 0))
+
+pulse_pal <- colorRampPalette(c("#FE1BE1", "#A300FF", "#57F7F5", "#57F7F5", "#57F7F5", "#57F7F5"))
+pulse_pal2 <- colorRampPalette(c("#FE1BE1", "#FE1BE1", "#57F7F5", "#57F7F5", "#57F7F5", "#57F7F5", "#57F7F5"))
+pulse_pal3 <- colorRampPalette(c("#FE1BE1", "#AA1BFE", "#57F7F5", "#57F7F5", "#57F7F5"))
+pulse_pal4 <- colorRampPalette(c("#AA1BFE", "#57F7F5", "#57F7F5"))
+
+#clip outer points
+xmax <- max(df3$x)/2.5
+xmin <- min(df3$x)/2.5
+ymax <- max(df3$y)/2.5
+ymin <- min(df3$y)/2.5
+
+df3_clip <- df3 |> 
+  filter(x > xmin & x < xmax) |> 
+    filter(y > ymin & y < ymax) |>
+    as_tibble() ## |>
+    ## mutate(dist=rescale(dist))
+
+print("df3clip ready")
+    
+#plot
+ggplot(df3_clip, aes(x, y)) + 
+    ## geom_point(aes(color = dist), shape=46, alpha=.01) +
+    geom_scattermore(aes(color = dist), shape=46, alpha=.03) +
+    scale_color_viridis(option="D")+
+    ## scale_color_gradientn(colors=pulse_pal4(500)) +
+    
+  opt
+
+ggsave("pulse.png", height=5000, width=5000, unit="px"
+       )
+
+toc()
+
+sessionInfo()