← Back to branch summary

~unicorn-core/unicorn/hpc

« back to all changes in this revision

Viewing changes to ucsolver/icns/unicorn/NSEContinuity3D.form

  • Committer: Niclas Jansson
  • Date: 2012-03-01 08:20:27 UTC
  • mfrom: (129.5.32 0.1.3-hpc)
  • Revision ID: njansson@csc.kth.se-20120301082027-snah9zahbb5d3t8n
Fixed F77 calling in elastic smoother

Show diffs side-by-side

added added

removed removed

Lines of Context:
ucsolver/icns/unicorn/NSEContinuity3D.form
26
26
 
27
27
delta1 = Function(constant_scalar) # stabilization parameter
28
28
 
 
29
P0  = Function(scalar) # trial basis function
 
30
k  = Constant("tetrahedron") # time step
 
31
 
29
32
i0 = Index() # index for tensor notation
30
33
i1 = Index() # index for tensor notation
31
34
 
32
35
def ugradu(u, v):
33
36
    return [dot(u, grad(v[i])) for i in range(d)]
34
37
 
 
38
alpha1 = 2.0 * k
 
39
 
35
40
# Bilinear and linear forms
36
 
a = delta1*dot(grad(q), grad(P))*dx + 0.001*P*q*dx
37
 
#L = - q*div(uc)*dx - delta1*dot(dot(uc,grad(uc)),grad(q))
38
 
L = - q*div(uc)*dx - mult(delta1, dot( ugradu(uc, uc), grad(q))) *dx
39
 
 
40
 
 
41
 
 
 
41
a = alpha1*dot(grad(q), grad(P))*dx + delta1*dot(grad(q), grad(P))*dx + 0.001*P*q*dx
 
42
L = alpha1*dot(grad(q), grad(P0))*dx - q*div(uc)*dx - mult(delta1, dot( ugradu(uc, uc), grad(q))) *dx
 
43
#L = alpha1*dot(grad(q), grad(P0))*dx - q*div(uc)*dx