def check_penaro(problemname='cylinderwake', N=2, Re=1.0e2, t0=0.0, tE=2.0, Nts=1e3+1, plot=False): femp = dnsps.cyl_fems(N, bccontrol=True, verbose=True) V = femp['V'] nu = femp['charlen']/Re bcsd = femp['contrbcssubdomains'] bcssfuns = femp['contrbcsshapefuns'] stokesmats = dts.get_stokessysmats(V, femp['Q'], nu=nu, bccontrol=True, cbclist=bcsd, cbshapefuns=bcssfuns) if plot: plt.figure(2) plt.spy(stokesmats['amatrob']) print 'Number of nonzeros in amatrob:', stokesmats['amatrob'].nnz plt.show()
def check_penaro(problemname='cylinderwake', N=2, Re=1.0e2, t0=0.0, tE=2.0, Nts=1e3+1, plot=False): femp = dnsps.cyl_fems(N, bccontrol=True, verbose=True) V = femp['V'] nu = femp['charlen']/Re bcsd = femp['contrbcssubdomains'] bcssfuns = femp['contrbcsshapefuns'] stokesmats = dts.get_stokessysmats(V, femp['Q'], nu=nu, bccontrol=True, cbclist=bcsd, cbshapefuns=bcssfuns) if plot: plt.figure(2) plt.spy(stokesmats['amatrob']) print('Number of nonzeros in amatrob:', stokesmats['amatrob'].nnz) plt.show()
from dolfin_navier_scipy.problem_setups import cyl_fems import optical_test_outpop as too import optical_test_inpop as tio N = 2 NU = 5 NY = 5 femp = cyl_fems(N) too.check_output_opa(NY, femp) tio.check_input_opa(NU, femp)