Python Ones Examples

Example #1

File: operators.py Project: majosm/grudge

def integral(arg, dd=None):
    import grudge.symbolic.primitives as prim

    if dd is None:
        dd = prim.DD_VOLUME
    dd = prim.as_dofdesc(dd)

    return NodalSum(dd)(arg * prim.cse(
        MassOperator(dd_in=dd)(prim.Ones(dd)), "mass_quad_weights",
        prim.cse_scope.DISCRETIZATION))

Example #2

File: operators.py Project: majosm/grudge

def h_max_from_volume(ambient_dim, dim=None, dd=None):
    """Defines a characteristic length based on the volume of the elements.
    This length may not be representative if the elements have very high
    aspect ratios.
    """

    import grudge.symbolic.primitives as prim
    if dd is None:
        dd = prim.DD_VOLUME
    dd = prim.as_dofdesc(dd)

    if dim is None:
        dim = ambient_dim

    return NodalMax(dd_in=dd)(ElementwiseSumOperator(dd)(
        MassOperator(dd_in=dd)(prim.Ones(dd))))**(1.0 / dim)