Update

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]] )