tinytemp

commit f38ebf41170d1c1b78cef3ce7f59e0ece394978a
Author: Martin Brien <martin@mbrien.net>
Date:   Tue, 30 Apr 2019 13:22:30 +0000

uiser

Diffstat:
tinytemp.r
 | 
82
+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

1 file changed, 82 insertions(+), 0 deletions(-)
diff --git a/tinytemp.r b/tinytemp.r
@@ -0,0 +1,82 @@
+#d0T0 =-(1A*1h01*(0T0-1T1)                    )/1K0
+#d1T1 =-(1A*1h12*(1T1-1T2) - 1A*1h01*(1T1-0T0))/1K1
+#d1T2 =-(1A*1h23*(1T2-1T3) + 1A*1h12*(1T2-1T1))/1K2
+
+K = 10000
+ref = c('01','02','03','04')
+A = c(100, 100, 100,100)
+U = c(5.0,5.0,5.0,5.0)
+K = c(1000,1000,1000,1000)
+
+df <- data.frame(ref, A, U, K)
+
+print(df)
+
+
+y = c(rep(20, nrow(df)*5+2))
+y = rbind(y)
+vars = vector()
+
+vars = append(vars, "NTN")
+
+# initial node temps
+for(j in 1:5) {
+    for(i in 1:nrow(df)) { 
+        vars=append(vars,paste0("S",i,"T",j))
+
+    }
+}
+
+vars = append(vars, "XTX")
+
+colnames(y) = vars
+
+#print("f")
+print(y)
+ls(all.names = TRUE)
+#print(state[1:nrow(df)])
+
+
+#XTX <- read.csv("iso.csv", header=FALSE,skip = 5)
+#XTX <- approxfun(
+#	x = XTX$V1,
+#	y = XTX$V6, 
+#	method = "linear", 
+#	rule = 2)
+
+#NTN=20
+#XTX=20
+Xh=5
+Nh=5
+
+
+dNTN = (-(sum(df$A,Nh,y[,'NTN'])-sum(df$A,Nh,y[2:nrow(df)]))/K)
+
+dS1T1 = (-df$A[1]*((1/df$U[1]/6)*(y[,()]
+assign(paste0("dS",i,"T1"), (-df$A[i]*((1/df$U[i])/6)*((get("S1T1")-get("S1T2")) +(Nh*(get("S1T1"))-get("NTN"))/(df$K[1]/3))))
+
+#for(i in 1:nrow(df)) {
+#        assign(paste0("dS",i,"T2"), (-df$A[i]*((2/df$U[i])/6)*(get(paste0("S",i,"T2"))-get(paste0("S",i,"T3"))) +
+#        (df$A[i]*((1/df$U[i])/6)*(get(paste0("S",i,"T2"))-get(paste0("S",i,"T1")))/(df$K[i]/3))))
+#}
+
+#for(i in 1:nrow(df)) {
+#        assign(paste0("dS",i,"T3"), (-df$A[i]*((2/df$U[i])/6)*(get(paste0("S",i,"T3"))-get(paste0("S",i,"T4"))) +
+#        (df$A[i]*((2/df$U[i])/6)*(get(paste0("S",i,"T3"))-get(paste0("S",i,"T2")))/(df$K[i]/3))))
+#}
+
+#for(i in 1:nrow(df)) {
+#        assign(paste0("dS",i,"T4"), (-df$A[i]*((2/df$U[i])/6)*(get(paste0("S",i,"T4"))-get(paste0("S",i,"T5"))) +
+#        (df$A[i]*((1/df$U[i])/6)*(get(paste0("S",i,"T4"))-get(paste0("S",i,"T3")))/(df$K[i]/3))))
+#}
+
+#for(i in 1:nrow(df)) {
+#        assign(paste0("dS",i,"T5"), (-df$A[i]*(Xh*(get(paste0("S",i,"T5"))-get(paste0("XTX"))) +
+#        (df$A[i]*((1/df$U[i])/6)*(get(paste0("S",i,"T5"))-get(paste0("S",i,"T4")))))))
+#}
+
+
+
+
+#ls(all.names = TRUE)
+#print(dS1T2)