コミット: c10ac77be5f614fec8c42094facd38606743edd8 - myprojects-hg-reborn (hg)

コミットメタ情報

• リビジョン: c10ac77be5f614fec8c42094facd38606743edd8 (tree)
• 日時: 2023-03-08 19:55:42
• 作者: Lorenzo Isella <lorenzo.isella@gmai...>
• コミッター: Lorenzo Isella

ログメッセージ

Code to parallelize the generation of Clifford attractors.

変更サマリ

• add: R-codes/gen_clifford_parallel.R

差分

diff -r 1774b95e0f05 -r c10ac77be5f6 R-codes/gen_clifford_parallel.R
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/R-codes/gen_clifford_parallel.R	Wed Mar 08 11:55:42 2023 +0100
@@ -0,0 +1,131 @@
+rm(list=ls())
+
+library(Rcpp)
+
+library(ggplot2)
+library(scattermore)
+library(dplyr)
+library(furrr)
+
+
+return_cores <- function(x){
+
+return(x)
+    
+}
+
+
+
+opt = theme(legend.position  = "none",
+
+            panel.background = element_rect(fill="white"),
+
+            axis.ticks       = element_blank(),
+
+            panel.grid       = element_blank(),
+
+            axis.title       = element_blank(),
+
+            axis.text        = element_blank())
+
+cppFunction('DataFrame createTrajectory(int n, double x0, double y0,
+
+            double a1, double a2, double a3, double a4) {
+
+            // 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(a1*y[i-1])+a3*cos(a1*x[i-1]);
+
+            y[i] = sin(a2*x[i-1])+a4*cos(a2*y[i-1]);
+
+            }
+
+            // return a new data frame
+
+            return DataFrame::create(_["x"]= x, _["y"]= y);
+
+            }
+
+            ')
+
+
+#########################################################################
+#########################################################################
+#########################################################################
+
+n_cores <- 8
+
+plan(multicore(workers=return_cores(n_cores)))
+
+
+set.seed(1253)
+
+n_trial <- 50
+
+npoints <- 10000000
+
+
+trial_seq <- 1:n_trial
+
+
+## a= -1. ## -1.24458046630025
+
+## b= -2.  ## -1.25191834103316
+
+## c= -1.  ## -1.81590817030519
+
+## d= -1.1 ## -1.90866735205054
+
+
+a1_seq <- runif( n_trial,-3, 3)
+a2_seq <-  runif( n_trial,-3, 3)
+a3_seq <- runif( n_trial,-3, 3)
+a4_seq <- runif( n_trial,-3, 3)
+
+
+df_seq <- tibble(a1=a1_seq, a2=a2_seq, a3=a3_seq,
+             a4=a4_seq)
+
+saveRDS(df_seq, "sequence_data.RDS")
+
+
+
+
+
+
+clifford_gen <- function( i){
+
+a1 <- df_seq$a1[i]
+a2 <- df_seq$a2[i]
+a3 <- df_seq$a3[i]
+a4 <- df_seq$a4[i]
+
+# Number of points to draw
+
+df=createTrajectory(npoints, 0, 0, a1, a2, a3, a4)
+
+fname <- paste("clifford_new_", i, ".png", sep="")
+
+
+gpl <- ggplot(df, aes(x, y)) +
+geom_point(color="black", shape=46, alpha=.01) + opt
+
+ggsave(plot=gpl, filename=fname, dpi=300)
+
+}
+
+future_map( trial_seq, clifford_gen )
+
+
+
+print("So far so good")