diff --git a/buzz_scripts/graphformGPS.bzz b/buzz_scripts/graphformGPS.bzz index 324c52e..f3b0bd5 100644 --- a/buzz_scripts/graphformGPS.bzz +++ b/buzz_scripts/graphformGPS.bzz @@ -729,7 +729,7 @@ function init() { #v_tag = stigmergy.create(m_lockstig) #uav_initstig() # go to diff. height since no collision avoidance implemented yet - TARGET_ALTITUDE = 3.0 + id * 1.25 + #TARGET_ALTITUDE = 3.0 + id * 1.25 statef=turnedoff UAVSTATE = "TURNEDOFF" } diff --git a/buzz_scripts/include/rttstar.bzz b/buzz_scripts/include/rttstar.bzz new file mode 100644 index 0000000..cfd4924 --- /dev/null +++ b/buzz_scripts/include/rttstar.bzz @@ -0,0 +1,233 @@ +##### +# RRT* Path Planing +# START, GOAL, TOL vec2 +# map table-based matrix +##### + +map = {.nb_col=100, .nb_row=100, .mat={}} + +function RRTSTAR(START,GOAL,TOL) { + # RRT_STAR + HEIGHT = map.nb_col + WIDTH = map.nb_row + RADIUS = 2 + + rng.setseed(11) + + goalBoundary = {.xmin=GOAL.x-TOL.x, .xmax=GOAL.x+TOL.x, .ymin=GOAL.y-TOL.y, .ymax=GOAL.y+TOL.y} + arrayOfPoints = {.nb_col=2, .nb_row=1, .mat={.0=START.x,.1=START.y}} + numberOfPoints = 1 + Q = {.nb_col=5,.nb_row=1,.mat={.0=START.x,.1=START.y,.2=0,.3=1,.4=0}} + log("Created ",Q) + + # Open space or obstacle + # DISPLAY_patchwork(map); + + goalReached = 0; + while(goalReached == 0) { + # Point generation + x = -HEIGHT*rng.uniform(1.0)-1; + y = WIDTH*rng.uniform(1.0)+1; + log("First trial point (", rng.uniform(1.0), "): ", x, " ", y) + + pointList = findPointsInRadius(x,y,Q,RADIUS); + + log(pointList.nb_col,pointList.nb_row,pointList.mat) + + # Find connection that provides the least cost to come + nbCount = 0; + minCounted = inf; + minCounter = 0; + + if(pointList.nb_col!=0) { + i=0 + while(i 0) { + arrayOfPoints.nb_row=arrayOfPoints.nb_row+1 + arrayOfPoints.mat[arrayOfPoints.nb_row]=x + arrayOfPoints.mat[arrayOfPoints.nb_row+1]=y + numberOfPoints = numberOfPoints + 1; + + wmat(Q,numberOfPoints,0, x) + wmat(Q,numberOfPoints,1, y) + wmat(Q,numberOfPoints,2, rmat(pointList,minCounter,4)); + wmat(Q,numberOfPoints,3, numberOfPoints) + wmat(Q,numberOfPoints,4, minCounted) + + # Now check to see if any of the other points can be redirected + nbCount = 0; + i = 0 + while(i goalBoundary.xmin and y > goalBoundary.ymin and y < goalBoundary.ymax) + goalReached = 1; + } + } else { + # Associate with the closest point + pointNum = findClosestPoint(x,y,arrayOfPoints); + + # Follow the line to see if it intersects anything + intersects = doesItIntersect(x,y,math.vec2.new(rmat(arrayOfPoints,1,pointNum),rmat(arrayOfPoints,2,pointNum)),map); + + # If there is no intersection we need to add to the tree + if(intersects != 1) { + arrayOfPoints.nb_row=arrayOfPoints.nb_row+1 + arrayOfPoints.mat[arrayOfPoints.nb_row]=x + arrayOfPoints.mat[arrayOfPoints.nb_row+1]=y + numberOfPoints = numberOfPoints + 1; + + wmat(Q,numberOfPoints,0, x) + wmat(Q,numberOfPoints,1, y) + wmat(Q,numberOfPoints,2, pointNum); + wmat(Q,numberOfPoints,3, numberOfPoints) + wmat(Q,numberOfPoints,4, rmat(Q,pointNum,5)+math.vec2.length(rmat(Q,pointNum,1)-x,rmat(Q,pointNum,2)-y)) + + # Check to see if this new point is within the goal + if(x < goalBoundary.xmax and x > goalBoundary.xmin and y > goalBoundary.ymin and y < goalBoundary.ymax) + goalReached = 1; + } + } + if(numberOfPoints % 100 == 0) { + log(numberOfPoints, " points processed. Still looking for goal."); + } + } + log("Goal found!"); +} + +function findClosestPoint(x,y,arrayOfPoints) { + # Go through each points and find the distances between them and the + # target point + distance = inf; + i=1 + while(i 0 and yi > 0) { + if(rmap(xi,yi) == 0) + return 1; + } else + return 1; + } + i = i + 1 + } + return 0 +} + +# Write to matrix +function wmat(mat, row, col, val) { + var index = (row-1)*mat.nb_col + col + mat.mat[index] = val +} + +# Read from matrix +function rmat(mat, row, col) { + var index = (row-1)*mat.nb_col + col + if (mat.mat[index] == nil) { + return -1 + } else { + return mat.mat[index] + } +} + +# copy a full matrix row +function mat_copyrow(out,ro,in,ri){ + var indexI = (ri-1)*in.nb_col + var indexO = (ro-1)*out.nb_col + i=0 + while(ilsyms->syms); ++i) { + for(uint32_t i = 1; i < buzzdarray_size(vm->lsyms->syms); ++i) { buzzvm_lload(vm, i); buzzobj_t o = buzzvm_stack_at(vm, 1); buzzvm_pop(vm);