Visualization of the refinement of a coarse grid

The coarse and the fine mesh can be visualized with simple plots. It is important to notice that the fine mesh works as a refinement of the coarse mesh.


In [1]:
import os
import sys
import numpy as np

%matplotlib notebook
import matplotlib.pyplot as plt
from matplotlib import cm

#coarse World
NWorldCoarse = np.array([11,11])
NpCoarse = np.prod(NWorldCoarse+1)

A = np.zeros(NWorldCoarse)
ABase = A.flatten()   

aCube = ABase.reshape(NWorldCoarse)

Coarse mesh


In [2]:
fig = plt.figure('Coarse Mesh')                                                               
ax = fig.add_subplot(1,1,1)     
major_ticks = np.arange(0, 11, 1)
ax.imshow(aCube, extent=[0, 11, 0, 11], cmap='Greys')
ax.axis([0, 11, 0, 11])
ax.set_xticks(major_ticks)                                                       
ax.set_yticks(major_ticks)  
ax.tick_params(axis='both', which='both', bottom='off', top='off', labelbottom='off', right='off', left='off', labelleft='off')
fig.subplots_adjust(left=0.00,bottom=0.02,right=1,top=0.95,wspace=0.2,hspace=0.2)
ax.grid(which='major', linestyle="-", color="black", alpha=0.8)
plt.title('Coarse Grid')
plt.show()


Refinement


In [3]:
fig = plt.figure('Refinement')                                                               
ax = fig.add_subplot(1,1,1)     
major_ticks = np.arange(0, 11, 1)
minor_ticks = np.arange(0, 11, 0.2)                                                
ax.imshow(aCube, extent=[0, 11, 0, 11], cmap='Greys')
ax.axis([0, 11, 0, 11])
ax.set_xticks(minor_ticks)
ax.set_xticks(minor_ticks)
ax.set_xticks(major_ticks)                                                       
ax.set_yticks(major_ticks)  
ax.minorticks_on()                                                     
ax.tick_params(axis='both', which='both', bottom='off', top='off', labelbottom='off', right='off', left='off', labelleft='off')
fig.subplots_adjust(left=0.00,bottom=0.02,right=1,top=0.95,wspace=0.2,hspace=0.2)
ax.grid(which='major', linestyle="-", color="black", alpha=0.8)                                                
ax.grid(which='minor', linestyle="-",linewidth=0.5 , color="black", alpha=0.4)
plt.title('Refinement')
plt.show()