[scipts] update PoU visualization

parent d5a397af
......@@ -35,14 +35,20 @@ def visualize(dd, pou, filename, local_space_type):
if __name__ == '__main__':
local_space_type = 'cg_p1'
num_local_refinements = 4
num_subdomains_per_dim = 3
num_local_refinements = 3
num_subdomains_per_dim = 1
num_elements_per_subdomain_and_dim = 2**num_local_refinements
dd = DomainDecomposition([num_subdomains_per_dim, num_subdomains_per_dim], num_local_refinements, [0, 0], [1, 1])
cg_PoU = ContinuousLagrangePartitionOfUnity(dd)
visualize(dd, cg_PoU, 'cg_PoU', local_space_type)
ft_PoU = ContinuousFlatTopPartitionOfUnity(dd, overlap=1 - 2/num_elements_per_subdomain_and_dim)
visualize(dd, ft_PoU, 'ft_PoU_Hminus2hOvlp', local_space_type)
ft_PoU = ContinuousFlatTopPartitionOfUnity(dd, overlap=0.5)
visualize(dd, ft_PoU, 'ft_PoU', local_space_type)
visualize(dd, ft_PoU, 'ft_PoU_halfHOvlp', local_space_type)
ft_PoU = ContinuousFlatTopPartitionOfUnity(dd, overlap=2/num_elements_per_subdomain_and_dim)
visualize(dd, ft_PoU, 'ft_PoU_2hOvlp', local_space_type)
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment