Skip to content
Snippets Groups Projects
Commit d00fe7c9 authored by David Seus's avatar David Seus
Browse files

uiea

parent aa4991f6
Branches
Tags
No related merge requests found
......@@ -352,6 +352,10 @@ cutoff = gaussian/(gaussian + zero_on_shrinking)
# 'nonwetting': (-1 -t*(1+y + x**2)**2)*cutoff},
# }
p_e_sym = {
0: {'wetting': -6 -(1 + x*x + y*y) + 0*t}
# 'nonwetting': -1 -(sym.sin(3*(1+y)/2*sym.pi)*sym.sin(5*(1+x)/2*sym.pi))**2},
}
print(f"\n\n\nsymbolic type is {type(p_e_sym[0]['wetting'])}\n\n\n")
# # pw_sym_x*pw_sym_y
# p_e_sym = {
......@@ -364,11 +368,6 @@ print(f"\n\n\nsymbolic type is {type(p_e_sym[0]['wetting'])}\n\n\n")
# 'nonwetting': -1*zero_on_shrinking+ 0*t},
# }
p_e_sym = {
0: {'wetting': -6 -(1 + x*x + y*y) + 0*t}
# 'nonwetting': -1 -(sym.sin(3*(1+y)/2*sym.pi)*sym.sin(5*(1+x)/2*sym.pi))**2},
}
pc_e_sym = dict()
for subdomain, isR in isRichards.items():
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment