class ElementTriArgyris(ElementH2):

    nodal_dofs = 6

    facet_dofs = 1

    dim = 2

    maxdeg = 5

    def gdof(self, u, du, ddu, v, e, n, i):

        if i < 18:

            j = i % 6

            k = int(i/6)

            if j == 0:

                return u(*v[k])

            elif j == 1:

                return du[0](*v[k])

            elif j == 2:

                return du[1](*v[k])

            elif j == 3:

                return ddu[0](*v[k])

            elif j == 4:

                return ddu[1](*v[k])

            elif j == 5:

                return ddu[2](*v[k])

        elif i == 18:

            return du[0](*e[0])*n[0, 0] + du[1](*e[0])*n[0, 1]

        elif i == 19:

            return du[0](*e[1])*n[1, 0] + du[1](*e[1])*n[1, 1]

        elif i == 20:

            return du[0](*e[2])*n[2, 0] + du[1](*e[2])*n[2, 1]