with open('./wu_et_al_data.csv', 'r') as f:
reader = csv.reader(f)
data = np.asarray(list(reader), dtype=np.float64)
####################
# plotting results #
####################
if not os.path.isfile(
'./batch_run.txt'): # We only visualize for runs of specific examples
from visualization import *
# plotting fracture footprint
Fr_list, properties = load_fractures(
address='./Data/Wu_et_al', time_srs=simulProp.get_solTimeSeries())
Fig = plot_fracture_list(Fr_list,
variable='mesh',
backGround_param='sigma0',
mat_properties=Solid)
plot_prop = PlotProperties(line_color='darkmagenta')
Fig = plot_fracture_list(Fr_list,
variable='footprint',
fig=Fig,
plot_prop=plot_prop)
# plotting footprint from experiment
ax = Fig.get_axes()[0]
ax.plot(data[:, 0] * 1e-3, -1e-3 * data[:, 1], 'k')
ax.plot(data[:, 2] * 1e-3, -1e-3 * data[:, 3], 'k')
ax.plot(data[:, 4] * 1e-3, -1e-3 * data[:, 5], 'k')
# loading the experiment data file
import csv
with open('./wu_et_al_data.csv', 'r') as f:
reader = csv.reader(f)
data = np.asarray(list(reader), dtype=np.float64)
####################
# plotting results #
####################
from visualization import *
# plotting fracture footprint
Fr_list, properties = load_fractures(address='./Data/Wu_et_al',
time_srs=simulProp.get_solTimeSeries())
Fig = plot_fracture_list(Fr_list,
variable='mesh',
backGround_param='sigma0',
mat_properties=Solid)
plot_prop = PlotProperties(line_color='darkmagenta')
Fig = plot_fracture_list(Fr_list,
variable='footprint',
fig=Fig,
plot_prop=plot_prop)
# plotting footprint from experiment
ax = Fig.get_axes()[0]
ax.plot(data[:, 0] * 1e-3, -1e-3 * data[:, 1], 'k')
ax.plot(data[:, 2] * 1e-3, -1e-3 * data[:, 3], 'k')
ax.plot(data[:, 4] * 1e-3, -1e-3 * data[:, 5], 'k')
