Exemplo n.º 1
0
def random_problem_2D(numpoints, radius=10.0, roundoff=0.0, angleratio=0.5):
    """Generate a random problem with given number of points, a roundoff
       value for the prototype points, a radius for the cloud of prototype points
       and a ratio of angle constraints over distance constraints"""
    group = {}
    problem = GeometricProblem(dimension=2)
    i = 0
    while i < numpoints:
        aname = 'p'+str(i)
        apoint = vector([
            _round(random.uniform(-radius,radius),roundoff),
            _round(random.uniform(-radius,radius),roundoff)
        ])
        unique = True
        for v in group:
            p = group[v]
            if tol_eq(apoint[0],p[0]) and tol_eq(apoint[1],p[1]):
                unique = False
                break
        if unique:
                problem.add_point(aname, apoint)
                group[aname] = apoint
                i = i + 1
    #next
    _constraint_group(problem, group, None, angleratio)
    return problem
Exemplo n.º 2
0
def random_problem_2D(numpoints, radius=10.0, roundoff=0.0, angleratio=0.5):
    """Generate a random problem with given number of points, a roundoff
       value for the prototype points, a radius for the cloud of prototype points
       and a ratio of angle constraints over distance constraints"""
    group = {}
    problem = GeometricProblem(dimension=2)
    i = 0
    while i < numpoints:
        aname = 'p' + str(i)
        apoint = vector([
            _round(random.uniform(-radius, radius), roundoff),
            _round(random.uniform(-radius, radius), roundoff)
        ])
        unique = True
        for v in group:
            p = group[v]
            if tol_eq(apoint[0], p[0]) and tol_eq(apoint[1], p[1]):
                unique = False
                break
        if unique:
            problem.add_point(aname, apoint)
            group[aname] = apoint
            i = i + 1
    #next
    _constraint_group(problem, group, None, angleratio)
    return problem
Exemplo n.º 3
0
def random_distance_problem_3D(npoints, radius, roundoff):
    """creates a 3D problem with random distances"""
    problem = GeometricProblem(dimension=3)
    for i in range(npoints):
        # add point
        newvar = 'v' + str(i)
        newpoint = vector([
            _round(random.uniform(-radius, radius), roundoff),
            _round(random.uniform(-radius, radius), roundoff),
            _round(random.uniform(-radius, radius), roundoff)
        ])
        sellist = list(problem.cg.variables())
        problem.add_point(newvar, newpoint)
        # add distance constraints
        for j in range(min(3, len(sellist))):
            index = random.randint(0, len(sellist) - 1)
            var = sellist.pop(index)
            point = problem.get_point(var)
            dist = distance_2p(point, newpoint)
            problem.add_constraint(DistanceConstraint(var, newvar, dist))
    return problem
Exemplo n.º 4
0
def random_distance_problem_3D(npoints, radius, roundoff):
	"""creates a 3D problem with random distances"""
	problem = GeometricProblem(dimension=3)
	for i in range(npoints):
		# add point
		newvar = 'v'+str(i)
		newpoint = vector([
            _round(random.uniform(-radius,radius),roundoff),
            _round(random.uniform(-radius,radius),roundoff),
            _round(random.uniform(-radius,radius),roundoff)
        ])
		sellist = list(problem.cg.variables())
		problem.add_point(newvar, newpoint)
		# add distance constraints    
		for j in range(min(3,len(sellist))):
			index = random.randint(0,len(sellist)-1)
			var = sellist.pop(index)
			point = problem.get_point(var)
			dist = distance_2p(point,newpoint)
			problem.add_constraint(DistanceConstraint(var,newvar,dist))
	return problem