2 This file is part of I am car.
4 I am car is nree software: you can redistribute it and/or modify
5 it under the terms of the GNU General Public License as published by
6 the Free Software Foundation, either version 3 of the License, or
7 (at your option) any later version.
9 I am car is distributed in the hope that it will be useful,
10 but WITHOUT ANY WARRANTY; without even the implied warranty of
11 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12 GNU General Public License for more details.
14 You should have received a copy of the GNU General Public License
15 along with I am car. If not, see <http://www.gnu.org/licenses/>.
23 #include <jsoncpp/json/json.h>
29 #include "rrtplanner.h"
30 #include "slotplanner.h"
40 //#define JSONLOGEDGES
41 //#define JSONLOGSAMPLES
47 //#define USE_INTERRUPT
57 #define SCREEN_WIDTH 1000
58 #define SCREEN_HEIGHT 1000
60 std::chrono::high_resolution_clock::time_point TSTART_;
61 std::chrono::high_resolution_clock::time_point TEND_;
64 void TSTART() {TSTART_ = std::chrono::high_resolution_clock::now();}
66 std::chrono::duration<float> DT_;
67 TEND_ = std::chrono::high_resolution_clock::now();
68 DT_ = std::chrono::duration_cast<std::chrono::duration<float>>(
71 TELAPSED += DT_.count();
72 ELAPSED = DT_.count();
74 void TPRINT(const char *what) {
75 std::chrono::duration<float> DT_;
76 DT_ = std::chrono::duration_cast<std::chrono::duration<float>>(
79 std::cerr << what << ": " << DT_.count() << std::endl;
82 bool run_planner = true;
85 SDL_Window* gw = NULL;
103 void *next_run(void *arg)
105 struct next_arg *na = (struct next_arg *) arg;
106 T2 *lp = (T2 *) na->p;
108 while (!*gf && lp->elapsed() < TMAX) {
120 float new_x = 1 + static_cast<float>(rand()) /
121 static_cast<float>(RAND_MAX / (6.6 - 1 - 1));
123 float new_h = M_PI / 2;
124 return new RRTNode(new_x, new_y, new_h);
129 srand(static_cast<unsigned>(time(0)));
134 Json::Value jvi; // JSON input
135 Json::Value jvo; // JSON output
139 std::string encoding = jvi.get("encoding", "UTF-8" ).asString();
141 if (jvi["init"] == Json::nullValue) {
142 std::cerr << "I need `init` in JSON input scenario";
143 std::cerr << std::endl;
147 if (jvi["goal"] == Json::nullValue) {
148 std::cerr << "I need `goal` in JSON input scenario";
149 std::cerr << std::endl;
155 jvi["init"][0].asFloat(),
156 jvi["init"][1].asFloat(),
157 jvi["init"][2].asFloat()
160 jvi["goal"][0].asFloat(),
161 jvi["goal"][1].asFloat(),
162 jvi["goal"][2].asFloat()
165 jvo["init"][0] = p.root()->x();
166 jvo["init"][1] = p.root()->y();
167 jvo["init"][2] = p.root()->h();
168 std::vector<CircleObstacle> co;
169 std::vector<SegmentObstacle> so;
170 p.HMIN = p.root()->x();
171 p.HMAX = p.root()->x();
172 p.VMIN = p.root()->y();
173 p.VMAX = p.root()->y();
175 for (auto o: jvi["obst"]) {
180 for (i = 0; i < o.size() - 1; i++) {
181 tmpx = o[i][0].asFloat();
182 tmpy = o[i][1].asFloat();
183 tmpr = o[i + 1][0].asFloat();
184 tmps = o[i + 1][1].asFloat();
185 so.push_back(SegmentObstacle(
186 new RRTNode(tmpx, tmpy, 0),
187 new RRTNode(tmpr, tmps, 0)
189 p.frame().add_bnode(so.back().init());
190 if (tmpx < p.HMIN) p.HMIN = tmpx;
191 if (tmpx > p.HMAX) p.HMAX = tmpx;
192 if (tmpy < p.VMIN) p.VMIN = tmpy;
193 if (tmpy > p.VMAX) p.VMAX = tmpy;
194 if (tmpr < p.HMIN) p.HMIN = tmpr;
195 if (tmpr > p.HMAX) p.HMAX = tmpr;
196 if (tmps < p.VMIN) p.VMIN = tmps;
197 if (tmps > p.VMAX) p.VMAX = tmps;
200 jvo["obst"][j][i][0] = tmpx;
201 jvo["obst"][j][i][1] = tmpy;
203 jvo["obst"][j][i][0] = tmpr;
204 jvo["obst"][j][i][1] = tmps;
207 p.defaultSamplingInfo();
208 p.link_obstacles(&co, &so);
212 ParallelSlot ps = ParallelSlot();
214 jvi["slot"] != Json::nullValue &&
215 jvi["slot"]["polygon"] != Json::nullValue
217 for (auto xy: jvi["slot"]["polygon"]) {
218 ps.slot().add_bnode(new RRTNode(
223 for (auto e: ps.slot().frame())
224 so.push_back(SegmentObstacle(e->init(), e->goal()));
226 #ifdef USE_SLOTPLANNER
228 if (ps.slot().bnodes().size() > 0)
231 jvo["ppse"] = ELAPSED;
232 TPRINT("ParallelSlot");
234 if (ps.slot().bnodes().size() > 0) {
236 //if (ps.getMidd() != nullptr)
237 // p.setSamplingInfo(ps.getSamplingInfo());
239 if (ps.cusp().size() > 0) {
240 p.goal(ps.getMidd());
241 p.slot_cusp(ps.cusp().front()); // use first found solution
243 jvo["midd"][0] = p.goal()->x();
244 jvo["midd"][1] = p.goal()->y();
245 jvo["midd"][2] = p.goal()->h();
246 jvo["goal"][0] = p.slot_cusp().back()->x();
247 jvo["goal"][1] = p.slot_cusp().back()->y();
248 jvo["goal"][2] = p.slot_cusp().back()->h();
250 jvo["goal"][0] = p.goal()->x();
251 jvo["goal"][1] = p.goal()->y();
252 jvo["goal"][2] = p.goal()->h();
255 std::cerr << "Slot Info:" << std::endl;
256 if (ps.slotSide() == LEFT)
257 std::cerr << "- LEFT" << std::endl;
259 std::cerr << "- RIGHT" << std::endl;
260 if (ps.slotType() == PARALLEL)
261 std::cerr << "- PARALLEL" << std::endl;
263 std::cerr << "- PERPENDICULAR" << std::endl;
265 std::vector<RRTNode *> steered;
266 for (auto jn: jvi["traj"][0]) {
267 steered.push_back(new RRTNode(
274 std::reverse(steered.begin(), steered.end());
275 RRTNode *pn = p.root();
276 for (auto n: steered) {
279 pn->add_child(n, p.cost(pn, n));
283 pn->add_child(p.goal(), p.cost(pn, p.goal()));
292 #elif defined USE_INTERRUPT
293 signal(SIGINT, hint);
294 signal(SIGTERM, hint);
296 while (run_planner) {
303 #elif defined USE_TMAX
306 while (!p.goal_found() && p.elapsed() < TMAX) {
313 if (p.goal_found() && ps.slotType() == PARALLEL)
314 p.tlog(p.findt(p.slot_cusp().back()));
315 #elif defined USE_PTHREAD
317 RRTNode *ron = nullptr;
318 RRTNode *gon = nullptr;
320 pthread_t rt; // root thread
321 pthread_t gt; // goal thread
322 pthread_t ct; // connect thread
335 pthread_create(&rt, NULL, &next_run, (void *) &ra);
336 pthread_create(>, NULL, &next_run, (void *) &ga);
340 while (!gf && p.elapsed() < TMAX &&
341 p.p_root_.nodes().size() < NOFNODES &&
342 p.p_goal_.nodes().size() < NOFNODES) {
345 for (int i = 0; i < IXSIZE; i++) {
346 for (int j = 0; j < IYSIZE; j++) {
347 if (p.p_root_.ixy_[i][j].changed() &&
348 p.p_goal_.ixy_[i][j].changed()) {
350 for (auto rn: p.p_root_.ixy_[i][j].nodes()) {
351 for (auto gn: p.p_goal_.ixy_[i][j].nodes()) {
352 if (rn->ccost() + gn->ccost() < mc &&
358 mc = rn->ccost() + gn->ccost();
366 if (p.elapsed() >= TMAX)
369 // end of overlap trees
373 pthread_join(rt, NULL);
374 pthread_join(gt, NULL);
375 float nodo = ((float) ndl / (float) tol);
376 std::cerr << "nothing done is " << 100.0 * nodo;
377 std::cerr << "%" << std::endl;
378 //std::cerr << "rgf is " << p.p_root_.goal_found() << std::endl;
379 //std::cerr << "ggf is " << p.p_goal_.goal_found() << std::endl;
380 //std::cerr << "cgf is " << p.goal_found() << std::endl;
381 if (p.p_root_.goal_found() && p.p_root_.goal()->ccost() < mc) {
382 ron = p.p_root_.goal()->parent();
383 gon = p.p_root_.goal();
384 mc = p.p_root_.goal()->ccost();
386 if (p.p_goal_.goal_found() && p.p_goal_.goal()->ccost() < mc) {
387 ron = p.p_goal_.goal();
388 gon = p.p_goal_.goal()->parent();
389 mc = p.p_goal_.goal()->ccost();
391 p.root()->remove_parent(); // needed if p.p_goal_.goal_found()
393 p.goal()->children().clear();
396 while (gon != p.goal()) {
397 p.p_root_.nodes().push_back(new RRTNode(
402 p.p_root_.nodes().back(),
405 p.p_root_.nodes().back()));
406 ron = p.p_root_.nodes().back();
409 ron->add_child(p.goal(), p.p_root_.cost(ron, p.goal()));
411 // end of connect trees
419 jvo["rrte"] = ELAPSED;
424 // statistics to error output
425 std::cerr << "TELAPSED is " << TELAPSED << std::endl;
426 std::cerr << "Elapsed is " << p.elapsed() << std::endl;
427 std::cerr << "Goal found is " << p.goal_found() << std::endl;
428 std::cerr << "#nodes is " << p.nodes().size() << std::endl;
429 std::cerr << "#samples is " << p.samples().size() << std::endl;
430 std::cerr << "`tlog` size is " << p.tlog().size() << std::endl;
431 std::cerr << "trajectories costs:" << std::endl;
432 for (j = 0; j < p.clog().size(); j++)
433 std::cerr << "- " << p.clog()[j] << std::endl;
434 std::cerr << "RRT #nodes:" << std::endl;
435 for (j = 0; j < p.nlog().size(); j++)
436 std::cerr << "- " << p.nlog()[j] << std::endl;
437 std::cerr << "trajectories seconds:" << std::endl;
438 for (j = 0; j < p.slog().size(); j++)
439 std::cerr << "- " << p.slog()[j] << std::endl;
440 std::cerr << "RRT edges (from root) log size: " << p.rlog().size();
441 std::cerr << std::endl;
442 for (auto edges: p.rlog())
443 std::cerr << "- " << edges.size() << std::endl;
446 jvo["elap"] = TELAPSED;
448 jvo["nodo"][0] = nodo;
451 for (j = 0; j < p.clog().size(); j++)
452 jvo["cost"][j] = p.clog()[j];
454 for (j = 0; j < p.nlog().size(); j++)
455 jvo["node"][j] = p.nlog()[j];
457 for (j = 0; j < p.slog().size(); j++)
458 jvo["secs"][j] = p.slog()[j];
462 for (auto traj: p.tlog()) {
465 jvo["traj"][j][i][0] = n->x();
466 jvo["traj"][j][i][1] = n->y();
467 jvo["traj"][j][i][2] = n->h();
468 jvo["traj"][j][i][3] = n->t();
469 jvo["traj"][j][i][4] = n->s();
477 for (auto edges: p.rlog()) {
479 for (auto e: edges) {
480 jvo["edge"][i][j][0][0] = e->init()->x();
481 jvo["edge"][i][j][0][1] = e->init()->y();
482 jvo["edge"][i][j][0][2] = e->init()->h();
483 jvo["edge"][i][j][1][0] = e->goal()->x();
484 jvo["edge"][i][j][1][1] = e->goal()->y();
485 jvo["edge"][i][j][1][2] = e->goal()->h();
491 #ifdef JSONLOGSAMPLES
494 for (auto s: p.samples()) {
495 jvo["samp"][j][0] = s->x();
496 jvo["samp"][j][1] = s->y();
497 jvo["samp"][j][2] = s->h();
502 std::cout << jvo << std::endl;
505 SDL_DestroyWindow(gw);
520 if (SDL_Init(SDL_INIT_VIDEO) < 0) {
521 std::cerr << "SDL could not initialize! SDL_Error: ";
522 std::cerr << SDL_GetError();
523 std::cerr << std::endl;
526 SDL_GL_SetAttribute(SDL_GL_CONTEXT_MAJOR_VERSION, 2);
527 SDL_GL_SetAttribute(SDL_GL_CONTEXT_MINOR_VERSION, 1);
528 gw = SDL_CreateWindow(
530 SDL_WINDOWPOS_UNDEFINED,
531 SDL_WINDOWPOS_UNDEFINED,
534 SDL_WINDOW_OPENGL | SDL_WINDOW_SHOWN);
536 std::cerr << "Window could not be created! SDL_Error: ";
537 std::cerr << SDL_GetError();
538 std::cerr << std::endl;
541 gc = SDL_GL_CreateContext(gw);
543 std::cerr << "OpenGL context couldn't be created! SDL Error: ";
544 std::cerr << SDL_GetError();
545 std::cerr << std::endl;
548 if (SDL_GL_SetSwapInterval(1) < 0) {
549 std::cerr << "Warning: Unable to set VSync! SDL Error: ";
550 std::cerr << SDL_GetError();
551 std::cerr << std::endl;
555 std::cerr << "Unable to initialize OpenGL!";
556 std::cerr << std::endl;
564 GLenum error = GL_NO_ERROR;
565 glMatrixMode(GL_PROJECTION);
567 error = glGetError();
568 if (error != GL_NO_ERROR) {
569 std::cerr << "Error initializing OpenGL! ";
570 std::cerr << gluErrorString(error);
571 std::cerr << std::endl;
574 glMatrixMode(GL_MODELVIEW);
576 error = glGetError();
577 if (error != GL_NO_ERROR) {
578 std::cerr << "Error initializing OpenGL! ";
579 std::cerr << gluErrorString(error);
580 std::cerr << std::endl;
583 glClearColor(1, 1, 1, 1);
584 error = glGetError();
585 if (error != GL_NO_ERROR) {
586 std::cerr << "Error initializing OpenGL! ";
587 std::cerr << gluErrorString(error);
588 std::cerr << std::endl;