    Uf = []
    for gdl in modeloL.geometria.gdls:
        u = None
        for e in modeloH.elementos:
            u = e.solucionInterpoladaLocal(gdl[0],gdl[1])
            if u:
                break
        if not u:
            raise Exception(gdl)
            u = 10
        Uf.append(u)
    Uf = np.array(Uf).reshape([len(Uf),1])
    error = lambda u,u_gt: np.abs((u-u_gt)/u_gt)
    errores = error(Ug,Uf)
