Beispiel #1
0
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)