__author__ = 'benjamin'
from dc2D import tworesolution_dual_contour
from dcSample import sample_data, sphere_f, doubletorus_f, torus_f
import numpy as np
dimensions = {'xmin': 0.0, 'xmax': 6.0, 'ymin': 2.0, 'ymax': 6.0}
res_fine = 1.0/16.0
res_coarse = res_fine * 8.0
resolutions = {'fine': res_fine,'coarse': res_coarse}
fine_data = sample_data(doubletorus_f, resolutions['fine'], dimensions)
[verts_out_dc, quads_out_dc] = tworesolution_dual_contour(fine_data, resolutions, dimensions)
import matplotlib.pyplot as plt
fig = plt.figure()
for q in quads_out_dc['fine']:
vtx = verts_out_dc['fine'][q]
x = vtx[:,0]
y = vtx[:,1]
#plt.plot(x, y, 'b')
#plt.plot(x, y, 'bo')
for q in quads_out_dc['coarse']:
vtx = verts_out_dc['coarse'][q]
x = vtx[:,0]
y = vtx[:,1]
from redis.utils import dummy
__author__ = 'benjamin'
from dc2D import dual_contour
from dcSample import sample_data, sphere_f, doubletorus_f, torus_f
from dcPlotting import plot_edges, plot_vertices, plot_non_manifold_vertices
from dcManifolds import detectManifolds2d
dimensions = {'xmin': 0.0, 'xmax': 8.0, 'ymin': 0.0, 'ymax': 8.0}
resolution = 1.0/4.0
data = sample_data(doubletorus_f, resolution, dimensions)
[dc_verts, dc_edges] = dual_contour(data, resolution, dimensions)
non_manifold_verts = detectManifolds2d(dc_edges)
import matplotlib.pyplot as plt
fig = plt.figure()
#plot_vertices(dc_verts)
plot_edges(dc_edges, dc_verts)
plot_non_manifold_vertices(dc_verts, non_manifold_verts)
plt.xlim(xmin = 0, xmax = 8)
plt.ylim(ymin = 0, ymax = 8)
plt.show()
