def rsolve(xcenter, ycenter, num, nred):
hcase = get_case(num)
rcase = get_reduced(num, nred)
mu = rcase.parameter(xcenter=xcenter, ycenter=ycenter)
with util.time():
lhs = solvers.poisson(rcase, mu)
@click.option('--fast/--no-fast', default=False)
@util.common_args
def disp(fast):
case = get_case(fast)
print(case)
@main.command()
@click.option('--velocity', default=1.0)
@click.option('--height', default=1.0)
@click.option('--fast/--no-fast', default=False)
@util.common_args
def solve(velocity, height, fast):
case = get_case(fast)
mu = case.parameter(velocity=velocity, height=height)
with util.time():
lhs = solvers.stokes(case, mu)
visualization.velocity(case,
mu,
lhs,
figsize=(10, 10),
name='full',
colorbar=True)
visualization.pressure(case,
mu,
lhs,
figsize=(10, 10),
name='full',
colorbar=True)