Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

problem with simulation from user specified location #27

Open
marysalvana opened this issue Feb 26, 2025 · 0 comments
Open

problem with simulation from user specified location #27

marysalvana opened this issue Feb 26, 2025 · 0 comments

Comments

@marysalvana
Copy link

Hi,

I get this error when trying to simulate from my own .txt file

[mls23007@cn559 codes]$ Rscript generating_spatial_from_data.R 
[starpu][check_bus_config_file] No performance model for the bus, calibrating...
[starpu][check_bus_config_file] ... done
		 ** Initialize ExaGeoStat hardware **
		 ** ExaGeoStat data generation/loading **
		 	Data is read from ./synthetic_data.txt successfully.

 *** caught segfault ***
address (nil), cause 'unknown'

Traceback:
 1: .External(list(name = "InternalFunction_invoke", address = <pointer: 0x2cfce40>,     dll = list(name = "Rcpp", path = "/home/mls23007/R/x86_64-redhat-linux-gnu-library/4.3/Rcpp/libs/Rcpp.so",         dynamicLookup = TRUE, handle = <pointer: 0x39c0a50>,         info = <pointer: 0x2d71e20>), numParameters = -1L), <pointer: 0x3d55390>,     kernel, initial_theta, distance_matrix, problem_size, seed,     dts, lts, dimension, log_path, data_path, observations_file,     recovery_file)
 2: simulate_data(kernel = kernel, initial_theta = initial_theta,     problem_size = 10000, dts = dts, dimension = dimension, data_path = data_path)
An irrecoverable exception occurred. R is aborting now ...

Here is my code:

library(assertthat)
library("ExaGeoStatCPP")

# Variables
computation <- "exact"
ncores <- 63
ngpus <- 0
p <- 1
q <- 1

hardware <- new(Hardware, computation, ncores, ngpus, p, q)

kernel <- "univariate_matern_stationary"
initial_theta <- c(1,0.1,0.5)
dts <- 360   #360 -560 -760 -960
dimension = "2D"

data_path <- "./synthetic_data.txt"

# Time the simulation with small size

start_time <- Sys.time()
exageostat_data <- simulate_data(kernel=kernel, initial_theta=initial_theta, problem_size=10000, dts=dts, dimension=dimension, data_path=data_path)
end_time <- Sys.time()

time_taken <- end_time - start_time
seconds_taken <- as.numeric(time_taken, units="secs")
print(paste("Time taken for 10000 locations:", seconds_taken, "seconds, using ", ncores, " cores."))

I am attaching synthetic_data.txt below.

Thank you!

synthetic_data.txt

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant