Repository navigation
Expand file tree
/
Copy pathstochastic_geo_display_dump.py
More file actions
71 lines (56 loc) · 1.54 KB
/
Copy pathstochastic_geo_display_dump.py
File metadata and controls
71 lines (56 loc) · 1.54 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
import numpy as np
import matplotlib.pylab as plt
#import scipy.stats as stats
#import pandas as pd
import useful_things as ut
##################
# h_gap and r_vals
##################
h_gap = 0.1;
r_vals = np.arange(h_gap, 40, h_gap);
r_vals = np.insert(r_vals, 0, 0);
r_vals = np.insert(r_vals, 0, 0);
########
# figure
########
fig1 = plt.figure(1);
fig1.patch.set_facecolor('white');
plt.clf();
ax = plt.gca();
ax.set_xlim(0, 10);
"""
# 1000
num_folder = "";
case_name = "1000"; # case
num_names = ["0", "1", "2", "3", "4"];
"""
for num_name in num_names:
str_cmd = "data_array = " "scf_vals_" + case_name + "_" + num_name;
print str_cmd;
exec(str_cmd);
data_array = [x for x in data_array if x != 0];
idx_array = r_vals[:len(data_array)];
plt.plot(idx_array, np.asarray(data_array), 'r', linewidth=1);
"""
"""
# 250
case_name = "250"; # case
for num_name in num_names:
str_cmd = "data_array = " "scf_vals_" + case_name + "_" + num_name;
print str_cmd;
exec(str_cmd);
data_array = [x for x in data_array if x != 0];
idx_array = r_vals[:len(data_array)];
plt.plot(idx_array, np.asarray(data_array), 'g', linewidth=1);
"""
"""
case_name = "5"; # case
for num_name in num_names:
str_cmd = "data_array = " "scf_vals_" + case_name + "_" + num_name;
print str_cmd;
exec(str_cmd);
data_array = [x for x in data_array if x != 0];
idx_array = r_vals[:len(data_array)];
plt.plot(idx_array, np.asarray(data_array), 'b', linewidth=1);
plt.title(r"SCPs of all");
plt.show();