Exemplo n.º 1
0
 def bilinf(u, p, v, q, w):
     from skfem.helpers import grad, ddot, div
     return (ddot(grad(u), grad(v)) - div(u) * q - div(v) * p -
             1e-2 * p * q)
Exemplo n.º 2
0
def bilinf(sigma, u, tau, v, w):
    return dot(sigma, tau) + div(sigma) * v + div(tau) * u
Exemplo n.º 3
0
def phipsi_load3(u, v, w):
    '''
    for 5.7b $(div_phi, p)$
    '''
    return div(u) * v
Exemplo n.º 4
0
def divu(u, v, w):
    return div(u) * v
Exemplo n.º 5
0
def phiq_load(u, v, w):
    '''
    for 5.7c $(div_phi, q)$
    '''
    return div(u) * v