X-Git-Url: https://git.auder.net/?a=blobdiff_plain;f=old_C_code%2Fwrapper.R;fp=old_C_code%2Fwrapper.R;h=ee15a7988b2a77c00179bc518871a7cd9dd9cfcb;hb=14cb6cf8266c0e1299f16a4b2352f54dbae26f44;hp=0000000000000000000000000000000000000000;hpb=5edda192f1a634ba284acdea7b6302f60df7b000;p=epclust.git

diff --git a/old_C_code/wrapper.R b/old_C_code/wrapper.R
new file mode 100644
index 0000000..ee15a79
--- /dev/null
+++ b/old_C_code/wrapper.R
@@ -0,0 +1,56 @@
+ppam_exe = function(path=".", np=parallel::detectCores(), data=NULL,
+	nbSeriesPerChunk, nbClusters, randomize=1, p_dissims=2)
+{
+	args = paste(nbSeriesPerChunk," ",nbClusters," ",randomize," ",p_dissims,sep="")
+
+	command_line = paste("mpirun -np ",np," ",path,"/ppam.exe cluster",sep="")
+
+	#if data provided (as data.frame or matrix...): binarize it, and add it as first argument
+	if (!is.null(data))
+	{
+		if (!is.character(data))
+		{
+			#assuming matrix or data.frame, WITH row names
+			#( identifiers; could be line number... e.g. data <- cbind(1:nrow(data),data) )
+			write.table(data, "/tmp/data_csv", sep=",", row.names=FALSE, col.names=FALSE)
+			system(paste(path,"/ppam.exe serialize /tmp/data_csv /tmp/data_bin 0 0",sep=""))
+		} else
+		{
+			system(paste(path,"/ppam.exe serialize ",data," /tmp/data_bin 0 0",sep=""))
+		}
+		command_line = paste(command_line," /tmp/data_bin",sep="")
+	}
+
+	command_line = paste(command_line," ",args,sep="")
+	system(command_line)
+}
+
+#NOTE: identifiers in first column
+getMedoids = function(path=".", xmlResult = "ppamResult.xml",
+	finalSeries = "ppamFinalSeries.bin")
+{
+	system(paste(path,"/ppam.exe deserialize ",finalSeries," ppamFinalSeries.csv -1",sep=""))
+	curves = read.table("ppamFinalSeries.csv", sep=",")
+	library(XML)
+	ranks = as.integer( xmlToList( xmlParse(xmlResult) )$ranks )
+	return ( curves[ranks,] ) # == medoids
+}
+
+getDistor = function(path=".", xmlResult = "ppamResult.xml",
+	finalSeries = "ppamFinalSeries.bin")
+{
+	system(paste(path,"/ppam.exe classif ",finalSeries," ",xmlResult,sep=""))
+}
+
+serialize = function(path=".", csvSeries, binSeries, byCols=0, nbSeries=0)
+{
+	system(paste(path,"/ppam.exe serialize ",csvSeries," ",binSeries," ",byCols," ",nbSeries,
+		sep=""))
+}
+
+deserialize = function(path=".", binSeries, csvSeries, ranks="-1", return=TRUE)
+{
+	system(paste(path,"/ppam.exe deserialize ",binSeries," ",csvSeries," ",ranks,sep=""))
+	if (return)
+		return ( read.table(csvSeries, sep=",") )
+}