diff --git a/.idea/.gitignore b/.idea/.gitignore
new file mode 100644
index 0000000..73f69e0
--- /dev/null
+++ b/.idea/.gitignore
@@ -0,0 +1,8 @@
+# Default ignored files
+/shelf/
+/workspace.xml
+# Datasource local storage ignored files
+/dataSources/
+/dataSources.local.xml
+# Editor-based HTTP Client requests
+/httpRequests/
diff --git a/.idea/misc.xml b/.idea/misc.xml
new file mode 100644
index 0000000..339cee9
--- /dev/null
+++ b/.idea/misc.xml
@@ -0,0 +1,6 @@
+
+
+
+
+
+
\ No newline at end of file
diff --git a/.idea/modules.xml b/.idea/modules.xml
new file mode 100644
index 0000000..1ef27e7
--- /dev/null
+++ b/.idea/modules.xml
@@ -0,0 +1,8 @@
+
+
+
+
+
+
+
+
\ No newline at end of file
diff --git a/.idea/rpg_svo.iml b/.idea/rpg_svo.iml
new file mode 100644
index 0000000..f08604b
--- /dev/null
+++ b/.idea/rpg_svo.iml
@@ -0,0 +1,2 @@
+
+
\ No newline at end of file
diff --git a/.idea/vcs.xml b/.idea/vcs.xml
new file mode 100644
index 0000000..94a25f7
--- /dev/null
+++ b/.idea/vcs.xml
@@ -0,0 +1,6 @@
+
+
+
+
+
+
\ No newline at end of file
diff --git a/CPU_version/Forster2014ICRA.pdf b/CPU_version/Forster2014ICRA.pdf
new file mode 100755
index 0000000..9567178
Binary files /dev/null and b/CPU_version/Forster2014ICRA.pdf differ
diff --git a/CPU_version/LICENCE b/CPU_version/LICENCE
new file mode 100755
index 0000000..70fc6da
--- /dev/null
+++ b/CPU_version/LICENCE
@@ -0,0 +1,675 @@
+GNU GENERAL PUBLIC LICENSE
+ Version 3, 29 June 2007
+
+ Copyright (C) 2007 Free Software Foundation, Inc. {http://fsf.org/}
+ Everyone is permitted to copy and distribute verbatim copies
+ of this license document, but changing it is not allowed.
+
+ Preamble
+
+ The GNU General Public License is a free, copyleft license for
+software and other kinds of works.
+
+ The licenses for most software and other practical works are designed
+to take away your freedom to share and change the works. By contrast,
+the GNU General Public License is intended to guarantee your freedom to
+share and change all versions of a program--to make sure it remains free
+software for all its users. We, the Free Software Foundation, use the
+GNU General Public License for most of our software; it applies also to
+any other work released this way by its authors. You can apply it to
+your programs, too.
+
+ When we speak of free software, we are referring to freedom, not
+price. Our General Public Licenses are designed to make sure that you
+have the freedom to distribute copies of free software (and charge for
+them if you wish), that you receive source code or can get it if you
+want it, that you can change the software or use pieces of it in new
+free programs, and that you know you can do these things.
+
+ To protect your rights, we need to prevent others from denying you
+these rights or asking you to surrender the rights. Therefore, you have
+certain responsibilities if you distribute copies of the software, or if
+you modify it: responsibilities to respect the freedom of others.
+
+ For example, if you distribute copies of such a program, whether
+gratis or for a fee, you must pass on to the recipients the same
+freedoms that you received. You must make sure that they, too, receive
+or can get the source code. And you must show them these terms so they
+know their rights.
+
+ Developers that use the GNU GPL protect your rights with two steps:
+(1) assert copyright on the software, and (2) offer you this License
+giving you legal permission to copy, distribute and/or modify it.
+
+ For the developers' and authors' protection, the GPL clearly explains
+that there is no warranty for this free software. For both users' and
+authors' sake, the GPL requires that modified versions be marked as
+changed, so that their problems will not be attributed erroneously to
+authors of previous versions.
+
+ Some devices are designed to deny users access to install or run
+modified versions of the software inside them, although the manufacturer
+can do so. This is fundamentally incompatible with the aim of
+protecting users' freedom to change the software. The systematic
+pattern of such abuse occurs in the area of products for individuals to
+use, which is precisely where it is most unacceptable. Therefore, we
+have designed this version of the GPL to prohibit the practice for those
+products. If such problems arise substantially in other domains, we
+stand ready to extend this provision to those domains in future versions
+of the GPL, as needed to protect the freedom of users.
+
+ Finally, every program is threatened constantly by software patents.
+States should not allow patents to restrict development and use of
+software on general-purpose computers, but in those that do, we wish to
+avoid the special danger that patents applied to a free program could
+make it effectively proprietary. To prevent this, the GPL assures that
+patents cannot be used to render the program non-free.
+
+ The precise terms and conditions for copying, distribution and
+modification follow.
+
+ TERMS AND CONDITIONS
+
+ 0. Definitions.
+
+ "This License" refers to version 3 of the GNU General Public License.
+
+ "Copyright" also means copyright-like laws that apply to other kinds of
+works, such as semiconductor masks.
+
+ "The Program" refers to any copyrightable work licensed under this
+License. Each licensee is addressed as "you". "Licensees" and
+"recipients" may be individuals or organizations.
+
+ To "modify" a work means to copy from or adapt all or part of the work
+in a fashion requiring copyright permission, other than the making of an
+exact copy. The resulting work is called a "modified version" of the
+earlier work or a work "based on" the earlier work.
+
+ A "covered work" means either the unmodified Program or a work based
+on the Program.
+
+ To "propagate" a work means to do anything with it that, without
+permission, would make you directly or secondarily liable for
+infringement under applicable copyright law, except executing it on a
+computer or modifying a private copy. Propagation includes copying,
+distribution (with or without modification), making available to the
+public, and in some countries other activities as well.
+
+ To "convey" a work means any kind of propagation that enables other
+parties to make or receive copies. Mere interaction with a user through
+a computer network, with no transfer of a copy, is not conveying.
+
+ An interactive user interface displays "Appropriate Legal Notices"
+to the extent that it includes a convenient and prominently visible
+feature that (1) displays an appropriate copyright notice, and (2)
+tells the user that there is no warranty for the work (except to the
+extent that warranties are provided), that licensees may convey the
+work under this License, and how to view a copy of this License. If
+the interface presents a list of user commands or options, such as a
+menu, a prominent item in the list meets this criterion.
+
+ 1. Source Code.
+
+ The "source code" for a work means the preferred form of the work
+for making modifications to it. "Object code" means any non-source
+form of a work.
+
+ A "Standard Interface" means an interface that either is an official
+standard defined by a recognized standards body, or, in the case of
+interfaces specified for a particular programming language, one that
+is widely used among developers working in that language.
+
+ The "System Libraries" of an executable work include anything, other
+than the work as a whole, that (a) is included in the normal form of
+packaging a Major Component, but which is not part of that Major
+Component, and (b) serves only to enable use of the work with that
+Major Component, or to implement a Standard Interface for which an
+implementation is available to the public in source code form. A
+"Major Component", in this context, means a major essential component
+(kernel, window system, and so on) of the specific operating system
+(if any) on which the executable work runs, or a compiler used to
+produce the work, or an object code interpreter used to run it.
+
+ The "Corresponding Source" for a work in object code form means all
+the source code needed to generate, install, and (for an executable
+work) run the object code and to modify the work, including scripts to
+control those activities. However, it does not include the work's
+System Libraries, or general-purpose tools or generally available free
+programs which are used unmodified in performing those activities but
+which are not part of the work. For example, Corresponding Source
+includes interface definition files associated with source files for
+the work, and the source code for shared libraries and dynamically
+linked subprograms that the work is specifically designed to require,
+such as by intimate data communication or control flow between those
+subprograms and other parts of the work.
+
+ The Corresponding Source need not include anything that users
+can regenerate automatically from other parts of the Corresponding
+Source.
+
+ The Corresponding Source for a work in source code form is that
+same work.
+
+ 2. Basic Permissions.
+
+ All rights granted under this License are granted for the term of
+copyright on the Program, and are irrevocable provided the stated
+conditions are met. This License explicitly affirms your unlimited
+permission to run the unmodified Program. The output from running a
+covered work is covered by this License only if the output, given its
+content, constitutes a covered work. This License acknowledges your
+rights of fair use or other equivalent, as provided by copyright law.
+
+ You may make, run and propagate covered works that you do not
+convey, without conditions so long as your license otherwise remains
+in force. You may convey covered works to others for the sole purpose
+of having them make modifications exclusively for you, or provide you
+with facilities for running those works, provided that you comply with
+the terms of this License in conveying all material for which you do
+not control copyright. Those thus making or running the covered works
+for you must do so exclusively on your behalf, under your direction
+and control, on terms that prohibit them from making any copies of
+your copyrighted material outside their relationship with you.
+
+ Conveying under any other circumstances is permitted solely under
+the conditions stated below. Sublicensing is not allowed; section 10
+makes it unnecessary.
+
+ 3. Protecting Users' Legal Rights From Anti-Circumvention Law.
+
+ No covered work shall be deemed part of an effective technological
+measure under any applicable law fulfilling obligations under article
+11 of the WIPO copyright treaty adopted on 20 December 1996, or
+similar laws prohibiting or restricting circumvention of such
+measures.
+
+ When you convey a covered work, you waive any legal power to forbid
+circumvention of technological measures to the extent such circumvention
+is effected by exercising rights under this License with respect to
+the covered work, and you disclaim any intention to limit operation or
+modification of the work as a means of enforcing, against the work's
+users, your or third parties' legal rights to forbid circumvention of
+technological measures.
+
+ 4. Conveying Verbatim Copies.
+
+ You may convey verbatim copies of the Program's source code as you
+receive it, in any medium, provided that you conspicuously and
+appropriately publish on each copy an appropriate copyright notice;
+keep intact all notices stating that this License and any
+non-permissive terms added in accord with section 7 apply to the code;
+keep intact all notices of the absence of any warranty; and give all
+recipients a copy of this License along with the Program.
+
+ You may charge any price or no price for each copy that you convey,
+and you may offer support or warranty protection for a fee.
+
+ 5. Conveying Modified Source Versions.
+
+ You may convey a work based on the Program, or the modifications to
+produce it from the Program, in the form of source code under the
+terms of section 4, provided that you also meet all of these conditions:
+
+ a) The work must carry prominent notices stating that you modified
+ it, and giving a relevant date.
+
+ b) The work must carry prominent notices stating that it is
+ released under this License and any conditions added under section
+ 7. This requirement modifies the requirement in section 4 to
+ "keep intact all notices".
+
+ c) You must license the entire work, as a whole, under this
+ License to anyone who comes into possession of a copy. This
+ License will therefore apply, along with any applicable section 7
+ additional terms, to the whole of the work, and all its parts,
+ regardless of how they are packaged. This License gives no
+ permission to license the work in any other way, but it does not
+ invalidate such permission if you have separately received it.
+
+ d) If the work has interactive user interfaces, each must display
+ Appropriate Legal Notices; however, if the Program has interactive
+ interfaces that do not display Appropriate Legal Notices, your
+ work need not make them do so.
+
+ A compilation of a covered work with other separate and independent
+works, which are not by their nature extensions of the covered work,
+and which are not combined with it such as to form a larger program,
+in or on a volume of a storage or distribution medium, is called an
+"aggregate" if the compilation and its resulting copyright are not
+used to limit the access or legal rights of the compilation's users
+beyond what the individual works permit. Inclusion of a covered work
+in an aggregate does not cause this License to apply to the other
+parts of the aggregate.
+
+ 6. Conveying Non-Source Forms.
+
+ You may convey a covered work in object code form under the terms
+of sections 4 and 5, provided that you also convey the
+machine-readable Corresponding Source under the terms of this License,
+in one of these ways:
+
+ a) Convey the object code in, or embodied in, a physical product
+ (including a physical distribution medium), accompanied by the
+ Corresponding Source fixed on a durable physical medium
+ customarily used for software interchange.
+
+ b) Convey the object code in, or embodied in, a physical product
+ (including a physical distribution medium), accompanied by a
+ written offer, valid for at least three years and valid for as
+ long as you offer spare parts or customer support for that product
+ model, to give anyone who possesses the object code either (1) a
+ copy of the Corresponding Source for all the software in the
+ product that is covered by this License, on a durable physical
+ medium customarily used for software interchange, for a price no
+ more than your reasonable cost of physically performing this
+ conveying of source, or (2) access to copy the
+ Corresponding Source from a network server at no charge.
+
+ c) Convey individual copies of the object code with a copy of the
+ written offer to provide the Corresponding Source. This
+ alternative is allowed only occasionally and noncommercially, and
+ only if you received the object code with such an offer, in accord
+ with subsection 6b.
+
+ d) Convey the object code by offering access from a designated
+ place (gratis or for a charge), and offer equivalent access to the
+ Corresponding Source in the same way through the same place at no
+ further charge. You need not require recipients to copy the
+ Corresponding Source along with the object code. If the place to
+ copy the object code is a network server, the Corresponding Source
+ may be on a different server (operated by you or a third party)
+ that supports equivalent copying facilities, provided you maintain
+ clear directions next to the object code saying where to find the
+ Corresponding Source. Regardless of what server hosts the
+ Corresponding Source, you remain obligated to ensure that it is
+ available for as long as needed to satisfy these requirements.
+
+ e) Convey the object code using peer-to-peer transmission, provided
+ you inform other peers where the object code and Corresponding
+ Source of the work are being offered to the general public at no
+ charge under subsection 6d.
+
+ A separable portion of the object code, whose source code is excluded
+from the Corresponding Source as a System Library, need not be
+included in conveying the object code work.
+
+ A "User Product" is either (1) a "consumer product", which means any
+tangible personal property which is normally used for personal, family,
+or household purposes, or (2) anything designed or sold for incorporation
+into a dwelling. In determining whether a product is a consumer product,
+doubtful cases shall be resolved in favor of coverage. For a particular
+product received by a particular user, "normally used" refers to a
+typical or common use of that class of product, regardless of the status
+of the particular user or of the way in which the particular user
+actually uses, or expects or is expected to use, the product. A product
+is a consumer product regardless of whether the product has substantial
+commercial, industrial or non-consumer uses, unless such uses represent
+the only significant mode of use of the product.
+
+ "Installation Information" for a User Product means any methods,
+procedures, authorization keys, or other information required to install
+and execute modified versions of a covered work in that User Product from
+a modified version of its Corresponding Source. The information must
+suffice to ensure that the continued functioning of the modified object
+code is in no case prevented or interfered with solely because
+modification has been made.
+
+ If you convey an object code work under this section in, or with, or
+specifically for use in, a User Product, and the conveying occurs as
+part of a transaction in which the right of possession and use of the
+User Product is transferred to the recipient in perpetuity or for a
+fixed term (regardless of how the transaction is characterized), the
+Corresponding Source conveyed under this section must be accompanied
+by the Installation Information. But this requirement does not apply
+if neither you nor any third party retains the ability to install
+modified object code on the User Product (for example, the work has
+been installed in ROM).
+
+ The requirement to provide Installation Information does not include a
+requirement to continue to provide support service, warranty, or updates
+for a work that has been modified or installed by the recipient, or for
+the User Product in which it has been modified or installed. Access to a
+network may be denied when the modification itself materially and
+adversely affects the operation of the network or violates the rules and
+protocols for communication across the network.
+
+ Corresponding Source conveyed, and Installation Information provided,
+in accord with this section must be in a format that is publicly
+documented (and with an implementation available to the public in
+source code form), and must require no special password or key for
+unpacking, reading or copying.
+
+ 7. Additional Terms.
+
+ "Additional permissions" are terms that supplement the terms of this
+License by making exceptions from one or more of its conditions.
+Additional permissions that are applicable to the entire Program shall
+be treated as though they were included in this License, to the extent
+that they are valid under applicable law. If additional permissions
+apply only to part of the Program, that part may be used separately
+under those permissions, but the entire Program remains governed by
+this License without regard to the additional permissions.
+
+ When you convey a copy of a covered work, you may at your option
+remove any additional permissions from that copy, or from any part of
+it. (Additional permissions may be written to require their own
+removal in certain cases when you modify the work.) You may place
+additional permissions on material, added by you to a covered work,
+for which you have or can give appropriate copyright permission.
+
+ Notwithstanding any other provision of this License, for material you
+add to a covered work, you may (if authorized by the copyright holders of
+that material) supplement the terms of this License with terms:
+
+ a) Disclaiming warranty or limiting liability differently from the
+ terms of sections 15 and 16 of this License; or
+
+ b) Requiring preservation of specified reasonable legal notices or
+ author attributions in that material or in the Appropriate Legal
+ Notices displayed by works containing it; or
+
+ c) Prohibiting misrepresentation of the origin of that material, or
+ requiring that modified versions of such material be marked in
+ reasonable ways as different from the original version; or
+
+ d) Limiting the use for publicity purposes of names of licensors or
+ authors of the material; or
+
+ e) Declining to grant rights under trademark law for use of some
+ trade names, trademarks, or service marks; or
+
+ f) Requiring indemnification of licensors and authors of that
+ material by anyone who conveys the material (or modified versions of
+ it) with contractual assumptions of liability to the recipient, for
+ any liability that these contractual assumptions directly impose on
+ those licensors and authors.
+
+ All other non-permissive additional terms are considered "further
+restrictions" within the meaning of section 10. If the Program as you
+received it, or any part of it, contains a notice stating that it is
+governed by this License along with a term that is a further
+restriction, you may remove that term. If a license document contains
+a further restriction but permits relicensing or conveying under this
+License, you may add to a covered work material governed by the terms
+of that license document, provided that the further restriction does
+not survive such relicensing or conveying.
+
+ If you add terms to a covered work in accord with this section, you
+must place, in the relevant source files, a statement of the
+additional terms that apply to those files, or a notice indicating
+where to find the applicable terms.
+
+ Additional terms, permissive or non-permissive, may be stated in the
+form of a separately written license, or stated as exceptions;
+the above requirements apply either way.
+
+ 8. Termination.
+
+ You may not propagate or modify a covered work except as expressly
+provided under this License. Any attempt otherwise to propagate or
+modify it is void, and will automatically terminate your rights under
+this License (including any patent licenses granted under the third
+paragraph of section 11).
+
+ However, if you cease all violation of this License, then your
+license from a particular copyright holder is reinstated (a)
+provisionally, unless and until the copyright holder explicitly and
+finally terminates your license, and (b) permanently, if the copyright
+holder fails to notify you of the violation by some reasonable means
+prior to 60 days after the cessation.
+
+ Moreover, your license from a particular copyright holder is
+reinstated permanently if the copyright holder notifies you of the
+violation by some reasonable means, this is the first time you have
+received notice of violation of this License (for any work) from that
+copyright holder, and you cure the violation prior to 30 days after
+your receipt of the notice.
+
+ Termination of your rights under this section does not terminate the
+licenses of parties who have received copies or rights from you under
+this License. If your rights have been terminated and not permanently
+reinstated, you do not qualify to receive new licenses for the same
+material under section 10.
+
+ 9. Acceptance Not Required for Having Copies.
+
+ You are not required to accept this License in order to receive or
+run a copy of the Program. Ancillary propagation of a covered work
+occurring solely as a consequence of using peer-to-peer transmission
+to receive a copy likewise does not require acceptance. However,
+nothing other than this License grants you permission to propagate or
+modify any covered work. These actions infringe copyright if you do
+not accept this License. Therefore, by modifying or propagating a
+covered work, you indicate your acceptance of this License to do so.
+
+ 10. Automatic Licensing of Downstream Recipients.
+
+ Each time you convey a covered work, the recipient automatically
+receives a license from the original licensors, to run, modify and
+propagate that work, subject to this License. You are not responsible
+for enforcing compliance by third parties with this License.
+
+ An "entity transaction" is a transaction transferring control of an
+organization, or substantially all assets of one, or subdividing an
+organization, or merging organizations. If propagation of a covered
+work results from an entity transaction, each party to that
+transaction who receives a copy of the work also receives whatever
+licenses to the work the party's predecessor in interest had or could
+give under the previous paragraph, plus a right to possession of the
+Corresponding Source of the work from the predecessor in interest, if
+the predecessor has it or can get it with reasonable efforts.
+
+ You may not impose any further restrictions on the exercise of the
+rights granted or affirmed under this License. For example, you may
+not impose a license fee, royalty, or other charge for exercise of
+rights granted under this License, and you may not initiate litigation
+(including a cross-claim or counterclaim in a lawsuit) alleging that
+any patent claim is infringed by making, using, selling, offering for
+sale, or importing the Program or any portion of it.
+
+ 11. Patents.
+
+ A "contributor" is a copyright holder who authorizes use under this
+License of the Program or a work on which the Program is based. The
+work thus licensed is called the contributor's "contributor version".
+
+ A contributor's "essential patent claims" are all patent claims
+owned or controlled by the contributor, whether already acquired or
+hereafter acquired, that would be infringed by some manner, permitted
+by this License, of making, using, or selling its contributor version,
+but do not include claims that would be infringed only as a
+consequence of further modification of the contributor version. For
+purposes of this definition, "control" includes the right to grant
+patent sublicenses in a manner consistent with the requirements of
+this License.
+
+ Each contributor grants you a non-exclusive, worldwide, royalty-free
+patent license under the contributor's essential patent claims, to
+make, use, sell, offer for sale, import and otherwise run, modify and
+propagate the contents of its contributor version.
+
+ In the following three paragraphs, a "patent license" is any express
+agreement or commitment, however denominated, not to enforce a patent
+(such as an express permission to practice a patent or covenant not to
+sue for patent infringement). To "grant" such a patent license to a
+party means to make such an agreement or commitment not to enforce a
+patent against the party.
+
+ If you convey a covered work, knowingly relying on a patent license,
+and the Corresponding Source of the work is not available for anyone
+to copy, free of charge and under the terms of this License, through a
+publicly available network server or other readily accessible means,
+then you must either (1) cause the Corresponding Source to be so
+available, or (2) arrange to deprive yourself of the benefit of the
+patent license for this particular work, or (3) arrange, in a manner
+consistent with the requirements of this License, to extend the patent
+license to downstream recipients. "Knowingly relying" means you have
+actual knowledge that, but for the patent license, your conveying the
+covered work in a country, or your recipient's use of the covered work
+in a country, would infringe one or more identifiable patents in that
+country that you have reason to believe are valid.
+
+ If, pursuant to or in connection with a single transaction or
+arrangement, you convey, or propagate by procuring conveyance of, a
+covered work, and grant a patent license to some of the parties
+receiving the covered work authorizing them to use, propagate, modify
+or convey a specific copy of the covered work, then the patent license
+you grant is automatically extended to all recipients of the covered
+work and works based on it.
+
+ A patent license is "discriminatory" if it does not include within
+the scope of its coverage, prohibits the exercise of, or is
+conditioned on the non-exercise of one or more of the rights that are
+specifically granted under this License. You may not convey a covered
+work if you are a party to an arrangement with a third party that is
+in the business of distributing software, under which you make payment
+to the third party based on the extent of your activity of conveying
+the work, and under which the third party grants, to any of the
+parties who would receive the covered work from you, a discriminatory
+patent license (a) in connection with copies of the covered work
+conveyed by you (or copies made from those copies), or (b) primarily
+for and in connection with specific products or compilations that
+contain the covered work, unless you entered into that arrangement,
+or that patent license was granted, prior to 28 March 2007.
+
+ Nothing in this License shall be construed as excluding or limiting
+any implied license or other defenses to infringement that may
+otherwise be available to you under applicable patent law.
+
+ 12. No Surrender of Others' Freedom.
+
+ If conditions are imposed on you (whether by court order, agreement or
+otherwise) that contradict the conditions of this License, they do not
+excuse you from the conditions of this License. If you cannot convey a
+covered work so as to satisfy simultaneously your obligations under this
+License and any other pertinent obligations, then as a consequence you may
+not convey it at all. For example, if you agree to terms that obligate you
+to collect a royalty for further conveying from those to whom you convey
+the Program, the only way you could satisfy both those terms and this
+License would be to refrain entirely from conveying the Program.
+
+ 13. Use with the GNU Affero General Public License.
+
+ Notwithstanding any other provision of this License, you have
+permission to link or combine any covered work with a work licensed
+under version 3 of the GNU Affero General Public License into a single
+combined work, and to convey the resulting work. The terms of this
+License will continue to apply to the part which is the covered work,
+but the special requirements of the GNU Affero General Public License,
+section 13, concerning interaction through a network will apply to the
+combination as such.
+
+ 14. Revised Versions of this License.
+
+ The Free Software Foundation may publish revised and/or new versions of
+the GNU General Public License from time to time. Such new versions will
+be similar in spirit to the present version, but may differ in detail to
+address new problems or concerns.
+
+ Each version is given a distinguishing version number. If the
+Program specifies that a certain numbered version of the GNU General
+Public License "or any later version" applies to it, you have the
+option of following the terms and conditions either of that numbered
+version or of any later version published by the Free Software
+Foundation. If the Program does not specify a version number of the
+GNU General Public License, you may choose any version ever published
+by the Free Software Foundation.
+
+ If the Program specifies that a proxy can decide which future
+versions of the GNU General Public License can be used, that proxy's
+public statement of acceptance of a version permanently authorizes you
+to choose that version for the Program.
+
+ Later license versions may give you additional or different
+permissions. However, no additional obligations are imposed on any
+author or copyright holder as a result of your choosing to follow a
+later version.
+
+ 15. Disclaimer of Warranty.
+
+ THERE IS NO WARRANTY FOR THE PROGRAM, TO THE EXTENT PERMITTED BY
+APPLICABLE LAW. EXCEPT WHEN OTHERWISE STATED IN WRITING THE COPYRIGHT
+HOLDERS AND/OR OTHER PARTIES PROVIDE THE PROGRAM "AS IS" WITHOUT WARRANTY
+OF ANY KIND, EITHER EXPRESSED OR IMPLIED, INCLUDING, BUT NOT LIMITED TO,
+THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
+PURPOSE. THE ENTIRE RISK AS TO THE QUALITY AND PERFORMANCE OF THE PROGRAM
+IS WITH YOU. SHOULD THE PROGRAM PROVE DEFECTIVE, YOU ASSUME THE COST OF
+ALL NECESSARY SERVICING, REPAIR OR CORRECTION.
+
+ 16. Limitation of Liability.
+
+ IN NO EVENT UNLESS REQUIRED BY APPLICABLE LAW OR AGREED TO IN WRITING
+WILL ANY COPYRIGHT HOLDER, OR ANY OTHER PARTY WHO MODIFIES AND/OR CONVEYS
+THE PROGRAM AS PERMITTED ABOVE, BE LIABLE TO YOU FOR DAMAGES, INCLUDING ANY
+GENERAL, SPECIAL, INCIDENTAL OR CONSEQUENTIAL DAMAGES ARISING OUT OF THE
+USE OR INABILITY TO USE THE PROGRAM (INCLUDING BUT NOT LIMITED TO LOSS OF
+DATA OR DATA BEING RENDERED INACCURATE OR LOSSES SUSTAINED BY YOU OR THIRD
+PARTIES OR A FAILURE OF THE PROGRAM TO OPERATE WITH ANY OTHER PROGRAMS),
+EVEN IF SUCH HOLDER OR OTHER PARTY HAS BEEN ADVISED OF THE POSSIBILITY OF
+SUCH DAMAGES.
+
+ 17. Interpretation of Sections 15 and 16.
+
+ If the disclaimer of warranty and limitation of liability provided
+above cannot be given local legal effect according to their terms,
+reviewing courts shall apply local law that most closely approximates
+an absolute waiver of all civil liability in connection with the
+Program, unless a warranty or assumption of liability accompanies a
+copy of the Program in return for a fee.
+
+ END OF TERMS AND CONDITIONS
+
+ How to Apply These Terms to Your New Programs
+
+ If you develop a new program, and you want it to be of the greatest
+possible use to the public, the best way to achieve this is to make it
+free software which everyone can redistribute and change under these terms.
+
+ To do so, attach the following notices to the program. It is safest
+to attach them to the start of each source file to most effectively
+state the exclusion of warranty; and each file should have at least
+the "copyright" line and a pointer to where the full notice is found.
+
+ {one line to give the program's name and a brief idea of what it does.}
+ Copyright (C) {year} {name of author}
+
+ This program is free software: you can redistribute it and/or modify
+ it under the terms of the GNU General Public License as published by
+ the Free Software Foundation, either version 3 of the License, or
+ (at your option) any later version.
+
+ This program is distributed in the hope that it will be useful,
+ but WITHOUT ANY WARRANTY; without even the implied warranty of
+ MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+ GNU General Public License for more details.
+
+ You should have received a copy of the GNU General Public License
+ along with this program. If not, see {http://www.gnu.org/licenses/}.
+
+Also add information on how to contact you by electronic and paper mail.
+
+ If the program does terminal interaction, make it output a short
+notice like this when it starts in an interactive mode:
+
+ PTAM-GPL Copyright (C) 2013 Oxford-PTAM
+ This program comes with ABSOLUTELY NO WARRANTY; for details type `show w'.
+ This is free software, and you are welcome to redistribute it
+ under certain conditions; type `show c' for details.
+
+The hypothetical commands `show w' and `show c' should show the appropriate
+parts of the General Public License. Of course, your program's commands
+might be different; for a GUI interface, you would use an "about box".
+
+ You should also get your employer (if you work as a programmer) or school,
+if any, to sign a "copyright disclaimer" for the program, if necessary.
+For more information on this, and how to apply and follow the GNU GPL, see
+{http://www.gnu.org/licenses/}.
+
+ The GNU General Public License does not permit incorporating your program
+into proprietary programs. If your program is a subroutine library, you
+may consider it more useful to permit linking proprietary applications with
+the library. If this is what you want to do, use the GNU Lesser General
+Public License instead of this License. But first, please read
+{http://www.gnu.org/philosophy/why-not-lgpl.html}.
+
diff --git a/CPU_version/README.md b/CPU_version/README.md
new file mode 100755
index 0000000..8387b0e
--- /dev/null
+++ b/CPU_version/README.md
@@ -0,0 +1,64 @@
+SVO
+===
+
+This code implements a semi-direct monocular visual odometry pipeline which has been integrated with GTSAM to infused IMU as well.
+
+
+Paper svo: http://rpg.ifi.uzh.ch/docs/ICRA14_Forster.pdf
+Paper Svo+GTSAM: https://arxiv.org/abs/1512.02363
+
+
+#### Disclaimer
+
+SVO has been tested under ROS Groovy, Hydro and Indigo with Ubuntu 12.04, 13.04 and 14.04. This is research code, any fitness for a particular purpose is disclaimed.
+
+
+#### Licence
+
+The source code is released under a GPLv3 licence. A closed-source professional edition is available for commercial purposes. In this case, please contact the authors for further info.
+
+
+#### Citing
+
+If you use SVO in an academic context, please cite the following publication:
+
+ @inproceedings{Forster2014ICRA,
+ author = {Forster, Christian and Pizzoli, Matia and Scaramuzza, Davide},
+ title = {{SVO}: Fast Semi-Direct Monocular Visual Odometry},
+ booktitle = {IEEE International Conference on Robotics and Automation (ICRA)},
+ year = {2014}
+ }
+
+
+#### Documentation
+
+The API is documented here: http://uzh-rpg.github.io/rpg_svo/doc/
+The GTSAM documentation https://gtsam.org/tutorials/
+
+#### Instructions
+
+First install gtsam library => https://github.com/borglab/gtsam
+
+Then follow the SVO instalation
+
+See the Wiki for more instructions. https://github.com/uzh-rpg/rpg_svo/wiki and https://gtsam.org/get_started/
+
+#### Contributing
+
+You are very welcome to contribute to SVO by opening a pull request via Github.
+I try to follow the ROS C++ style guide http://wiki.ros.org/CppStyleGuide
+
+## Setting
+
+You can edit the camera position in the launch file as (The IMU frame will rotate to align with camera frame based on the follwoing Eular angles):
+CPU_version/vio_svo/launch/Px30.launch
+
+
+
+
+
+
+You can also edit the algorithm parameters in the follwoing file, please look at the config file to underestand the meaning of the parameters (CPU_version/svo/include/svo/config.h)
+
+CPU_version/vio_svo/param/vo_fast.yaml
+
diff --git a/CPU_version/__init__.py b/CPU_version/__init__.py
new file mode 100644
index 0000000..e69de29
diff --git a/CPU_version/svo/.idea/.gitignore b/CPU_version/svo/.idea/.gitignore
new file mode 100644
index 0000000..73f69e0
--- /dev/null
+++ b/CPU_version/svo/.idea/.gitignore
@@ -0,0 +1,8 @@
+# Default ignored files
+/shelf/
+/workspace.xml
+# Datasource local storage ignored files
+/dataSources/
+/dataSources.local.xml
+# Editor-based HTTP Client requests
+/httpRequests/
diff --git a/CPU_version/svo/.idea/misc.xml b/CPU_version/svo/.idea/misc.xml
new file mode 100644
index 0000000..79b3c94
--- /dev/null
+++ b/CPU_version/svo/.idea/misc.xml
@@ -0,0 +1,4 @@
+
+
+
+
\ No newline at end of file
diff --git a/CPU_version/svo/.idea/modules.xml b/CPU_version/svo/.idea/modules.xml
new file mode 100644
index 0000000..4c8cf2b
--- /dev/null
+++ b/CPU_version/svo/.idea/modules.xml
@@ -0,0 +1,8 @@
+
+
+
+
+
+
+
+
\ No newline at end of file
diff --git a/CPU_version/svo/.idea/svo.iml b/CPU_version/svo/.idea/svo.iml
new file mode 100644
index 0000000..f08604b
--- /dev/null
+++ b/CPU_version/svo/.idea/svo.iml
@@ -0,0 +1,2 @@
+
+
\ No newline at end of file
diff --git a/CPU_version/svo/.idea/vcs.xml b/CPU_version/svo/.idea/vcs.xml
new file mode 100644
index 0000000..6c0b863
--- /dev/null
+++ b/CPU_version/svo/.idea/vcs.xml
@@ -0,0 +1,6 @@
+
+
+
+
+
+
\ No newline at end of file
diff --git a/CPU_version/svo/CMakeLists.txt b/CPU_version/svo/CMakeLists.txt
new file mode 100644
index 0000000..c3ace90
--- /dev/null
+++ b/CPU_version/svo/CMakeLists.txt
@@ -0,0 +1,156 @@
+################################################################################
+# user build settings
+
+SET(TRACE TRUE)
+SET(HAVE_G2O FALSE)
+SET(USE_ROS TRUE) # Set FALSE if you want to build this package without Catkin
+SET(DEBUG_OUTPUT FALSE) # Only relevant if build without ROS
+
+################################################################################
+
+SET(PROJECT_NAME svo)
+PROJECT(${PROJECT_NAME})
+CMAKE_MINIMUM_REQUIRED (VERSION 2.8.3)
+SET(CMAKE_BUILD_TYPE Release) # Release, RelWithDebInfo
+SET(CMAKE_VERBOSE_MAKEFILE OFF)
+SET(CMAKE_MODULE_PATH ${CMAKE_MODULE_PATH} "${PROJECT_SOURCE_DIR}/CMakeModules/")
+set(CMAKE_C_STANDARD 14)
+set(CMAKE_CXX_STANDARD 14)
+
+# Set definitions
+IF(TRACE)
+ ADD_DEFINITIONS(-DSVO_TRACE)
+ENDIF()
+IF(HAVE_G2O)
+ ADD_DEFINITIONS(-DUSE_BUNDLE_ADJUSTMENT)
+ENDIF()
+IF(USE_ROS)
+ ADD_DEFINITIONS(-DSVO_USE_ROS)
+ENDIF()
+IF(DEBUG_OUTPUT)
+ ADD_DEFINITIONS(-DSVO_DEBUG_OUTPUT)
+ENDIF()
+
+# Set build flags, set ARM_ARCHITECTURE environment variable on Odroid
+# Set build flags. Set IS_ARM on odroid board as environment variable
+SET(CMAKE_CXX_FLAGS "-Wall -D_LINUX -D_REENTRANT -march=native -Wno-unused-variable -Wno-unused-but-set-variable -Wno-unknown-pragmas")
+IF(DEFINED ENV{ARM_ARCHITECTURE})
+ SET(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} ")
+ELSE()
+ SET(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} -mmmx -msse -msse -msse2 -msse3 -mssse3")
+ENDIF()
+IF(CMAKE_COMPILER_IS_GNUCC)
+ SET(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} -std=c++0x")
+ELSE()
+ SET(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} -std=c++11")
+ENDIF()
+SET(CMAKE_CXX_FLAGS_RELEASE "${CMAKE_CXX_FLAGS} -O3 -fsee -fomit-frame-pointer -fno-signed-zeros -fno-math-errno -funroll-loops")
+
+# Add plain cmake packages
+FIND_PACKAGE(OpenCV 3 REQUIRED)
+FIND_PACKAGE(Eigen REQUIRED)
+FIND_PACKAGE(Sophus REQUIRED)
+FIND_PACKAGE(fast REQUIRED)
+FIND_PACKAGE(Boost REQUIRED COMPONENTS thread system)
+FIND_PACKAGE(GTSAM REQUIRED)
+
+set(GTSAM_INCLUDE_DIRS "${GTSAM_INCLUDE_DIR}")
+IF(USE_ROS)
+ FIND_PACKAGE(catkin REQUIRED COMPONENTS roscpp roslib cmake_modules vikit_common vikit_ros)
+ catkin_package(
+ DEPENDS Eigen OpenCV Sophus Boost fast GTSAM
+ CATKIN_DEPENDS roscpp roslib vikit_common vikit_ros
+ INCLUDE_DIRS include
+ LIBRARIES svo
+ )
+ELSE()
+ FIND_PACKAGE(vikit_common REQUIRED)
+ SET(EXECUTABLE_OUTPUT_PATH ${PROJECT_SOURCE_DIR}/bin)
+ SET(LIBRARY_OUTPUT_PATH ${PROJECT_SOURCE_DIR}/lib)
+ENDIF()
+
+# Include dirs
+INCLUDE_DIRECTORIES(
+ include
+ ${Eigen_INCLUDE_DIRS}
+ ${OpenCV_INCLUDE_DIRS}
+ ${Sophus_INCLUDE_DIRS}
+ ${Boost_INCLUDE_DIRS}
+ ${fast_INCLUDE_DIRS} ${GTSAM_INCLUDE_DIR}
+)
+
+# Set link libraries
+LIST(APPEND LINK_LIBS
+ ${OpenCV_LIBS}
+ ${Sophus_LIBRARIES}
+ ${Boost_LIBRARIES}
+ ${fast_LIBRARIES}
+)
+IF(USE_ROS)
+ INCLUDE_DIRECTORIES(${catkin_INCLUDE_DIRS})
+ LIST(APPEND LINK_LIBS ${catkin_LIBRARIES})
+ELSE()
+ INCLUDE_DIRECTORIES(${vikit_common_INCLUDE_DIRS})
+ LIST(APPEND LINK_LIBS ${vikit_common_LIBRARIES})
+ENDIF()
+
+# Set sourcefiles
+LIST(APPEND SOURCEFILES
+ src/frame_handler_mono.cpp
+ src/frame_handler_base.cpp
+ src/frame.cpp
+ src/point.cpp
+ src/map.cpp
+ src/pose_optimizer.cpp
+ src/initialization.cpp
+ src/matcher.cpp
+ src/reprojector.cpp
+ src/feature_alignment.cpp
+ src/feature_detection.cpp
+ src/depth_filter.cpp
+ src/config.cpp
+ src/sparse_img_align.cpp src/imu_integration.cpp)
+
+# Add g2o if available
+IF(HAVE_G2O)
+ FIND_PACKAGE(G2O REQUIRED)
+ INCLUDE_DIRECTORIES(
+ ${G2O_INCLUDE_DIR}
+ /usr/include/suitesparse # for cholmod
+ )
+ LIST(APPEND LINK_LIBS
+ ${G2O_CORE_LIBRARY} ${G2O_STUFF_LIBRARY} ${G2O_SOLVER_CHOLMOD} ${G2O_SOLVER_CSPARSE} ${G2O_SOLVER_DENSE}
+ ${G2O_SOLVER_PCG} ${G2O_TYPES_SBA} cholmod cxsparse
+ )
+ LIST(APPEND SOURCEFILES src/bundle_adjustment.cpp)
+ENDIF(HAVE_G2O)
+
+# Create svo library
+ADD_LIBRARY(svo SHARED ${SOURCEFILES})
+TARGET_LINK_LIBRARIES(svo ${LINK_LIBS} gtsam)
+set_property(TARGET ${PROJECT_NAME} PROPERTY CXX_STANDARD 14)
+set_property(TARGET ${PROJECT_NAME} PROPERTY CXX_STANDARD_REQUIRED ON)
+
+################################################################################
+# TESTS
+
+#ADD_EXECUTABLE(test_feature_align test/test_feature_alignment.cpp)
+#TARGET_LINK_LIBRARIES(test_feature_align svo)
+
+#ADD_EXECUTABLE(test_pipeline test/test_pipeline.cpp)
+#TARGET_LINK_LIBRARIES(test_pipeline svo)
+
+#ADD_EXECUTABLE(test_matcher test/test_matcher.cpp)
+#TARGET_LINK_LIBRARIES(test_matcher svo)
+
+#ADD_EXECUTABLE(test_feature_detection test/test_feature_detection.cpp)
+#TARGET_LINK_LIBRARIES(test_feature_detection svo)
+
+#ADD_EXECUTABLE(test_depth_filter test/test_depth_filter.cpp)
+#TARGET_LINK_LIBRARIES(test_depth_filter svo)
+
+#ADD_EXECUTABLE(test_sparse_img_align test/test_sparse_img_align.cpp)
+#TARGET_LINK_LIBRARIES(test_sparse_img_align svo)
+
+#ADD_EXECUTABLE(test_pose_optimizer test/test_pose_optimizer.cpp)
+#TARGET_LINK_LIBRARIES(test_pose_optimizer svo)
diff --git a/CPU_version/svo/CMakeModules/FindEigen.cmake b/CPU_version/svo/CMakeModules/FindEigen.cmake
new file mode 100644
index 0000000..8587367
--- /dev/null
+++ b/CPU_version/svo/CMakeModules/FindEigen.cmake
@@ -0,0 +1,81 @@
+###############################################################################
+#
+# CMake script for finding the Eigen library.
+#
+# http://eigen.tuxfamily.org/index.php?title=Main_Page
+#
+# Copyright (c) 2006, 2007 Montel Laurent,
+# Copyright (c) 2008, 2009 Gael Guennebaud,
+# Copyright (c) 2009 Benoit Jacob
+# Redistribution and use is allowed according to the terms of the 2-clause BSD
+# license.
+#
+#
+# Input variables:
+#
+# - Eigen_ROOT_DIR (optional): When specified, header files and libraries
+# will be searched for in `${Eigen_ROOT_DIR}/include` and
+# `${Eigen_ROOT_DIR}/libs` respectively, and the default CMake search order
+# will be ignored. When unspecified, the default CMake search order is used.
+# This variable can be specified either as a CMake or environment variable.
+# If both are set, preference is given to the CMake variable.
+# Use this variable for finding packages installed in a nonstandard location,
+# or for enforcing that one of multiple package installations is picked up.
+#
+# Cache variables (not intended to be used in CMakeLists.txt files)
+#
+# - Eigen_INCLUDE_DIR: Absolute path to package headers.
+#
+#
+# Output variables:
+#
+# - Eigen_FOUND: Boolean that indicates if the package was found
+# - Eigen_INCLUDE_DIRS: Paths to the necessary header files
+# - Eigen_VERSION: Version of Eigen library found
+# - Eigen_DEFINITIONS: Definitions to be passed on behalf of eigen
+#
+#
+# Example usage:
+#
+# # Passing the version means Eigen_FOUND will only be TRUE if a
+# # version >= the provided version is found.
+# find_package(Eigen 3.1.2)
+# if(NOT Eigen_FOUND)
+# # Error handling
+# endif()
+# ...
+# add_definitions(${Eigen_DEFINITIONS})
+# ...
+# include_directories(${Eigen_INCLUDE_DIRS} ...)
+#
+###############################################################################
+
+find_package(PkgConfig)
+pkg_check_modules(PC_EIGEN eigen3)
+set(EIGEN_DEFINITIONS ${PC_EIGEN_CFLAGS_OTHER})
+
+
+find_path(EIGEN_INCLUDE_DIR Eigen/Core
+ HINTS ${PC_EIGEN_INCLUDEDIR} ${PC_EIGEN_INCLUDE_DIRS}
+ "${Eigen_ROOT_DIR}" "$ENV{EIGEN_ROOT_DIR}"
+ "${EIGEN_ROOT}" "$ENV{EIGEN_ROOT}" # Backwards Compatibility
+ PATHS "$ENV{PROGRAMFILES}/Eigen" "$ENV{PROGRAMW6432}/Eigen"
+ "$ENV{PROGRAMFILES}/Eigen 3.0.0" "$ENV{PROGRAMW6432}/Eigen 3.0.0"
+ PATH_SUFFIXES eigen3 include/eigen3 include)
+
+set(EIGEN_INCLUDE_DIRS ${EIGEN_INCLUDE_DIR})
+
+include(FindPackageHandleStandardArgs)
+find_package_handle_standard_args(Eigen DEFAULT_MSG EIGEN_INCLUDE_DIR)
+
+mark_as_advanced(EIGEN_INCLUDE_DIR)
+
+if(EIGEN_FOUND)
+ message(STATUS "Eigen found (include: ${EIGEN_INCLUDE_DIRS})")
+endif(EIGEN_FOUND)
+
+
+set(Eigen_INCLUDE_DIRS ${EIGEN_INCLUDE_DIRS})
+set(Eigen_FOUND ${EIGEN_FOUND})
+set(Eigen_VERSION ${EIGEN_VERSION})
+set(Eigen_DEFINITIONS ${EIGEN_DEFINITIONS})
diff --git a/CPU_version/svo/CMakeModules/FindG2O.cmake b/CPU_version/svo/CMakeModules/FindG2O.cmake
new file mode 100644
index 0000000..655600f
--- /dev/null
+++ b/CPU_version/svo/CMakeModules/FindG2O.cmake
@@ -0,0 +1,113 @@
+# Find the header files
+
+FIND_PATH(G2O_INCLUDE_DIR g2o/core/base_vertex.h
+ ${G2O_ROOT}/include
+ $ENV{G2O_ROOT}/include
+ $ENV{G2O_ROOT}
+ /usr/local/include
+ /usr/include
+ /opt/local/include
+ /sw/local/include
+ /sw/include
+ NO_DEFAULT_PATH
+ )
+
+# Macro to unify finding both the debug and release versions of the
+# libraries; this is adapted from the OpenSceneGraph FIND_LIBRARY
+# macro.
+
+MACRO(FIND_G2O_LIBRARY MYLIBRARY MYLIBRARYNAME)
+
+ FIND_LIBRARY("${MYLIBRARY}_DEBUG"
+ NAMES "g2o_${MYLIBRARYNAME}_d"
+ PATHS
+ ${G2O_ROOT}/lib/Debug
+ ${G2O_ROOT}/lib
+ $ENV{G2O_ROOT}/lib/Debug
+ $ENV{G2O_ROOT}/lib
+ NO_DEFAULT_PATH
+ )
+
+ FIND_LIBRARY("${MYLIBRARY}_DEBUG"
+ NAMES "g2o_${MYLIBRARYNAME}_d"
+ PATHS
+ ~/Library/Frameworks
+ /Library/Frameworks
+ /usr/local/lib
+ /usr/local/lib64
+ /usr/lib
+ /usr/lib64
+ /opt/local/lib
+ /sw/local/lib
+ /sw/lib
+ )
+
+ FIND_LIBRARY(${MYLIBRARY}
+ NAMES "g2o_${MYLIBRARYNAME}"
+ PATHS
+ ${G2O_ROOT}/lib/Release
+ ${G2O_ROOT}/lib
+ $ENV{G2O_ROOT}/lib/Release
+ $ENV{G2O_ROOT}/lib
+ NO_DEFAULT_PATH
+ )
+
+ FIND_LIBRARY(${MYLIBRARY}
+ NAMES "g2o_${MYLIBRARYNAME}"
+ PATHS
+ ~/Library/Frameworks
+ /Library/Frameworks
+ /usr/local/lib
+ /usr/local/lib64
+ /usr/lib
+ /usr/lib64
+ /opt/local/lib
+ /sw/local/lib
+ /sw/lib
+ )
+
+ IF(NOT ${MYLIBRARY}_DEBUG)
+ IF(MYLIBRARY)
+ SET(${MYLIBRARY}_DEBUG ${MYLIBRARY})
+ ENDIF(MYLIBRARY)
+ ENDIF( NOT ${MYLIBRARY}_DEBUG)
+
+ENDMACRO(FIND_G2O_LIBRARY LIBRARY LIBRARYNAME)
+
+# Find the core elements
+FIND_G2O_LIBRARY(G2O_STUFF_LIBRARY stuff)
+FIND_G2O_LIBRARY(G2O_CORE_LIBRARY core)
+
+# Find the CLI library
+FIND_G2O_LIBRARY(G2O_CLI_LIBRARY cli)
+
+# Find the pluggable solvers
+FIND_G2O_LIBRARY(G2O_SOLVER_CHOLMOD solver_cholmod)
+FIND_G2O_LIBRARY(G2O_SOLVER_CSPARSE solver_csparse)
+FIND_G2O_LIBRARY(G2O_SOLVER_CSPARSE_EXTENSION csparse_extension)
+FIND_G2O_LIBRARY(G2O_SOLVER_DENSE solver_dense)
+FIND_G2O_LIBRARY(G2O_SOLVER_PCG solver_pcg)
+FIND_G2O_LIBRARY(G2O_SOLVER_SLAM2D_LINEAR solver_slam2d_linear)
+FIND_G2O_LIBRARY(G2O_SOLVER_STRUCTURE_ONLY solver_structure_only)
+FIND_G2O_LIBRARY(G2O_SOLVER_EIGEN solver_eigen)
+
+# Find the predefined types
+FIND_G2O_LIBRARY(G2O_TYPES_DATA types_data)
+FIND_G2O_LIBRARY(G2O_TYPES_ICP types_icp)
+FIND_G2O_LIBRARY(G2O_TYPES_SBA types_sba)
+FIND_G2O_LIBRARY(G2O_TYPES_SCLAM2D types_sclam2d)
+FIND_G2O_LIBRARY(G2O_TYPES_SIM3 types_sim3)
+FIND_G2O_LIBRARY(G2O_TYPES_SLAM2D types_slam2d)
+FIND_G2O_LIBRARY(G2O_TYPES_SLAM3D types_slam3d)
+
+# G2O solvers declared found if we found at least one solver
+SET(G2O_SOLVERS_FOUND "NO")
+IF(G2O_SOLVER_CHOLMOD OR G2O_SOLVER_CSPARSE OR G2O_SOLVER_DENSE OR G2O_SOLVER_PCG OR G2O_SOLVER_SLAM2D_LINEAR OR G2O_SOLVER_STRUCTURE_ONLY OR G2O_SOLVER_EIGEN)
+ SET(G2O_SOLVERS_FOUND "YES")
+ENDIF(G2O_SOLVER_CHOLMOD OR G2O_SOLVER_CSPARSE OR G2O_SOLVER_DENSE OR G2O_SOLVER_PCG OR G2O_SOLVER_SLAM2D_LINEAR OR G2O_SOLVER_STRUCTURE_ONLY OR G2O_SOLVER_EIGEN)
+
+# G2O itself declared found if we found the core libraries and at least one solver
+SET(G2O_FOUND "NO")
+IF(G2O_STUFF_LIBRARY AND G2O_CORE_LIBRARY AND G2O_INCLUDE_DIR AND G2O_SOLVERS_FOUND)
+ SET(G2O_FOUND "YES")
+ENDIF(G2O_STUFF_LIBRARY AND G2O_CORE_LIBRARY AND G2O_INCLUDE_DIR AND G2O_SOLVERS_FOUND)
diff --git a/CPU_version/svo/Doxyfile b/CPU_version/svo/Doxyfile
new file mode 100644
index 0000000..b1b4256
--- /dev/null
+++ b/CPU_version/svo/Doxyfile
@@ -0,0 +1,1782 @@
+# Doxyfile 1.7.6.1
+
+# This file describes the settings to be used by the documentation system
+# doxygen (www.doxygen.org) for a project.
+#
+# All text after a hash (#) is considered a comment and will be ignored.
+# The format is:
+# TAG = value [value, ...]
+# For lists items can also be appended using:
+# TAG += value [value, ...]
+# Values that contain spaces should be placed between quotes (" ").
+
+#---------------------------------------------------------------------------
+# Project related configuration options
+#---------------------------------------------------------------------------
+
+# This tag specifies the encoding used for all characters in the config file
+# that follow. The default is UTF-8 which is also the encoding used for all
+# text before the first occurrence of this tag. Doxygen uses libiconv (or the
+# iconv built into libc) for the transcoding. See
+# http://www.gnu.org/software/libiconv for the list of possible encodings.
+
+DOXYFILE_ENCODING = UTF-8
+
+# The PROJECT_NAME tag is a single word (or sequence of words) that should
+# identify the project. Note that if you do not use Doxywizard you need
+# to put quotes around the project name if it contains spaces.
+
+PROJECT_NAME = "svo"
+
+# The PROJECT_NUMBER tag can be used to enter a project or revision number.
+# This could be handy for archiving the generated documentation or
+# if some version control system is used.
+
+PROJECT_NUMBER =
+
+# Using the PROJECT_BRIEF tag one can provide an optional one line description
+# for a project that appears at the top of each page and should give viewer
+# a quick idea about the purpose of the project. Keep the description short.
+
+PROJECT_BRIEF = "Semi-Direct Visual Odometry"
+
+# With the PROJECT_LOGO tag one can specify an logo or icon that is
+# included in the documentation. The maximum height of the logo should not
+# exceed 55 pixels and the maximum width should not exceed 200 pixels.
+# Doxygen will copy the logo to the output directory.
+
+PROJECT_LOGO =
+
+# The OUTPUT_DIRECTORY tag is used to specify the (relative or absolute)
+# base path where the generated documentation will be put.
+# If a relative path is entered, it will be relative to the location
+# where doxygen was started. If left blank the current directory will be used.
+
+OUTPUT_DIRECTORY = "doc"
+
+# If the CREATE_SUBDIRS tag is set to YES, then doxygen will create
+# 4096 sub-directories (in 2 levels) under the output directory of each output
+# format and will distribute the generated files over these directories.
+# Enabling this option can be useful when feeding doxygen a huge amount of
+# source files, where putting all generated files in the same directory would
+# otherwise cause performance problems for the file system.
+
+CREATE_SUBDIRS = NO
+
+# The OUTPUT_LANGUAGE tag is used to specify the language in which all
+# documentation generated by doxygen is written. Doxygen will use this
+# information to generate all constant output in the proper language.
+# The default language is English, other supported languages are:
+# Afrikaans, Arabic, Brazilian, Catalan, Chinese, Chinese-Traditional,
+# Croatian, Czech, Danish, Dutch, Esperanto, Farsi, Finnish, French, German,
+# Greek, Hungarian, Italian, Japanese, Japanese-en (Japanese with English
+# messages), Korean, Korean-en, Lithuanian, Norwegian, Macedonian, Persian,
+# Polish, Portuguese, Romanian, Russian, Serbian, Serbian-Cyrillic, Slovak,
+# Slovene, Spanish, Swedish, Ukrainian, and Vietnamese.
+
+OUTPUT_LANGUAGE = English
+
+# If the BRIEF_MEMBER_DESC tag is set to YES (the default) Doxygen will
+# include brief member descriptions after the members that are listed in
+# the file and class documentation (similar to JavaDoc).
+# Set to NO to disable this.
+
+BRIEF_MEMBER_DESC = YES
+
+# If the REPEAT_BRIEF tag is set to YES (the default) Doxygen will prepend
+# the brief description of a member or function before the detailed description.
+# Note: if both HIDE_UNDOC_MEMBERS and BRIEF_MEMBER_DESC are set to NO, the
+# brief descriptions will be completely suppressed.
+
+REPEAT_BRIEF = YES
+
+# This tag implements a quasi-intelligent brief description abbreviator
+# that is used to form the text in various listings. Each string
+# in this list, if found as the leading text of the brief description, will be
+# stripped from the text and the result after processing the whole list, is
+# used as the annotated text. Otherwise, the brief description is used as-is.
+# If left blank, the following values are used ("$name" is automatically
+# replaced with the name of the entity): "The $name class" "The $name widget"
+# "The $name file" "is" "provides" "specifies" "contains"
+# "represents" "a" "an" "the"
+
+ABBREVIATE_BRIEF =
+
+# If the ALWAYS_DETAILED_SEC and REPEAT_BRIEF tags are both set to YES then
+# Doxygen will generate a detailed section even if there is only a brief
+# description.
+
+ALWAYS_DETAILED_SEC = NO
+
+# If the INLINE_INHERITED_MEMB tag is set to YES, doxygen will show all
+# inherited members of a class in the documentation of that class as if those
+# members were ordinary class members. Constructors, destructors and assignment
+# operators of the base classes will not be shown.
+
+INLINE_INHERITED_MEMB = NO
+
+# If the FULL_PATH_NAMES tag is set to YES then Doxygen will prepend the full
+# path before files name in the file list and in the header files. If set
+# to NO the shortest path that makes the file name unique will be used.
+
+FULL_PATH_NAMES = YES
+
+# If the FULL_PATH_NAMES tag is set to YES then the STRIP_FROM_PATH tag
+# can be used to strip a user-defined part of the path. Stripping is
+# only done if one of the specified strings matches the left-hand part of
+# the path. The tag can be used to show relative paths in the file list.
+# If left blank the directory from which doxygen is run is used as the
+# path to strip.
+
+STRIP_FROM_PATH =
+
+# The STRIP_FROM_INC_PATH tag can be used to strip a user-defined part of
+# the path mentioned in the documentation of a class, which tells
+# the reader which header file to include in order to use a class.
+# If left blank only the name of the header file containing the class
+# definition is used. Otherwise one should specify the include paths that
+# are normally passed to the compiler using the -I flag.
+
+STRIP_FROM_INC_PATH =
+
+# If the SHORT_NAMES tag is set to YES, doxygen will generate much shorter
+# (but less readable) file names. This can be useful if your file system
+# doesn't support long names like on DOS, Mac, or CD-ROM.
+
+SHORT_NAMES = NO
+
+# If the JAVADOC_AUTOBRIEF tag is set to YES then Doxygen
+# will interpret the first line (until the first dot) of a JavaDoc-style
+# comment as the brief description. If set to NO, the JavaDoc
+# comments will behave just like regular Qt-style comments
+# (thus requiring an explicit @brief command for a brief description.)
+
+JAVADOC_AUTOBRIEF = YES
+
+# If the QT_AUTOBRIEF tag is set to YES then Doxygen will
+# interpret the first line (until the first dot) of a Qt-style
+# comment as the brief description. If set to NO, the comments
+# will behave just like regular Qt-style comments (thus requiring
+# an explicit \brief command for a brief description.)
+
+QT_AUTOBRIEF = NO
+
+# The MULTILINE_CPP_IS_BRIEF tag can be set to YES to make Doxygen
+# treat a multi-line C++ special comment block (i.e. a block of //! or ///
+# comments) as a brief description. This used to be the default behaviour.
+# The new default is to treat a multi-line C++ comment block as a detailed
+# description. Set this tag to YES if you prefer the old behaviour instead.
+
+MULTILINE_CPP_IS_BRIEF = NO
+
+# If the INHERIT_DOCS tag is set to YES (the default) then an undocumented
+# member inherits the documentation from any documented member that it
+# re-implements.
+
+INHERIT_DOCS = YES
+
+# If the SEPARATE_MEMBER_PAGES tag is set to YES, then doxygen will produce
+# a new page for each member. If set to NO, the documentation of a member will
+# be part of the file/class/namespace that contains it.
+
+SEPARATE_MEMBER_PAGES = NO
+
+# The TAB_SIZE tag can be used to set the number of spaces in a tab.
+# Doxygen uses this value to replace tabs by spaces in code fragments.
+
+TAB_SIZE = 8
+
+# This tag can be used to specify a number of aliases that acts
+# as commands in the documentation. An alias has the form "name=value".
+# For example adding "sideeffect=\par Side Effects:\n" will allow you to
+# put the command \sideeffect (or @sideeffect) in the documentation, which
+# will result in a user-defined paragraph with heading "Side Effects:".
+# You can put \n's in the value part of an alias to insert newlines.
+
+ALIASES =
+
+# This tag can be used to specify a number of word-keyword mappings (TCL only).
+# A mapping has the form "name=value". For example adding
+# "class=itcl::class" will allow you to use the command class in the
+# itcl::class meaning.
+
+TCL_SUBST =
+
+# Set the OPTIMIZE_OUTPUT_FOR_C tag to YES if your project consists of C
+# sources only. Doxygen will then generate output that is more tailored for C.
+# For instance, some of the names that are used will be different. The list
+# of all members will be omitted, etc.
+
+OPTIMIZE_OUTPUT_FOR_C = NO
+
+# Set the OPTIMIZE_OUTPUT_JAVA tag to YES if your project consists of Java
+# sources only. Doxygen will then generate output that is more tailored for
+# Java. For instance, namespaces will be presented as packages, qualified
+# scopes will look different, etc.
+
+OPTIMIZE_OUTPUT_JAVA = NO
+
+# Set the OPTIMIZE_FOR_FORTRAN tag to YES if your project consists of Fortran
+# sources only. Doxygen will then generate output that is more tailored for
+# Fortran.
+
+OPTIMIZE_FOR_FORTRAN = NO
+
+# Set the OPTIMIZE_OUTPUT_VHDL tag to YES if your project consists of VHDL
+# sources. Doxygen will then generate output that is tailored for
+# VHDL.
+
+OPTIMIZE_OUTPUT_VHDL = NO
+
+# Doxygen selects the parser to use depending on the extension of the files it
+# parses. With this tag you can assign which parser to use for a given extension.
+# Doxygen has a built-in mapping, but you can override or extend it using this
+# tag. The format is ext=language, where ext is a file extension, and language
+# is one of the parsers supported by doxygen: IDL, Java, Javascript, CSharp, C,
+# C++, D, PHP, Objective-C, Python, Fortran, VHDL, C, C++. For instance to make
+# doxygen treat .inc files as Fortran files (default is PHP), and .f files as C
+# (default is Fortran), use: inc=Fortran f=C. Note that for custom extensions
+# you also need to set FILE_PATTERNS otherwise the files are not read by doxygen.
+
+EXTENSION_MAPPING =
+
+# If you use STL classes (i.e. std::string, std::vector, etc.) but do not want
+# to include (a tag file for) the STL sources as input, then you should
+# set this tag to YES in order to let doxygen match functions declarations and
+# definitions whose arguments contain STL classes (e.g. func(std::string); v.s.
+# func(std::string) {}). This also makes the inheritance and collaboration
+# diagrams that involve STL classes more complete and accurate.
+
+BUILTIN_STL_SUPPORT = NO
+
+# If you use Microsoft's C++/CLI language, you should set this option to YES to
+# enable parsing support.
+
+CPP_CLI_SUPPORT = NO
+
+# Set the SIP_SUPPORT tag to YES if your project consists of sip sources only.
+# Doxygen will parse them like normal C++ but will assume all classes use public
+# instead of private inheritance when no explicit protection keyword is present.
+
+SIP_SUPPORT = NO
+
+# For Microsoft's IDL there are propget and propput attributes to indicate getter
+# and setter methods for a property. Setting this option to YES (the default)
+# will make doxygen replace the get and set methods by a property in the
+# documentation. This will only work if the methods are indeed getting or
+# setting a simple type. If this is not the case, or you want to show the
+# methods anyway, you should set this option to NO.
+
+IDL_PROPERTY_SUPPORT = YES
+
+# If member grouping is used in the documentation and the DISTRIBUTE_GROUP_DOC
+# tag is set to YES, then doxygen will reuse the documentation of the first
+# member in the group (if any) for the other members of the group. By default
+# all members of a group must be documented explicitly.
+
+DISTRIBUTE_GROUP_DOC = NO
+
+# Set the SUBGROUPING tag to YES (the default) to allow class member groups of
+# the same type (for instance a group of public functions) to be put as a
+# subgroup of that type (e.g. under the Public Functions section). Set it to
+# NO to prevent subgrouping. Alternatively, this can be done per class using
+# the \nosubgrouping command.
+
+SUBGROUPING = YES
+
+# When the INLINE_GROUPED_CLASSES tag is set to YES, classes, structs and
+# unions are shown inside the group in which they are included (e.g. using
+# @ingroup) instead of on a separate page (for HTML and Man pages) or
+# section (for LaTeX and RTF).
+
+INLINE_GROUPED_CLASSES = NO
+
+# When the INLINE_SIMPLE_STRUCTS tag is set to YES, structs, classes, and
+# unions with only public data fields will be shown inline in the documentation
+# of the scope in which they are defined (i.e. file, namespace, or group
+# documentation), provided this scope is documented. If set to NO (the default),
+# structs, classes, and unions are shown on a separate page (for HTML and Man
+# pages) or section (for LaTeX and RTF).
+
+INLINE_SIMPLE_STRUCTS = NO
+
+# When TYPEDEF_HIDES_STRUCT is enabled, a typedef of a struct, union, or enum
+# is documented as struct, union, or enum with the name of the typedef. So
+# typedef struct TypeS {} TypeT, will appear in the documentation as a struct
+# with name TypeT. When disabled the typedef will appear as a member of a file,
+# namespace, or class. And the struct will be named TypeS. This can typically
+# be useful for C code in case the coding convention dictates that all compound
+# types are typedef'ed and only the typedef is referenced, never the tag name.
+
+TYPEDEF_HIDES_STRUCT = NO
+
+# The SYMBOL_CACHE_SIZE determines the size of the internal cache use to
+# determine which symbols to keep in memory and which to flush to disk.
+# When the cache is full, less often used symbols will be written to disk.
+# For small to medium size projects (<1000 input files) the default value is
+# probably good enough. For larger projects a too small cache size can cause
+# doxygen to be busy swapping symbols to and from disk most of the time
+# causing a significant performance penalty.
+# If the system has enough physical memory increasing the cache will improve the
+# performance by keeping more symbols in memory. Note that the value works on
+# a logarithmic scale so increasing the size by one will roughly double the
+# memory usage. The cache size is given by this formula:
+# 2^(16+SYMBOL_CACHE_SIZE). The valid range is 0..9, the default is 0,
+# corresponding to a cache size of 2^16 = 65536 symbols.
+
+SYMBOL_CACHE_SIZE = 0
+
+# Similar to the SYMBOL_CACHE_SIZE the size of the symbol lookup cache can be
+# set using LOOKUP_CACHE_SIZE. This cache is used to resolve symbols given
+# their name and scope. Since this can be an expensive process and often the
+# same symbol appear multiple times in the code, doxygen keeps a cache of
+# pre-resolved symbols. If the cache is too small doxygen will become slower.
+# If the cache is too large, memory is wasted. The cache size is given by this
+# formula: 2^(16+LOOKUP_CACHE_SIZE). The valid range is 0..9, the default is 0,
+# corresponding to a cache size of 2^16 = 65536 symbols.
+
+LOOKUP_CACHE_SIZE = 0
+
+#---------------------------------------------------------------------------
+# Build related configuration options
+#---------------------------------------------------------------------------
+
+# If the EXTRACT_ALL tag is set to YES doxygen will assume all entities in
+# documentation are documented, even if no documentation was available.
+# Private class members and static file members will be hidden unless
+# the EXTRACT_PRIVATE and EXTRACT_STATIC tags are set to YES
+
+EXTRACT_ALL = YES
+
+# If the EXTRACT_PRIVATE tag is set to YES all private members of a class
+# will be included in the documentation.
+
+EXTRACT_PRIVATE = YES
+
+# If the EXTRACT_STATIC tag is set to YES all static members of a file
+# will be included in the documentation.
+
+EXTRACT_STATIC = YES
+
+# If the EXTRACT_LOCAL_CLASSES tag is set to YES classes (and structs)
+# defined locally in source files will be included in the documentation.
+# If set to NO only classes defined in header files are included.
+
+EXTRACT_LOCAL_CLASSES = YES
+
+# This flag is only useful for Objective-C code. When set to YES local
+# methods, which are defined in the implementation section but not in
+# the interface are included in the documentation.
+# If set to NO (the default) only methods in the interface are included.
+
+EXTRACT_LOCAL_METHODS = NO
+
+# If this flag is set to YES, the members of anonymous namespaces will be
+# extracted and appear in the documentation as a namespace called
+# 'anonymous_namespace{file}', where file will be replaced with the base
+# name of the file that contains the anonymous namespace. By default
+# anonymous namespaces are hidden.
+
+EXTRACT_ANON_NSPACES = NO
+
+# If the HIDE_UNDOC_MEMBERS tag is set to YES, Doxygen will hide all
+# undocumented members of documented classes, files or namespaces.
+# If set to NO (the default) these members will be included in the
+# various overviews, but no documentation section is generated.
+# This option has no effect if EXTRACT_ALL is enabled.
+
+HIDE_UNDOC_MEMBERS = NO
+
+# If the HIDE_UNDOC_CLASSES tag is set to YES, Doxygen will hide all
+# undocumented classes that are normally visible in the class hierarchy.
+# If set to NO (the default) these classes will be included in the various
+# overviews. This option has no effect if EXTRACT_ALL is enabled.
+
+HIDE_UNDOC_CLASSES = NO
+
+# If the HIDE_FRIEND_COMPOUNDS tag is set to YES, Doxygen will hide all
+# friend (class|struct|union) declarations.
+# If set to NO (the default) these declarations will be included in the
+# documentation.
+
+HIDE_FRIEND_COMPOUNDS = NO
+
+# If the HIDE_IN_BODY_DOCS tag is set to YES, Doxygen will hide any
+# documentation blocks found inside the body of a function.
+# If set to NO (the default) these blocks will be appended to the
+# function's detailed documentation block.
+
+HIDE_IN_BODY_DOCS = NO
+
+# The INTERNAL_DOCS tag determines if documentation
+# that is typed after a \internal command is included. If the tag is set
+# to NO (the default) then the documentation will be excluded.
+# Set it to YES to include the internal documentation.
+
+INTERNAL_DOCS = NO
+
+# If the CASE_SENSE_NAMES tag is set to NO then Doxygen will only generate
+# file names in lower-case letters. If set to YES upper-case letters are also
+# allowed. This is useful if you have classes or files whose names only differ
+# in case and if your file system supports case sensitive file names. Windows
+# and Mac users are advised to set this option to NO.
+
+CASE_SENSE_NAMES = YES
+
+# If the HIDE_SCOPE_NAMES tag is set to NO (the default) then Doxygen
+# will show members with their full class and namespace scopes in the
+# documentation. If set to YES the scope will be hidden.
+
+HIDE_SCOPE_NAMES = NO
+
+# If the SHOW_INCLUDE_FILES tag is set to YES (the default) then Doxygen
+# will put a list of the files that are included by a file in the documentation
+# of that file.
+
+SHOW_INCLUDE_FILES = YES
+
+# If the FORCE_LOCAL_INCLUDES tag is set to YES then Doxygen
+# will list include files with double quotes in the documentation
+# rather than with sharp brackets.
+
+FORCE_LOCAL_INCLUDES = NO
+
+# If the INLINE_INFO tag is set to YES (the default) then a tag [inline]
+# is inserted in the documentation for inline members.
+
+INLINE_INFO = YES
+
+# If the SORT_MEMBER_DOCS tag is set to YES (the default) then doxygen
+# will sort the (detailed) documentation of file and class members
+# alphabetically by member name. If set to NO the members will appear in
+# declaration order.
+
+SORT_MEMBER_DOCS = YES
+
+# If the SORT_BRIEF_DOCS tag is set to YES then doxygen will sort the
+# brief documentation of file, namespace and class members alphabetically
+# by member name. If set to NO (the default) the members will appear in
+# declaration order.
+
+SORT_BRIEF_DOCS = NO
+
+# If the SORT_MEMBERS_CTORS_1ST tag is set to YES then doxygen
+# will sort the (brief and detailed) documentation of class members so that
+# constructors and destructors are listed first. If set to NO (the default)
+# the constructors will appear in the respective orders defined by
+# SORT_MEMBER_DOCS and SORT_BRIEF_DOCS.
+# This tag will be ignored for brief docs if SORT_BRIEF_DOCS is set to NO
+# and ignored for detailed docs if SORT_MEMBER_DOCS is set to NO.
+
+SORT_MEMBERS_CTORS_1ST = NO
+
+# If the SORT_GROUP_NAMES tag is set to YES then doxygen will sort the
+# hierarchy of group names into alphabetical order. If set to NO (the default)
+# the group names will appear in their defined order.
+
+SORT_GROUP_NAMES = NO
+
+# If the SORT_BY_SCOPE_NAME tag is set to YES, the class list will be
+# sorted by fully-qualified names, including namespaces. If set to
+# NO (the default), the class list will be sorted only by class name,
+# not including the namespace part.
+# Note: This option is not very useful if HIDE_SCOPE_NAMES is set to YES.
+# Note: This option applies only to the class list, not to the
+# alphabetical list.
+
+SORT_BY_SCOPE_NAME = NO
+
+# If the STRICT_PROTO_MATCHING option is enabled and doxygen fails to
+# do proper type resolution of all parameters of a function it will reject a
+# match between the prototype and the implementation of a member function even
+# if there is only one candidate or it is obvious which candidate to choose
+# by doing a simple string match. By disabling STRICT_PROTO_MATCHING doxygen
+# will still accept a match between prototype and implementation in such cases.
+
+STRICT_PROTO_MATCHING = NO
+
+# The GENERATE_TODOLIST tag can be used to enable (YES) or
+# disable (NO) the todo list. This list is created by putting \todo
+# commands in the documentation.
+
+GENERATE_TODOLIST = YES
+
+# The GENERATE_TESTLIST tag can be used to enable (YES) or
+# disable (NO) the test list. This list is created by putting \test
+# commands in the documentation.
+
+GENERATE_TESTLIST = YES
+
+# The GENERATE_BUGLIST tag can be used to enable (YES) or
+# disable (NO) the bug list. This list is created by putting \bug
+# commands in the documentation.
+
+GENERATE_BUGLIST = YES
+
+# The GENERATE_DEPRECATEDLIST tag can be used to enable (YES) or
+# disable (NO) the deprecated list. This list is created by putting
+# \deprecated commands in the documentation.
+
+GENERATE_DEPRECATEDLIST= YES
+
+# The ENABLED_SECTIONS tag can be used to enable conditional
+# documentation sections, marked by \if sectionname ... \endif.
+
+ENABLED_SECTIONS =
+
+# The MAX_INITIALIZER_LINES tag determines the maximum number of lines
+# the initial value of a variable or macro consists of for it to appear in
+# the documentation. If the initializer consists of more lines than specified
+# here it will be hidden. Use a value of 0 to hide initializers completely.
+# The appearance of the initializer of individual variables and macros in the
+# documentation can be controlled using \showinitializer or \hideinitializer
+# command in the documentation regardless of this setting.
+
+MAX_INITIALIZER_LINES = 30
+
+# Set the SHOW_USED_FILES tag to NO to disable the list of files generated
+# at the bottom of the documentation of classes and structs. If set to YES the
+# list will mention the files that were used to generate the documentation.
+
+SHOW_USED_FILES = YES
+
+# If the sources in your project are distributed over multiple directories
+# then setting the SHOW_DIRECTORIES tag to YES will show the directory hierarchy
+# in the documentation. The default is NO.
+
+SHOW_DIRECTORIES = NO
+
+# Set the SHOW_FILES tag to NO to disable the generation of the Files page.
+# This will remove the Files entry from the Quick Index and from the
+# Folder Tree View (if specified). The default is YES.
+
+SHOW_FILES = YES
+
+# Set the SHOW_NAMESPACES tag to NO to disable the generation of the
+# Namespaces page.
+# This will remove the Namespaces entry from the Quick Index
+# and from the Folder Tree View (if specified). The default is YES.
+
+SHOW_NAMESPACES = YES
+
+# The FILE_VERSION_FILTER tag can be used to specify a program or script that
+# doxygen should invoke to get the current version for each file (typically from
+# the version control system). Doxygen will invoke the program by executing (via
+# popen()) the command , where is the value of
+# the FILE_VERSION_FILTER tag, and is the name of an input file
+# provided by doxygen. Whatever the program writes to standard output
+# is used as the file version. See the manual for examples.
+
+FILE_VERSION_FILTER =
+
+# The LAYOUT_FILE tag can be used to specify a layout file which will be parsed
+# by doxygen. The layout file controls the global structure of the generated
+# output files in an output format independent way. The create the layout file
+# that represents doxygen's defaults, run doxygen with the -l option.
+# You can optionally specify a file name after the option, if omitted
+# DoxygenLayout.xml will be used as the name of the layout file.
+
+LAYOUT_FILE =
+
+# The CITE_BIB_FILES tag can be used to specify one or more bib files
+# containing the references data. This must be a list of .bib files. The
+# .bib extension is automatically appended if omitted. Using this command
+# requires the bibtex tool to be installed. See also
+# http://en.wikipedia.org/wiki/BibTeX for more info. For LaTeX the style
+# of the bibliography can be controlled using LATEX_BIB_STYLE. To use this
+# feature you need bibtex and perl available in the search path.
+
+CITE_BIB_FILES =
+
+#---------------------------------------------------------------------------
+# configuration options related to warning and progress messages
+#---------------------------------------------------------------------------
+
+# The QUIET tag can be used to turn on/off the messages that are generated
+# by doxygen. Possible values are YES and NO. If left blank NO is used.
+
+QUIET = NO
+
+# The WARNINGS tag can be used to turn on/off the warning messages that are
+# generated by doxygen. Possible values are YES and NO. If left blank
+# NO is used.
+
+WARNINGS = YES
+
+# If WARN_IF_UNDOCUMENTED is set to YES, then doxygen will generate warnings
+# for undocumented members. If EXTRACT_ALL is set to YES then this flag will
+# automatically be disabled.
+
+WARN_IF_UNDOCUMENTED = YES
+
+# If WARN_IF_DOC_ERROR is set to YES, doxygen will generate warnings for
+# potential errors in the documentation, such as not documenting some
+# parameters in a documented function, or documenting parameters that
+# don't exist or using markup commands wrongly.
+
+WARN_IF_DOC_ERROR = YES
+
+# The WARN_NO_PARAMDOC option can be enabled to get warnings for
+# functions that are documented, but have no documentation for their parameters
+# or return value. If set to NO (the default) doxygen will only warn about
+# wrong or incomplete parameter documentation, but not about the absence of
+# documentation.
+
+WARN_NO_PARAMDOC = NO
+
+# The WARN_FORMAT tag determines the format of the warning messages that
+# doxygen can produce. The string should contain the $file, $line, and $text
+# tags, which will be replaced by the file and line number from which the
+# warning originated and the warning text. Optionally the format may contain
+# $version, which will be replaced by the version of the file (if it could
+# be obtained via FILE_VERSION_FILTER)
+
+WARN_FORMAT = "$file:$line: $text"
+
+# The WARN_LOGFILE tag can be used to specify a file to which warning
+# and error messages should be written. If left blank the output is written
+# to stderr.
+
+WARN_LOGFILE =
+
+#---------------------------------------------------------------------------
+# configuration options related to the input files
+#---------------------------------------------------------------------------
+
+# The INPUT tag can be used to specify the files and/or directories that contain
+# documented source files. You may enter file names like "myfile.cpp" or
+# directories like "/usr/src/myproject". Separate the files or directories
+# with spaces.
+
+INPUT = "src", "include/svo"
+
+# This tag can be used to specify the character encoding of the source files
+# that doxygen parses. Internally doxygen uses the UTF-8 encoding, which is
+# also the default input encoding. Doxygen uses libiconv (or the iconv built
+# into libc) for the transcoding. See http://www.gnu.org/software/libiconv for
+# the list of possible encodings.
+
+INPUT_ENCODING = UTF-8
+
+# If the value of the INPUT tag contains directories, you can use the
+# FILE_PATTERNS tag to specify one or more wildcard pattern (like *.cpp
+# and *.h) to filter out the source-files in the directories. If left
+# blank the following patterns are tested:
+# *.c *.cc *.cxx *.cpp *.c++ *.d *.java *.ii *.ixx *.ipp *.i++ *.inl *.h *.hh
+# *.hxx *.hpp *.h++ *.idl *.odl *.cs *.php *.php3 *.inc *.m *.mm *.dox *.py
+# *.f90 *.f *.for *.vhd *.vhdl
+
+FILE_PATTERNS =
+
+# The RECURSIVE tag can be used to turn specify whether or not subdirectories
+# should be searched for input files as well. Possible values are YES and NO.
+# If left blank NO is used.
+
+RECURSIVE = NO
+
+# The EXCLUDE tag can be used to specify files and/or directories that should be
+# excluded from the INPUT source files. This way you can easily exclude a
+# subdirectory from a directory tree whose root is specified with the INPUT tag.
+# Note that relative paths are relative to the directory from which doxygen is
+# run.
+
+EXCLUDE =
+
+# The EXCLUDE_SYMLINKS tag can be used to select whether or not files or
+# directories that are symbolic links (a Unix file system feature) are excluded
+# from the input.
+
+EXCLUDE_SYMLINKS = NO
+
+# If the value of the INPUT tag contains directories, you can use the
+# EXCLUDE_PATTERNS tag to specify one or more wildcard patterns to exclude
+# certain files from those directories. Note that the wildcards are matched
+# against the file with absolute path, so to exclude all test directories
+# for example use the pattern */test/*
+
+EXCLUDE_PATTERNS =
+
+# The EXCLUDE_SYMBOLS tag can be used to specify one or more symbol names
+# (namespaces, classes, functions, etc.) that should be excluded from the
+# output. The symbol name can be a fully qualified name, a word, or if the
+# wildcard * is used, a substring. Examples: ANamespace, AClass,
+# AClass::ANamespace, ANamespace::*Test
+
+EXCLUDE_SYMBOLS = g2o
+EXCLUDE_SYMBOLS += vk
+
+# The EXAMPLE_PATH tag can be used to specify one or more files or
+# directories that contain example code fragments that are included (see
+# the \include command).
+
+EXAMPLE_PATH =
+
+# If the value of the EXAMPLE_PATH tag contains directories, you can use the
+# EXAMPLE_PATTERNS tag to specify one or more wildcard pattern (like *.cpp
+# and *.h) to filter out the source-files in the directories. If left
+# blank all files are included.
+
+EXAMPLE_PATTERNS =
+
+# If the EXAMPLE_RECURSIVE tag is set to YES then subdirectories will be
+# searched for input files to be used with the \include or \dontinclude
+# commands irrespective of the value of the RECURSIVE tag.
+# Possible values are YES and NO. If left blank NO is used.
+
+EXAMPLE_RECURSIVE = NO
+
+# The IMAGE_PATH tag can be used to specify one or more files or
+# directories that contain image that are included in the documentation (see
+# the \image command).
+
+IMAGE_PATH =
+
+# The INPUT_FILTER tag can be used to specify a program that doxygen should
+# invoke to filter for each input file. Doxygen will invoke the filter program
+# by executing (via popen()) the command , where
+# is the value of the INPUT_FILTER tag, and is the name of an
+# input file. Doxygen will then use the output that the filter program writes
+# to standard output.
+# If FILTER_PATTERNS is specified, this tag will be
+# ignored.
+
+INPUT_FILTER =
+
+# The FILTER_PATTERNS tag can be used to specify filters on a per file pattern
+# basis.
+# Doxygen will compare the file name with each pattern and apply the
+# filter if there is a match.
+# The filters are a list of the form:
+# pattern=filter (like *.cpp=my_cpp_filter). See INPUT_FILTER for further
+# info on how filters are used. If FILTER_PATTERNS is empty or if
+# non of the patterns match the file name, INPUT_FILTER is applied.
+
+FILTER_PATTERNS =
+
+# If the FILTER_SOURCE_FILES tag is set to YES, the input filter (if set using
+# INPUT_FILTER) will be used to filter the input files when producing source
+# files to browse (i.e. when SOURCE_BROWSER is set to YES).
+
+FILTER_SOURCE_FILES = NO
+
+# The FILTER_SOURCE_PATTERNS tag can be used to specify source filters per file
+# pattern. A pattern will override the setting for FILTER_PATTERN (if any)
+# and it is also possible to disable source filtering for a specific pattern
+# using *.ext= (so without naming a filter). This option only has effect when
+# FILTER_SOURCE_FILES is enabled.
+
+FILTER_SOURCE_PATTERNS =
+
+#---------------------------------------------------------------------------
+# configuration options related to source browsing
+#---------------------------------------------------------------------------
+
+# If the SOURCE_BROWSER tag is set to YES then a list of source files will
+# be generated. Documented entities will be cross-referenced with these sources.
+# Note: To get rid of all source code in the generated output, make sure also
+# VERBATIM_HEADERS is set to NO.
+
+SOURCE_BROWSER = NO
+
+# Setting the INLINE_SOURCES tag to YES will include the body
+# of functions and classes directly in the documentation.
+
+INLINE_SOURCES = NO
+
+# Setting the STRIP_CODE_COMMENTS tag to YES (the default) will instruct
+# doxygen to hide any special comment blocks from generated source code
+# fragments. Normal C and C++ comments will always remain visible.
+
+STRIP_CODE_COMMENTS = YES
+
+# If the REFERENCED_BY_RELATION tag is set to YES
+# then for each documented function all documented
+# functions referencing it will be listed.
+
+REFERENCED_BY_RELATION = NO
+
+# If the REFERENCES_RELATION tag is set to YES
+# then for each documented function all documented entities
+# called/used by that function will be listed.
+
+REFERENCES_RELATION = NO
+
+# If the REFERENCES_LINK_SOURCE tag is set to YES (the default)
+# and SOURCE_BROWSER tag is set to YES, then the hyperlinks from
+# functions in REFERENCES_RELATION and REFERENCED_BY_RELATION lists will
+# link to the source code.
+# Otherwise they will link to the documentation.
+
+REFERENCES_LINK_SOURCE = YES
+
+# If the USE_HTAGS tag is set to YES then the references to source code
+# will point to the HTML generated by the htags(1) tool instead of doxygen
+# built-in source browser. The htags tool is part of GNU's global source
+# tagging system (see http://www.gnu.org/software/global/global.html). You
+# will need version 4.8.6 or higher.
+
+USE_HTAGS = NO
+
+# If the VERBATIM_HEADERS tag is set to YES (the default) then Doxygen
+# will generate a verbatim copy of the header file for each class for
+# which an include is specified. Set to NO to disable this.
+
+VERBATIM_HEADERS = YES
+
+#---------------------------------------------------------------------------
+# configuration options related to the alphabetical class index
+#---------------------------------------------------------------------------
+
+# If the ALPHABETICAL_INDEX tag is set to YES, an alphabetical index
+# of all compounds will be generated. Enable this if the project
+# contains a lot of classes, structs, unions or interfaces.
+
+ALPHABETICAL_INDEX = YES
+
+# If the alphabetical index is enabled (see ALPHABETICAL_INDEX) then
+# the COLS_IN_ALPHA_INDEX tag can be used to specify the number of columns
+# in which this list will be split (can be a number in the range [1..20])
+
+COLS_IN_ALPHA_INDEX = 5
+
+# In case all classes in a project start with a common prefix, all
+# classes will be put under the same header in the alphabetical index.
+# The IGNORE_PREFIX tag can be used to specify one or more prefixes that
+# should be ignored while generating the index headers.
+
+IGNORE_PREFIX =
+
+#---------------------------------------------------------------------------
+# configuration options related to the HTML output
+#---------------------------------------------------------------------------
+
+# If the GENERATE_HTML tag is set to YES (the default) Doxygen will
+# generate HTML output.
+
+GENERATE_HTML = YES
+
+# The HTML_OUTPUT tag is used to specify where the HTML docs will be put.
+# If a relative path is entered the value of OUTPUT_DIRECTORY will be
+# put in front of it. If left blank `html' will be used as the default path.
+
+HTML_OUTPUT = html
+
+# The HTML_FILE_EXTENSION tag can be used to specify the file extension for
+# each generated HTML page (for example: .htm,.php,.asp). If it is left blank
+# doxygen will generate files with .html extension.
+
+HTML_FILE_EXTENSION = .html
+
+# The HTML_HEADER tag can be used to specify a personal HTML header for
+# each generated HTML page. If it is left blank doxygen will generate a
+# standard header. Note that when using a custom header you are responsible
+# for the proper inclusion of any scripts and style sheets that doxygen
+# needs, which is dependent on the configuration options used.
+# It is advised to generate a default header using "doxygen -w html
+# header.html footer.html stylesheet.css YourConfigFile" and then modify
+# that header. Note that the header is subject to change so you typically
+# have to redo this when upgrading to a newer version of doxygen or when
+# changing the value of configuration settings such as GENERATE_TREEVIEW!
+
+HTML_HEADER =
+
+# The HTML_FOOTER tag can be used to specify a personal HTML footer for
+# each generated HTML page. If it is left blank doxygen will generate a
+# standard footer.
+
+HTML_FOOTER =
+
+# The HTML_STYLESHEET tag can be used to specify a user-defined cascading
+# style sheet that is used by each HTML page. It can be used to
+# fine-tune the look of the HTML output. If the tag is left blank doxygen
+# will generate a default style sheet. Note that doxygen will try to copy
+# the style sheet file to the HTML output directory, so don't put your own
+# style sheet in the HTML output directory as well, or it will be erased!
+
+HTML_STYLESHEET =
+
+# The HTML_EXTRA_FILES tag can be used to specify one or more extra images or
+# other source files which should be copied to the HTML output directory. Note
+# that these files will be copied to the base HTML output directory. Use the
+# $relpath$ marker in the HTML_HEADER and/or HTML_FOOTER files to load these
+# files. In the HTML_STYLESHEET file, use the file name only. Also note that
+# the files will be copied as-is; there are no commands or markers available.
+
+HTML_EXTRA_FILES =
+
+# The HTML_COLORSTYLE_HUE tag controls the color of the HTML output.
+# Doxygen will adjust the colors in the style sheet and background images
+# according to this color. Hue is specified as an angle on a colorwheel,
+# see http://en.wikipedia.org/wiki/Hue for more information.
+# For instance the value 0 represents red, 60 is yellow, 120 is green,
+# 180 is cyan, 240 is blue, 300 purple, and 360 is red again.
+# The allowed range is 0 to 359.
+
+HTML_COLORSTYLE_HUE = 220
+
+# The HTML_COLORSTYLE_SAT tag controls the purity (or saturation) of
+# the colors in the HTML output. For a value of 0 the output will use
+# grayscales only. A value of 255 will produce the most vivid colors.
+
+HTML_COLORSTYLE_SAT = 100
+
+# The HTML_COLORSTYLE_GAMMA tag controls the gamma correction applied to
+# the luminance component of the colors in the HTML output. Values below
+# 100 gradually make the output lighter, whereas values above 100 make
+# the output darker. The value divided by 100 is the actual gamma applied,
+# so 80 represents a gamma of 0.8, The value 220 represents a gamma of 2.2,
+# and 100 does not change the gamma.
+
+HTML_COLORSTYLE_GAMMA = 80
+
+# If the HTML_TIMESTAMP tag is set to YES then the footer of each generated HTML
+# page will contain the date and time when the page was generated. Setting
+# this to NO can help when comparing the output of multiple runs.
+
+HTML_TIMESTAMP = YES
+
+# If the HTML_ALIGN_MEMBERS tag is set to YES, the members of classes,
+# files or namespaces will be aligned in HTML using tables. If set to
+# NO a bullet list will be used.
+
+HTML_ALIGN_MEMBERS = YES
+
+# If the HTML_DYNAMIC_SECTIONS tag is set to YES then the generated HTML
+# documentation will contain sections that can be hidden and shown after the
+# page has loaded. For this to work a browser that supports
+# JavaScript and DHTML is required (for instance Mozilla 1.0+, Firefox
+# Netscape 6.0+, Internet explorer 5.0+, Konqueror, or Safari).
+
+HTML_DYNAMIC_SECTIONS = NO
+
+# If the GENERATE_DOCSET tag is set to YES, additional index files
+# will be generated that can be used as input for Apple's Xcode 3
+# integrated development environment, introduced with OSX 10.5 (Leopard).
+# To create a documentation set, doxygen will generate a Makefile in the
+# HTML output directory. Running make will produce the docset in that
+# directory and running "make install" will install the docset in
+# ~/Library/Developer/Shared/Documentation/DocSets so that Xcode will find
+# it at startup.
+# See http://developer.apple.com/tools/creatingdocsetswithdoxygen.html
+# for more information.
+
+GENERATE_DOCSET = NO
+
+# When GENERATE_DOCSET tag is set to YES, this tag determines the name of the
+# feed. A documentation feed provides an umbrella under which multiple
+# documentation sets from a single provider (such as a company or product suite)
+# can be grouped.
+
+DOCSET_FEEDNAME = "Doxygen generated docs"
+
+# When GENERATE_DOCSET tag is set to YES, this tag specifies a string that
+# should uniquely identify the documentation set bundle. This should be a
+# reverse domain-name style string, e.g. com.mycompany.MyDocSet. Doxygen
+# will append .docset to the name.
+
+DOCSET_BUNDLE_ID = org.doxygen.Project
+
+# When GENERATE_PUBLISHER_ID tag specifies a string that should uniquely identify
+# the documentation publisher. This should be a reverse domain-name style
+# string, e.g. com.mycompany.MyDocSet.documentation.
+
+DOCSET_PUBLISHER_ID = org.doxygen.Publisher
+
+# The GENERATE_PUBLISHER_NAME tag identifies the documentation publisher.
+
+DOCSET_PUBLISHER_NAME = Publisher
+
+# If the GENERATE_HTMLHELP tag is set to YES, additional index files
+# will be generated that can be used as input for tools like the
+# Microsoft HTML help workshop to generate a compiled HTML help file (.chm)
+# of the generated HTML documentation.
+
+GENERATE_HTMLHELP = NO
+
+# If the GENERATE_HTMLHELP tag is set to YES, the CHM_FILE tag can
+# be used to specify the file name of the resulting .chm file. You
+# can add a path in front of the file if the result should not be
+# written to the html output directory.
+
+CHM_FILE =
+
+# If the GENERATE_HTMLHELP tag is set to YES, the HHC_LOCATION tag can
+# be used to specify the location (absolute path including file name) of
+# the HTML help compiler (hhc.exe). If non-empty doxygen will try to run
+# the HTML help compiler on the generated index.hhp.
+
+HHC_LOCATION =
+
+# If the GENERATE_HTMLHELP tag is set to YES, the GENERATE_CHI flag
+# controls if a separate .chi index file is generated (YES) or that
+# it should be included in the master .chm file (NO).
+
+GENERATE_CHI = NO
+
+# If the GENERATE_HTMLHELP tag is set to YES, the CHM_INDEX_ENCODING
+# is used to encode HtmlHelp index (hhk), content (hhc) and project file
+# content.
+
+CHM_INDEX_ENCODING =
+
+# If the GENERATE_HTMLHELP tag is set to YES, the BINARY_TOC flag
+# controls whether a binary table of contents is generated (YES) or a
+# normal table of contents (NO) in the .chm file.
+
+BINARY_TOC = NO
+
+# The TOC_EXPAND flag can be set to YES to add extra items for group members
+# to the contents of the HTML help documentation and to the tree view.
+
+TOC_EXPAND = NO
+
+# If the GENERATE_QHP tag is set to YES and both QHP_NAMESPACE and
+# QHP_VIRTUAL_FOLDER are set, an additional index file will be generated
+# that can be used as input for Qt's qhelpgenerator to generate a
+# Qt Compressed Help (.qch) of the generated HTML documentation.
+
+GENERATE_QHP = NO
+
+# If the QHG_LOCATION tag is specified, the QCH_FILE tag can
+# be used to specify the file name of the resulting .qch file.
+# The path specified is relative to the HTML output folder.
+
+QCH_FILE =
+
+# The QHP_NAMESPACE tag specifies the namespace to use when generating
+# Qt Help Project output. For more information please see
+# http://doc.trolltech.com/qthelpproject.html#namespace
+
+QHP_NAMESPACE = org.doxygen.Project
+
+# The QHP_VIRTUAL_FOLDER tag specifies the namespace to use when generating
+# Qt Help Project output. For more information please see
+# http://doc.trolltech.com/qthelpproject.html#virtual-folders
+
+QHP_VIRTUAL_FOLDER = doc
+
+# If QHP_CUST_FILTER_NAME is set, it specifies the name of a custom filter to
+# add. For more information please see
+# http://doc.trolltech.com/qthelpproject.html#custom-filters
+
+QHP_CUST_FILTER_NAME =
+
+# The QHP_CUST_FILT_ATTRS tag specifies the list of the attributes of the
+# custom filter to add. For more information please see
+#
+# Qt Help Project / Custom Filters.
+
+QHP_CUST_FILTER_ATTRS =
+
+# The QHP_SECT_FILTER_ATTRS tag specifies the list of the attributes this
+# project's
+# filter section matches.
+#
+# Qt Help Project / Filter Attributes.
+
+QHP_SECT_FILTER_ATTRS =
+
+# If the GENERATE_QHP tag is set to YES, the QHG_LOCATION tag can
+# be used to specify the location of Qt's qhelpgenerator.
+# If non-empty doxygen will try to run qhelpgenerator on the generated
+# .qhp file.
+
+QHG_LOCATION =
+
+# If the GENERATE_ECLIPSEHELP tag is set to YES, additional index files
+# will be generated, which together with the HTML files, form an Eclipse help
+# plugin. To install this plugin and make it available under the help contents
+# menu in Eclipse, the contents of the directory containing the HTML and XML
+# files needs to be copied into the plugins directory of eclipse. The name of
+# the directory within the plugins directory should be the same as
+# the ECLIPSE_DOC_ID value. After copying Eclipse needs to be restarted before
+# the help appears.
+
+GENERATE_ECLIPSEHELP = NO
+
+# A unique identifier for the eclipse help plugin. When installing the plugin
+# the directory name containing the HTML and XML files should also have
+# this name.
+
+ECLIPSE_DOC_ID = org.doxygen.Project
+
+# The DISABLE_INDEX tag can be used to turn on/off the condensed index (tabs)
+# at top of each HTML page. The value NO (the default) enables the index and
+# the value YES disables it. Since the tabs have the same information as the
+# navigation tree you can set this option to NO if you already set
+# GENERATE_TREEVIEW to YES.
+
+DISABLE_INDEX = NO
+
+# The GENERATE_TREEVIEW tag is used to specify whether a tree-like index
+# structure should be generated to display hierarchical information.
+# If the tag value is set to YES, a side panel will be generated
+# containing a tree-like index structure (just like the one that
+# is generated for HTML Help). For this to work a browser that supports
+# JavaScript, DHTML, CSS and frames is required (i.e. any modern browser).
+# Windows users are probably better off using the HTML help feature.
+# Since the tree basically has the same information as the tab index you
+# could consider to set DISABLE_INDEX to NO when enabling this option.
+
+GENERATE_TREEVIEW = NO
+
+# The ENUM_VALUES_PER_LINE tag can be used to set the number of enum values
+# (range [0,1..20]) that doxygen will group on one line in the generated HTML
+# documentation. Note that a value of 0 will completely suppress the enum
+# values from appearing in the overview section.
+
+ENUM_VALUES_PER_LINE = 4
+
+# By enabling USE_INLINE_TREES, doxygen will generate the Groups, Directories,
+# and Class Hierarchy pages using a tree view instead of an ordered list.
+
+USE_INLINE_TREES = NO
+
+# If the treeview is enabled (see GENERATE_TREEVIEW) then this tag can be
+# used to set the initial width (in pixels) of the frame in which the tree
+# is shown.
+
+TREEVIEW_WIDTH = 250
+
+# When the EXT_LINKS_IN_WINDOW option is set to YES doxygen will open
+# links to external symbols imported via tag files in a separate window.
+
+EXT_LINKS_IN_WINDOW = NO
+
+# Use this tag to change the font size of Latex formulas included
+# as images in the HTML documentation. The default is 10. Note that
+# when you change the font size after a successful doxygen run you need
+# to manually remove any form_*.png images from the HTML output directory
+# to force them to be regenerated.
+
+FORMULA_FONTSIZE = 10
+
+# Use the FORMULA_TRANPARENT tag to determine whether or not the images
+# generated for formulas are transparent PNGs. Transparent PNGs are
+# not supported properly for IE 6.0, but are supported on all modern browsers.
+# Note that when changing this option you need to delete any form_*.png files
+# in the HTML output before the changes have effect.
+
+FORMULA_TRANSPARENT = YES
+
+# Enable the USE_MATHJAX option to render LaTeX formulas using MathJax
+# (see http://www.mathjax.org) which uses client side Javascript for the
+# rendering instead of using prerendered bitmaps. Use this if you do not
+# have LaTeX installed or if you want to formulas look prettier in the HTML
+# output. When enabled you also need to install MathJax separately and
+# configure the path to it using the MATHJAX_RELPATH option.
+
+USE_MATHJAX = NO
+
+# When MathJax is enabled you need to specify the location relative to the
+# HTML output directory using the MATHJAX_RELPATH option. The destination
+# directory should contain the MathJax.js script. For instance, if the mathjax
+# directory is located at the same level as the HTML output directory, then
+# MATHJAX_RELPATH should be ../mathjax. The default value points to the
+# mathjax.org site, so you can quickly see the result without installing
+# MathJax, but it is strongly recommended to install a local copy of MathJax
+# before deployment.
+
+MATHJAX_RELPATH = http://www.mathjax.org/mathjax
+
+# The MATHJAX_EXTENSIONS tag can be used to specify one or MathJax extension
+# names that should be enabled during MathJax rendering.
+
+MATHJAX_EXTENSIONS =
+
+# When the SEARCHENGINE tag is enabled doxygen will generate a search box
+# for the HTML output. The underlying search engine uses javascript
+# and DHTML and should work on any modern browser. Note that when using
+# HTML help (GENERATE_HTMLHELP), Qt help (GENERATE_QHP), or docsets
+# (GENERATE_DOCSET) there is already a search function so this one should
+# typically be disabled. For large projects the javascript based search engine
+# can be slow, then enabling SERVER_BASED_SEARCH may provide a better solution.
+
+SEARCHENGINE = YES
+
+# When the SERVER_BASED_SEARCH tag is enabled the search engine will be
+# implemented using a PHP enabled web server instead of at the web client
+# using Javascript. Doxygen will generate the search PHP script and index
+# file to put on the web server. The advantage of the server
+# based approach is that it scales better to large projects and allows
+# full text search. The disadvantages are that it is more difficult to setup
+# and does not have live searching capabilities.
+
+SERVER_BASED_SEARCH = NO
+
+#---------------------------------------------------------------------------
+# configuration options related to the LaTeX output
+#---------------------------------------------------------------------------
+
+# If the GENERATE_LATEX tag is set to YES (the default) Doxygen will
+# generate Latex output.
+
+GENERATE_LATEX = YES
+
+# The LATEX_OUTPUT tag is used to specify where the LaTeX docs will be put.
+# If a relative path is entered the value of OUTPUT_DIRECTORY will be
+# put in front of it. If left blank `latex' will be used as the default path.
+
+LATEX_OUTPUT = latex
+
+# The LATEX_CMD_NAME tag can be used to specify the LaTeX command name to be
+# invoked. If left blank `latex' will be used as the default command name.
+# Note that when enabling USE_PDFLATEX this option is only used for
+# generating bitmaps for formulas in the HTML output, but not in the
+# Makefile that is written to the output directory.
+
+LATEX_CMD_NAME = latex
+
+# The MAKEINDEX_CMD_NAME tag can be used to specify the command name to
+# generate index for LaTeX. If left blank `makeindex' will be used as the
+# default command name.
+
+MAKEINDEX_CMD_NAME = makeindex
+
+# If the COMPACT_LATEX tag is set to YES Doxygen generates more compact
+# LaTeX documents. This may be useful for small projects and may help to
+# save some trees in general.
+
+COMPACT_LATEX = NO
+
+# The PAPER_TYPE tag can be used to set the paper type that is used
+# by the printer. Possible values are: a4, letter, legal and
+# executive. If left blank a4wide will be used.
+
+PAPER_TYPE = a4
+
+# The EXTRA_PACKAGES tag can be to specify one or more names of LaTeX
+# packages that should be included in the LaTeX output.
+
+EXTRA_PACKAGES =
+
+# The LATEX_HEADER tag can be used to specify a personal LaTeX header for
+# the generated latex document. The header should contain everything until
+# the first chapter. If it is left blank doxygen will generate a
+# standard header. Notice: only use this tag if you know what you are doing!
+
+LATEX_HEADER =
+
+# The LATEX_FOOTER tag can be used to specify a personal LaTeX footer for
+# the generated latex document. The footer should contain everything after
+# the last chapter. If it is left blank doxygen will generate a
+# standard footer. Notice: only use this tag if you know what you are doing!
+
+LATEX_FOOTER =
+
+# If the PDF_HYPERLINKS tag is set to YES, the LaTeX that is generated
+# is prepared for conversion to pdf (using ps2pdf). The pdf file will
+# contain links (just like the HTML output) instead of page references
+# This makes the output suitable for online browsing using a pdf viewer.
+
+PDF_HYPERLINKS = YES
+
+# If the USE_PDFLATEX tag is set to YES, pdflatex will be used instead of
+# plain latex in the generated Makefile. Set this option to YES to get a
+# higher quality PDF documentation.
+
+USE_PDFLATEX = YES
+
+# If the LATEX_BATCHMODE tag is set to YES, doxygen will add the \\batchmode.
+# command to the generated LaTeX files. This will instruct LaTeX to keep
+# running if errors occur, instead of asking the user for help.
+# This option is also used when generating formulas in HTML.
+
+LATEX_BATCHMODE = NO
+
+# If LATEX_HIDE_INDICES is set to YES then doxygen will not
+# include the index chapters (such as File Index, Compound Index, etc.)
+# in the output.
+
+LATEX_HIDE_INDICES = NO
+
+# If LATEX_SOURCE_CODE is set to YES then doxygen will include
+# source code with syntax highlighting in the LaTeX output.
+# Note that which sources are shown also depends on other settings
+# such as SOURCE_BROWSER.
+
+LATEX_SOURCE_CODE = NO
+
+# The LATEX_BIB_STYLE tag can be used to specify the style to use for the
+# bibliography, e.g. plainnat, or ieeetr. The default style is "plain". See
+# http://en.wikipedia.org/wiki/BibTeX for more info.
+
+LATEX_BIB_STYLE = plain
+
+#---------------------------------------------------------------------------
+# configuration options related to the RTF output
+#---------------------------------------------------------------------------
+
+# If the GENERATE_RTF tag is set to YES Doxygen will generate RTF output
+# The RTF output is optimized for Word 97 and may not look very pretty with
+# other RTF readers or editors.
+
+GENERATE_RTF = NO
+
+# The RTF_OUTPUT tag is used to specify where the RTF docs will be put.
+# If a relative path is entered the value of OUTPUT_DIRECTORY will be
+# put in front of it. If left blank `rtf' will be used as the default path.
+
+RTF_OUTPUT = rtf
+
+# If the COMPACT_RTF tag is set to YES Doxygen generates more compact
+# RTF documents. This may be useful for small projects and may help to
+# save some trees in general.
+
+COMPACT_RTF = NO
+
+# If the RTF_HYPERLINKS tag is set to YES, the RTF that is generated
+# will contain hyperlink fields. The RTF file will
+# contain links (just like the HTML output) instead of page references.
+# This makes the output suitable for online browsing using WORD or other
+# programs which support those fields.
+# Note: wordpad (write) and others do not support links.
+
+RTF_HYPERLINKS = NO
+
+# Load style sheet definitions from file. Syntax is similar to doxygen's
+# config file, i.e. a series of assignments. You only have to provide
+# replacements, missing definitions are set to their default value.
+
+RTF_STYLESHEET_FILE =
+
+# Set optional variables used in the generation of an rtf document.
+# Syntax is similar to doxygen's config file.
+
+RTF_EXTENSIONS_FILE =
+
+#---------------------------------------------------------------------------
+# configuration options related to the man page output
+#---------------------------------------------------------------------------
+
+# If the GENERATE_MAN tag is set to YES (the default) Doxygen will
+# generate man pages
+
+GENERATE_MAN = NO
+
+# The MAN_OUTPUT tag is used to specify where the man pages will be put.
+# If a relative path is entered the value of OUTPUT_DIRECTORY will be
+# put in front of it. If left blank `man' will be used as the default path.
+
+MAN_OUTPUT = man
+
+# The MAN_EXTENSION tag determines the extension that is added to
+# the generated man pages (default is the subroutine's section .3)
+
+MAN_EXTENSION = .3
+
+# If the MAN_LINKS tag is set to YES and Doxygen generates man output,
+# then it will generate one additional man file for each entity
+# documented in the real man page(s). These additional files
+# only source the real man page, but without them the man command
+# would be unable to find the correct page. The default is NO.
+
+MAN_LINKS = NO
+
+#---------------------------------------------------------------------------
+# configuration options related to the XML output
+#---------------------------------------------------------------------------
+
+# If the GENERATE_XML tag is set to YES Doxygen will
+# generate an XML file that captures the structure of
+# the code including all documentation.
+
+GENERATE_XML = NO
+
+# The XML_OUTPUT tag is used to specify where the XML pages will be put.
+# If a relative path is entered the value of OUTPUT_DIRECTORY will be
+# put in front of it. If left blank `xml' will be used as the default path.
+
+XML_OUTPUT = xml
+
+# The XML_SCHEMA tag can be used to specify an XML schema,
+# which can be used by a validating XML parser to check the
+# syntax of the XML files.
+
+XML_SCHEMA =
+
+# The XML_DTD tag can be used to specify an XML DTD,
+# which can be used by a validating XML parser to check the
+# syntax of the XML files.
+
+XML_DTD =
+
+# If the XML_PROGRAMLISTING tag is set to YES Doxygen will
+# dump the program listings (including syntax highlighting
+# and cross-referencing information) to the XML output. Note that
+# enabling this will significantly increase the size of the XML output.
+
+XML_PROGRAMLISTING = YES
+
+#---------------------------------------------------------------------------
+# configuration options for the AutoGen Definitions output
+#---------------------------------------------------------------------------
+
+# If the GENERATE_AUTOGEN_DEF tag is set to YES Doxygen will
+# generate an AutoGen Definitions (see autogen.sf.net) file
+# that captures the structure of the code including all
+# documentation. Note that this feature is still experimental
+# and incomplete at the moment.
+
+GENERATE_AUTOGEN_DEF = NO
+
+#---------------------------------------------------------------------------
+# configuration options related to the Perl module output
+#---------------------------------------------------------------------------
+
+# If the GENERATE_PERLMOD tag is set to YES Doxygen will
+# generate a Perl module file that captures the structure of
+# the code including all documentation. Note that this
+# feature is still experimental and incomplete at the
+# moment.
+
+GENERATE_PERLMOD = NO
+
+# If the PERLMOD_LATEX tag is set to YES Doxygen will generate
+# the necessary Makefile rules, Perl scripts and LaTeX code to be able
+# to generate PDF and DVI output from the Perl module output.
+
+PERLMOD_LATEX = NO
+
+# If the PERLMOD_PRETTY tag is set to YES the Perl module output will be
+# nicely formatted so it can be parsed by a human reader.
+# This is useful
+# if you want to understand what is going on.
+# On the other hand, if this
+# tag is set to NO the size of the Perl module output will be much smaller
+# and Perl will parse it just the same.
+
+PERLMOD_PRETTY = YES
+
+# The names of the make variables in the generated doxyrules.make file
+# are prefixed with the string contained in PERLMOD_MAKEVAR_PREFIX.
+# This is useful so different doxyrules.make files included by the same
+# Makefile don't overwrite each other's variables.
+
+PERLMOD_MAKEVAR_PREFIX =
+
+#---------------------------------------------------------------------------
+# Configuration options related to the preprocessor
+#---------------------------------------------------------------------------
+
+# If the ENABLE_PREPROCESSING tag is set to YES (the default) Doxygen will
+# evaluate all C-preprocessor directives found in the sources and include
+# files.
+
+ENABLE_PREPROCESSING = YES
+
+# If the MACRO_EXPANSION tag is set to YES Doxygen will expand all macro
+# names in the source code. If set to NO (the default) only conditional
+# compilation will be performed. Macro expansion can be done in a controlled
+# way by setting EXPAND_ONLY_PREDEF to YES.
+
+MACRO_EXPANSION = NO
+
+# If the EXPAND_ONLY_PREDEF and MACRO_EXPANSION tags are both set to YES
+# then the macro expansion is limited to the macros specified with the
+# PREDEFINED and EXPAND_AS_DEFINED tags.
+
+EXPAND_ONLY_PREDEF = NO
+
+# If the SEARCH_INCLUDES tag is set to YES (the default) the includes files
+# pointed to by INCLUDE_PATH will be searched when a #include is found.
+
+SEARCH_INCLUDES = YES
+
+# The INCLUDE_PATH tag can be used to specify one or more directories that
+# contain include files that are not input files but should be processed by
+# the preprocessor.
+
+INCLUDE_PATH =
+
+# You can use the INCLUDE_FILE_PATTERNS tag to specify one or more wildcard
+# patterns (like *.h and *.hpp) to filter out the header-files in the
+# directories. If left blank, the patterns specified with FILE_PATTERNS will
+# be used.
+
+INCLUDE_FILE_PATTERNS =
+
+# The PREDEFINED tag can be used to specify one or more macro names that
+# are defined before the preprocessor is started (similar to the -D option of
+# gcc). The argument of the tag is a list of macros of the form: name
+# or name=definition (no spaces). If the definition and the = are
+# omitted =1 is assumed. To prevent a macro definition from being
+# undefined via #undef or recursively expanded use the := operator
+# instead of the = operator.
+
+PREDEFINED =
+
+# If the MACRO_EXPANSION and EXPAND_ONLY_PREDEF tags are set to YES then
+# this tag can be used to specify a list of macro names that should be expanded.
+# The macro definition that is found in the sources will be used.
+# Use the PREDEFINED tag if you want to use a different macro definition that
+# overrules the definition found in the source code.
+
+EXPAND_AS_DEFINED =
+
+# If the SKIP_FUNCTION_MACROS tag is set to YES (the default) then
+# doxygen's preprocessor will remove all references to function-like macros
+# that are alone on a line, have an all uppercase name, and do not end with a
+# semicolon, because these will confuse the parser if not removed.
+
+SKIP_FUNCTION_MACROS = YES
+
+#---------------------------------------------------------------------------
+# Configuration::additions related to external references
+#---------------------------------------------------------------------------
+
+# The TAGFILES option can be used to specify one or more tagfiles.
+# Optionally an initial location of the external documentation
+# can be added for each tagfile. The format of a tag file without
+# this location is as follows:
+#
+# TAGFILES = file1 file2 ...
+# Adding location for the tag files is done as follows:
+#
+# TAGFILES = file1=loc1 "file2 = loc2" ...
+# where "loc1" and "loc2" can be relative or absolute paths or
+# URLs. If a location is present for each tag, the installdox tool
+# does not have to be run to correct the links.
+# Note that each tag file must have a unique name
+# (where the name does NOT include the path)
+# If a tag file is not located in the directory in which doxygen
+# is run, you must also specify the path to the tagfile here.
+
+TAGFILES =
+
+# When a file name is specified after GENERATE_TAGFILE, doxygen will create
+# a tag file that is based on the input files it reads.
+
+GENERATE_TAGFILE =
+
+# If the ALLEXTERNALS tag is set to YES all external classes will be listed
+# in the class index. If set to NO only the inherited external classes
+# will be listed.
+
+ALLEXTERNALS = NO
+
+# If the EXTERNAL_GROUPS tag is set to YES all external groups will be listed
+# in the modules index. If set to NO, only the current project's groups will
+# be listed.
+
+EXTERNAL_GROUPS = YES
+
+# The PERL_PATH should be the absolute path and name of the perl script
+# interpreter (i.e. the result of `which perl').
+
+PERL_PATH = /usr/bin/perl
+
+#---------------------------------------------------------------------------
+# Configuration options related to the dot tool
+#---------------------------------------------------------------------------
+
+# If the CLASS_DIAGRAMS tag is set to YES (the default) Doxygen will
+# generate a inheritance diagram (in HTML, RTF and LaTeX) for classes with base
+# or super classes. Setting the tag to NO turns the diagrams off. Note that
+# this option also works with HAVE_DOT disabled, but it is recommended to
+# install and use dot, since it yields more powerful graphs.
+
+CLASS_DIAGRAMS = YES
+
+# You can define message sequence charts within doxygen comments using the \msc
+# command. Doxygen will then run the mscgen tool (see
+# http://www.mcternan.me.uk/mscgen/) to produce the chart and insert it in the
+# documentation. The MSCGEN_PATH tag allows you to specify the directory where
+# the mscgen tool resides. If left empty the tool is assumed to be found in the
+# default search path.
+
+MSCGEN_PATH =
+
+# If set to YES, the inheritance and collaboration graphs will hide
+# inheritance and usage relations if the target is undocumented
+# or is not a class.
+
+HIDE_UNDOC_RELATIONS = YES
+
+# If you set the HAVE_DOT tag to YES then doxygen will assume the dot tool is
+# available from the path. This tool is part of Graphviz, a graph visualization
+# toolkit from AT&T and Lucent Bell Labs. The other options in this section
+# have no effect if this option is set to NO (the default)
+
+HAVE_DOT = NO
+
+# The DOT_NUM_THREADS specifies the number of dot invocations doxygen is
+# allowed to run in parallel. When set to 0 (the default) doxygen will
+# base this on the number of processors available in the system. You can set it
+# explicitly to a value larger than 0 to get control over the balance
+# between CPU load and processing speed.
+
+DOT_NUM_THREADS = 0
+
+# By default doxygen will use the Helvetica font for all dot files that
+# doxygen generates. When you want a differently looking font you can specify
+# the font name using DOT_FONTNAME. You need to make sure dot is able to find
+# the font, which can be done by putting it in a standard location or by setting
+# the DOTFONTPATH environment variable or by setting DOT_FONTPATH to the
+# directory containing the font.
+
+DOT_FONTNAME = Helvetica
+
+# The DOT_FONTSIZE tag can be used to set the size of the font of dot graphs.
+# The default size is 10pt.
+
+DOT_FONTSIZE = 10
+
+# By default doxygen will tell dot to use the Helvetica font.
+# If you specify a different font using DOT_FONTNAME you can use DOT_FONTPATH to
+# set the path where dot can find it.
+
+DOT_FONTPATH =
+
+# If the CLASS_GRAPH and HAVE_DOT tags are set to YES then doxygen
+# will generate a graph for each documented class showing the direct and
+# indirect inheritance relations. Setting this tag to YES will force the
+# CLASS_DIAGRAMS tag to NO.
+
+CLASS_GRAPH = YES
+
+# If the COLLABORATION_GRAPH and HAVE_DOT tags are set to YES then doxygen
+# will generate a graph for each documented class showing the direct and
+# indirect implementation dependencies (inheritance, containment, and
+# class references variables) of the class with other documented classes.
+
+COLLABORATION_GRAPH = YES
+
+# If the GROUP_GRAPHS and HAVE_DOT tags are set to YES then doxygen
+# will generate a graph for groups, showing the direct groups dependencies
+
+GROUP_GRAPHS = YES
+
+# If the UML_LOOK tag is set to YES doxygen will generate inheritance and
+# collaboration diagrams in a style similar to the OMG's Unified Modeling
+# Language.
+
+UML_LOOK = NO
+
+# If set to YES, the inheritance and collaboration graphs will show the
+# relations between templates and their instances.
+
+TEMPLATE_RELATIONS = NO
+
+# If the ENABLE_PREPROCESSING, SEARCH_INCLUDES, INCLUDE_GRAPH, and HAVE_DOT
+# tags are set to YES then doxygen will generate a graph for each documented
+# file showing the direct and indirect include dependencies of the file with
+# other documented files.
+
+INCLUDE_GRAPH = YES
+
+# If the ENABLE_PREPROCESSING, SEARCH_INCLUDES, INCLUDED_BY_GRAPH, and
+# HAVE_DOT tags are set to YES then doxygen will generate a graph for each
+# documented header file showing the documented files that directly or
+# indirectly include this file.
+
+INCLUDED_BY_GRAPH = YES
+
+# If the CALL_GRAPH and HAVE_DOT options are set to YES then
+# doxygen will generate a call dependency graph for every global function
+# or class method. Note that enabling this option will significantly increase
+# the time of a run. So in most cases it will be better to enable call graphs
+# for selected functions only using the \callgraph command.
+
+CALL_GRAPH = NO
+
+# If the CALLER_GRAPH and HAVE_DOT tags are set to YES then
+# doxygen will generate a caller dependency graph for every global function
+# or class method. Note that enabling this option will significantly increase
+# the time of a run. So in most cases it will be better to enable caller
+# graphs for selected functions only using the \callergraph command.
+
+CALLER_GRAPH = NO
+
+# If the GRAPHICAL_HIERARCHY and HAVE_DOT tags are set to YES then doxygen
+# will generate a graphical hierarchy of all classes instead of a textual one.
+
+GRAPHICAL_HIERARCHY = YES
+
+# If the DIRECTORY_GRAPH, SHOW_DIRECTORIES and HAVE_DOT tags are set to YES
+# then doxygen will show the dependencies a directory has on other directories
+# in a graphical way. The dependency relations are determined by the #include
+# relations between the files in the directories.
+
+DIRECTORY_GRAPH = YES
+
+# The DOT_IMAGE_FORMAT tag can be used to set the image format of the images
+# generated by dot. Possible values are svg, png, jpg, or gif.
+# If left blank png will be used. If you choose svg you need to set
+# HTML_FILE_EXTENSION to xhtml in order to make the SVG files
+# visible in IE 9+ (other browsers do not have this requirement).
+
+DOT_IMAGE_FORMAT = png
+
+# If DOT_IMAGE_FORMAT is set to svg, then this option can be set to YES to
+# enable generation of interactive SVG images that allow zooming and panning.
+# Note that this requires a modern browser other than Internet Explorer.
+# Tested and working are Firefox, Chrome, Safari, and Opera. For IE 9+ you
+# need to set HTML_FILE_EXTENSION to xhtml in order to make the SVG files
+# visible. Older versions of IE do not have SVG support.
+
+INTERACTIVE_SVG = NO
+
+# The tag DOT_PATH can be used to specify the path where the dot tool can be
+# found. If left blank, it is assumed the dot tool can be found in the path.
+
+DOT_PATH =
+
+# The DOTFILE_DIRS tag can be used to specify one or more directories that
+# contain dot files that are included in the documentation (see the
+# \dotfile command).
+
+DOTFILE_DIRS =
+
+# The MSCFILE_DIRS tag can be used to specify one or more directories that
+# contain msc files that are included in the documentation (see the
+# \mscfile command).
+
+MSCFILE_DIRS =
+
+# The DOT_GRAPH_MAX_NODES tag can be used to set the maximum number of
+# nodes that will be shown in the graph. If the number of nodes in a graph
+# becomes larger than this value, doxygen will truncate the graph, which is
+# visualized by representing a node as a red box. Note that doxygen if the
+# number of direct children of the root node in a graph is already larger than
+# DOT_GRAPH_MAX_NODES then the graph will not be shown at all. Also note
+# that the size of a graph can be further restricted by MAX_DOT_GRAPH_DEPTH.
+
+DOT_GRAPH_MAX_NODES = 50
+
+# The MAX_DOT_GRAPH_DEPTH tag can be used to set the maximum depth of the
+# graphs generated by dot. A depth value of 3 means that only nodes reachable
+# from the root by following a path via at most 3 edges will be shown. Nodes
+# that lay further from the root node will be omitted. Note that setting this
+# option to 1 or 2 may greatly reduce the computation time needed for large
+# code bases. Also note that the size of a graph can be further restricted by
+# DOT_GRAPH_MAX_NODES. Using a depth of 0 means no depth restriction.
+
+MAX_DOT_GRAPH_DEPTH = 0
+
+# Set the DOT_TRANSPARENT tag to YES to generate images with a transparent
+# background. This is disabled by default, because dot on Windows does not
+# seem to support this out of the box. Warning: Depending on the platform used,
+# enabling this option may lead to badly anti-aliased labels on the edges of
+# a graph (i.e. they become hard to read).
+
+DOT_TRANSPARENT = NO
+
+# Set the DOT_MULTI_TARGETS tag to YES allow dot to generate multiple output
+# files in one run (i.e. multiple -o and -T options on the command line). This
+# makes dot run faster, but since only newer versions of dot (>1.8.10)
+# support this, this feature is disabled by default.
+
+DOT_MULTI_TARGETS = YES
+
+# If the GENERATE_LEGEND tag is set to YES (the default) Doxygen will
+# generate a legend page explaining the meaning of the various boxes and
+# arrows in the dot generated graphs.
+
+GENERATE_LEGEND = YES
+
+# If the DOT_CLEANUP tag is set to YES (the default) Doxygen will
+# remove the intermediate dot files that are used to generate
+# the various graphs.
+
+DOT_CLEANUP = YES
diff --git a/CPU_version/svo/doc/.gitignore b/CPU_version/svo/doc/.gitignore
new file mode 100644
index 0000000..cf57c3b
--- /dev/null
+++ b/CPU_version/svo/doc/.gitignore
@@ -0,0 +1,2 @@
+.gitignore
+!.gitignore
\ No newline at end of file
diff --git a/CPU_version/svo/doc/notation.ipe b/CPU_version/svo/doc/notation.ipe
new file mode 100644
index 0000000..f8c1efe
--- /dev/null
+++ b/CPU_version/svo/doc/notation.ipe
@@ -0,0 +1,312 @@
+
+
+
+
+
+
+
+0 0 m
+-1 0.333 l
+-1 -0.333 l
+h
+
+
+
+
+0 0 m
+-1 0.333 l
+-1 -0.333 l
+h
+
+
+
+
+0.6 0 0 0.6 0 0 e
+0.4 0 0 0.4 0 0 e
+
+
+
+
+0.6 0 0 0.6 0 0 e
+
+
+
+
+
+0.6 0 0 0.6 0 0 e
+0.4 0 0 0.4 0 0 e
+
+
+0.4 0 0 0.4 0 0 e
+
+
+
+
+
+-0.6 -0.6 m
+0.6 -0.6 l
+0.6 0.6 l
+-0.6 0.6 l
+h
+-0.4 -0.4 m
+0.4 -0.4 l
+0.4 0.4 l
+-0.4 0.4 l
+h
+
+
+
+
+-0.6 -0.6 m
+0.6 -0.6 l
+0.6 0.6 l
+-0.6 0.6 l
+h
+
+
+
+
+
+-0.6 -0.6 m
+0.6 -0.6 l
+0.6 0.6 l
+-0.6 0.6 l
+h
+-0.4 -0.4 m
+0.4 -0.4 l
+0.4 0.4 l
+-0.4 0.4 l
+h
+
+
+-0.4 -0.4 m
+0.4 -0.4 l
+0.4 0.4 l
+-0.4 0.4 l
+h
+
+
+
+
+
+
+-0.43 -0.57 m
+0.57 0.43 l
+0.43 0.57 l
+-0.57 -0.43 l
+h
+
+
+-0.43 0.57 m
+0.57 -0.43 l
+0.43 -0.57 l
+-0.57 0.43 l
+h
+
+
+
+
+
+0 0 m
+-1 0.333 l
+-1 -0.333 l
+h
+
+
+
+
+0 0 m
+-1 0.333 l
+-0.8 0 l
+-1 -0.333 l
+h
+
+
+
+
+0 0 m
+-1 0.333 l
+-0.8 0 l
+-1 -0.333 l
+h
+
+
+
+
+-1 0.333 m
+0 0 l
+-1 -0.333 l
+
+
+
+
+0 0 m
+-1 0.333 l
+-1 -0.333 l
+h
+-1 0 m
+-2 0.333 l
+-2 -0.333 l
+h
+
+
+
+
+0 0 m
+-1 0.333 l
+-1 -0.333 l
+h
+-1 0 m
+-2 0.333 l
+-2 -0.333 l
+h
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+$_w\mathtt{p}$
+$\overrightarrow{\mathcal{F}}_w$
+
+64 576 m
+96 576 l
+
+
+64 576 m
+64 608 l
+
+
+64 576 m
+80 592 l
+
+
+144 688 m
+176 688 l
+
+
+144 688 m
+144 720 l
+
+
+144 688 m
+160 704 l
+
+
+285.262 759.615 m
+162.017 735.456
+129.232 628.801 s
+
+
+$\mathtt{T}_{f,w}$
+
+128.441 626.745 m
+320.625 647.224 l
+
+
+331.681 587.399 m
+331.681 587.399 l
+
+$x$
+$z$
+$y$
+
+256 736 m
+224 704 l
+320 704 l
+352 736 l
+h
+
+
+285.262 759.265 m
+322.176 651.119 l
+
+$\mathtt{px}$
+
+285.612 758.915 m
+315.723 669.632 l
+
+$\mathtt{f}$
+$\overrightarrow{\mathcal{F}}_f$
+
+255.668 719.918 m
+274.728 719.957 l
+
+
+256.107 720.028 m
+246.655 710.527 l
+
+$v$
+$u$
+
+$\mathbf{p}$
+$_f\mathtt{p}$
+
+
diff --git a/CPU_version/svo/doc/notation.pdf b/CPU_version/svo/doc/notation.pdf
new file mode 100644
index 0000000..3bc3fcd
Binary files /dev/null and b/CPU_version/svo/doc/notation.pdf differ
diff --git a/CPU_version/svo/doc/notation.png b/CPU_version/svo/doc/notation.png
new file mode 100644
index 0000000..d613b70
Binary files /dev/null and b/CPU_version/svo/doc/notation.png differ
diff --git a/CPU_version/svo/include/svo/bundle_adjustment.h b/CPU_version/svo/include/svo/bundle_adjustment.h
new file mode 100644
index 0000000..f3c3501
--- /dev/null
+++ b/CPU_version/svo/include/svo/bundle_adjustment.h
@@ -0,0 +1,109 @@
+// This file is part of SVO - Semi-direct Visual Odometry.
+//
+// Copyright (C) 2014 Christian Forster
+// (Robotics and Perception Group, University of Zurich, Switzerland).
+//
+// SVO is free software: you can redistribute it and/or modify it under the
+// terms of the GNU General Public License as published by the Free Software
+// Foundation, either version 3 of the License, or any later version.
+//
+// SVO is distributed in the hope that it will be useful, but WITHOUT ANY
+// WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
+// FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.
+//
+// You should have received a copy of the GNU General Public License
+// along with this program. If not, see .
+
+#ifndef SVO_BUNDLE_ADJUSTMENT_H_
+#define SVO_BUNDLE_ADJUSTMENT_H_
+
+#include
+
+namespace g2o {
+class EdgeProjectXYZ2UV;
+class SparseOptimizer;
+class VertexSE3Expmap;
+class VertexSBAPointXYZ;
+}
+
+namespace svo {
+
+typedef g2o::EdgeProjectXYZ2UV g2oEdgeSE3;
+typedef g2o::VertexSE3Expmap g2oFrameSE3;
+typedef g2o::VertexSBAPointXYZ g2oPoint;
+
+class Frame;
+class Point;
+class Feature;
+class Map;
+
+/// Local, global and 2-view bundle adjustment with g2o
+namespace ba {
+
+/// Temporary container to hold the g2o edge with reference to frame and point.
+struct EdgeContainerSE3
+{
+ g2oEdgeSE3* edge;
+ Frame* frame;
+ Feature* feature;
+ bool is_deleted;
+ EdgeContainerSE3(g2oEdgeSE3* e, Frame* frame, Feature* feature) :
+ edge(e), frame(frame), feature(feature), is_deleted(false)
+ {}
+};
+
+/// Optimize two camera frames and their observed 3D points.
+/// Is used after initialization.
+void twoViewBA(Frame* frame1, Frame* frame2, double reproj_thresh, Map* map);
+
+/// Local bundle adjustment.
+/// Optimizes core_kfs and their observed map points while keeping the
+/// neighbourhood fixed.
+void localBA(
+ Frame* center_kf,
+ set* core_kfs,
+ Map* map,
+ size_t& n_incorrect_edges_1,
+ size_t& n_incorrect_edges_2,
+ double& init_error,
+ double& final_error);
+
+/// Global bundle adjustment.
+/// Optimizes the whole map. Is currently not used in SVO.
+void globalBA(Map* map);
+
+/// Initialize g2o with solver type, optimization strategy and camera model.
+void setupG2o(g2o::SparseOptimizer * optimizer);
+
+/// Run the optimization on the provided graph.
+void runSparseBAOptimizer(
+ g2o::SparseOptimizer* optimizer,
+ unsigned int num_iter,
+ double& init_error,
+ double& final_error);
+
+/// Create a g2o vertice from a keyframe object.
+g2oFrameSE3* createG2oFrameSE3(
+ Frame* kf,
+ size_t id,
+ bool fixed);
+
+/// Creates a g2o vertice from a mappoint object.
+g2oPoint* createG2oPoint(
+ Vector3d pos,
+ size_t id,
+ bool fixed);
+
+/// Creates a g2o edge between a g2o keyframe and mappoint vertice with the provided measurement.
+g2oEdgeSE3* createG2oEdgeSE3(
+ g2oFrameSE3* v_kf,
+ g2oPoint* v_mp,
+ const Vector2d& f_up,
+ bool robust_kernel,
+ double huber_width,
+ double weight = 1);
+
+} // namespace ba
+} // namespace svo
+
+#endif // SVO_BUNDLE_ADJUSTMENT_H_
diff --git a/CPU_version/svo/include/svo/config.h b/CPU_version/svo/include/svo/config.h
new file mode 100644
index 0000000..e1c97b1
--- /dev/null
+++ b/CPU_version/svo/include/svo/config.h
@@ -0,0 +1,177 @@
+// This file is part of SVO - Semi-direct Visual Odometry.
+//
+// Copyright (C) 2014 Christian Forster
+// (Robotics and Perception Group, University of Zurich, Switzerland).
+//
+// SVO is free software: you can redistribute it and/or modify it under the
+// terms of the GNU General Public License as published by the Free Software
+// Foundation, either version 3 of the License, or any later version.
+//
+// SVO is distributed in the hope that it will be useful, but WITHOUT ANY
+// WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
+// FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.
+//
+// You should have received a copy of the GNU General Public License
+// along with this program. If not, see .
+
+#ifndef SVO_CONFIG_H_
+#define SVO_CONFIG_H_
+
+#include
+#include
+#include
+
+namespace svo {
+
+using std::string;
+
+/// Global configuration file of SVO.
+/// Implements the Singleton design pattern to allow global access and to ensure
+/// that only one instance exists.
+class Config
+{
+public:
+ static Config& getInstance();
+
+ /// Base-name of the tracefiles.
+ static string& traceName() { return getInstance().trace_name; }
+
+ /// Directory where the tracefiles are saved.
+ static string& traceDir() { return getInstance().trace_dir; }
+
+ /// Number of pyramid levels used for features.
+ static size_t& nPyrLevels() { return getInstance().n_pyr_levels; }
+
+ /// Use the IMU to get relative rotations.
+ static bool& useImu() { return getInstance().use_imu; }
+
+ /// Number of keyframes in the core. The core-kfs are optimized through bundle adjustment.
+ static size_t& coreNKfs() { return getInstance().core_n_kfs; }
+
+ /// Initial scale of the map. Depends on the distance the camera is moved for the initialization.
+ static double& mapScale() { return getInstance().map_scale; }
+
+ /// Feature grid size of a cell in [px].
+ static size_t& gridSize() { return getInstance().grid_size; }
+
+ /// Initialization: Minimum required disparity between the first two frames.
+ static double& initMinDisparity() { return getInstance().init_min_disparity; }
+
+ /// Initialization: Minimum number of tracked features.
+ static size_t& initMinTracked() { return getInstance().init_min_tracked; }
+
+ /// Initialization: Minimum number of inliers after RANSAC.
+ static size_t& initMinInliers() { return getInstance().init_min_inliers; }
+
+ /// Maximum level of the Lucas Kanade tracker.
+ static size_t& kltMaxLevel() { return getInstance().klt_max_level; }
+
+ /// Minimum level of the Lucas Kanade tracker.
+ static size_t& kltMinLevel() { return getInstance().klt_min_level; }
+
+ /// Reprojection threshold [px].
+ static double& reprojThresh() { return getInstance().reproj_thresh; }
+
+ /// Reprojection threshold after pose optimization.
+ static double& poseOptimThresh() { return getInstance().poseoptim_thresh; }
+
+ /// Number of iterations in local bundle adjustment.
+ static size_t& poseOptimNumIter() { return getInstance().poseoptim_num_iter; }
+
+ /// Maximum number of points to optimize at every iteration.
+ static size_t& structureOptimMaxPts() { return getInstance().structureoptim_max_pts; }
+
+ /// Number of iterations in structure optimization.
+ static size_t& structureOptimNumIter() { return getInstance().structureoptim_num_iter; }
+
+ /// Reprojection threshold after bundle adjustment.
+ static double& lobaThresh() { return getInstance().loba_thresh; }
+
+ /// Threshold for the robust Huber kernel of the local bundle adjustment.
+ static double& lobaRobustHuberWidth() { return getInstance().loba_robust_huber_width; }
+
+ /// Number of iterations in the local bundle adjustment.
+ static size_t& lobaNumIter() { return getInstance().loba_num_iter; }
+
+ /// Minimum distance between two keyframes. Relative to the average height in the map.
+ static double& kfSelectMinDist() { return getInstance().kfselect_mindist; }
+
+ /// Select only features with a minimum Harris corner score for triangulation.
+ static double& triangMinCornerScore() { return getInstance().triang_min_corner_score; }
+
+ /// Subpixel refinement of reprojection and triangulation. Set to 0 if no subpix refinement required!
+ static size_t& subpixNIter() { return getInstance().subpix_n_iter; }
+
+ /// Limit the number of keyframes in the map. This makes nslam essentially.
+ /// a Visual Odometry. Set to 0 if unlimited number of keyframes are allowed.
+ /// Minimum number of keyframes is 3.
+ static size_t& maxNKfs() { return getInstance().max_n_kfs; }
+
+ /// How much (in milliseconds) is the camera delayed with respect to the imu.
+ static double& imgImuDelay() { return getInstance().img_imu_delay; }
+
+ /// Maximum number of features that should be tracked.
+ static size_t& maxFts() { return getInstance().max_fts; }
+
+ /// If the number of tracked features drops below this threshold. Tracking quality is bad.
+ static size_t& qualityMinFts() { return getInstance().quality_min_fts; }
+
+ /// If within one frame, this amount of features are dropped. Tracking quality is bad.
+ static int& qualityMaxFtsDrop() { return getInstance().quality_max_drop_fts; }
+
+ /// acc white noise in continuous.
+ static double& ACC_Noise() { return getInstance().ACC_noise; }
+
+ /// acc white noise in continuous.
+ static double& GYO_Noise() { return getInstance().GYO_noise; }
+ /// integration uncertainty continuous
+ static double& IUC() { return getInstance().iuc; }
+ // acc bias in continuous
+ static double& ABC() { return getInstance().abc; }
+ // gyro bias in continuous
+ static double& GBC() { return getInstance().gbc; }
+ // error in the bias used for preintegration
+ static double& EBP() { return getInstance().ebp; }
+ // the noise used for projection
+ static double& NRP() { return getInstance().nrp; }
+
+
+private:
+ Config();
+ Config(Config const&);
+ void operator=(Config const&);
+ string trace_name;
+ string trace_dir;
+ size_t n_pyr_levels;
+ bool use_imu;
+ size_t core_n_kfs;
+ double map_scale;
+ size_t grid_size;
+ double init_min_disparity;
+ size_t init_min_tracked;
+ size_t init_min_inliers;
+ size_t klt_max_level;
+ size_t klt_min_level;
+ double reproj_thresh;
+ double poseoptim_thresh;
+ size_t poseoptim_num_iter;
+ size_t structureoptim_max_pts;
+ size_t structureoptim_num_iter;
+ double loba_thresh;
+ double loba_robust_huber_width;
+ size_t loba_num_iter;
+ double kfselect_mindist;
+ double triang_min_corner_score;
+ size_t triang_half_patch_size;
+ size_t subpix_n_iter;
+ size_t max_n_kfs;
+ double img_imu_delay;
+ size_t max_fts;
+ size_t quality_min_fts;
+ int quality_max_drop_fts;
+ double nrp,ebp,gbc,abc,iuc,ACC_noise,GYO_noise;
+};
+
+} // namespace svo
+
+#endif // SVO_CONFIG_H_
diff --git a/CPU_version/svo/include/svo/depth_filter.h b/CPU_version/svo/include/svo/depth_filter.h
new file mode 100644
index 0000000..1612d0e
--- /dev/null
+++ b/CPU_version/svo/include/svo/depth_filter.h
@@ -0,0 +1,171 @@
+// This file is part of SVO - Semi-direct Visual Odometry.
+//
+// Copyright (C) 2014 Christian Forster
+// (Robotics and Perception Group, University of Zurich, Switzerland).
+//
+// SVO is free software: you can redistribute it and/or modify it under the
+// terms of the GNU General Public License as published by the Free Software
+// Foundation, either version 3 of the License, or any later version.
+//
+// SVO is distributed in the hope that it will be useful, but WITHOUT ANY
+// WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
+// FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.
+//
+// You should have received a copy of the GNU General Public License
+// along with this program. If not, see .
+
+#ifndef SVO_DEPTH_FILTER_H_
+#define SVO_DEPTH_FILTER_H_
+
+#include
+#include
+#include
+#include
+#include
+#include
+#include
+
+namespace svo {
+
+class Frame;
+class Feature;
+class Point;
+
+/// A seed is a probabilistic depth estimate for a single pixel.
+struct Seed
+{
+ EIGEN_MAKE_ALIGNED_OPERATOR_NEW
+
+ static int batch_counter;
+ static int seed_counter;
+ int batch_id; //!< Batch id is the id of the keyframe for which the seed was created.
+ int id; //!< Seed ID, only used for visualization.
+ Feature* ftr; //!< Feature in the keyframe for which the depth should be computed.
+ float a; //!< a of Beta distribution: When high, probability of inlier is large.
+ float b; //!< b of Beta distribution: When high, probability of outlier is large.
+ float mu; //!< Mean of normal distribution.
+ float z_range; //!< Max range of the possible depth.
+ float sigma2; //!< Variance of normal distribution.
+ Matrix2d patch_cov; //!< Patch covariance in reference image.
+ Seed(Feature* ftr, float depth_mean, float depth_min);
+};
+
+/// Depth filter implements the Bayesian Update proposed in:
+/// "Video-based, Real-Time Multi View Stereo" by G. Vogiatzis and C. Hernández.
+/// In Image and Vision Computing, 29(7):434-441, 2011.
+///
+/// The class uses a callback mechanism such that it can be used also by other
+/// algorithms than nslam and for simplified testing.
+class DepthFilter
+{
+public:
+ EIGEN_MAKE_ALIGNED_OPERATOR_NEW
+
+ typedef boost::unique_lock lock_t;
+ typedef boost::function callback_t;
+
+ /// Depth-filter config parameters
+ struct Options
+ {
+ bool check_ftr_angle; //!< gradient features are only updated if the epipolar line is orthogonal to the gradient.
+ bool epi_search_1d; //!< restrict Gauss Newton in the epipolar search to the epipolar line.
+ bool verbose; //!< display output.
+ bool use_photometric_disparity_error; //!< use photometric disparity error instead of 1px error in tau computation.
+ int max_n_kfs; //!< maximum number of keyframes for which we maintain seeds.
+ double sigma_i_sq; //!< image noise.
+ double seed_convergence_sigma2_thresh; //!< threshold on depth uncertainty for convergence.
+ Options()
+ : check_ftr_angle(false),
+ epi_search_1d(false),
+ verbose(false),
+ use_photometric_disparity_error(false),
+ max_n_kfs(3),
+ sigma_i_sq(5e-4),
+ seed_convergence_sigma2_thresh(200.0)
+ {}
+ } options_;
+
+ DepthFilter(
+ feature_detection::DetectorPtr feature_detector,
+ callback_t seed_converged_cb);
+
+ virtual ~DepthFilter();
+
+ /// Start this thread when seed updating should be in a parallel thread.
+ void startThread();
+
+ /// Stop the parallel thread that is running.
+ void stopThread();
+
+ /// Add frame to the queue to be processed.
+ void addFrame(FramePtr frame);
+
+ /// Add new keyframe to the queue
+ void addKeyframe(FramePtr frame, double depth_mean, double depth_min);
+
+ /// Remove all seeds which are initialized from the specified keyframe. This
+ /// function is used to make sure that no seeds points to a non-existent frame
+ /// when a frame is removed from the map.
+ void removeKeyframe(FramePtr frame);
+
+ /// If the map is reset, call this function such that we don't have pointers
+ /// to old frames.
+ void reset();
+
+ /// Returns a copy of the seeds belonging to frame. Thread-safe.
+ /// Can be used to compute the Next-Best-View in parallel.
+ /// IMPORTANT! Make sure you hold a valid reference counting pointer to frame
+ /// so it is not being deleted while you use it.
+ void getSeedsCopy(const FramePtr& frame, std::list& seeds);
+
+ /// Return a reference to the seeds. This is NOT THREAD SAFE!
+ //std::list& getSeeds() { return seeds_; } Bug 238
+ std::list>& getSeeds() { return seeds_; }
+
+ /// Bayes update of the seed, x is the measurement, tau2 the measurement uncertainty
+ static void updateSeed(
+ const float x,
+ const float tau2,
+ Seed* seed);
+
+ /// Compute the uncertainty of the measurement.
+ static double computeTau(
+ const SE3& T_ref_cur,
+ const Vector3d& f,
+ const double z,
+ const double px_error_angle);
+
+protected:
+ feature_detection::DetectorPtr feature_detector_;
+ callback_t seed_converged_cb_;
+ //std::list seeds_;
+ std::list> seeds_;
+ boost::mutex seeds_mut_;
+ bool seeds_updating_halt_; //!< Set this value to true when seeds updating should be interrupted.
+ boost::thread* thread_;
+ std::queue frame_queue_;
+ boost::mutex frame_queue_mut_;
+ boost::condition_variable frame_queue_cond_;
+ FramePtr new_keyframe_; //!< Next keyframe to extract new seeds.
+ bool new_keyframe_set_; //!< Do we have a new keyframe to process?.
+ double new_keyframe_min_depth_; //!< Minimum depth in the new keyframe. Used for range in new seeds.
+ double new_keyframe_mean_depth_; //!< Maximum depth in the new keyframe. Used for range in new seeds.
+ vk::PerformanceMonitor permon_; //!< Separate performance monitor since the DepthFilter runs in a parallel thread.
+ Matcher matcher_;
+
+ /// Initialize new seeds from a frame.
+ void initializeSeeds(FramePtr frame);
+
+ /// Update all seeds with a new measurement frame.
+ virtual void updateSeeds(FramePtr frame);
+
+ /// When a new keyframe arrives, the frame queue should be cleared.
+ void clearFrameQueue();
+
+ /// A thread that is continuously updating the seeds.
+ void updateSeedsLoop();
+};
+
+} // namespace svo
+
+#endif // SVO_DEPTH_FILTER_H_
diff --git a/CPU_version/svo/include/svo/feature.h b/CPU_version/svo/include/svo/feature.h
new file mode 100644
index 0000000..6f90a9a
--- /dev/null
+++ b/CPU_version/svo/include/svo/feature.h
@@ -0,0 +1,75 @@
+// This file is part of SVO - Semi-direct Visual Odometry.
+//
+// Copyright (C) 2014 Christian Forster
+// (Robotics and Perception Group, University of Zurich, Switzerland).
+//
+// SVO is free software: you can redistribute it and/or modify it under the
+// terms of the GNU General Public License as published by the Free Software
+// Foundation, either version 3 of the License, or any later version.
+//
+// SVO is distributed in the hope that it will be useful, but WITHOUT ANY
+// WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
+// FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.
+//
+// You should have received a copy of the GNU General Public License
+// along with this program. If not, see .
+
+#ifndef SVO_FEATURE_H_
+#define SVO_FEATURE_H_
+
+#include
+
+namespace svo {
+
+/// A salient image region that is tracked across frames.
+struct Feature
+{
+ EIGEN_MAKE_ALIGNED_OPERATOR_NEW
+
+ enum FeatureType {
+ CORNER,
+ EDGELET
+ };
+
+ FeatureType type; //!< Type can be corner or edgelet.
+ Frame* frame; //!< Pointer to frame in which the feature was detected.
+ Vector2d px; //!< Coordinates in pixels on pyramid level 0.
+ Vector3d f; //!< Unit-bearing vector of the feature.
+ int level; //!< Image pyramid level where feature was extracted.
+ Point* point; //!< Pointer to 3D point which corresponds to the feature.
+ Vector2d grad; //!< Dominant gradient direction for edglets, normalized.
+
+ Feature(Frame* _frame, const Vector2d& _px, int _level) :
+ type(CORNER),
+ frame(_frame),
+ px(_px),
+ f(frame->cam_->cam2world(px)),
+ level(_level),
+ point(NULL),
+ grad(1.0,0.0)
+ {}
+
+ Feature(Frame* _frame, const Vector2d& _px, const Vector3d& _f, int _level) :
+ type(CORNER),
+ frame(_frame),
+ px(_px),
+ f(_f),
+ level(_level),
+ point(NULL),
+ grad(1.0,0.0)
+ {}
+
+ Feature(Frame* _frame, Point* _point, const Vector2d& _px, const Vector3d& _f, int _level) :
+ type(CORNER),
+ frame(_frame),
+ px(_px),
+ f(_f),
+ level(_level),
+ point(_point),
+ grad(1.0,0.0)
+ {}
+};
+
+} // namespace svo
+
+#endif // SVO_FEATURE_H_
diff --git a/CPU_version/svo/include/svo/feature_alignment.h b/CPU_version/svo/include/svo/feature_alignment.h
new file mode 100644
index 0000000..12523e3
--- /dev/null
+++ b/CPU_version/svo/include/svo/feature_alignment.h
@@ -0,0 +1,63 @@
+// This file is part of SVO - Semi-direct Visual Odometry.
+//
+// Copyright (C) 2014 Christian Forster
+// (Robotics and Perception Group, University of Zurich, Switzerland).
+//
+// SVO is free software: you can redistribute it and/or modify it under the
+// terms of the GNU General Public License as published by the Free Software
+// Foundation, either version 3 of the License, or any later version.
+//
+// SVO is distributed in the hope that it will be useful, but WITHOUT ANY
+// WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
+// FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.
+//
+// You should have received a copy of the GNU General Public License
+// along with this program. If not, see .
+
+#ifndef SVO_FEATURE_ALIGNMENT_H_
+#define SVO_FEATURE_ALIGNMENT_H_
+
+#include
+
+namespace svo {
+
+/// Subpixel refinement of a reference feature patch with the current image.
+/// Implements the inverse-compositional approach (see "Lucas-Kanade 20 Years on"
+/// paper by Baker.
+namespace feature_alignment {
+
+bool align1D(
+ const cv::Mat& cur_img,
+ const Vector2f& dir, // direction in which the patch is allowed to move
+ uint8_t* ref_patch_with_border,
+ uint8_t* ref_patch,
+ const int n_iter,
+ Vector2d& cur_px_estimate,
+ double& h_inv);
+
+bool align2D(
+ const cv::Mat& cur_img,
+ uint8_t* ref_patch_with_border,
+ uint8_t* ref_patch,
+ const int n_iter,
+ Vector2d& cur_px_estimate,
+ bool no_simd = false);
+
+bool align2D_SSE2(
+ const cv::Mat& cur_img,
+ uint8_t* ref_patch_with_border,
+ uint8_t* ref_patch,
+ const int n_iter,
+ Vector2d& cur_px_estimate);
+
+bool align2D_NEON(
+ const cv::Mat& cur_img,
+ uint8_t* ref_patch_with_border,
+ uint8_t* ref_patch,
+ const int n_iter,
+ Vector2d& cur_px_estimate);
+
+} // namespace feature_alignment
+} // namespace svo
+
+#endif // SVO_FEATURE_ALIGNMENT_H_
diff --git a/CPU_version/svo/include/svo/feature_detection.h b/CPU_version/svo/include/svo/feature_detection.h
new file mode 100644
index 0000000..d0990a8
--- /dev/null
+++ b/CPU_version/svo/include/svo/feature_detection.h
@@ -0,0 +1,107 @@
+// This file is part of SVO - Semi-direct Visual Odometry.
+//
+// Copyright (C) 2014 Christian Forster
+// (Robotics and Perception Group, University of Zurich, Switzerland).
+//
+// SVO is free software: you can redistribute it and/or modify it under the
+// terms of the GNU General Public License as published by the Free Software
+// Foundation, either version 3 of the License, or any later version.
+//
+// SVO is distributed in the hope that it will be useful, but WITHOUT ANY
+// WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
+// FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.
+//
+// You should have received a copy of the GNU General Public License
+// along with this program. If not, see .
+
+#ifndef SVO_FEATURE_DETECTION_H_
+#define SVO_FEATURE_DETECTION_H_
+
+#include
+#include
+
+namespace svo {
+
+/// Implementation of various feature detectors.
+namespace feature_detection {
+
+/// Temporary container used for corner detection. Features are initialized from these.
+struct Corner
+{
+ int x; //!< x-coordinate of corner in the image.
+ int y; //!< y-coordinate of corner in the image.
+ int level; //!< pyramid level of the corner.
+ float score; //!< shi-tomasi score of the corner.
+ float angle; //!< for gradient-features: dominant gradient angle.
+ Corner(int x, int y, float score, int level, float angle) :
+ x(x), y(y), level(level), score(score), angle(angle)
+ {}
+};
+typedef vector Corners;
+
+/// All detectors should derive from this abstract class.
+class AbstractDetector
+{
+public:
+ AbstractDetector(
+ const int img_width,
+ const int img_height,
+ const int cell_size,
+ const int n_pyr_levels);
+
+ virtual ~AbstractDetector() {};
+
+ virtual void detect(
+ Frame* frame,
+ const ImgPyr& img_pyr,
+ const double detection_threshold,
+ Features& fts) = 0;
+
+ /// Flag the grid cell as occupied
+ void setGridOccpuancy(const Vector2d& px);
+
+ /// Set grid cells of existing features as occupied
+ void setExistingFeatures(const Features& fts);
+
+protected:
+
+ static const int border_ = 8; //!< no feature should be within 8px of border.
+ const int cell_size_;
+ const int n_pyr_levels_;
+ const int grid_n_cols_;
+ const int grid_n_rows_;
+ vector grid_occupancy_;
+
+ void resetGrid();
+
+ inline int getCellIndex(int x, int y, int level)
+ {
+ const int scale = (1< DetectorPtr;
+
+/// FAST detector by Edward Rosten.
+class FastDetector : public AbstractDetector
+{
+public:
+ FastDetector(
+ const int img_width,
+ const int img_height,
+ const int cell_size,
+ const int n_pyr_levels);
+
+ virtual ~FastDetector() {}
+
+ virtual void detect(
+ Frame* frame,
+ const ImgPyr& img_pyr,
+ const double detection_threshold,
+ Features& fts);
+};
+
+} // namespace feature_detection
+} // namespace svo
+
+#endif // SVO_FEATURE_DETECTION_H_
diff --git a/CPU_version/svo/include/svo/frame.h b/CPU_version/svo/include/svo/frame.h
new file mode 100644
index 0000000..c926277
--- /dev/null
+++ b/CPU_version/svo/include/svo/frame.h
@@ -0,0 +1,154 @@
+// This file is part of SVO - Semi-direct Visual Odometry.
+//
+// Copyright (C) 2014 Christian Forster
+// (Robotics and Perception Group, University of Zurich, Switzerland).
+//
+// SVO is free software: you can redistribute it and/or modify it under the
+// terms of the GNU General Public License as published by the Free Software
+// Foundation, either version 3 of the License, or any later version.
+//
+// SVO is distributed in the hope that it will be useful, but WITHOUT ANY
+// WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
+// FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.
+//
+// You should have received a copy of the GNU General Public License
+// along with this program. If not, see .
+
+#ifndef SVO_FRAME_H_
+#define SVO_FRAME_H_
+
+#include
+#include
+#include
+#include
+#include
+
+namespace g2o {
+class VertexSE3Expmap;
+}
+typedef g2o::VertexSE3Expmap g2oFrameSE3;
+
+namespace svo {
+
+class Point;
+struct Feature;
+
+typedef list Features;
+typedef vector ImgPyr;
+
+/// A frame saves the image, the associated features and the estimated pose.
+class Frame : boost::noncopyable
+{
+public:
+ EIGEN_MAKE_ALIGNED_OPERATOR_NEW
+
+ static int frame_counter_; //!< Counts the number of created frames. Used to set the unique id.
+ int id_; //!< Unique id of the frame.
+ double timestamp_; //!< Timestamp of when the image was recorded.
+ vk::AbstractCamera* cam_; //!< Camera model.
+ Sophus::SE3 T_f_w_; //!< Transform (f)rame from (w)orld.
+ Matrix Cov_; //!< Covariance.
+ ImgPyr img_pyr_; //!< Image Pyramid.
+ Features fts_; //!< List of features in the image.
+ vector key_pts_; //!< Five features and associated 3D points which are used to detect if two frames have overlapping field of view.
+ bool is_keyframe_; //!< Was this frames selected as keyframe?
+ g2oFrameSE3* v_kf_; //!< Temporary pointer to the g2o node object of the keyframe.
+ int last_published_ts_; //!< Timestamp of last publishing.
+
+ Frame(vk::AbstractCamera* cam, const cv::Mat& img, double timestamp);
+ ~Frame();
+
+ /// Initialize new frame and create image pyramid.
+ void initFrame(const cv::Mat& img);
+
+ /// Select this frame as keyframe.
+ void setKeyframe();
+
+ /// Add a feature to the image
+ void addFeature(Feature* ftr);
+
+ /// The KeyPoints are those five features which are closest to the 4 image corners
+ /// and to the center and which have a 3D point assigned. These points are used
+ /// to quickly check whether two frames have overlapping field of view.
+ void setKeyPoints();
+
+ /// Check if we can select five better key-points.
+ void checkKeyPoints(Feature* ftr);
+
+ /// If a point is deleted, we must remove the corresponding key-point.
+ void removeKeyPoint(Feature* ftr);
+
+ /// Return number of point observations.
+ inline size_t nObs() const { return fts_.size(); }
+
+ /// Check if a point in (w)orld coordinate frame is visible in the image.
+ bool isVisible(const Vector3d& xyz_w) const;
+
+ /// Full resolution image stored in the frame.
+ inline const cv::Mat& img() const { return img_pyr_[0]; }
+
+ /// Was this frame selected as keyframe?
+ inline bool isKeyframe() const { return is_keyframe_; }
+
+ /// Transforms point coordinates in world-frame (w) to camera pixel coordinates (c).
+ inline Vector2d w2c(const Vector3d& xyz_w) const { return cam_->world2cam( T_f_w_ * xyz_w ); }
+
+ /// Transforms pixel coordinates (c) to frame unit sphere coordinates (f).
+ inline Vector3d c2f(const Vector2d& px) const { return cam_->cam2world(px[0], px[1]); }
+
+ /// Transforms pixel coordinates (c) to frame unit sphere coordinates (f).
+ inline Vector3d c2f(const double x, const double y) const { return cam_->cam2world(x, y); }
+
+ /// Transforms point coordinates in world-frame (w) to camera-frams (f).
+ inline Vector3d w2f(const Vector3d& xyz_w) const { return T_f_w_ * xyz_w; }
+
+ /// Transforms point from frame unit sphere (f) frame to world coordinate frame (w).
+ inline Vector3d f2w(const Vector3d& f) const { return T_f_w_.inverse() * f; }
+
+ /// Projects Point from unit sphere (f) in camera pixels (c).
+ inline Vector2d f2c(const Vector3d& f) const { return cam_->world2cam( f ); }
+
+ /// Return the pose of the frame in the (w)orld coordinate frame.
+ inline Vector3d pos() const { return T_f_w_.inverse().translation(); }
+
+ /// Frame jacobian for projection of 3D point in (f)rame coordinate to
+ /// unit plane coordinates uv (focal length = 1).
+ inline static void jacobian_xyz2uv(
+ const Vector3d& xyz_in_f,
+ Matrix& J)
+ {
+ const double x = xyz_in_f[0];
+ const double y = xyz_in_f[1];
+ const double z_inv = 1./xyz_in_f[2];
+ const double z_inv_2 = z_inv*z_inv;
+
+ J(0,0) = -z_inv; // -1/z
+ J(0,1) = 0.0; // 0
+ J(0,2) = x*z_inv_2; // x/z^2
+ J(0,3) = y*J(0,2); // x*y/z^2
+ J(0,4) = -(1.0 + x*J(0,2)); // -(1.0 + x^2/z^2)
+ J(0,5) = y*z_inv; // y/z
+
+ J(1,0) = 0.0; // 0
+ J(1,1) = -z_inv; // -1/z
+ J(1,2) = y*z_inv_2; // y/z^2
+ J(1,3) = 1.0 + y*J(1,2); // 1.0 + y^2/z^2
+ J(1,4) = -J(0,3); // -x*y/z^2
+ J(1,5) = -x*z_inv; // x/z
+ }
+};
+
+
+/// Some helper functions for the frame object.
+namespace frame_utils {
+
+/// Creates an image pyramid of half-sampled images.
+void createImgPyramid(const cv::Mat& img_level_0, int n_levels, ImgPyr& pyr);
+
+/// Get the average depth of the features in the image.
+bool getSceneDepth(const Frame& frame, double& depth_mean, double& depth_min);
+
+} // namespace frame_utils
+} // namespace svo
+
+#endif // SVO_FRAME_H_
diff --git a/CPU_version/svo/include/svo/frame_handler_base.h b/CPU_version/svo/include/svo/frame_handler_base.h
new file mode 100644
index 0000000..9a78639
--- /dev/null
+++ b/CPU_version/svo/include/svo/frame_handler_base.h
@@ -0,0 +1,123 @@
+// This file is part of SVO - Semi-direct Visual Odometry.
+//
+// Copyright (C) 2014 Christian Forster
+// (Robotics and Perception Group, University of Zurich, Switzerland).
+//
+// SVO is free software: you can redistribute it and/or modify it under the
+// terms of the GNU General Public License as published by the Free Software
+// Foundation, either version 3 of the License, or any later version.
+//
+// SVO is distributed in the hope that it will be useful, but WITHOUT ANY
+// WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
+// FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.
+//
+// You should have received a copy of the GNU General Public License
+// along with this program. If not, see .
+
+#ifndef SVO_FRAME_HANDLER_BASE_H_
+#define SVO_FRAME_HANDLER_BASE_H_
+
+#include
+#include
+#include
+#include
+#include
+#include
+#include
+#include
+
+namespace vk
+{
+class AbstractCamera;
+class PerformanceMonitor;
+}
+
+namespace svo
+{
+class Point;
+class Matcher;
+class DepthFilter;
+
+/// Base class for various VO pipelines. Manages the map and the state machine.
+class FrameHandlerBase : boost::noncopyable
+{
+public:
+ enum Stage {
+ STAGE_PAUSED,
+ STAGE_FIRST_FRAME,
+ STAGE_SECOND_FRAME,
+ STAGE_DEFAULT_FRAME,
+ STAGE_RELOCALIZING
+ };
+ enum TrackingQuality {
+ TRACKING_INSUFFICIENT,
+ TRACKING_BAD,
+ TRACKING_GOOD
+ };
+ enum UpdateResult {
+ RESULT_NO_KEYFRAME,
+ RESULT_IS_KEYFRAME,
+ RESULT_FAILURE
+ };
+
+ FrameHandlerBase();
+
+ virtual ~FrameHandlerBase();
+
+ /// Get the current map.
+ const Map& map() const { return map_; }
+
+ /// Will reset the map as soon as the current frame is finished processing.
+ void reset() { set_reset_ = true; }
+
+ /// Start processing.
+ void start() { set_start_ = true; }
+
+ /// Get the current stage of the algorithm.
+ Stage stage() const { return stage_; }
+
+ /// Get tracking quality.
+ TrackingQuality trackingQuality() const { return tracking_quality_; }
+
+ /// Get the processing time of the previous iteration.
+ double lastProcessingTime() const { return timer_.getTime(); }
+
+ /// Get the number of feature observations of the last frame.
+ size_t lastNumObservations() const { return num_obs_last_; }
+
+protected:
+ Stage stage_; //!< Current stage of the algorithm.
+ bool set_reset_; //!< Flag that the user can set. Will reset the system before the next iteration.
+ bool set_start_; //!< Flag the user can set to start the system when the next image is received.
+ Map map_; //!< Map of keyframes created by the slam system.
+ vk::Timer timer_; //!< Stopwatch to measure time to process frame.
+ vk::RingBuffer acc_frame_timings_; //!< Total processing time of the last 10 frames, used to give some user feedback on the performance.
+ vk::RingBuffer acc_num_obs_; //!< Number of observed features of the last 10 frames, used to give some user feedback on the tracking performance.
+ size_t num_obs_last_; //!< Number of observations in the previous frame.
+ TrackingQuality tracking_quality_; //!< An estimate of the tracking quality based on the number of tracked features.
+
+ /// Before a frame is processed, this function is called.
+ bool startFrameProcessingCommon(const double timestamp);
+
+ /// When a frame is finished processing, this function is called.
+ int finishFrameProcessingCommon(
+ const size_t update_id,
+ const UpdateResult dropout,
+ const size_t num_observations);
+
+ /// Reset the map and frame handler to start from scratch.
+ void resetCommon();
+
+ /// Reset the frame handler. Implement in derived class.
+ virtual void resetAll() { resetCommon(); }
+
+ /// Set the tracking quality based on the number of tracked features.
+ virtual void setTrackingQuality(const size_t num_observations);
+
+ /// Optimize some of the observed 3D points.
+ virtual void optimizeStructure(FramePtr frame, size_t max_n_pts, int max_iter);
+};
+
+} // namespace nslam
+
+#endif // SVO_FRAME_HANDLER_BASE_H_
diff --git a/CPU_version/svo/include/svo/frame_handler_mono.h b/CPU_version/svo/include/svo/frame_handler_mono.h
new file mode 100644
index 0000000..58789c8
--- /dev/null
+++ b/CPU_version/svo/include/svo/frame_handler_mono.h
@@ -0,0 +1,104 @@
+// This file is part of SVO - Semi-direct Visual Odometry.
+//
+// Copyright (C) 2014 Christian Forster
+// (Robotics and Perception Group, University of Zurich, Switzerland).
+//
+// SVO is free software: you can redistribute it and/or modify it under the
+// terms of the GNU General Public License as published by the Free Software
+// Foundation, either version 3 of the License, or any later version.
+//
+// SVO is distributed in the hope that it will be useful, but WITHOUT ANY
+// WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
+// FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.
+//
+// You should have received a copy of the GNU General Public License
+// along with this program. If not, see .
+
+#ifndef SVO_FRAME_HANDLER_H_
+#define SVO_FRAME_HANDLER_H_
+
+#include
+#include
+#include
+#include
+#include
+#include
+
+namespace svo {
+
+/// Monocular Visual Odometry Pipeline as described in the SVO paper.
+class FrameHandlerMono : public FrameHandlerBase
+{
+public:
+ EIGEN_MAKE_ALIGNED_OPERATOR_NEW
+
+ FrameHandlerMono(vk::AbstractCamera* cam,Sophus::SE3& SE_init);
+ virtual ~FrameHandlerMono();
+
+ /// Provide an image.
+ void addImage(const cv::Mat& img, double timestamp);
+
+ /// Set the first frame (used for synthetic datasets in benchmark node)
+ void setFirstFrame(const FramePtr& first_frame);
+
+ /// Get the last frame that has been processed.
+ FramePtr lastFrame() { return last_frame_; }
+
+ /// Get the set of spatially closest keyframes of the last frame.
+ const set& coreKeyframes() { return core_kfs_; }
+
+ /// Return the feature track to visualize the KLT tracking during initialization.
+ const vector& initFeatureTrackRefPx() const { return klt_homography_init_.px_ref_; }
+ const vector& initFeatureTrackCurPx() const { return klt_homography_init_.px_cur_; }
+
+ /// Access the depth filter.
+ DepthFilter* depthFilter() const { return depth_filter_; }
+
+ /// An external place recognition module may know where to relocalize.
+ bool relocalizeFrameAtPose(
+ const int keyframe_id,
+ const SE3& T_kf_f,
+ const cv::Mat& img,
+ const double timestamp);
+ std::unique_ptr imu_integPtr_; //Integration of IMU
+
+protected:
+ vk::AbstractCamera* cam_; //!< Camera model, can be ATAN, Pinhole or Ocam (see vikit).
+ Reprojector reprojector_; //!< Projects points from other keyframes into the current frame
+ FramePtr new_frame_; //!< Current frame.
+ FramePtr last_frame_; //!< Last frame, not necessarily a keyframe.
+ set core_kfs_; //!< Keyframes in the closer neighbourhood.
+ vector< pair > overlap_kfs_; //!< All keyframes with overlapping field of view. the paired number specifies how many common mappoints are observed TODO: why vector!?
+ initialization::KltHomographyInit klt_homography_init_; //!< Used to estimate pose of the first two keyframes by estimating a homography.
+ DepthFilter* depth_filter_; //!< Depth estimation algorithm runs in a parallel thread and is used to initialize new 3D points.
+
+
+ /// Initialize the visual odometry algorithm.
+ virtual void initialize();
+
+ /// Processes the first frame and sets it as a keyframe.
+ virtual UpdateResult processFirstFrame();
+
+ /// Processes all frames after the first frame until a keyframe is selected.
+ virtual UpdateResult processSecondFrame();
+
+ /// Processes all frames after the first two keyframes.
+ virtual UpdateResult processFrame();
+
+ /// Try relocalizing the frame at relative position to provided keyframe.
+ virtual UpdateResult relocalizeFrame(
+ const SE3& T_cur_ref,
+ FramePtr ref_keyframe);
+
+ /// Reset the frame handler. Implement in derived class.
+ virtual void resetAll();
+
+ /// Keyframe selection criterion.
+ virtual bool needNewKf(double scene_depth_mean);
+
+ void setCoreKfs(size_t n_closest);
+};
+
+} // namespace svo
+
+#endif // SVO_FRAME_HANDLER_H_
diff --git a/CPU_version/svo/include/svo/global.h b/CPU_version/svo/include/svo/global.h
new file mode 100644
index 0000000..03aff4e
--- /dev/null
+++ b/CPU_version/svo/include/svo/global.h
@@ -0,0 +1,101 @@
+// This file is part of SVO - Semi-direct Visual Odometry.
+//
+// Copyright (C) 2014 Christian Forster
+// (Robotics and Perception Group, University of Zurich, Switzerland).
+//
+// SVO is free software: you can redistribute it and/or modify it under the
+// terms of the GNU General Public License as published by the Free Software
+// Foundation, either version 3 of the License, or any later version.
+//
+// SVO is distributed in the hope that it will be useful, but WITHOUT ANY
+// WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
+// FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.
+//
+// You should have received a copy of the GNU General Public License
+// along with this program. If not, see .
+
+#ifndef SVO_GLOBAL_H_
+#define SVO_GLOBAL_H_
+
+#include
+#include
+#include
+#include
+#include
+#include
+
+#include
+#include
+#include
+#include
+#include
+#include
+#ifndef RPG_SVO_VIKIT_IS_VECTOR_SPECIALIZED //Guard for rpg_vikit
+#define RPG_SVO_VIKIT_IS_VECTOR_SPECIALIZED
+EIGEN_DEFINE_STL_VECTOR_SPECIALIZATION(Eigen::Vector3d)
+EIGEN_DEFINE_STL_VECTOR_SPECIALIZATION(Eigen::Vector2d)
+#endif
+
+#ifdef SVO_USE_ROS
+ #include
+ #define SVO_DEBUG_STREAM(x) ROS_DEBUG_STREAM(x)
+ #define SVO_INFO_STREAM(x) ROS_INFO_STREAM(x)
+ #define SVO_WARN_STREAM(x) ROS_WARN_STREAM(x)
+ #define SVO_WARN_STREAM_THROTTLE(rate, x) ROS_WARN_STREAM_THROTTLE(rate, x)
+ #define SVO_ERROR_STREAM(x) ROS_ERROR_STREAM(x)
+#else
+ #define SVO_INFO_STREAM(x) std::cerr<<"\033[0;0m[INFO] "< // Adapted from rosconsole. Copyright (c) 2008, Willow Garage, Inc.
+ #define SVO_WARN_STREAM_THROTTLE(rate, x) \
+ do { \
+ static double __log_stream_throttle__last_hit__ = 0.0; \
+ std::chrono::time_point __log_stream_throttle__now__ = \
+ std::chrono::system_clock::now(); \
+ if (__log_stream_throttle__last_hit__ + rate <= \
+ std::chrono::duration_cast( \
+ __log_stream_throttle__now__.time_since_epoch()).count()) { \
+ __log_stream_throttle__last_hit__ = \
+ std::chrono::duration_cast( \
+ __log_stream_throttle__now__.time_since_epoch()).count(); \
+ SVO_WARN_STREAM(x); \
+ } \
+ } while(0)
+#endif
+
+namespace svo
+{
+ using namespace Eigen;
+ using namespace Sophus;
+
+ const double EPS = 0.0000000001;
+ const double PI = 3.14159265;
+
+#ifdef SVO_TRACE
+ extern vk::PerformanceMonitor* g_permon;
+ #define SVO_LOG(value) g_permon->log(std::string((#value)),(value))
+ #define SVO_LOG2(value1, value2) SVO_LOG(value1); SVO_LOG(value2)
+ #define SVO_LOG3(value1, value2, value3) SVO_LOG2(value1, value2); SVO_LOG(value3)
+ #define SVO_LOG4(value1, value2, value3, value4) SVO_LOG2(value1, value2); SVO_LOG2(value3, value4)
+ #define SVO_START_TIMER(name) g_permon->startTimer((name))
+ #define SVO_STOP_TIMER(name) g_permon->stopTimer((name))
+#else
+ #define SVO_LOG(v)
+ #define SVO_LOG2(v1, v2)
+ #define SVO_LOG3(v1, v2, v3)
+ #define SVO_LOG4(v1, v2, v3, v4)
+ #define SVO_START_TIMER(name)
+ #define SVO_STOP_TIMER(name)
+#endif
+
+ class Frame;
+ typedef boost::shared_ptr FramePtr;
+} // namespace svo
+
+#endif // SVO_GLOBAL_H_
diff --git a/CPU_version/svo/include/svo/imu_integration.h b/CPU_version/svo/include/svo/imu_integration.h
new file mode 100644
index 0000000..11b9055
--- /dev/null
+++ b/CPU_version/svo/include/svo/imu_integration.h
@@ -0,0 +1,121 @@
+//
+// Created by gufe on 4/9/21.
+//
+
+#ifndef SVO_IMU_INTEGRATION_H
+#define SVO_IMU_INTEGRATION_H
+#include
+#include
+#include
+#include
+#include
+#include
+#include
+#include
+#include
+#include
+
+using gtsam::symbol_shorthand::B; // Bias (ax,ay,az,gx,gy,gz)
+using gtsam::symbol_shorthand::V; // Vel (xdot,ydot,zdot)
+using gtsam::symbol_shorthand::P; // Pose3 (x,y,z,r,p,y)
+using gtsam::symbol_shorthand::C; // Camera pose
+
+//https://stackoverflow.com/questions/24881799/get-index-in-c11-foreach-loop
+// Wrapper class
+template
+class enumerate_impl
+{
+public:
+ // The return value of the operator* of the iterator, this
+ // is what you will get inside of the for loop
+ struct item
+ {
+ size_t index;
+ typename T::value_type & item;
+ };
+ typedef item value_type;
+
+ // Custom iterator with minimal interface
+ struct iterator
+ {
+ iterator(typename T::iterator _it, size_t counter=0) :
+ it(_it), counter(counter)
+ {}
+
+ iterator operator++()
+ {
+ return iterator(++it, ++counter);
+ }
+
+ bool operator!=(iterator other)
+ {
+ return it != other.it;
+ }
+
+ typename T::iterator::value_type item()
+ {
+ return *it;
+ }
+
+ value_type operator*()
+ {
+ return value_type{counter, *it};
+ }
+
+ size_t index()
+ {
+ return counter;
+ }
+
+ private:
+ typename T::iterator it;
+ size_t counter;
+ };
+
+ enumerate_impl(T & t) : container(t) {}
+
+ iterator begin()
+ {
+ return iterator(container.begin());
+ }
+
+ iterator end()
+ {
+ return iterator(container.end());
+ }
+
+private:
+ T & container;
+};
+
+// A templated free function allows you to create the wrapper class
+// conveniently
+template
+enumerate_impl enumerate(T & t)
+{
+ return enumerate_impl(t);
+}
+
+class Imu_Integration{
+public:
+ Imu_Integration(Sophus::SE3& SE_init);
+ ~Imu_Integration();
+ bool reset(gtsam::ISAM2& optimizer,boost::shared_ptr& new_frame);
+ bool reset(gtsam::LevenbergMarquardtOptimizer& optimizer,boost::shared_ptr& new_frame);
+ bool update(double* imu= nullptr);
+ bool predict(boost::shared_ptr&,std::size_t&,const double reproj_thresh);
+private:
+ std::shared_ptr statePtr;
+ std::shared_ptr imu_biasPtr;
+ std::shared_ptr valuesPtr;
+ std::shared_ptr graphPtr;
+ boost::shared_ptr parameterPtr;
+ std::shared_ptr preintegratedPtr;
+ boost::shared_ptr ProjectNoisePtr;
+ std::chrono::steady_clock::time_point t_1;
+ std::uint32_t imu_factor_id=0;
+ uint imu_n=0;
+ std::shared_ptr optimizerParamPtr;
+ bool syn=false;
+};
+#endif //SVO_IMU_INTEGRATION_H
diff --git a/CPU_version/svo/include/svo/initialization.h b/CPU_version/svo/include/svo/initialization.h
new file mode 100644
index 0000000..cb562e5
--- /dev/null
+++ b/CPU_version/svo/include/svo/initialization.h
@@ -0,0 +1,83 @@
+// This file is part of SVO - Semi-direct Visual Odometry.
+//
+// Copyright (C) 2014 Christian Forster
+// (Robotics and Perception Group, University of Zurich, Switzerland).
+//
+// SVO is free software: you can redistribute it and/or modify it under the
+// terms of the GNU General Public License as published by the Free Software
+// Foundation, either version 3 of the License, or any later version.
+//
+// SVO is distributed in the hope that it will be useful, but WITHOUT ANY
+// WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
+// FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.
+//
+// You should have received a copy of the GNU General Public License
+// along with this program. If not, see .
+
+#ifndef SVO_INITIALIZATION_H
+#define SVO_INITIALIZATION_H
+
+#include
+
+namespace svo {
+
+class FrameHandlerMono;
+
+/// Bootstrapping the map from the first two views.
+namespace initialization {
+
+enum InitResult { FAILURE, NO_KEYFRAME, SUCCESS };
+
+/// Tracks features using Lucas-Kanade tracker and then estimates a homography.
+class KltHomographyInit {
+ friend class svo::FrameHandlerMono;
+public:
+ EIGEN_MAKE_ALIGNED_OPERATOR_NEW
+
+ FramePtr frame_ref_;
+ KltHomographyInit() {};
+ ~KltHomographyInit() {};
+ InitResult addFirstFrame(FramePtr frame_ref);
+ InitResult addSecondFrame(FramePtr frame_ref);
+ void reset();
+
+protected:
+ vector px_ref_; //!< keypoints to be tracked in reference frame.
+ vector px_cur_; //!< tracked keypoints in current frame.
+ vector f_ref_; //!< bearing vectors corresponding to the keypoints in the reference image.
+ vector f_cur_; //!< bearing vectors corresponding to the keypoints in the current image.
+ vector disparities_; //!< disparity between first and second frame.
+ vector inliers_; //!< inliers after the geometric check (e.g., Homography).
+ vector xyz_in_cur_; //!< 3D points computed during the geometric check.
+ SE3 T_cur_from_ref_; //!< computed transformation between the first two frames.
+};
+
+/// Detect Fast corners in the image.
+void detectFeatures(
+ FramePtr frame,
+ vector& px_vec,
+ vector& f_vec);
+
+/// Compute optical flow (Lucas Kanade) for selected keypoints.
+void trackKlt(
+ FramePtr frame_ref,
+ FramePtr frame_cur,
+ vector& px_ref,
+ vector& px_cur,
+ vector& f_ref,
+ vector& f_cur,
+ vector& disparities);
+
+void computeHomography(
+ const vector& f_ref,
+ const vector& f_cur,
+ double focal_length,
+ double reprojection_threshold,
+ vector& inliers,
+ vector& xyz_in_cur,
+ SE3& T_cur_from_ref);
+
+} // namespace initialization
+} // namespace svo
+
+#endif // SVO_INITIALIZATION_H
diff --git a/CPU_version/svo/include/svo/map.h b/CPU_version/svo/include/svo/map.h
new file mode 100644
index 0000000..458622b
--- /dev/null
+++ b/CPU_version/svo/include/svo/map.h
@@ -0,0 +1,136 @@
+// This file is part of SVO - Semi-direct Visual Odometry.
+//
+// Copyright (C) 2014 Christian Forster
+// (Robotics and Perception Group, University of Zurich, Switzerland).
+//
+// SVO is free software: you can redistribute it and/or modify it under the
+// terms of the GNU General Public License as published by the Free Software
+// Foundation, either version 3 of the License, or any later version.
+//
+// SVO is distributed in the hope that it will be useful, but WITHOUT ANY
+// WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
+// FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.
+//
+// You should have received a copy of the GNU General Public License
+// along with this program. If not, see .
+
+#ifndef SVO_MAP_H_
+#define SVO_MAP_H_
+
+#include
+#include
+#include
+#include
+
+namespace svo {
+
+class Point;
+class Feature;
+class Seed;
+
+/// Container for converged 3D points that are not already assigned to two keyframes.
+class MapPointCandidates
+{
+public:
+ typedef pair PointCandidate;
+ typedef list PointCandidateList;
+
+ /// The depth-filter is running in a parallel thread and fills the canidate list.
+ /// This mutex controls concurrent access to point_candidates.
+ boost::mutex mut_;
+
+ /// Candidate points are created from converged seeds.
+ /// Until the next keyframe, these points can be used for reprojection and pose optimization.
+ PointCandidateList candidates_;
+ list< Point* > trash_points_;
+
+ MapPointCandidates();
+ ~MapPointCandidates();
+
+ /// Add a candidate point.
+ void newCandidatePoint(Point* point, double depth_sigma2);
+
+ /// Adds the feature to the frame and deletes candidate from list.
+ void addCandidatePointToFrame(FramePtr frame);
+
+ /// Remove a candidate point from the list of candidates.
+ bool deleteCandidatePoint(Point* point);
+
+ /// Remove all candidates that belong to a frame.
+ void removeFrameCandidates(FramePtr frame);
+
+ /// Reset the candidate list, remove and delete all points.
+ void reset();
+
+ void deleteCandidate(PointCandidate& c);
+
+ void emptyTrash();
+};
+
+/// Map object which saves all keyframes which are in a map.
+class Map : boost::noncopyable
+{
+public:
+ list< FramePtr > keyframes_; //!< List of keyframes in the map.
+ list< Point* > trash_points_; //!< A deleted point is moved to the trash bin. Now and then this is cleaned. One reason is that the visualizer must remove the points also.
+ MapPointCandidates point_candidates_;
+
+ Map();
+ ~Map();
+
+ /// Reset the map. Delete all keyframes and reset the frame and point counters.
+ void reset();
+
+ /// Delete a point in the map and remove all references in keyframes to it.
+ void safeDeletePoint(Point* pt);
+
+ /// Moves the point to the trash queue which is cleaned now and then.
+ void deletePoint(Point* pt);
+
+ /// Moves the frame to the trash queue which is cleaned now and then.
+ bool safeDeleteFrame(FramePtr frame);
+
+ /// Remove the references between a point and a frame.
+ void removePtFrameRef(Frame* frame, Feature* ftr);
+
+ /// Add a new keyframe to the map.
+ void addKeyframe(FramePtr new_keyframe);
+
+ /// Given a frame, return all keyframes which have an overlapping field of view.
+ void getCloseKeyframes(const FramePtr& frame, list< pair >& close_kfs) const;
+
+ /// Return the keyframe which is spatially closest and has overlapping field of view.
+ FramePtr getClosestKeyframe(const FramePtr& frame) const;
+
+ /// Return the keyframe which is furthest apart from pos.
+ FramePtr getFurthestKeyframe(const Vector3d& pos) const;
+
+ bool getKeyframeById(const int id, FramePtr& frame) const;
+
+ /// Transform the whole map with rotation R, translation t and scale s.
+ void transform(const Matrix3d& R, const Vector3d& t, const double& s);
+
+ /// Empty trash bin of deleted keyframes and map points. We don't delete the
+ /// points immediately to ensure proper cleanup and to provide the visualizer
+ /// a list of objects which must be removed.
+ void emptyTrash();
+
+ /// Return the keyframe which was last inserted in the map.
+ inline FramePtr lastKeyframe() { return keyframes_.back(); }
+
+ /// Return the number of keyframes in the map
+ inline size_t size() const { return keyframes_.size(); }
+};
+
+/// A collection of debug functions to check the data consistency.
+namespace map_debug {
+
+void mapStatistics(Map* map);
+void mapValidation(Map* map, int id);
+void frameValidation(Frame* frame, int id);
+void pointValidation(Point* point, int id);
+
+} // namespace map_debug
+} // namespace svo
+
+#endif // SVO_MAP_H_
diff --git a/CPU_version/svo/include/svo/matcher.h b/CPU_version/svo/include/svo/matcher.h
new file mode 100644
index 0000000..b53340e
--- /dev/null
+++ b/CPU_version/svo/include/svo/matcher.h
@@ -0,0 +1,130 @@
+// This file is part of SVO - Semi-direct Visual Odometry.
+//
+// Copyright (C) 2014 Christian Forster
+// (Robotics and Perception Group, University of Zurich, Switzerland).
+//
+// SVO is free software: you can redistribute it and/or modify it under the
+// terms of the GNU General Public License as published by the Free Software
+// Foundation, either version 3 of the License, or any later version.
+//
+// SVO is distributed in the hope that it will be useful, but WITHOUT ANY
+// WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
+// FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.
+//
+// You should have received a copy of the GNU General Public License
+// along with this program. If not, see .
+
+#ifndef SVO_MATCHER_H_
+#define SVO_MATCHER_H_
+
+#include
+
+namespace vk {
+ class AbstractCamera;
+ namespace patch_score {
+ template class ZMSSD;
+ }
+}
+
+namespace svo {
+
+class Point;
+class Frame;
+class Feature;
+
+/// Warp a patch from the reference view to the current view.
+namespace warp {
+
+void getWarpMatrixAffine(
+ const vk::AbstractCamera& cam_ref,
+ const vk::AbstractCamera& cam_cur,
+ const Vector2d& px_ref,
+ const Vector3d& f_ref,
+ const double depth_ref,
+ const SE3& T_cur_ref,
+ const int level_ref,
+ Matrix2d& A_cur_ref);
+
+int getBestSearchLevel(
+ const Matrix2d& A_cur_ref,
+ const int max_level);
+
+void warpAffine(
+ const Matrix2d& A_cur_ref,
+ const cv::Mat& img_ref,
+ const Vector2d& px_ref,
+ const int level_ref,
+ const int level_cur,
+ const int halfpatch_size,
+ uint8_t* patch);
+
+} // namespace warp
+
+/// Patch-matcher for reprojection-matching and epipolar search in triangulation.
+class Matcher
+{
+public:
+ EIGEN_MAKE_ALIGNED_OPERATOR_NEW
+
+ static const int halfpatch_size_ = 4;
+ static const int patch_size_ = 8;
+
+ typedef vk::patch_score::ZMSSD PatchScore;
+
+ struct Options
+ {
+ bool align_1d; //!< in epipolar search: align patch 1D along epipolar line
+ int align_max_iter; //!< number of iterations for aligning the feature patches in gauss newton
+ double max_epi_length_optim;//!< max length of epipolar line to skip epipolar search and directly go to img align
+ size_t max_epi_search_steps;//!< max number of evaluations along epipolar line
+ bool subpix_refinement; //!< do gauss newton feature patch alignment after epipolar search
+ bool epi_search_edgelet_filtering;
+ double epi_search_edgelet_max_angle;
+ Options() :
+ align_1d(false),
+ align_max_iter(10),
+ max_epi_length_optim(2.0),
+ max_epi_search_steps(1000),
+ subpix_refinement(true),
+ epi_search_edgelet_filtering(true),
+ epi_search_edgelet_max_angle(0.7)
+ {}
+ } options_;
+
+ uint8_t patch_[patch_size_*patch_size_] __attribute__ ((aligned (16)));
+ uint8_t patch_with_border_[(patch_size_+2)*(patch_size_+2)] __attribute__ ((aligned (16)));
+ Matrix2d A_cur_ref_; //!< affine warp matrix
+ Vector2d epi_dir_;
+ double epi_length_; //!< length of epipolar line segment in pixels (only used for epipolar search)
+ double h_inv_; //!< hessian of 1d image alignment along epipolar line
+ int search_level_;
+ bool reject_;
+ Feature* ref_ftr_;
+ Vector2d px_cur_;
+
+ Matcher() = default;
+ ~Matcher() = default;
+
+ /// Find a match by directly applying subpix refinement.
+ /// IMPORTANT! This function assumes that px_cur is already set to an estimate that is within ~2-3 pixel of the final result!
+ bool findMatchDirect(
+ const Point& pt,
+ const Frame& frame,
+ Vector2d& px_cur);
+
+ /// Find a match by searching along the epipolar line without using any features.
+ bool findEpipolarMatchDirect(
+ const Frame& ref_frame,
+ const Frame& cur_frame,
+ const Feature& ref_ftr,
+ const double d_estimate,
+ const double d_min,
+ const double d_max,
+ double& depth);
+
+ void createPatchFromPatchWithBorder();
+};
+
+} // namespace svo
+
+#endif // SVO_MATCHER_H_
diff --git a/CPU_version/svo/include/svo/point.h b/CPU_version/svo/include/svo/point.h
new file mode 100644
index 0000000..e1bb81b
--- /dev/null
+++ b/CPU_version/svo/include/svo/point.h
@@ -0,0 +1,108 @@
+// This file is part of SVO - Semi-direct Visual Odometry.
+//
+// Copyright (C) 2014 Christian Forster
+// (Robotics and Perception Group, University of Zurich, Switzerland).
+//
+// SVO is free software: you can redistribute it and/or modify it under the
+// terms of the GNU General Public License as published by the Free Software
+// Foundation, either version 3 of the License, or any later version.
+//
+// SVO is distributed in the hope that it will be useful, but WITHOUT ANY
+// WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
+// FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.
+//
+// You should have received a copy of the GNU General Public License
+// along with this program. If not, see .
+
+#ifndef SVO_POINT_H_
+#define SVO_POINT_H_
+
+#include
+#include
+
+namespace g2o {
+ class VertexSBAPointXYZ;
+}
+typedef g2o::VertexSBAPointXYZ g2oPoint;
+
+namespace svo {
+
+class Feature;
+
+typedef Matrix Matrix23d;
+
+/// A 3D point on the surface of the scene.
+class Point : boost::noncopyable
+{
+public:
+ EIGEN_MAKE_ALIGNED_OPERATOR_NEW
+
+ enum PointType {
+ TYPE_DELETED,
+ TYPE_CANDIDATE,
+ TYPE_UNKNOWN,
+ TYPE_GOOD
+ };
+
+ static int point_counter_; //!< Counts the number of created points. Used to set the unique id.
+ int id_; //!< Unique ID of the point.
+ Vector3d pos_; //!< 3d pos of the point in the world coordinate frame.
+ Vector3d normal_; //!< Surface normal at point.
+ Matrix3d normal_information_; //!< Inverse covariance matrix of normal estimation.
+ bool normal_set_; //!< Flag whether the surface normal was estimated or not.
+ list obs_; //!< References to keyframes which observe the point.
+ size_t n_obs_; //!< Number of obervations: Keyframes AND successful reprojections in intermediate frames.
+ g2oPoint* v_pt_; //!< Temporary pointer to the point-vertex in g2o during bundle adjustment.
+ int last_published_ts_; //!< Timestamp of last publishing.
+ int last_projected_kf_id_; //!< Flag for the reprojection: don't reproject a pt twice.
+ PointType type_; //!< Quality of the point.
+ int n_failed_reproj_; //!< Number of failed reprojections. Used to assess the quality of the point.
+ int n_succeeded_reproj_; //!< Number of succeeded reprojections. Used to assess the quality of the point.
+ int last_structure_optim_; //!< Timestamp of last point optimization
+
+ Point(const Vector3d& pos);
+ Point(const Vector3d& pos, Feature* ftr);
+ ~Point();
+
+ /// Add a reference to a frame.
+ void addFrameRef(Feature* ftr);
+
+ /// Remove reference to a frame.
+ bool deleteFrameRef(Frame* frame);
+
+ /// Initialize point normal. The inital estimate will point towards the frame.
+ void initNormal();
+
+ /// Check whether mappoint has reference to a frame.
+ Feature* findFrameRef(Frame* frame);
+
+ /// Get Frame with similar viewpoint.
+ bool getCloseViewObs(const Vector3d& pos, Feature*& obs) const;
+
+ /// Get number of observations.
+ inline size_t nRefs() const { return obs_.size(); }
+
+ /// Optimize point position through minimizing the reprojection error.
+ void optimize(const size_t n_iter);
+
+ /// Jacobian of point projection on unit plane (focal length = 1) in frame (f).
+ inline static void jacobian_xyz2uv(
+ const Vector3d& p_in_f,
+ const Matrix3d& R_f_w,
+ Matrix23d& point_jac)
+ {
+ const double z_inv = 1.0/p_in_f[2];
+ const double z_inv_sq = z_inv*z_inv;
+ point_jac(0, 0) = z_inv;
+ point_jac(0, 1) = 0.0;
+ point_jac(0, 2) = -p_in_f[0] * z_inv_sq;
+ point_jac(1, 0) = 0.0;
+ point_jac(1, 1) = z_inv;
+ point_jac(1, 2) = -p_in_f[1] * z_inv_sq;
+ point_jac = - point_jac * R_f_w;
+ }
+};
+
+} // namespace svo
+
+#endif // SVO_POINT_H_
diff --git a/CPU_version/svo/include/svo/pose_optimizer.h b/CPU_version/svo/include/svo/pose_optimizer.h
new file mode 100644
index 0000000..ca048ae
--- /dev/null
+++ b/CPU_version/svo/include/svo/pose_optimizer.h
@@ -0,0 +1,50 @@
+// This file is part of SVO - Semi-direct Visual Odometry.
+//
+// Copyright (C) 2014 Christian Forster
+// (Robotics and Perception Group, University of Zurich, Switzerland).
+//
+// SVO is free software: you can redistribute it and/or modify it under the
+// terms of the GNU General Public License as published by the Free Software
+// Foundation, either version 3 of the License, or any later version.
+//
+// SVO is distributed in the hope that it will be useful, but WITHOUT ANY
+// WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
+// FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.
+//
+// You should have received a copy of the GNU General Public License
+// along with this program. If not, see .
+
+#ifndef SVO_POSE_OPTIMIZER_H_
+#define SVO_POSE_OPTIMIZER_H_
+
+#include
+
+namespace svo {
+
+using namespace Eigen;
+using namespace Sophus;
+using namespace std;
+
+typedef Matrix Matrix6d;
+typedef Matrix Matrix26d;
+typedef Matrix Vector6d;
+
+class Point;
+
+/// Motion-only bundle adjustment. Minimize the reprojection error of a single frame.
+namespace pose_optimizer {
+
+void optimizeGaussNewton(
+ const double reproj_thresh,
+ const size_t n_iter,
+ const bool verbose,
+ FramePtr& frame,
+ double& estimated_scale,
+ double& error_init,
+ double& error_final,
+ size_t& num_obs);
+
+} // namespace pose_optimizer
+} // namespace svo
+
+#endif // SVO_POSE_OPTIMIZER_H_
diff --git a/CPU_version/svo/include/svo/reprojector.h b/CPU_version/svo/include/svo/reprojector.h
new file mode 100644
index 0000000..e345f0d
--- /dev/null
+++ b/CPU_version/svo/include/svo/reprojector.h
@@ -0,0 +1,101 @@
+// This file is part of SVO - Semi-direct Visual Odometry.
+//
+// Copyright (C) 2014 Christian Forster
+// (Robotics and Perception Group, University of Zurich, Switzerland).
+//
+// SVO is free software: you can redistribute it and/or modify it under the
+// terms of the GNU General Public License as published by the Free Software
+// Foundation, either version 3 of the License, or any later version.
+//
+// SVO is distributed in the hope that it will be useful, but WITHOUT ANY
+// WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
+// FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.
+//
+// You should have received a copy of the GNU General Public License
+// along with this program. If not, see .
+
+#ifndef SVO_REPROJECTION_H_
+#define SVO_REPROJECTION_H_
+
+#include
+#include
+
+namespace vk {
+class AbstractCamera;
+}
+
+namespace svo {
+
+class Map;
+class Point;
+
+/// Project points from the map into the image and find the corresponding
+/// feature (corner). We don't search a match for every point but only for one
+/// point per cell. Thereby, we achieve a homogeneously distributed set of
+/// matched features and at the same time we can save processing time by not
+/// projecting all points.
+class Reprojector
+{
+public:
+ EIGEN_MAKE_ALIGNED_OPERATOR_NEW
+
+ /// Reprojector config parameters
+ struct Options {
+ size_t max_n_kfs; //!< max number of keyframes to reproject from
+ bool find_match_direct;
+ Options()
+ : max_n_kfs(10),
+ find_match_direct(true)
+ {}
+ } options_;
+
+ size_t n_matches_;
+ size_t n_trials_;
+
+ Reprojector(vk::AbstractCamera* cam, Map& map);
+
+ ~Reprojector();
+
+ /// Project points from the map into the image. First finds keyframes with
+ /// overlapping field of view and projects only those map-points.
+ void reprojectMap(
+ FramePtr frame,
+ std::vector< std::pair >& overlap_kfs);
+
+private:
+
+ /// A candidate is a point that projects into the image plane and for which we
+ /// will search a maching feature in the image.
+ struct Candidate {
+ EIGEN_MAKE_ALIGNED_OPERATOR_NEW
+ Point* pt; //!< 3D point.
+ Vector2d px; //!< projected 2D pixel location.
+ Candidate(Point* pt, Vector2d& px) : pt(pt), px(px) {}
+ };
+ typedef std::list Cell;
+ typedef std::vector| CandidateGrid;
+
+ /// The grid stores a set of candidate matches. For every grid cell we try to find one match.
+ struct Grid
+ {
+ CandidateGrid cells;
+ vector cell_order;
+ int cell_size;
+ int grid_n_cols;
+ int grid_n_rows;
+ };
+
+ Grid grid_;
+ Matcher matcher_;
+ Map& map_;
+
+ static bool pointQualityComparator(Candidate& lhs, Candidate& rhs);
+ void initializeGrid(vk::AbstractCamera* cam);
+ void resetGrid();
+ bool reprojectCell(Cell& cell, FramePtr frame);
+ bool reprojectPoint(FramePtr frame, Point* point);
+};
+
+} // namespace svo
+
+#endif // SVO_REPROJECTION_H_
diff --git a/CPU_version/svo/include/svo/sparse_img_align.h b/CPU_version/svo/include/svo/sparse_img_align.h
new file mode 100644
index 0000000..1bedfed
--- /dev/null
+++ b/CPU_version/svo/include/svo/sparse_img_align.h
@@ -0,0 +1,83 @@
+// This file is part of SVO - Semi-direct Visual Odometry.
+//
+// Copyright (C) 2014 Christian Forster
+// (Robotics and Perception Group, University of Zurich, Switzerland).
+//
+// SVO is free software: you can redistribute it and/or modify it under the
+// terms of the GNU General Public License as published by the Free Software
+// Foundation, either version 3 of the License, or any later version.
+//
+// SVO is distributed in the hope that it will be useful, but WITHOUT ANY
+// WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
+// FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.
+//
+// You should have received a copy of the GNU General Public License
+// along with this program. If not, see .
+
+#ifndef SVO_SPARSE_IMG_ALIGN_H_
+#define SVO_SPARSE_IMG_ALIGN_H_
+
+#include
+#include
+#include
+
+namespace vk {
+class AbstractCamera;
+}
+
+namespace svo {
+
+class Feature;
+
+/// Optimize the pose of the frame by minimizing the photometric error of feature patches.
+class SparseImgAlign : public vk::NLLSSolver<6, SE3>
+{
+ static const int patch_halfsize_ = 2;
+ static const int patch_size_ = 2*patch_halfsize_;
+ static const int patch_area_ = patch_size_*patch_size_;
+public:
+ EIGEN_MAKE_ALIGNED_OPERATOR_NEW
+
+ cv::Mat resimg_;
+
+ SparseImgAlign(
+ int n_levels,
+ int min_level,
+ int n_iter,
+ Method method,
+ bool display,
+ bool verbose);
+
+ size_t run(
+ FramePtr ref_frame,
+ FramePtr cur_frame);
+
+ /// Return fisher information matrix, i.e. the Hessian of the log-likelihood
+ /// at the converged state.
+ Matrix getFisherInformation();
+
+protected:
+ FramePtr ref_frame_; //!< reference frame, has depth for gradient pixels.
+ FramePtr cur_frame_; //!< only the image is known!
+ int level_; //!< current pyramid level on which the optimization runs.
+ bool display_; //!< display residual image.
+ int max_level_; //!< coarsest pyramid level for the alignment.
+ int min_level_; //!< finest pyramid level for the alignment.
+
+ // cache:
+ Matrix jacobian_cache_;
+ bool have_ref_patch_cache_;
+ cv::Mat ref_patch_cache_;
+ std::vector visible_fts_;
+
+ void precomputeReferencePatches();
+ virtual double computeResiduals(const SE3& model, bool linearize_system, bool compute_weight_scale = false);
+ virtual int solve();
+ virtual void update (const ModelType& old_model, ModelType& new_model);
+ virtual void startIteration();
+ virtual void finishIteration();
+};
+
+} // namespace svo
+
+#endif // SVO_SPARSE_IMG_ALIGN_H_
diff --git a/CPU_version/svo/include/svo/svoVisionFactor.h b/CPU_version/svo/include/svo/svoVisionFactor.h
new file mode 100644
index 0000000..cd00816
--- /dev/null
+++ b/CPU_version/svo/include/svo/svoVisionFactor.h
@@ -0,0 +1,70 @@
+//
+// Created by root on 4/13/21.
+//
+
+#ifndef SVO_SVOVISIONFACTOR_H
+#define SVO_SVOVISIONFACTOR_H
+#pragma once
+
+#include
+#include
+#include
+#include
+#include
+#include
+#include
+
+namespace svo {
+
+ class VisionFactor: public gtsam::NoiseModelFactor1 {
+
+ private:
+ // measurement information
+ gtsam::Point3 feature; /// body_P_sensor_; ///< TODO The pose of the sensor in the body frame
+ int level; ///< Image pyramid level where feature was extracted.
+
+ public:
+ /// shorthand for a smart pointer to a factor
+ typedef boost::shared_ptr shared_ptr;
+
+ /**
+ * Constructor
+ * @param feature Unit-bearing vector of the feature.
+ * @param model is the standard deviation
+ * @param pointKey is the index of the landmark
+ * @param poseKey is the index of the camera pose
+ * @param landmark 3D position of the landmark
+ * @param body_P_sensor is the transform from body to sensor frame (default identity)
+ * @param level Image pyramid level where feature was extracted.
+ */
+ VisionFactor(const gtsam::Point3& feature, const gtsam::Point3& landmark,const int level, const gtsam::SharedNoiseModel& Hp,
+ gtsam::Key poseKey,boost::optional body_P_sensor = boost::none) :
+ gtsam::NoiseModelFactor1(Hp, poseKey), feature(feature), landmark(landmark),body_P_sensor_(body_P_sensor),level(level){}
+ virtual ~VisionFactor(){}
+
+ /// error function
+ /// @param p the Estimated pose in Pose3
+ /// @param H the optional Jacobian matrix, which use boost optional and has default null pointer
+ gtsam::Vector evaluateError(const gtsam::Pose3& p, boost::optional Hp = boost::none) const {
+ // note that use boost optional like a pointer
+ // only calculate jacobian matrix when non-null pointer exists
+ gtsam::Matrix26 J;
+ gtsam::Point3 point(p*landmark);
+ svo::Frame::jacobian_xyz2uv(point, J);
+ Vector2d e = (vk::project2d(feature) - vk::project2d(point))*1.0 / (1<(
+ gtsam::NonlinearFactor::shared_ptr(new VisionFactor(*this))); }
+
+ };
+
+} // namespace gtsamexamples
+#endif //SVO_SVOVISIONFACTOR_H
diff --git a/CPU_version/svo/package.xml b/CPU_version/svo/package.xml
new file mode 100644
index 0000000..4d71cb3
--- /dev/null
+++ b/CPU_version/svo/package.xml
@@ -0,0 +1,27 @@
+
+
+ svo
+ 0.1.0
+
+ Monocular Visual Odometry
+
+ Christian Forster
+ Christian Forster
+ GPLv3
+
+
+ catkin
+
+
+ cmake_modules
+ roscpp
+ roslib
+ vikit_common
+ vikit_ros
+
+
+ roscpp
+ roslib
+ vikit_common
+ vikit_ros
+
\ No newline at end of file
diff --git a/CPU_version/svo/src/bundle_adjustment.cpp b/CPU_version/svo/src/bundle_adjustment.cpp
new file mode 100644
index 0000000..eaba2b0
--- /dev/null
+++ b/CPU_version/svo/src/bundle_adjustment.cpp
@@ -0,0 +1,435 @@
+// This file is part of SVO - Semi-direct Visual Odometry.
+//
+// Copyright (C) 2014 Christian Forster
+// (Robotics and Perception Group, University of Zurich, Switzerland).
+//
+// SVO is free software: you can redistribute it and/or modify it under the
+// terms of the GNU General Public License as published by the Free Software
+// Foundation, either version 3 of the License, or any later version.
+//
+// SVO is distributed in the hope that it will be useful, but WITHOUT ANY
+// WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
+// FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.
+//
+// You should have received a copy of the GNU General Public License
+// along with this program. If not, see .
+
+#include
+#include
+#include
+#include
+#include
+#include
+#include
+#include
+#include
+#include
+#include
+#include
+#include
+#include
+#include
+#include
+#include
+#include
+
+#define SCHUR_TRICK 1
+
+namespace svo {
+namespace ba {
+
+void twoViewBA(
+ Frame* frame1,
+ Frame* frame2,
+ double reproj_thresh,
+ Map* map)
+{
+ // scale reprojection threshold in pixels to unit plane
+ reproj_thresh /= frame1->cam_->errorMultiplier2();
+
+ // init g2o
+ g2o::SparseOptimizer optimizer;
+ setupG2o(&optimizer);
+
+ list edges;
+ size_t v_id = 0;
+
+ // New Keyframe Vertex 1: This Keyframe is set to fixed!
+ g2oFrameSE3* v_frame1 = createG2oFrameSE3(frame1, v_id++, true);
+ optimizer.addVertex(v_frame1);
+
+ // New Keyframe Vertex 2
+ g2oFrameSE3* v_frame2 = createG2oFrameSE3(frame2, v_id++, false);
+ optimizer.addVertex(v_frame2);
+
+ // Create Point Vertices
+ for(Features::iterator it_ftr=frame1->fts_.begin(); it_ftr!=frame1->fts_.end(); ++it_ftr)
+ {
+ Point* pt = (*it_ftr)->point;
+ if(pt == NULL)
+ continue;
+ g2oPoint* v_pt = createG2oPoint(pt->pos_, v_id++, false);
+ optimizer.addVertex(v_pt);
+ pt->v_pt_ = v_pt;
+ g2oEdgeSE3* e = createG2oEdgeSE3(v_frame1, v_pt, vk::project2d((*it_ftr)->f), true, reproj_thresh*Config::lobaRobustHuberWidth());
+ optimizer.addEdge(e);
+ edges.push_back(EdgeContainerSE3(e, frame1, *it_ftr)); // TODO feature now links to frame, so we can simplify edge container!
+
+ // find at which index the second frame observes the point
+ Feature* ftr_frame2 = pt->findFrameRef(frame2);
+ e = createG2oEdgeSE3(v_frame2, v_pt, vk::project2d(ftr_frame2->f), true, reproj_thresh*Config::lobaRobustHuberWidth());
+ optimizer.addEdge(e);
+ edges.push_back(EdgeContainerSE3(e, frame2, ftr_frame2));
+ }
+
+ // Optimization
+ double init_error, final_error;
+ runSparseBAOptimizer(&optimizer, Config::lobaNumIter(), init_error, final_error);
+ printf("2-View BA: Error before/after = %f / %f\n", init_error, final_error);
+
+ // Update Keyframe Positions
+ frame1->T_f_w_.rotation_matrix() = v_frame1->estimate().rotation().toRotationMatrix();
+ frame1->T_f_w_.translation() = v_frame1->estimate().translation();
+ frame2->T_f_w_.rotation_matrix() = v_frame2->estimate().rotation().toRotationMatrix();
+ frame2->T_f_w_.translation() = v_frame2->estimate().translation();
+
+ // Update Mappoint Positions
+ for(Features::iterator it=frame1->fts_.begin(); it!=frame1->fts_.end(); ++it)
+ {
+ if((*it)->point == NULL)
+ continue;
+ (*it)->point->pos_ = (*it)->point->v_pt_->estimate();
+ (*it)->point->v_pt_ = NULL;
+ }
+
+ // Find Mappoints with too large reprojection error
+ const double reproj_thresh_squared = reproj_thresh*reproj_thresh;
+ size_t n_incorrect_edges = 0;
+ for(list::iterator it_e = edges.begin(); it_e != edges.end(); ++it_e)
+ if(it_e->edge->chi2() > reproj_thresh_squared)
+ {
+ if(it_e->feature->point != NULL)
+ {
+ map->safeDeletePoint(it_e->feature->point);
+ it_e->feature->point = NULL;
+ }
+ ++n_incorrect_edges;
+ }
+
+ printf("2-View BA: Wrong edges = %zu\n", n_incorrect_edges);
+}
+
+void localBA(
+ Frame* center_kf,
+ set* core_kfs,
+ Map* map,
+ size_t& n_incorrect_edges_1,
+ size_t& n_incorrect_edges_2,
+ double& init_error,
+ double& final_error)
+{
+
+ // init g2o
+ g2o::SparseOptimizer optimizer;
+ setupG2o(&optimizer);
+
+ list edges;
+ set mps;
+ list neib_kfs;
+ size_t v_id = 0;
+ size_t n_mps = 0;
+ size_t n_fix_kfs = 0;
+ size_t n_var_kfs = 1;
+ size_t n_edges = 0;
+ n_incorrect_edges_1 = 0;
+ n_incorrect_edges_2 = 0;
+
+ // Add all core keyframes
+ for(set::iterator it_kf = core_kfs->begin(); it_kf != core_kfs->end(); ++it_kf)
+ {
+ g2oFrameSE3* v_kf = createG2oFrameSE3(it_kf->get(), v_id++, false);
+ (*it_kf)->v_kf_ = v_kf;
+ ++n_var_kfs;
+ assert(optimizer.addVertex(v_kf));
+
+ // all points that the core keyframes observe are also optimized:
+ for(Features::iterator it_pt=(*it_kf)->fts_.begin(); it_pt!=(*it_kf)->fts_.end(); ++it_pt)
+ if((*it_pt)->point != NULL)
+ mps.insert((*it_pt)->point);
+ }
+
+ // Now go throug all the points and add a measurement. Add a fixed neighbour
+ // Keyframe if it is not in the set of core kfs
+ double reproj_thresh_2 = Config::lobaThresh() / center_kf->cam_->errorMultiplier2();
+ double reproj_thresh_1 = Config::poseOptimThresh() / center_kf->cam_->errorMultiplier2();
+ double reproj_thresh_1_squared = reproj_thresh_1*reproj_thresh_1;
+ for(set::iterator it_pt = mps.begin(); it_pt!=mps.end(); ++it_pt)
+ {
+ // Create point vertex
+ g2oPoint* v_pt = createG2oPoint((*it_pt)->pos_, v_id++, false);
+ (*it_pt)->v_pt_ = v_pt;
+ assert(optimizer.addVertex(v_pt));
+ ++n_mps;
+
+ // Add edges
+ list::iterator it_obs=(*it_pt)->obs_.begin();
+ while(it_obs!=(*it_pt)->obs_.end())
+ {
+ Vector2d error = vk::project2d((*it_obs)->f) - vk::project2d((*it_obs)->frame->w2f((*it_pt)->pos_));
+
+ if((*it_obs)->frame->v_kf_ == NULL)
+ {
+ // frame does not have a vertex yet -> it belongs to the neib kfs and
+ // is fixed. create one:
+ g2oFrameSE3* v_kf = createG2oFrameSE3((*it_obs)->frame, v_id++, true);
+ (*it_obs)->frame->v_kf_ = v_kf;
+ ++n_fix_kfs;
+ assert(optimizer.addVertex(v_kf));
+ neib_kfs.push_back((*it_obs)->frame);
+ }
+
+ // create edge
+ g2oEdgeSE3* e = createG2oEdgeSE3((*it_obs)->frame->v_kf_, v_pt,
+ vk::project2d((*it_obs)->f),
+ true,
+ reproj_thresh_2*Config::lobaRobustHuberWidth(),
+ 1.0 / (1<<(*it_obs)->level));
+ assert(optimizer.addEdge(e));
+ edges.push_back(EdgeContainerSE3(e, (*it_obs)->frame, *it_obs));
+ ++n_edges;
+ ++it_obs;
+ }
+ }
+
+ // structure only
+ g2o::StructureOnlySolver<3> structure_only_ba;
+ g2o::OptimizableGraph::VertexContainer points;
+ for (g2o::OptimizableGraph::VertexIDMap::const_iterator it = optimizer.vertices().begin(); it != optimizer.vertices().end(); ++it)
+ {
+ g2o::OptimizableGraph::Vertex* v = static_cast(it->second);
+ if (v->dimension() == 3 && v->edges().size() >= 2)
+ points.push_back(v);
+ }
+ structure_only_ba.calc(points, 10);
+
+ // Optimization
+ if(Config::lobaNumIter() > 0)
+ runSparseBAOptimizer(&optimizer, Config::lobaNumIter(), init_error, final_error);
+
+ // Update Keyframes
+ for(set::iterator it = core_kfs->begin(); it != core_kfs->end(); ++it)
+ {
+ (*it)->T_f_w_ = SE3( (*it)->v_kf_->estimate().rotation(),
+ (*it)->v_kf_->estimate().translation());
+ (*it)->v_kf_ = NULL;
+ }
+
+ for(list::iterator it = neib_kfs.begin(); it != neib_kfs.end(); ++it)
+ (*it)->v_kf_ = NULL;
+
+ // Update Mappoints
+ for(set::iterator it = mps.begin(); it != mps.end(); ++it)
+ {
+ (*it)->pos_ = (*it)->v_pt_->estimate();
+ (*it)->v_pt_ = NULL;
+ }
+
+ // Remove Measurements with too large reprojection error
+ double reproj_thresh_2_squared = reproj_thresh_2*reproj_thresh_2;
+ for(list::iterator it = edges.begin(); it != edges.end(); ++it)
+ {
+ if(it->edge->chi2() > reproj_thresh_2_squared) //*(1<feature_->level))
+ {
+ map->removePtFrameRef(it->frame, it->feature);
+ ++n_incorrect_edges_2;
+ }
+ }
+
+ // TODO: delete points and edges!
+ init_error = sqrt(init_error)*center_kf->cam_->errorMultiplier2();
+ final_error = sqrt(final_error)*center_kf->cam_->errorMultiplier2();
+}
+
+void globalBA(Map* map)
+{
+ // init g2o
+ g2o::SparseOptimizer optimizer;
+ setupG2o(&optimizer);
+
+ list edges;
+ list< pair > incorrect_edges;
+
+ // Go through all Keyframes
+ size_t v_id = 0;
+ double reproj_thresh_2 = Config::lobaThresh() / map->lastKeyframe()->cam_->errorMultiplier2();
+ double reproj_thresh_1_squared = Config::poseOptimThresh() / map->lastKeyframe()->cam_->errorMultiplier2();
+ reproj_thresh_1_squared *= reproj_thresh_1_squared;
+ for(list::iterator it_kf = map->keyframes_.begin();
+ it_kf != map->keyframes_.end(); ++it_kf)
+ {
+ // New Keyframe Vertex
+ g2oFrameSE3* v_kf = createG2oFrameSE3(it_kf->get(), v_id++, false);
+ (*it_kf)->v_kf_ = v_kf;
+ optimizer.addVertex(v_kf);
+ for(Features::iterator it_ftr=(*it_kf)->fts_.begin(); it_ftr!=(*it_kf)->fts_.end(); ++it_ftr)
+ {
+ // for each keyframe add edges to all observed mapoints
+ Point* mp = (*it_ftr)->point;
+ if(mp == NULL)
+ continue;
+ g2oPoint* v_mp = mp->v_pt_;
+ if(v_mp == NULL)
+ {
+ // mappoint-vertex doesn't exist yet. create a new one:
+ v_mp = createG2oPoint(mp->pos_, v_id++, false);
+ mp->v_pt_ = v_mp;
+ optimizer.addVertex(v_mp);
+ }
+
+ // Due to merging of mappoints it is possible that references in kfs suddenly
+ // have a very large reprojection error which may result in distorted results.
+ Vector2d error = vk::project2d((*it_ftr)->f) - vk::project2d((*it_kf)->w2f(mp->pos_));
+ if(error.squaredNorm() > reproj_thresh_1_squared)
+ incorrect_edges.push_back(pair(*it_kf, *it_ftr));
+ else
+ {
+ g2oEdgeSE3* e = createG2oEdgeSE3(v_kf, v_mp, vk::project2d((*it_ftr)->f),
+ true,
+ reproj_thresh_2*Config::lobaRobustHuberWidth());
+
+ edges.push_back(EdgeContainerSE3(e, it_kf->get(), *it_ftr));
+ optimizer.addEdge(e);
+ }
+ }
+ }
+
+ // Optimization
+ double init_error=0.0, final_error=0.0;
+ if(Config::lobaNumIter() > 0)
+ runSparseBAOptimizer(&optimizer, Config::lobaNumIter(), init_error, final_error);
+
+ // Update Keyframe and MapPoint Positions
+ for(list::iterator it_kf = map->keyframes_.begin();
+ it_kf != map->keyframes_.end(); ++it_kf)
+ {
+ (*it_kf)->T_f_w_ = SE3( (*it_kf)->v_kf_->estimate().rotation(),
+ (*it_kf)->v_kf_->estimate().translation());
+ (*it_kf)->v_kf_ = NULL;
+ for(Features::iterator it_ftr=(*it_kf)->fts_.begin(); it_ftr!=(*it_kf)->fts_.end(); ++it_ftr)
+ {
+ Point* mp = (*it_ftr)->point;
+ if(mp == NULL)
+ continue;
+ if(mp->v_pt_ == NULL)
+ continue; // mp was updated before
+ mp->pos_ = mp->v_pt_->estimate();
+ mp->v_pt_ = NULL;
+ }
+ }
+
+ // Remove Measurements with too large reprojection error
+ for(list< pair >::iterator it=incorrect_edges.begin();
+ it!=incorrect_edges.end(); ++it)
+ map->removePtFrameRef(it->first.get(), it->second);
+
+ double reproj_thresh_2_squared = reproj_thresh_2*reproj_thresh_2;
+ for(list::iterator it = edges.begin(); it != edges.end(); ++it)
+ {
+ if(it->edge->chi2() > reproj_thresh_2_squared)
+ {
+ map->removePtFrameRef(it->frame, it->feature);
+ }
+ }
+}
+
+void setupG2o(g2o::SparseOptimizer * optimizer)
+{
+ // TODO: What's happening with all this HEAP stuff? Memory Leak?
+ optimizer->setVerbose(false);
+
+#if SCHUR_TRICK
+ // solver
+ g2o::BlockSolver_6_3::LinearSolverType* linearSolver;
+ linearSolver = new g2o::LinearSolverCholmod();
+ //linearSolver = new g2o::LinearSolverCSparse();
+
+ g2o::BlockSolver_6_3 * solver_ptr = new g2o::BlockSolver_6_3(linearSolver);
+ g2o::OptimizationAlgorithmLevenberg * solver = new g2o::OptimizationAlgorithmLevenberg(solver_ptr);
+#else
+ g2o::BlockSolverX::LinearSolverType * linearSolver;
+ linearSolver = new g2o::LinearSolverCholmod();
+ //linearSolver = new g2o::LinearSolverCSparse();
+ g2o::BlockSolverX * solver_ptr = new g2o::BlockSolverX(linearSolver);
+ g2o::OptimizationAlgorithmLevenberg* solver = new g2o::OptimizationAlgorithmLevenberg(solver_ptr);
+#endif
+
+ solver->setMaxTrialsAfterFailure(5);
+ optimizer->setAlgorithm(solver);
+
+ // setup camera
+ g2o::CameraParameters * cam_params = new g2o::CameraParameters(1.0, Vector2d(0.,0.), 0.);
+ cam_params->setId(0);
+ if (!optimizer->addParameter(cam_params)) {
+ assert(false);
+ }
+}
+
+void
+runSparseBAOptimizer(g2o::SparseOptimizer* optimizer,
+ unsigned int num_iter,
+ double& init_error, double& final_error)
+{
+ optimizer->initializeOptimization();
+ optimizer->computeActiveErrors();
+ init_error = optimizer->activeChi2();
+ optimizer->optimize(num_iter);
+ final_error = optimizer->activeChi2();
+}
+
+g2oFrameSE3*
+createG2oFrameSE3(Frame* frame, size_t id, bool fixed)
+{
+ g2oFrameSE3* v = new g2oFrameSE3();
+ v->setId(id);
+ v->setFixed(fixed);
+ v->setEstimate(g2o::SE3Quat(frame->T_f_w_.unit_quaternion(), frame->T_f_w_.translation()));
+ return v;
+}
+
+g2oPoint*
+createG2oPoint(Vector3d pos,
+ size_t id,
+ bool fixed)
+{
+ g2oPoint* v = new g2oPoint();
+ v->setId(id);
+#if SCHUR_TRICK
+ v->setMarginalized(true);
+#endif
+ v->setFixed(fixed);
+ v->setEstimate(pos);
+ return v;
+}
+
+g2oEdgeSE3*
+createG2oEdgeSE3( g2oFrameSE3* v_frame,
+ g2oPoint* v_point,
+ const Vector2d& f_up,
+ bool robust_kernel,
+ double huber_width,
+ double weight)
+{
+ g2oEdgeSE3* e = new g2oEdgeSE3();
+ e->setVertex(0, dynamic_cast(v_point));
+ e->setVertex(1, dynamic_cast(v_frame));
+ e->setMeasurement(f_up);
+ e->information() = weight * Eigen::Matrix2d::Identity(2,2);
+ g2o::RobustKernelHuber* rk = new g2o::RobustKernelHuber(); // TODO: memory leak
+ rk->setDelta(huber_width);
+ e->setRobustKernel(rk);
+ e->setParameterId(0, 0); //old: e->setId(v_point->id());
+ return e;
+}
+
+} // namespace ba
+} // namespace svo
diff --git a/CPU_version/svo/src/config.cpp b/CPU_version/svo/src/config.cpp
new file mode 100644
index 0000000..e8475fe
--- /dev/null
+++ b/CPU_version/svo/src/config.cpp
@@ -0,0 +1,102 @@
+// This file is part of SVO - Semi-direct Visual Odometry.
+//
+// Copyright (C) 2014 Christian Forster
+// (Robotics and Perception Group, University of Zurich, Switzerland).
+//
+// SVO is free software: you can redistribute it and/or modify it under the
+// terms of the GNU General Public License as published by the Free Software
+// Foundation, either version 3 of the License, or any later version.
+//
+// SVO is distributed in the hope that it will be useful, but WITHOUT ANY
+// WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
+// FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.
+//
+// You should have received a copy of the GNU General Public License
+// along with this program. If not, see .
+
+#ifdef SVO_USE_ROS
+#include
+#endif
+#include
+
+namespace svo {
+
+Config::Config() :
+#ifdef SVO_USE_ROS
+ trace_name(vk::getParam("svo/trace_name", "svo")),
+ trace_dir(vk::getParam("svo/trace_dir", "/tmp")),
+ n_pyr_levels(vk::getParam("svo/n_pyr_levels", 3)),
+ use_imu(vk::getParam("svo/use_imu", false)),
+ core_n_kfs(vk::getParam("svo/core_n_kfs", 3)),
+ map_scale(vk::getParam("svo/map_scale", 1.0)),
+ grid_size(vk::getParam("svo/grid_size", 30)),
+ init_min_disparity(vk::getParam("svo/init_min_disparity", 50.0)),
+ init_min_tracked(vk::getParam("svo/init_min_tracked", 50)),
+ init_min_inliers(vk::getParam("svo/init_min_inliers", 40)),
+ klt_max_level(vk::getParam("svo/klt_max_level", 4)),
+ klt_min_level(vk::getParam("svo/klt_min_level", 2)),
+ reproj_thresh(vk::getParam("svo/reproj_thresh", 2.0)),
+ poseoptim_thresh(vk::getParam("svo/poseoptim_thresh", 2.0)),
+ poseoptim_num_iter(vk::getParam |