From ed75f41ed3623ddb8d616fb1a180d018deac5c38 Mon Sep 17 00:00:00 2001
From: Klaus <klaus.boehnlein@tu-dresden.de>
Date: Wed, 13 Oct 2021 23:23:07 +0200
Subject: [PATCH] Update
---
src/Plot_Angle_Lemma1.4.py | 43 ++++++++++++++++++++++----------------
1 file changed, 25 insertions(+), 18 deletions(-)
diff --git a/src/Plot_Angle_Lemma1.4.py b/src/Plot_Angle_Lemma1.4.py
index 7f2cb216..a6c39d3b 100644
--- a/src/Plot_Angle_Lemma1.4.py
+++ b/src/Plot_Angle_Lemma1.4.py
@@ -112,8 +112,8 @@ theta = 1.0/4.0
lambda1 = 0.0
gamma = 1.0/4.0
-gamma = 'infinity'
-# gamma = '0'
+gamma = 'infinity' #Elliptic Setting
+# gamma = '0' #Hyperbolic Setting
# gamma = 0.5
@@ -169,8 +169,8 @@ xmax = 0.41
# xmin = 0.01
# xmax = 3.0
-# numPoints = 200
-numPoints = 101
+numPoints = 200
+# numPoints = 101
X_Values = np.linspace(xmin, xmax, num=numPoints)
print(X_Values)
@@ -347,7 +347,7 @@ ax = plt.axes((0.1,0.1,0.5,0.8))
# ax.tick_params(axis='x',which='major',direction='out',length=10,width=5,color='red',pad=15,labelsize=15,labelcolor='green',
# labelrotation=15)
ax.tick_params(axis='x',which='major', direction='out',pad=5,labelsize=10)
-
+ax.tick_params(axis='y',which='major', length=5, width=1, direction='out',pad=5,labelsize=10)
ax.xaxis.set_major_locator(MultipleLocator(0.05))
ax.xaxis.set_minor_locator(MultipleLocator(0.025))
@@ -436,9 +436,14 @@ ax.grid(True,which='major',axis='both',alpha=0.3)
# --- leave out jumps:
# ax.scatter(X_Values, Y_Values)
-# leave out jumps:
-ax.plot(X_Values[X_Values>=jump_xValues[0]], Y_Values[X_Values>=jump_xValues[0]] , 'blue')
-ax.plot(X_Values[X_Values<jump_xValues[0]], Y_Values[X_Values<jump_xValues[0]], 'blue')
+
+
+
+
+# --- leave out jumps:
+if gamma == 'infinity':
+ ax.plot(X_Values[X_Values>=jump_xValues[0]], Y_Values[X_Values>=jump_xValues[0]] , 'royalblue')
+ ax.plot(X_Values[X_Values<jump_xValues[0]], Y_Values[X_Values<jump_xValues[0]], 'royalblue')
@@ -456,8 +461,8 @@ ax.plot(X_Values[X_Values<jump_xValues[0]], Y_Values[X_Values<jump_xValues[0]],
# plt.plot([X_Values[0],X_Values[-1]], [Y_Values[0],Y_Values[-1]])
# plt.axis([0, 6, 0, 20])
-ax.set_xlabel(r"$\theta$")
-ax.set_ylabel('angle')
+ax.set_xlabel(r"volume fraction $\theta$")
+ax.set_ylabel(r"angle $\angle$")
# plt.ylabel('$\kappa$')
@@ -468,14 +473,16 @@ ax.set_ylabel('angle')
# Plot every other line.. not the jumps...
-tmp = 1
-# for idx, x in enumerate(x_plotValues):
-# if idx > 0 and tmp == 1:
-# # plt.plot([x_plotValues[idx-1],x_plotValues[idx]] ,[y_plotValues[idx-1],y_plotValues[idx]] )
-# ax.plot([x_plotValues[idx-1],x_plotValues[idx]] ,[y_plotValues[idx-1],y_plotValues[idx]] )
-# tmp = 0
-# else:
-# tmp = 1
+
+if gamma == '0':
+ tmp = 1
+ for idx, x in enumerate(x_plotValues):
+ if idx > 0 and tmp == 1:
+ # plt.plot([x_plotValues[idx-1],x_plotValues[idx]] ,[y_plotValues[idx-1],y_plotValues[idx]] )
+ ax.plot([x_plotValues[idx-1],x_plotValues[idx]] ,[y_plotValues[idx-1],y_plotValues[idx]], 'royalblue' )
+ tmp = 0
+ else:
+ tmp = 1
# plt.plot([x_plotValues[0],x_plotValues[1]] ,[y_plotValues[0],y_plotValues[1]] )
# plt.plot([x_plotValues[2],x_plotValues[3]] ,[y_plotValues[2],y_plotValues[3]] )
--
GitLab