zhangyiss/gm3d-gui
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Compare commits

base: zhangyiss:master
Branches Tags
zhangyiss:master
zhangyiss:yzhang
zhangyiss:ssd_fm
zhangyiss:v1.0
..
compare: zhangyiss:ssd_fm
Branches Tags
zhangyiss:master
zhangyiss:yzhang
zhangyiss:ssd_fm
zhangyiss:v1.0

2 Commits
master ... ssd_fm

Author SHA1 Message Date
Yi Zhang
3c505c0cfe Update read_model_mst.cpp 2019-09-13 13:16:00 +08:00
kiciollo
716939c93a add function for forward mdling using meshtool3d mesh and phy files 2019-09-12 18:55:37 +08:00
67 changed files with 415 additions and 942748 deletions
Show Stats Expand all files Collapse all files

4
.gitignore vendored
View File

@ -1,4 +1,4 @@
# Prerequisites
# Prerequisites
*.d
# Compiled Object files
@ -33,4 +33,4 @@
# Mac folder
.DS_Store
build/
build_mac/

15
CMakeLists.txt
View File

@ -1,12 +1,7 @@
cmake_minimum_required(VERSION 3.15.2)
# 设置工程名称
project(GM3D_GUI VERSION 0.3 LANGUAGES CXX)
cmake_minimum_required(VERSION 3.15.2)
message(STATUS "Platform: " ${CMAKE_HOST_SYSTEM_NAME})
# CMake默认的安装路径 Windows下为C:/Program\ Files/${Project_Name} Linux/Unix下为/usr/local
message(STATUS "Install prefix: " ${CMAKE_INSTALL_PREFIX})
# CMake默认的变异类型为空
message(STATUS "Build type: " ${CMAKE_BUILD_TYPE})
set(CMAKE_CXX_COMPILER /usr/local/bin/g++-9)
# 添加源文件地址
add_subdirectory(src/)
project(GM3D_GUI VERSION 0.2 LANGUAGES CXX)
add_subdirectory(src)

665
LICENSE
View File

@ -1,665 +0,0 @@
gm3d License
--------------
gm3d is distributed under a dual licensing scheme. You can
redistribute it and/or modify it under the terms of the GNU Affero
General Public License as published by the Free Software Foundation,
either version 3 of the License, or (at your option) any later
version. A copy of the GNU Affero General Public License is reproduced
below.
If the terms and conditions of the AGPL v.3. would prevent you from
using gm3d, please consider the option to obtain a commercial
license for a fee. These licenses are offered by the Author. As a rule,
licenses are provided "as-is", unlimited in time for a one time
fee. Please send corresponding requests to:
zhangyiss@icloud.com. Please do not forget to include some
description of your company and the realm of its activities.
=====================================================================
GNU AFFERO GENERAL PUBLIC LICENSE
Version 3, 19 November 2007
Copyright © 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 Affero General Public License is a free, copyleft license for
software and other kinds of works, specifically designed to ensure
cooperation with the community in the case of network server software.
The licenses for most software and other practical works are designed
to take away your freedom to share and change the works. By contrast,
our General Public Licenses are 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.
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.
Developers that use our General Public Licenses protect your rights
with two steps: (1) assert copyright on the software, and (2) offer
you this License which gives you legal permission to copy, distribute
and/or modify the software.
A secondary benefit of defending all users' freedom is that
improvements made in alternate versions of the program, if they
receive widespread use, become available for other developers to
incorporate. Many developers of free software are heartened and
encouraged by the resulting cooperation. However, in the case of
software used on network servers, this result may fail to come
about. The GNU General Public License permits making a modified
version and letting the public access it on a server without ever
releasing its source code to the public.
The GNU Affero General Public License is designed specifically to
ensure that, in such cases, the modified source code becomes available
to the community. It requires the operator of a network server to
provide the source code of the modified version running there to the
users of that server. Therefore, public use of a modified version, on
a publicly accessible server, gives the public access to the source
code of the modified version.
An older license, called the Affero General Public License and
published by Affero, was designed to accomplish similar goals. This is
a different license, not a version of the Affero GPL, but Affero has
released a new version of the Affero GPL which permits relicensing
under this license.
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 Affero 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. Remote Network Interaction; Use with the GNU General Public
License.
Notwithstanding any other provision of this License, if you modify the
Program, your modified version must prominently offer all users
interacting with it remotely through a computer network (if your
version supports such interaction) an opportunity to receive the
Corresponding Source of your version by providing access to the
Corresponding Source from a network server at no charge, through some
standard or customary means of facilitating copying of software. This
Corresponding Source shall include the Corresponding Source for any
work covered by version 3 of the GNU General Public License that is
incorporated pursuant to the following paragraph.
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 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 work with which it is combined will remain governed by version
3 of the GNU General Public License.
14. Revised Versions of this License.
The Free Software Foundation may publish revised and/or new versions
of the GNU Affero 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 Affero 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 Affero 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 Affero 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 Affero 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
Affero General Public License for more details.
You should have received a copy of the GNU Affero 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 your software can interact with users remotely through a computer
network, you should also make sure that it provides a way for users to
get its source. For example, if your program is a web application, its
interface could display a "Source" link that leads users to an archive
of the code. There are many ways you could offer source, and different
solutions will be better for different programs; see section 13 for
the specific requirements.
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 AGPL, see <http://www.gnu.org/licenses/>.

539
archive/gm3d_gui.cpp
View File

@ -1,539 +0,0 @@
#include "gm3d_gui.h"
#include "gm3d.h"
Fl_Double_Window *main_window=(Fl_Double_Window *)0;
Fl_Tabs *main_tabs=(Fl_Tabs *)0;
Fl_Group *model_tab=(Fl_Group *)0;
Fl_Input *mesh_para_input=(Fl_Input *)0;
Fl_Button *mesh_file_btn=(Fl_Button *)0;
Fl_Button *mod_para_file_btn=(Fl_Button *)0;
Fl_Input *mod_ele_input_build=(Fl_Input *)0;
Fl_Button *build_mod_btn=(Fl_Button *)0;
Fl_Output *mesh_para_output=(Fl_Output *)0;
Fl_Button *add_mod_btn=(Fl_Button *)0;
Fl_Button *del_mod_btn=(Fl_Button *)0;
Fl_Check_Button *rm_emp_bok_check=(Fl_Check_Button *)0;
Fl_Browser *mod_para_brw=(Fl_Browser *)0;
Fl_Button *mod_file_out_btn=(Fl_Button *)0;
Fl_Output *mod_out_file_output=(Fl_Output *)0;
Fl_Group *forward_tab=(Fl_Group *)0;
Fl_Group *grav_group=(Fl_Group *)0;
Fl_Check_Button *Vz_check=(Fl_Check_Button *)0;
Fl_Check_Button *Vzx_check=(Fl_Check_Button *)0;
Fl_Check_Button *Vzy_check=(Fl_Check_Button *)0;
Fl_Check_Button *Vzz_check=(Fl_Check_Button *)0;
Fl_Button *mod_file_btn=(Fl_Button *)0;
Fl_Button *obs_file_btn=(Fl_Button *)0;
Fl_Input *mod_file_input=(Fl_Input *)0;
Fl_Output *mod_file_output=(Fl_Output *)0;
Fl_Input *obs_file_input=(Fl_Input *)0;
Fl_Output *obs_file_output=(Fl_Output *)0;
Fl_Input *mod_ele_input=(Fl_Input *)0;
Fl_Group *mag_group=(Fl_Group *)0;
Fl_Check_Button *DeltaT_check=(Fl_Check_Button *)0;
Fl_Check_Button *DeltaTx_check=(Fl_Check_Button *)0;
Fl_Check_Button *DeltaTy_check=(Fl_Check_Button *)0;
Fl_Check_Button *DeltaTz_check=(Fl_Check_Button *)0;
Fl_Check_Button *Hax_check=(Fl_Check_Button *)0;
Fl_Check_Button *Hay_check=(Fl_Check_Button *)0;
Fl_Check_Button *Za_check=(Fl_Check_Button *)0;
Fl_Input *mag_para_input=(Fl_Input *)0;
Fl_Check_Button *mag_data_check=(Fl_Check_Button *)0;
Fl_Button *cal_btn=(Fl_Button *)0;
Fl_Input *noise_para_input=(Fl_Input *)0;
Fl_Check_Button *noise_check=(Fl_Check_Button *)0;
Fl_Input *res_file_input=(Fl_Input *)0;
Fl_Button *res_file_btn=(Fl_Button *)0;
Fl_Output *res_file_output=(Fl_Output *)0;
char mesh_filename[1024];
char out_msh_filename[1024];
void cb_mesh_para_input(Fl_Input*, void*){
strcpy(mesh_filename,mesh_para_input->value());
mesh_para_output->value(mesh_para_input->value());
return;
}
void cb_mesh_file_btn(Fl_Widget*, void*){
// Create the file chooser, and show it
// directory, filter, chooser type, title
Fl_File_Chooser chooser(".", "*.txt, *.dat, *.xyz", Fl_File_Chooser::SINGLE, "Choose a mesh file.");
chooser.show();
// Block until user picks something.
while(chooser.shown()) { Fl::wait(); }
//从文件的绝对路径中分离文件名并从输出框输出
if (chooser.value() != NULL)
{
strcpy(mesh_filename,chooser.value());
std::string abs_filename = chooser.value();
int pos = abs_filename.find_last_of("/");
std::stringstream temp_ss;
temp_ss << abs_filename.substr(pos+1);
const char* disp_filename = temp_ss.str().c_str();
mesh_para_output->value(disp_filename);
}
else
{
mesh_para_output->value("Unset");
}
return;
}
void cb_mod_para_file_btn(Fl_Widget*, void*){
// Create the file chooser, and show it
Fl_File_Chooser chooser(".", "*.txt, *.dat, *.xyz", Fl_File_Chooser::SINGLE, "Choose a model list file.");
chooser.show();
// Block until user picks something.
while(chooser.shown()) { Fl::wait(); }
//下面我们读入文件并添加模型参数到列表显示
std::ifstream modpara_in;
modpara_in.open(chooser.value());
if (!modpara_in)
{
fl_message("Error! Fail to open file: %s",chooser.value());
return;
}
else{
std::string temp_str;
std::stringstream temp_ss;
while(getline(modpara_in,temp_str)){
if (*(temp_str.begin()) == '#') continue;
else{
temp_ss.clear(); temp_ss.str("");
temp_ss << temp_str;
mod_para_brw->add(temp_ss.str().c_str());
}
}
modpara_in.close();
}
return;
}
void cb_del_mod_btn(Fl_Widget*, void*){
int index = mod_para_brw->value();
mod_para_brw->remove(index);
del_mod_btn->deactivate();
}
void cb_mod_para_brw(Fl_Browser*, void*){
for ( int t=1; t<=mod_para_brw->size(); t++ ) {
if ( mod_para_brw->selected(t) ) {
//printf("%d) '%s'\n", t, mod_para_brw->text(t));
del_mod_btn->activate();
}
}
return;
}
bool remove_null = false;
void cb_rm_emp_bok_check(Fl_Check_Button*, void*){
if (rm_emp_bok_check->value() == 1) remove_null = true;
else if (rm_emp_bok_check->value() == 0) remove_null = false;
return;
}
void cb_mod_file_out_btn(Fl_Widget*, void*){
// Create the file chooser, and show it
// directory, filter, chooser type, title
Fl_File_Chooser chooser(".", "*.msh", Fl_File_Chooser::CREATE, "Create a new model file.");
chooser.show();
// Block until user picks something.
while(chooser.shown()) { Fl::wait(); }
//从文件的绝对路径中分离文件名并从输出框输出
if (chooser.value() != NULL)
{
strcpy(out_msh_filename,chooser.value());
std::string abs_filename = chooser.value();
int pos = abs_filename.find_last_of("/");
std::stringstream temp_ss;
temp_ss << abs_filename.substr(pos+1);
const char* disp_filename = temp_ss.str().c_str();
mod_out_file_output->value(disp_filename);
}
else
{
mod_out_file_output->value("Unset");
}
return;
}
void cb_build_mod_btn(Fl_Widget*, void*){
//首先构建一个GM3D实例
GM3D gm3d_instance;
//设置初始化变量
char dimension[1024] = "10/20/990/10/20/990/10/20/490";
char out_mshname[1024] = "Untitled.msh";
char elename[1024] = "Untitled";
//将mesh_para_output拷贝到dimension
strcpy(dimension, mesh_filename);
//std::cout << dimension << std::endl;
//将mod_para_brw中的模型参数列表拷贝到一个临时参数列表
std::stringstream temp_ss;
modelist temp_list;
modelistArray brw_model_list;
for ( int t=1; t<=mod_para_brw->size(); t++ ) {
temp_ss.clear(); temp_ss.str("");
temp_ss << mod_para_brw->text(t);
if(temp_ss >> temp_list.mod_type >> temp_list.val_type >> temp_list.mod_value >> temp_list.mod_para){
brw_model_list.push_back(temp_list);
}
else{
fl_message("Wrong model parameter: %s. Skipped !",mod_para_brw->text(t));
}
}
//将brw_model_list拷贝到GM3D中
gm3d_instance.get_model_list(brw_model_list);
//将mod_out_file_input拷贝到out_mshname
strcpy(out_mshname, out_msh_filename);
if (!strcmp(out_mshname,"")){
fl_message("Output file's name can't be empty !");
return;
}
//将mod_ele_input_build拷贝到elename
strcpy(elename, mod_ele_input_build->value());
if (!strcmp(elename,"")){
fl_message("Element data's name can't be empty !");
return;
}
//构建模型网络
if(gm3d_instance.BuildRegularGrid(dimension)){
fl_message("Mesh Parameters Load Error !");
return;
}
gm3d_instance.AddModels_GUI();
gm3d_instance.RegisteredOuput(remove_null);
if (gm3d_instance.OutMshFile(out_mshname,elename)) return;
//这里需要一个输出信息窗口
fl_message("Model Construction Completed !");
return;
}
/***********************************以下是正演计算GUI函数******************************/
char in_msh_filename[1024];
char in_obs_filename[1024];
char out_res_filename[1024];
void cb_mod_file_input(Fl_Input*, void*){
strcpy(in_msh_filename,mod_file_input->value());
mod_file_output->value(mod_file_input->value());
return;
}
void cb_mod_file_btn(Fl_Widget*, void*){
// Create the file chooser, and show it
// directory, filter, chooser type, title
Fl_File_Chooser chooser(".", "*.msh", Fl_File_Chooser::SINGLE, "Create a model file.");
chooser.show();
// Block until user picks something.
while(chooser.shown()) { Fl::wait(); }
//从文件的绝对路径中分离文件名并从输出框输出
if (chooser.value() != NULL)
{
strcpy(in_msh_filename,chooser.value());
std::string abs_filename = chooser.value();
int pos = abs_filename.find_last_of("/");
std::stringstream temp_ss;
temp_ss << abs_filename.substr(pos+1);
const char* disp_filename = temp_ss.str().c_str();
mod_file_output->value(disp_filename);
}
else
{
mod_file_output->value("Unset");
}
return;
}
void cb_obs_file_input(Fl_Input*, void*){
strcpy(in_obs_filename,obs_file_input->value());
obs_file_output->value(obs_file_input->value());
return;
}
void cb_obs_file_btn(Fl_Widget*, void*){
// Create the file chooser, and show it
// directory, filter, chooser type, title
Fl_File_Chooser chooser(".", "*.txt, *.dat, *.xyz", Fl_File_Chooser::SINGLE, "Create observe file.");
chooser.show();
// Block until user picks something.
while(chooser.shown()) { Fl::wait(); }
//从文件的绝对路径中分离文件名并从输出框输出
if (chooser.value() != NULL)
{
strcpy(in_obs_filename,chooser.value());
std::string abs_filename = chooser.value();
int pos = abs_filename.find_last_of("/");
std::stringstream temp_ss;
temp_ss << abs_filename.substr(pos+1);
const char* disp_filename = temp_ss.str().c_str();
obs_file_output->value(disp_filename);
}
else
{
obs_file_output->value("Unset");
}
return;
}
void cb_res_file_input(Fl_Input*, void*){
strcpy(out_res_filename,res_file_input->value());
res_file_output->value(res_file_input->value());
return;
}
void cb_res_file_btn(Fl_Widget*, void*){
// Create the file chooser, and show it
// directory, filter, chooser type, title
Fl_File_Chooser chooser(".", "*.txt, *.dat", Fl_File_Chooser::CREATE, "Create a new observe file.");
chooser.show();
// Block until user picks something.
while(chooser.shown()) { Fl::wait(); }
//从文件的绝对路径中分离文件名并从输出框输出
if (chooser.value() != NULL)
{
strcpy(out_res_filename,chooser.value());
std::string abs_filename = chooser.value();
int pos = abs_filename.find_last_of("/");
std::stringstream temp_ss;
temp_ss << abs_filename.substr(pos+1);
const char* disp_filename = temp_ss.str().c_str();
res_file_output->value(disp_filename);
}
else
{
res_file_output->value("Unset");
}
return;
}
void cb_noise_check(Fl_Check_Button*, void*){
if (noise_check->value() == 1){
noise_para_input->activate();
}
else if (noise_check->value() == 0){
noise_para_input->deactivate();
}
return;
}
void cb_mag_data_check(Fl_Check_Button*, void*){
if (mag_data_check->value() == 1){
mag_para_input->activate();
DeltaT_check->activate();
DeltaTx_check->activate();
DeltaTy_check->activate();
DeltaTz_check->activate();
Hax_check->activate();
Hay_check->activate();
Za_check->activate();
}
else if (mag_data_check->value() == 0){
mag_para_input->deactivate();
DeltaT_check->deactivate();
DeltaTx_check->deactivate();
DeltaTy_check->deactivate();
DeltaTz_check->deactivate();
Hax_check->deactivate();
Hay_check->deactivate();
Za_check->deactivate();
}
return;
}
void cb_cal_btn(Fl_Button*, void*){
//首先构建一个GM3D实例
GM3D gm3d_instance;
//设置初始化变量
char in_mshname[1024] = "Untitled.msh";
char in_obspara[1024] = "Untitled.txt";
char res_outfile[1024] = "Untitled";
char res_outfile_full[1024];
char ele_name[1024] = "Untitled";
char noise_para[1024] = "0.0/0.0";
char mag_para[1024] = "90.0/0.0/90.0/0.0";
strcpy(in_mshname,in_msh_filename);
strcpy(in_obspara,in_obs_filename);
strcpy(res_outfile,out_res_filename);
strcpy(ele_name,mod_ele_input->value());
strcpy(noise_para,noise_para_input->value());
strcpy(mag_para,mag_para_input->value());
if (gm3d_instance.ReadModel(in_mshname,ele_name))
{
fl_message("Model Load Error !");
return;
}
if (gm3d_instance.InitObs(in_obspara))
{
fl_message("Observe File Load Error !");
return;
}
if (Vz_check->value() == 1)
{
gm3d_instance.ForwardVz(noise_para);
strcpy(res_outfile_full,res_outfile);
strcat(res_outfile_full,"_Vz.txt");
if (gm3d_instance.OutObs(res_outfile_full))
{
fl_message("Observe File Output Error !");
return;
}
}
if (Vzx_check->value() == 1)
{
gm3d_instance.ForwardVzx(noise_para);
strcpy(res_outfile_full,res_outfile);
strcat(res_outfile_full,"_Vzx.txt");
if (gm3d_instance.OutObs(res_outfile_full))
{
fl_message("Observe File Output Error !");
return;
}
}
if (Vzy_check->value() == 1)
{
gm3d_instance.ForwardVzy(noise_para);
strcpy(res_outfile_full,res_outfile);
strcat(res_outfile_full,"_Vzy.txt");
if (gm3d_instance.OutObs(res_outfile_full))
{
fl_message("Observe File Output Error !");
return;
}
}
if (Vzz_check->value() == 1)
{
gm3d_instance.ForwardVzz(noise_para);
strcpy(res_outfile_full,res_outfile);
strcat(res_outfile_full,"_Vzz.txt");
if (gm3d_instance.OutObs(res_outfile_full))
{
fl_message("Observe File Output Error !");
return;
}
}
if (mag_data_check->value() == 1)
{
if (DeltaT_check->value() == 1)
{
gm3d_instance.ForwardDeltaT(noise_para,mag_para);
strcpy(res_outfile_full,res_outfile);
strcat(res_outfile_full,"_DeltaT.txt");
if (gm3d_instance.OutObs(res_outfile_full))
{
fl_message("Observe File Output Error !");
return;
}
}
if (DeltaTx_check->value() == 1)
{
gm3d_instance.ForwardDeltaTx(noise_para,mag_para);
strcpy(res_outfile_full,res_outfile);
strcat(res_outfile_full,"_DeltaTx.txt");
if (gm3d_instance.OutObs(res_outfile_full))
{
fl_message("Observe File Output Error !");
return;
}
}
if (DeltaTy_check->value() == 1)
{
gm3d_instance.ForwardDeltaTy(noise_para,mag_para);
strcpy(res_outfile_full,res_outfile);
strcat(res_outfile_full,"_DeltaTy.txt");
if (gm3d_instance.OutObs(res_outfile_full))
{
fl_message("Observe File Output Error !");
return;
}
}
if (DeltaTz_check->value() == 1)
{
gm3d_instance.ForwardDeltaTz(noise_para,mag_para);
strcpy(res_outfile_full,res_outfile);
strcat(res_outfile_full,"_DeltaTz.txt");
if (gm3d_instance.OutObs(res_outfile_full))
{
fl_message("Observe File Output Error !");
return;
}
}
if (Hax_check->value() == 1)
{
gm3d_instance.ForwardHax(noise_para,mag_para);
strcpy(res_outfile_full,res_outfile);
strcat(res_outfile_full,"_Hax.txt");
if (gm3d_instance.OutObs(res_outfile_full))
{
fl_message("Observe File Output Error !");
return;
}
}
if (Hay_check->value() == 1)
{
gm3d_instance.ForwardHay(noise_para,mag_para);
strcpy(res_outfile_full,res_outfile);
strcat(res_outfile_full,"_Hay.txt");
if (gm3d_instance.OutObs(res_outfile_full))
{
fl_message("Observe File Output Error !");
return;
}
}
if (Za_check->value() == 1)
{
gm3d_instance.ForwardZa(noise_para,mag_para);
strcpy(res_outfile_full,res_outfile);
strcat(res_outfile_full,"_Za.txt");
if (gm3d_instance.OutObs(res_outfile_full))
{
fl_message("Observe File Output Error !");
return;
}
}
}
fl_message("Forward calculation completed !");
return;
}

115
archive/gm3d_gui.h
View File

@ -1,115 +0,0 @@
// generated by Fast Light User Interface Designer (fluid) version 1.0305
#ifndef _GM3D_GUI_H
#define _GM3D_GUI_H
#include <FL/Fl.H>
#include <FL/Fl_Double_Window.H>
#include <FL/Fl_Tabs.H>
#include <FL/Fl_Group.H>
#include <FL/Fl_Input.H>
#include <FL/Fl_Button.H>
#include <FL/Fl_Output.H>
#include <FL/Fl_Check_Button.H>
#include <FL/Fl_Radio_Round_Button.H>
#include <FL/Fl_Browser.H>
#include <FL/Fl_File_Chooser.H>
#include <iostream>
#include <string.h>
#include <sstream>
#include <fstream>
extern void cb_mesh_para_input(Fl_Input*, void*);
extern void cb_mesh_file_btn(Fl_Widget*, void*);
extern void cb_add_mod_btn(Fl_Widget*, void*);
extern void cb_mod_para_file_btn(Fl_Widget*, void*);
extern void cb_del_mod_btn(Fl_Widget*, void*);
extern void cb_mod_para_brw(Fl_Browser*, void*);
extern void cb_rm_emp_bok_check(Fl_Check_Button*, void*);
extern void cb_mod_file_out_btn(Fl_Widget*, void*);
extern void cb_build_mod_btn(Fl_Widget*, void*);
extern void cb_mod_file_input(Fl_Input*, void*);
extern void cb_mod_file_btn(Fl_Widget*, void*);
extern void cb_obs_file_input(Fl_Input*, void*);
extern void cb_obs_file_btn(Fl_Widget*, void*);
extern void cb_res_file_input(Fl_Input*, void*);
extern void cb_res_file_btn(Fl_Widget*, void*);
extern void cb_noise_check(Fl_Check_Button*, void*);
extern void cb_mag_data_check(Fl_Check_Button*, void*);
extern void cb_cal_btn(Fl_Button*, void*);
extern void cb_reg_bok_rbtn(Fl_Button*, void*);
extern void cb_til_bok_rbtn(Fl_Button*, void*);
extern void cb_sph_rbtn(Fl_Button*, void*);
extern void cb_int_face_rbtn(Fl_Button*, void*);
extern void cb_app_val_rbtn(Fl_Button*, void*);
extern void cb_rep_val_rbtn(Fl_Button*, void*);
extern void cb_era_val_rbtn(Fl_Button*, void*);
extern void cb_top_val_rbtn(Fl_Button*, void*);
extern void cb_bot_val_rbtn(Fl_Button*, void*);
extern void cb_sig_add_btn(Fl_Widget*, void*);
extern void cb_can_add_btn(Fl_Widget*, void*);
extern Fl_Double_Window *main_window;
extern Fl_Tabs *main_tabs;
extern Fl_Group *model_tab;
extern Fl_Input *mesh_para_input;
extern Fl_Button *mesh_file_btn;
extern Fl_Button *mod_para_file_btn;
extern Fl_Input *mod_ele_input_build;
extern Fl_Button *build_mod_btn;
extern Fl_Output *mesh_para_output;
extern Fl_Button *add_mod_btn;
extern Fl_Button *del_mod_btn;
extern Fl_Check_Button *rm_emp_bok_check;
extern Fl_Browser *mod_para_brw;
extern Fl_Button *mod_file_out_btn;
extern Fl_Output *mod_out_file_output;
extern Fl_Group *forward_tab;
extern Fl_Group *grav_group;
extern Fl_Check_Button *Vz_check;
extern Fl_Check_Button *Vzx_check;
extern Fl_Check_Button *Vzy_check;
extern Fl_Check_Button *Vzz_check;
extern Fl_Button *mod_file_btn;
extern Fl_Button *obs_file_btn;
extern Fl_Input *mod_file_input;
extern Fl_Output *mod_file_output;
extern Fl_Input *obs_file_input;
extern Fl_Output *obs_file_output;
extern Fl_Input *mod_ele_input;
extern Fl_Group *mag_group;
extern Fl_Check_Button *DeltaT_check;
extern Fl_Check_Button *DeltaTx_check;
extern Fl_Check_Button *DeltaTy_check;
extern Fl_Check_Button *DeltaTz_check;
extern Fl_Check_Button *Hax_check;
extern Fl_Check_Button *Hay_check;
extern Fl_Check_Button *Za_check;
extern Fl_Input *mag_para_input;
extern Fl_Check_Button *mag_data_check;
extern Fl_Button *cal_btn;
extern Fl_Input *noise_para_input;
extern Fl_Check_Button *noise_check;
extern Fl_Input *res_file_input;
extern Fl_Button *res_file_btn;
extern Fl_Output *res_file_output;
extern Fl_Double_Window *add_mod_win;
extern Fl_Group *mod_type_group;
extern Fl_Radio_Round_Button *reg_bok_rbtn;
extern Fl_Radio_Round_Button *til_bok_rbtn;
extern Fl_Radio_Round_Button *sph_rbtn;
extern Fl_Radio_Round_Button *int_face_rbtn;
extern Fl_Group *val_type_group;
extern Fl_Radio_Round_Button *app_val_rbtn;
extern Fl_Radio_Round_Button *rep_val_rbtn;
extern Fl_Radio_Round_Button *era_val_rbtn;
extern Fl_Group *agn_part_group;
extern Fl_Radio_Round_Button *top_val_btn;
extern Fl_Radio_Round_Button *bot_val_btn;
extern Fl_Input *sig_mod_para_input;
extern Fl_Return_Button *can_add_btn;
extern Fl_Button *sig_add_btn;
extern Fl_Input *mod_val_input;
#endif

BIN
assert/logo/Logo_ZY.icns
View File

Binary file not shown.

38
fluid_project/add_model.cxx
View File

@ -2,7 +2,7 @@
#include "add_model.h"
Fl_Double_Window *edit_mod_win=(Fl_Double_Window *)0;
Fl_Double_Window *add_mod_win=(Fl_Double_Window *)0;
Fl_Group *mod_type_group=(Fl_Group *)0;
@ -30,18 +30,14 @@ Fl_Round_Button *bot_val_btn=(Fl_Round_Button *)0;
Fl_Input *sig_mod_para_input=(Fl_Input *)0;
Fl_Input *mod_val_input=(Fl_Input *)0;
Fl_Return_Button *can_add_btn=(Fl_Return_Button *)0;
Fl_Button *sig_add_btn=(Fl_Button *)0;
Fl_Button *sig_replace_btn=(Fl_Button *)0;
Fl_Button *sig_delete_btn=(Fl_Button *)0;
Fl_Return_Button *can_add_btn=(Fl_Return_Button *)0;
Fl_Input *mod_val_input=(Fl_Input *)0;
void cb_add_mod_btn(Fl_Button*, void*) {
{ edit_mod_win = new Fl_Double_Window(314, 320, "Model editor (gm3d)");
{ add_mod_win = new Fl_Double_Window(315, 300, "Add model (gm3d)");
{ mod_type_group = new Fl_Group(20, 25, 250, 53, "Model Type :");
mod_type_group->align(Fl_Align(FL_ALIGN_TOP_LEFT));
{ reg_bok_rbtn = new Fl_Round_Button(20, 25, 110, 28, "Regular Block");
@ -50,7 +46,7 @@ void cb_add_mod_btn(Fl_Button*, void*) {
{ til_bok_rbtn = new Fl_Round_Button(160, 25, 100, 28, "Tilted Block");
til_bok_rbtn->down_box(FL_ROUND_DOWN_BOX);
} // Fl_Round_Button* til_bok_rbtn
{ sph_rbtn = new Fl_Round_Button(20, 50, 70, 28, "Ellipsoid");
{ sph_rbtn = new Fl_Round_Button(20, 50, 70, 28, "Sphere");
sph_rbtn->down_box(FL_ROUND_DOWN_BOX);
} // Fl_Round_Button* sph_rbtn
{ int_face_rbtn = new Fl_Round_Button(160, 50, 80, 28, "Interface");
@ -71,13 +67,13 @@ void cb_add_mod_btn(Fl_Button*, void*) {
} // Fl_Round_Button* era_val_rbtn
val_type_group->end();
} // Fl_Group* val_type_group
{ agn_part_group = new Fl_Group(20, 150, 160, 28, "Assign Part :");
{ agn_part_group = new Fl_Group(20, 150, 160, 28, "Assgin Part :");
agn_part_group->align(Fl_Align(FL_ALIGN_TOP_LEFT));
{ top_val_btn = new Fl_Round_Button(20, 150, 50, 28, "Top");
top_val_btn->down_box(FL_ROUND_DOWN_BOX);
top_val_btn->deactivate();
} // Fl_Round_Button* top_val_btn
{ bot_val_btn = new Fl_Round_Button(90, 150, 70, 28, "Bottom");
{ bot_val_btn = new Fl_Round_Button(110, 150, 70, 28, "Bottom");
bot_val_btn->down_box(FL_ROUND_DOWN_BOX);
bot_val_btn->deactivate();
} // Fl_Round_Button* bot_val_btn
@ -86,19 +82,13 @@ void cb_add_mod_btn(Fl_Button*, void*) {
{ sig_mod_para_input = new Fl_Input(20, 200, 275, 28, "Model Parameter :");
sig_mod_para_input->align(Fl_Align(FL_ALIGN_TOP_LEFT));
} // Fl_Input* sig_mod_para_input
{ mod_val_input = new Fl_Input(180, 150, 115, 28, "Model Value :");
{ can_add_btn = new Fl_Return_Button(210, 255, 85, 28, "Cancel");
} // Fl_Return_Button* can_add_btn
{ sig_add_btn = new Fl_Button(140, 255, 60, 28, "Add");
} // Fl_Button* sig_add_btn
{ mod_val_input = new Fl_Input(20, 255, 110, 28, "Model Value :");
mod_val_input->align(Fl_Align(FL_ALIGN_TOP_LEFT));
} // Fl_Input* mod_val_input
{ sig_add_btn = new Fl_Button(19, 240, 80, 28, "Add");
} // Fl_Button* sig_add_btn
{ sig_replace_btn = new Fl_Button(118, 240, 80, 28, "Replace");
sig_replace_btn->deactivate();
} // Fl_Button* sig_replace_btn
{ sig_delete_btn = new Fl_Button(215, 240, 80, 28, "Remove");
sig_delete_btn->deactivate();
} // Fl_Button* sig_delete_btn
{ can_add_btn = new Fl_Return_Button(215, 280, 80, 28, "Cancel");
} // Fl_Return_Button* can_add_btn
edit_mod_win->end();
} // Fl_Double_Window* edit_mod_win
add_mod_win->end();
} // Fl_Double_Window* add_mod_win
}

44
fluid_project/add_model.fl
View File

@ -4,12 +4,12 @@ header_name {.h}
code_name {.cxx}
Function {cb_add_mod_btn(Fl_Button*, void*)} {open return_type void
} {
Fl_Window edit_mod_win {
label {gm3d editor} open selected
xywh {645 271 314 320} type Double visible
Fl_Window add_mod_win {
label {Add model (gm3d)} open
xywh {735 331 315 300} type Double visible
} {
Fl_Group mod_type_group {
label {Entity Type :} open
label {Model Type :} open
xywh {20 25 250 53} align 5
} {
Fl_Round_Button reg_bok_rbtn {
@ -30,7 +30,7 @@ Function {cb_add_mod_btn(Fl_Button*, void*)} {open return_type void
}
}
Fl_Group val_type_group {
label {Value Type:} open
label {Value Type :} open
xywh {20 100 250 30} align 5
} {
Fl_Round_Button app_val_rbtn {
@ -47,41 +47,33 @@ Function {cb_add_mod_btn(Fl_Button*, void*)} {open return_type void
}
}
Fl_Group agn_part_group {
label {Solid Side:} open
label {Assgin Part :} open
xywh {20 150 160 28} align 5
} {
Fl_Round_Button top_val_btn {
label Above
xywh {20 150 70 28} down_box ROUND_DOWN_BOX deactivate
label Top
xywh {20 150 50 28} down_box ROUND_DOWN_BOX deactivate
}
Fl_Round_Button bot_val_btn {
label Below
xywh {95 150 70 28} down_box ROUND_DOWN_BOX deactivate
label Bottom
xywh {110 150 70 28} down_box ROUND_DOWN_BOX deactivate
}
}
Fl_Input sig_mod_para_input {
label {Geometric Parameters:}
label {Model Parameter :}
xywh {20 200 275 28} align 5
}
Fl_Input mod_val_input {
label {Physical Value:}
xywh {180 150 115 28} align 5
Fl_Return_Button can_add_btn {
label Cancel
xywh {210 255 85 28}
}
Fl_Button sig_add_btn {
label Add
xywh {19 240 80 28}
xywh {140 255 60 28}
}
Fl_Button sig_replace_btn {
label Replace
xywh {118 240 80 28} deactivate
}
Fl_Button sig_delete_btn {
label Remove
xywh {215 240 80 28} deactivate
}
Fl_Return_Button can_add_btn {
label Quit
xywh {215 280 80 28}
Fl_Input mod_val_input {
label {Model Value :}
xywh {20 255 110 28} align 5
}
}
}

10
fluid_project/add_model.h
View File

@ -4,7 +4,7 @@
#define add_model_h
#include <FL/Fl.H>
#include <FL/Fl_Double_Window.H>
extern Fl_Double_Window *edit_mod_win;
extern Fl_Double_Window *add_mod_win;
#include <FL/Fl_Group.H>
extern Fl_Group *mod_type_group;
#include <FL/Fl_Round_Button.H>
@ -21,12 +21,10 @@ extern Fl_Round_Button *top_val_btn;
extern Fl_Round_Button *bot_val_btn;
#include <FL/Fl_Input.H>
extern Fl_Input *sig_mod_para_input;
extern Fl_Input *mod_val_input;
#include <FL/Fl_Button.H>
extern Fl_Button *sig_add_btn;
extern Fl_Button *sig_replace_btn;
extern Fl_Button *sig_delete_btn;
#include <FL/Fl_Return_Button.H>
extern Fl_Return_Button *can_add_btn;
#include <FL/Fl_Button.H>
extern Fl_Button *sig_add_btn;
extern Fl_Input *mod_val_input;
void cb_add_mod_btn(Fl_Button*, void*);
#endif

185
fluid_project/add_model2.fl
View File

@ -1,185 +0,0 @@
# data file for the Fltk User Interface Designer (fluid)
version 1.0305
header_name {.h}
code_name {.cxx}
Function {cb_add_mod_btn(Fl_Button*, void*)} {open return_type void
} {
Fl_Window edit_mod_win {
label {Physical Entity (gm3d)} open
xywh {551 245 405 493} type Double visible
} {
Fl_Group val_type_group {
label {Value Type:}
xywh {20 340 250 30} align 5
} {
Fl_Round_Button app_val_rbtn {
label Append
xywh {20 340 80 28} down_box ROUND_DOWN_BOX
}
Fl_Round_Button rep_val_rbtn {
label Replace
xywh {100 340 80 28} down_box ROUND_DOWN_BOX
}
Fl_Round_Button era_val_rbtn {
label Erase
xywh {210 340 60 28} down_box ROUND_DOWN_BOX
}
}
Fl_Group agn_part_group {
label {Solid Side:}
xywh {20 400 160 28} align 5
} {
Fl_Round_Button top_val_btn {
label Above
xywh {20 400 70 28} down_box ROUND_DOWN_BOX deactivate
}
Fl_Round_Button bot_val_btn {
label Below
xywh {100 400 70 28} down_box ROUND_DOWN_BOX deactivate
}
}
Fl_Input mod_val_input {
label {Physical Value:}
xywh {210 400 175 28} align 5
}
Fl_Button sig_add_btn {
label Add
xywh {20 447 80 28}
}
Fl_Button sig_replace_btn {
label Replace
xywh {115 447 80 28} deactivate
}
Fl_Button sig_delete_btn {
label Remove
xywh {210 447 80 28} deactivate
}
Fl_Return_Button can_add_btn {
label Quit
xywh {305 447 80 28}
}
Fl_Tabs {} {open
xywh {20 15 365 290} box UP_BOX align 5
} {
Fl_Group block_geo_tab {
label Block
xywh {20 15 365 260} hide
} {
Fl_Input {} {
label {xmin:}
xywh {75 35 110 29}
}
Fl_Input {} {
label {xmax:}
xywh {250 35 110 29}
}
Fl_Input {} {
label {xmin:}
xywh {75 80 110 29}
}
Fl_Input {} {
label {xmax:}
xywh {250 80 110 29}
}
Fl_Input {} {
label {xmin:}
xywh {75 125 110 29}
}
Fl_Input {} {
label {xmax:}
xywh {250 125 110 29}
}
}
Fl_Group tilt_block_geo_tab {
label {Tilted Block}
xywh {20 15 365 260} hide
} {
Fl_Input {} {
label {xmin:}
xywh {75 35 110 29}
}
Fl_Input {} {
label {xmin:}
xywh {250 35 110 29}
}
Fl_Input {} {
label {xmin:}
xywh {75 80 110 29}
}
Fl_Input {} {
label {xmin:}
xywh {250 80 110 29}
}
Fl_Input {} {
label {xmin:}
xywh {75 125 110 29}
}
Fl_Input {} {
label {xmin:}
xywh {250 125 110 29}
}
Fl_Input {} {
label {xmin:}
xywh {75 170 110 29}
}
Fl_Input {} {
label {xmin:}
xywh {250 170 110 29}
}
Fl_Input {} {
label {xmin:}
xywh {75 215 110 29}
}
Fl_Input {} {
label {xmin:}
xywh {250 215 110 29}
}
}
Fl_Group sphere_geo_tab {
label Sphere
xywh {20 15 365 260} hide
} {
Fl_Input {} {
label {xmin:}
xywh {75 35 110 29}
}
Fl_Input {} {
label {xmax:}
xywh {250 35 110 29}
}
Fl_Input {} {
label {xmin:}
xywh {75 80 110 29}
}
Fl_Input {} {
label {xmax:}
xywh {250 80 110 29}
}
Fl_Input {} {
label {xmin:}
xywh {75 125 110 29}
}
Fl_Input {} {
label {xmax:}
xywh {250 125 110 29}
}
}
Fl_Group interface_geo_tab {
label Interface
xywh {20 15 365 260}
} {
Fl_Button {} {
label {Import From}
xywh {40 35 100 28}
}
Fl_Browser {} {
xywh {40 75 325 185}
}
Fl_Output {} {
label {File:}
xywh {185 36 180 28} box FLAT_BOX
}
}
}
}
}

223
fluid_project/gm3d_gui2.cxx
View File

@ -1,223 +0,0 @@
// generated by Fast Light User Interface Designer (fluid) version 1.0305
#include "gm3d_gui2.h"
Fl_Double_Window *main_window=(Fl_Double_Window *)0;
Fl_Tabs *main_tabs=(Fl_Tabs *)0;
Fl_Group *model_tab=(Fl_Group *)0;
Fl_Input *mesh_para_input=(Fl_Input *)0;
Fl_Button *mesh_file_btn=(Fl_Button *)0;
Fl_Button *mod_para_file_btn=(Fl_Button *)0;
Fl_Input *mod_ele_input_build=(Fl_Input *)0;
Fl_Button *build_mod_btn=(Fl_Button *)0;
Fl_Output *mesh_para_output=(Fl_Output *)0;
Fl_Button *edit_mod_btn=(Fl_Button *)0;
Fl_Check_Button *rm_emp_bok_check=(Fl_Check_Button *)0;
Fl_Browser *mod_para_brw=(Fl_Browser *)0;
Fl_Button *mod_file_out_btn=(Fl_Button *)0;
Fl_Output *mod_out_file_output=(Fl_Output *)0;
Fl_Group *forward_tab=(Fl_Group *)0;
Fl_Group *grav_group=(Fl_Group *)0;
Fl_Check_Button *Vz_check=(Fl_Check_Button *)0;
Fl_Check_Button *Vzx_check=(Fl_Check_Button *)0;
Fl_Check_Button *Vzy_check=(Fl_Check_Button *)0;
Fl_Check_Button *Vzz_check=(Fl_Check_Button *)0;
Fl_Button *mod_file_btn=(Fl_Button *)0;
Fl_Button *obs_file_btn=(Fl_Button *)0;
Fl_Input *mod_file_input=(Fl_Input *)0;
Fl_Output *mod_file_output=(Fl_Output *)0;
Fl_Input *obs_file_input=(Fl_Input *)0;
Fl_Output *obs_file_output=(Fl_Output *)0;
Fl_Input *mod_ele_input=(Fl_Input *)0;
Fl_Group *mag_group=(Fl_Group *)0;
Fl_Check_Button *DeltaT_check=(Fl_Check_Button *)0;
Fl_Check_Button *DeltaTx_check=(Fl_Check_Button *)0;
Fl_Check_Button *DeltaTy_check=(Fl_Check_Button *)0;
Fl_Check_Button *DeltaTz_check=(Fl_Check_Button *)0;
Fl_Check_Button *Hax_check=(Fl_Check_Button *)0;
Fl_Check_Button *Hay_check=(Fl_Check_Button *)0;
Fl_Check_Button *Za_check=(Fl_Check_Button *)0;
Fl_Input *mag_para_input=(Fl_Input *)0;
Fl_Check_Button *mag_data_check=(Fl_Check_Button *)0;
Fl_Button *cal_btn=(Fl_Button *)0;
Fl_Input *noise_para_input=(Fl_Input *)0;
Fl_Check_Button *noise_check=(Fl_Check_Button *)0;
Fl_Input *res_file_input=(Fl_Input *)0;
Fl_Button *res_file_btn=(Fl_Button *)0;
Fl_Output *res_out_file_output=(Fl_Output *)0;
int main(int argc, char **argv) {
{ main_window = new Fl_Double_Window(500, 600, "gm3d");
{ main_tabs = new Fl_Tabs(10, 10, 480, 580);
{ model_tab = new Fl_Group(10, 40, 480, 550, "Build Model");
{ mesh_para_input = new Fl_Input(40, 70, 300, 28, "Input Mesh Parameters :");
mesh_para_input->callback((Fl_Callback*)cb_mesh_para_input);
mesh_para_input->align(Fl_Align(FL_ALIGN_TOP_LEFT));
mesh_para_input->tooltip("mesh parameters:\n<xmin>/<dx>/<xmax>/<ymin>/<dy>/<ymax>/<zmin>/<dz>/<zmax>");
} // Fl_Input* mesh_para_input
{ mesh_file_btn = new Fl_Button(360, 70, 100, 28, "Mesh File");
mesh_file_btn->callback((Fl_Callback*)cb_mesh_file_btn);
mesh_file_btn->tooltip("read mesh parameters from an input file.");
} // Fl_Button* mesh_file_btn
{ mod_para_file_btn = new Fl_Button(170, 160, 160, 28, "Add Model From File");
mod_para_file_btn->callback((Fl_Callback*)cb_mod_para_file_btn);
mod_para_file_btn->tooltip("read model parameters from an input file.");
} // Fl_Button* mod_para_file_btn
{ mod_ele_input_build = new Fl_Input(40, 445, 215, 28, "Input Model Element Data Name:");
mod_ele_input_build->align(Fl_Align(FL_ALIGN_TOP_LEFT));
mod_ele_input_build->tooltip("a string. e.g. my_model");
} // Fl_Input* mod_ele_input_build
{ build_mod_btn = new Fl_Button(40, 545, 420, 28, "Build Model !");
build_mod_btn->callback((Fl_Callback*)cb_build_mod_btn);
build_mod_btn->tooltip("callback a gm3d function to build the model.");
} // Fl_Button* build_mod_btn
{ mesh_para_output = new Fl_Output(40, 120, 300, 28, "Mesh Parameters :");
mesh_para_output->value("Unset");
mesh_para_output->box(FL_FLAT_BOX);
mesh_para_output->color(FL_BACKGROUND_COLOR);
mesh_para_output->align(Fl_Align(FL_ALIGN_TOP_LEFT));
} // Fl_Output* mesh_para_output
{ edit_mod_btn = new Fl_Button(40, 160, 100, 28, "Model Editor");
edit_mod_btn->callback((Fl_Callback*)cb_edit_mod_btn);
edit_mod_btn->tooltip("add or edit model in a gui interface.");
} // Fl_Button* edit_mod_btn
{ rm_emp_bok_check = new Fl_Check_Button(280, 447, 180, 28, "Remove Empty Blocks");
rm_emp_bok_check->down_box(FL_DOWN_BOX);
rm_emp_bok_check->callback((Fl_Callback*)cb_rm_emp_bok_check);
rm_emp_bok_check->tooltip("suppress the output of blocks with zero values.");
} // Fl_Check_Button* rm_emp_bok_check
{ mod_para_brw = new Fl_Browser(40, 205, 420, 210);
mod_para_brw->callback((Fl_Callback*)cb_mod_para_brw);
} // Fl_Browser* mod_para_brw
{ mod_file_out_btn = new Fl_Button(360, 500, 100, 28, "Model File");
mod_file_out_btn->callback((Fl_Callback*)cb_mod_file_out_btn);
mod_file_out_btn->tooltip("create or choose a filename for outputting.");
} // Fl_Button* mod_file_out_btn
{ mod_out_file_output = new Fl_Output(40, 500, 300, 28, "Output File Name :");
mod_out_file_output->box(FL_FLAT_BOX);
mod_out_file_output->color(FL_BACKGROUND_COLOR);
mod_out_file_output->align(Fl_Align(FL_ALIGN_TOP_LEFT));
} // Fl_Output* mod_out_file_output
model_tab->end();
} // Fl_Group* model_tab
{ forward_tab = new Fl_Group(10, 40, 480, 550, "Forward Modeling");
forward_tab->hide();
{ grav_group = new Fl_Group(40, 340, 190, 58, "Forward Gravitational Data :");
grav_group->align(Fl_Align(FL_ALIGN_TOP_LEFT));
{ Vz_check = new Fl_Check_Button(40, 340, 50, 28, "Vz");
Vz_check->down_box(FL_DOWN_BOX);
} // Fl_Check_Button* Vz_check
{ Vzx_check = new Fl_Check_Button(105, 340, 50, 28, "Vzx");
Vzx_check->down_box(FL_DOWN_BOX);
} // Fl_Check_Button* Vzx_check
{ Vzy_check = new Fl_Check_Button(170, 340, 50, 28, "Vzy");
Vzy_check->down_box(FL_DOWN_BOX);
} // Fl_Check_Button* Vzy_check
{ Vzz_check = new Fl_Check_Button(40, 370, 50, 28, "Vzz");
Vzz_check->down_box(FL_DOWN_BOX);
} // Fl_Check_Button* Vzz_check
grav_group->end();
} // Fl_Group* grav_group
{ mod_file_btn = new Fl_Button(360, 70, 100, 28, "Model File");
mod_file_btn->callback((Fl_Callback*)cb_mod_file_btn);
} // Fl_Button* mod_file_btn
{ obs_file_btn = new Fl_Button(360, 130, 100, 28, "Observe File");
obs_file_btn->callback((Fl_Callback*)cb_obs_file_btn);
} // Fl_Button* obs_file_btn
{ mod_file_input = new Fl_Input(40, 70, 300, 28, "Input Model FIle :");
mod_file_input->callback((Fl_Callback*)cb_mod_file_input);
mod_file_input->align(Fl_Align(FL_ALIGN_TOP_LEFT));
} // Fl_Input* mod_file_input
{ mod_file_output = new Fl_Output(40, 240, 200, 28, "Chosen Model File :");
mod_file_output->box(FL_FLAT_BOX);
mod_file_output->color(FL_BACKGROUND_COLOR);
mod_file_output->align(Fl_Align(FL_ALIGN_TOP_LEFT));
mod_file_output->value("Unset");
} // Fl_Output* mod_file_output
{ obs_file_input = new Fl_Input(40, 130, 300, 28, "Input Observe FIle :");
obs_file_input->callback((Fl_Callback*)cb_obs_file_input);
obs_file_input->align(Fl_Align(FL_ALIGN_TOP_LEFT));
} // Fl_Input* obs_file_input
{ obs_file_output = new Fl_Output(260, 240, 200, 28, "Chosen Observe File :");
obs_file_output->box(FL_FLAT_BOX);
obs_file_output->color(FL_BACKGROUND_COLOR);
obs_file_output->align(Fl_Align(FL_ALIGN_TOP_LEFT));
obs_file_output->value("Unset");
} // Fl_Output* obs_file_output
{ mod_ele_input = new Fl_Input(260, 315, 200, 28, "Model Element Data Name :");
mod_ele_input->align(Fl_Align(FL_ALIGN_TOP_LEFT));
} // Fl_Input* mod_ele_input
{ mag_group = new Fl_Group(40, 440, 420, 90);
{ DeltaT_check = new Fl_Check_Button(40, 475, 65, 25, "DeltaT");
DeltaT_check->down_box(FL_DOWN_BOX);
DeltaT_check->deactivate();
} // Fl_Check_Button* DeltaT_check
{ DeltaTx_check = new Fl_Check_Button(120, 475, 70, 25, "DeltaTx");
DeltaTx_check->down_box(FL_DOWN_BOX);
DeltaTx_check->deactivate();
} // Fl_Check_Button* DeltaTx_check
{ DeltaTy_check = new Fl_Check_Button(200, 475, 70, 25, "DeltaTy");
DeltaTy_check->down_box(FL_DOWN_BOX);
DeltaTy_check->deactivate();
} // Fl_Check_Button* DeltaTy_check
{ DeltaTz_check = new Fl_Check_Button(280, 475, 70, 25, "DeltaTz");
DeltaTz_check->down_box(FL_DOWN_BOX);
DeltaTz_check->deactivate();
} // Fl_Check_Button* DeltaTz_check
{ Hax_check = new Fl_Check_Button(40, 505, 65, 25, "Hax");
Hax_check->down_box(FL_DOWN_BOX);
Hax_check->deactivate();
} // Fl_Check_Button* Hax_check
{ Hay_check = new Fl_Check_Button(120, 505, 65, 25, "Hay");
Hay_check->down_box(FL_DOWN_BOX);
Hay_check->deactivate();
} // Fl_Check_Button* Hay_check
{ Za_check = new Fl_Check_Button(200, 505, 65, 25, "Za");
Za_check->down_box(FL_DOWN_BOX);
Za_check->deactivate();
} // Fl_Check_Button* Za_check
{ mag_para_input = new Fl_Input(220, 440, 240, 28, "Magnetization Parameters : ");
mag_para_input->tooltip("<I0>/<D0>/<I>/<D>");
mag_para_input->callback((Fl_Callback*)cb_mag_para_input);
mag_para_input->deactivate();
} // Fl_Input* mag_para_input
mag_group->end();
} // Fl_Group* mag_group
{ mag_data_check = new Fl_Check_Button(40, 405, 170, 30, "Forward Magnetic Data");
mag_data_check->down_box(FL_DOWN_BOX);
mag_data_check->callback((Fl_Callback*)cb_mag_data_check);
} // Fl_Check_Button* mag_data_check
{ cal_btn = new Fl_Button(40, 545, 420, 28, "Calculate !");
cal_btn->callback((Fl_Callback*)cb_cal_btn);
} // Fl_Button* cal_btn
{ noise_para_input = new Fl_Input(260, 370, 200, 28);
noise_para_input->tooltip("<noise-mean>/<nosie-deviation>");
noise_para_input->align(Fl_Align(FL_ALIGN_TOP_LEFT));
noise_para_input->deactivate();
} // Fl_Input* noise_para_input
{ noise_check = new Fl_Check_Button(260, 345, 180, 30, "Input Noise Parameters :");
noise_check->down_box(FL_DOWN_BOX);
noise_check->callback((Fl_Callback*)cb_noise_check);
} // Fl_Check_Button* noise_check
{ res_file_input = new Fl_Input(40, 190, 300, 28, "Input Prefix of Output FIle :");
res_file_input->callback((Fl_Callback*)cb_res_file_input);
res_file_input->align(Fl_Align(FL_ALIGN_TOP_LEFT));
} // Fl_Input* res_file_input
{ res_file_btn = new Fl_Button(360, 190, 100, 28, "Result File");
res_file_btn->callback((Fl_Callback*)cb_res_file_btn);
} // Fl_Button* res_file_btn
{ res_out_file_output = new Fl_Output(40, 290, 200, 28, "Prefix of Output File :");
res_file_output->box(FL_FLAT_BOX);
res_file_output->color(FL_BACKGROUND_COLOR);
res_file_output->align(Fl_Align(FL_ALIGN_TOP_LEFT));
res_file_output->value("Unset");
} // Fl_Output* res_out_file_output
forward_tab->end();
} // Fl_Group* forward_tab
main_tabs->end();
} // Fl_Tabs* main_tabs
main_window->end();
} // Fl_Double_Window* main_window
main_window->show(argc, argv);
return Fl::run();
}

271
fluid_project/gm3d_gui2.fl
View File

@ -1,271 +0,0 @@
# data file for the Fltk User Interface Designer (fluid)
version 1.0305
header_name {.h}
code_name {.cxx}
Function {cb_edit_mod_btn(Fl_Widget*, void*)} {open
} {}
Function {cb_mod_para_file_btn(Fl_Widget*, void*)} {open
} {}
Function {cb_mod_para_save_btn(Fl_Widget*, void*)} {open
} {}
Function {cb_mod_para_brw(Fl_Browser*, void*)} {open
} {}
Function {cb_mesh_mod_savename(Fl_Input*, void*)} {open
} {}
Function {cb_rm_emp_bok_check(Fl_Check_Button*, void*)} {open
} {}
Function {cb_build_mod_btn(Fl_Widget*, void*)} {open
} {}
Function {cb_mod_file_out_btn(Fl_Widget*, void*)} {open
} {}
Function {cb_mod_file_btn(Fl_Widget*, void*)} {open
} {}
Function {cb_obs_file_input(Fl_Input*, void*)} {open
} {}
Function {cb_obs_file_btn(Fl_Widget*, void*)} {open
} {}
Function {cb_res_file_btn(Fl_Widget*, void*)} {open
} {}
Function {cb_nosie_check(Fl_Check_Button*, void*)} {open
} {}
Function {cb_noise_para_input(Fl_Input*, void*)} {open
} {}
Function {cb_mag_data_check(Fl_Check_Button*, void*)} {open
} {}
Function {cb_mag_para_input(Fl_Input*, void*)} {open
} {}
Function {} {open
} {
Fl_Window main_window {
label gm3d open
xywh {616 107 500 580} type Double visible
} {
Fl_Tabs main_tabs {open
xywh {10 10 480 560}
} {
Fl_Group model_tab {
label {Build 3D Mesh Model}
xywh {10 40 480 530} hide
} {
Fl_Group {} {
label {3D Mesh:} open
xywh {40 70 340 75} align 5
} {
Fl_Button mesh_para_file_btn {
label {Import From}
xywh {160 75 100 28}
}
Fl_Button edit_mesh_btn {
label Editor
xywh {40 75 100 28}
}
Fl_Button mesh_para_save_btn {
label {Save As}
xywh {280 75 100 28} deactivate
}
Fl_Output mesh_para_output {
label {Mesh:}
xywh {80 115 300 28} box FLAT_BOX
}
}
Fl_Group {} {
label {Physical Entities:}
xywh {40 175 369 40} align 5
} {
Fl_Button edit_mod_btn {
label Editor
callback cb_edit_mod_btn
tooltip {Open the model editor to add or edit a model.} xywh {40 180 100 28}
}
Fl_Button mod_para_file_btn {
label {Import From}
callback cb_mod_para_file_btn
tooltip {Import models from a text file.} xywh {160 180 100 28}
}
Fl_Button mod_para_save_btn {
label {Save As}
callback cb_mod_para_save_btn
tooltip {Save models to a text file.} xywh {280 180 100 28} deactivate
}
}
Fl_Browser mod_para_brw {
label {Added Physical Entities}
callback cb_mod_para_brw
tooltip {Select a model to edit or remove it using the model editor.} xywh {40 240 420 205} align 5
}
Fl_Input mesh_model_save_name {
label {Mesh Model's Name:}
callback cb_mesh_mod_savename
tooltip {Name of the meshed 3D model.} xywh {40 475 215 28} align 5
}
Fl_Check_Button rm_emp_bok_check {
label {Remove Empty Blocks}
callback cb_rm_emp_bok_check
tooltip {Remove blocks with zero value from output mesh model files.} xywh {280 475 180 28} down_box DOWN_BOX
}
Fl_Button build_mod_btn {
label {Build 3D Model}
callback cb_build_mod_btn
tooltip {Build 3D meshed model uisng entity parameters listed in the browser.} xywh {40 520 300 28} deactivate
}
Fl_Button mod_file_out_btn {
label {Save As}
callback cb_mod_file_out_btn
tooltip {Save 3D meshed model to a file.} xywh {360 520 100 28} deactivate
}
}
Fl_Group forward_tab {
label {Forward Modeling} selected
xywh {10 40 480 530} box ROUNDED_FRAME
} {
Fl_Group {} {
label {Mesh Model:} open
xywh {40 70 120 35} align 5
} {
Fl_Button mod_file_btn {
label {Import From}
callback cb_mod_file_btn
xywh {40 75 110 28}
}
}
Fl_Output mesh_mod_disp {
label {Mesh Model:}
xywh {120 110 340 28} box FLAT_BOX color 49
}
Fl_Output obs_file_disp {
label {Observations:}
xywh {125 135 335 28} box FLAT_BOX color 49
}
Fl_Group {} {
label {Observations:}
xywh {210 70 250 35} align 5
} {
Fl_Button obs_file_btn {
label {Import From}
callback cb_obs_file_btn
xywh {340 75 120 28}
}
Fl_Button {} {
label Editor
xywh {210 75 100 28}
}
}
Fl_Input mod_ele_input {
label {Model Data Name:}
callback cb_mod_ele_input
xywh {155 175 305 28}
}
Fl_Group grav_group {
label {Forward Gravitational Data:}
xywh {40 235 190 75} align 5
} {
Fl_Check_Button Vz_check {
label Vz
xywh {40 242 50 28} down_box DOWN_BOX
}
Fl_Check_Button Vzx_check {
label Vzx
xywh {105 242 50 28} down_box DOWN_BOX
}
Fl_Check_Button Vzy_check {
label Vzy
xywh {170 242 50 28} down_box DOWN_BOX
}
Fl_Check_Button Vzz_check {
label Vzz
xywh {40 275 50 28} down_box DOWN_BOX
}
}
Fl_Group mag_group {
xywh {40 346 425 114}
} {
Fl_Check_Button DeltaT_check {
label DeltaT
xywh {40 394 60 28} down_box DOWN_BOX deactivate
}
Fl_Check_Button DeltaTx_check {
label DeltaTx
xywh {115 394 70 28} down_box DOWN_BOX deactivate
}
Fl_Check_Button DeltaTy_check {
label DeltaTy
xywh {195 394 70 28} down_box DOWN_BOX deactivate
}
Fl_Check_Button DeltaTz_check {
label DeltaTz
xywh {275 394 70 28} down_box DOWN_BOX deactivate
}
Fl_Check_Button Hax_check {
label Hax
xywh {40 429 60 28} down_box DOWN_BOX deactivate
}
Fl_Check_Button Hay_check {
label Hay
xywh {115 429 65 28} down_box DOWN_BOX deactivate
}
Fl_Check_Button Za_check {
label Za
xywh {195 429 65 28} down_box DOWN_BOX deactivate
}
Fl_Button {} {
label Magnetization
xywh {40 351 110 28} deactivate
}
Fl_Output {} {
xywh {160 351 300 28} box FLAT_BOX deactivate
}
}
Fl_Check_Button mag_data_check {
label {Forward Magnetic Data}
callback cb_mag_data_check
xywh {40 312 170 28} down_box DOWN_BOX
}
Fl_Button cal_btn {
label Calculate
callback cb_cal_btn
xywh {40 470 300 28} deactivate
}
Fl_Input noise_para_input {
label {Mean:}
callback cb_noise_para_input
tooltip {<noise-mean>/<nosie-deviation>} xywh {300 242 160 28} deactivate
}
Fl_Check_Button noise_check {
label {Add Noises:}
callback cb_noise_check
xywh {255 213 105 28} down_box DOWN_BOX
}
Fl_Button res_file_btn {
label {Save As}
callback cb_res_file_btn
xywh {360 470 100 28} deactivate
}
Fl_Progress cal_progress_bar {
label {Calculating...}
xywh {40 520 420 28} selection_color 139
}
Fl_Input noise_para_input {
label {Deviation:}
callback cb_noise_para_input
tooltip {<noise-mean>/<nosie-deviation>} xywh {325 282 135 28} deactivate
}
}
}
}
}

83
fluid_project/gm3d_gui2.h
View File

@ -1,83 +0,0 @@
// generated by Fast Light User Interface Designer (fluid) version 1.0305
#ifndef _GM3D_GUI_H
#define _GM3D_GUI_H
#include <FL/Fl.H>
#include <FL/Fl_Double_Window.H>
#include <FL/Fl_Tabs.H>
#include <FL/Fl_Group.H>
#include <FL/Fl_Input.H>
#include <FL/Fl_Button.H>
#include <FL/Fl_Output.H>
#include <FL/Fl_Check_Button.H>
#include <FL/Fl_Browser.H>
#include <FL/Fl_File_Chooser.H>
#include <iostream>
#include <string.h>
#include <sstream>
#include <fstream>
extern void cb_mesh_para_input(Fl_Input*, void*);
extern void cb_mesh_file_btn(Fl_Widget*, void*);
extern void cb_edit_mod_btn(Fl_Widget*, void*);
extern void cb_mod_para_file_btn(Fl_Widget*, void*);
extern void cb_mod_para_brw(Fl_Browser*, void*);
extern void cb_rm_emp_bok_check(Fl_Check_Button*, void*);
extern void cb_mod_file_out_btn(Fl_Widget*, void*);
extern void cb_build_mod_btn(Fl_Widget*, void*);
extern void cb_mod_file_input(Fl_Input*, void*);
extern void cb_mod_file_btn(Fl_Widget*, void*);
extern void cb_obs_file_input(Fl_Input*, void*);
extern void cb_obs_file_btn(Fl_Widget*, void*);
extern void cb_res_file_input(Fl_Input*, void*);
extern void cb_res_file_btn(Fl_Widget*, void*);
extern void cb_nosie_check(Fl_Check_Button*, void*);
extern void cb_mag_data_check(Fl_Check_Button*, void*);
extern void cb_mag_para_input(Fl_Input*, void*);
extern void cb_cal_btn(Fl_Button*, void*);
extern Fl_Double_Window *main_window;
extern Fl_Tabs *main_tabs;
extern Fl_Group *model_tab;
extern Fl_Input *mesh_para_input;
extern Fl_Button *mesh_file_btn;
extern Fl_Button *mod_para_file_btn;
extern Fl_Input *mod_ele_input_build;
extern Fl_Button *build_mod_btn;
extern Fl_Output *mesh_para_output;
extern Fl_Button *edit_mod_btn;
extern Fl_Check_Button *rm_emp_bok_check;
extern Fl_Browser *mod_para_brw;
extern Fl_Button *mod_file_out_btn;
extern Fl_Output *mod_out_file_output;
extern Fl_Group *forward_tab;
extern Fl_Group *grav_group;
extern Fl_Check_Button *Vz_check;
extern Fl_Check_Button *Vzx_check;
extern Fl_Check_Button *Vzy_check;
extern Fl_Check_Button *Vzz_check;
extern Fl_Button *mod_file_btn;
extern Fl_Button *obs_file_btn;
extern Fl_Input *mod_file_input;
extern Fl_Output *mod_file_output;
extern Fl_Input *obs_file_input;
extern Fl_Output *obs_file_output;
extern Fl_Input *mod_ele_input;
extern Fl_Group *mag_group;
extern Fl_Check_Button *DeltaT_check;
extern Fl_Check_Button *DeltaTx_check;
extern Fl_Check_Button *DeltaTy_check;
extern Fl_Check_Button *DeltaTz_check;
extern Fl_Check_Button *Hax_check;
extern Fl_Check_Button *Hay_check;
extern Fl_Check_Button *Za_check;
extern Fl_Input *mag_para_input;
extern Fl_Check_Button *mag_data_check;
extern Fl_Button *cal_btn;
extern Fl_Input *noise_para_input;
extern Fl_Check_Button *noise_check;
extern Fl_Input *res_file_input;
extern Fl_Button *res_file_btn;
extern Fl_Output *res_file_output;
#endif //_GM3D_GUI_H

BIN
logo/Logo_ZY.icns
View File

Binary file not shown.

BIN
logo/Logo_ZY.iconset/icon_128x128.png
View File

Binary file not shown.
Side by Side

Before

Width:  |  Height:  |  Size: 28 KiB

BIN
logo/Logo_ZY.iconset/icon_128x128@2x.png
View File

Binary file not shown.
Side by Side

Before

Width:  |  Height:  |  Size: 79 KiB

BIN
logo/Logo_ZY.iconset/icon_16x16.png
View File

Binary file not shown.
Side by Side

Before

Width:  |  Height:  |  Size: 1008 B

BIN
logo/Logo_ZY.iconset/icon_16x16@2x.png
View File

Binary file not shown.
Side by Side

Before

Width:  |  Height:  |  Size: 3.1 KiB

BIN
logo/Logo_ZY.iconset/icon_256x256.png
View File

Binary file not shown.
Side by Side

Before

Width:  |  Height:  |  Size: 79 KiB

BIN
logo/Logo_ZY.iconset/icon_256x256@2x.png
View File

Binary file not shown.
Side by Side

Before

Width:  |  Height:  |  Size: 226 KiB

BIN
logo/Logo_ZY.iconset/icon_32x32.png
View File

Binary file not shown.
Side by Side

Before

Width:  |  Height:  |  Size: 3.1 KiB

BIN
logo/Logo_ZY.iconset/icon_32x32@2x.png
View File

Binary file not shown.
Side by Side

Before

Width:  |  Height:  |  Size: 9.6 KiB

BIN
logo/Logo_ZY.iconset/icon_512x512.png
View File

Binary file not shown.
Side by Side

Before

Width:  |  Height:  |  Size: 226 KiB

BIN
logo/Logo_ZY.iconset/icon_512x512@2x.png
View File

Binary file not shown.
Side by Side

Before

Width:  |  Height:  |  Size: 668 KiB

BIN
logo/Logo_ZY.iconset/icon_64x64.png
View File

Binary file not shown.
Side by Side

Before

Width:  |  Height:  |  Size: 9.6 KiB

BIN
logo/Logo_ZY.iconset/icon_64x64@2x.png
View File

Binary file not shown.
Side by Side

Before

Width:  |  Height:  |  Size: 28 KiB

BIN
logo/logo-zy.png
View File

Binary file not shown.
Side by Side

Before

Width:  |  Height:  |  Size: 159 KiB

21
logo/png2icns.sh
View File

@ -1,21 +0,0 @@
#/bin/bash
png_file=${1}
icns_name=${2}
mkdir ${icns_name}".iconset"
sips -z 16 16 ${png_file} --out ${icns_name}.iconset/icon_16x16.png
sips -z 32 32 ${png_file} --out ${icns_name}.iconset/icon_16x16@2x.png
sips -z 32 32 ${png_file} --out ${icns_name}.iconset/icon_32x32.png
sips -z 64 64 ${png_file} --out ${icns_name}.iconset/icon_32x32@2x.png
sips -z 64 64 ${png_file} --out ${icns_name}.iconset/icon_64x64.png
sips -z 128 128 ${png_file} --out ${icns_name}.iconset/icon_64x64@2x.png
sips -z 128 128 ${png_file} --out ${icns_name}.iconset/icon_128x128.png
sips -z 256 256 ${png_file} --out ${icns_name}.iconset/icon_128x128@2x.png
sips -z 256 256 ${png_file} --out ${icns_name}.iconset/icon_256x256.png
sips -z 512 512 ${png_file} --out ${icns_name}.iconset/icon_256x256@2x.png
sips -z 512 512 ${png_file} --out ${icns_name}.iconset/icon_512x512.png
sips -z 1024 1024 ${png_file} --out ${icns_name}.iconset/icon_512x512@2x.png
iconutil -c icns ${icns_name}.iconset

130183
logo/syn_model_out.ai
View File

File diff suppressed because one or more lines are too long

809657
logo/syn_model_out.eps
View File

File diff suppressed because it is too large Load Diff

11
routine.sh
View File

@ -1,11 +0,0 @@
#!/bin/bash
if [[ ${1} == "configure" && ! -d "build/" ]]; then
mkdir build && cd build && rm -rf * && cmake ..
elif [[ ${1} == "configure" ]]; then
cd build && rm -rf * && cmake ..
elif [[ ${1} == "build" ]]; then
cd build && make
elif [[ ${1} == "install" ]]; then
cd build && make install
fi

22
src/CMakeLists.txt
View File

@ -1,10 +1,20 @@
aux_source_directory(. DIR_SRC)
add_executable(gm3d ${DIR_SRC})
set(EXECUTABLE_OUTPUT_PATH ${PROJECT_BINARY_DIR}/bin)
# 添加头文件地址
include_directories(E:/Library/include)
# 添加动态库地址
find_library(FLTK_LIBRARY fltk HINTS E:/Library/lib)
find_package(OpenMP)
if (OPENMP_FOUND)
set (CMAKE_C_FLAGS "${CMAKE_C_FLAGS} ${OpenMP_C_FLAGS}")
set (CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} ${OpenMP_CXX_FLAGS}")
set (CMAKE_EXE_LINKER_FLAGS "${CMAKE_EXE_LINKER_FLAGS} ${OpenMP_EXE_LINKER_FLAGS}")
endif()
# 生成windows和mac的应用程序文件
add_executable(gm3d_gui WIN32 MACOSX_BUNDLE ${DIR_SRC})
target_link_libraries(gm3d ${FLTK_LIBRARY})
include_directories(/usr/local/include)
find_library(HINTS /usr/local/lib)
#find_library(FLTK_LIBRARY fltk HINTS /usr/local/lib)
target_link_libraries(gm3d_gui PUBLIC ${FLTK_LIBRARY})

39
src/CMakeLists.txt.backup
View File

@ -1,39 +0,0 @@
aux_source_directory(. DIR_SRC)
set(EXECUTABLE_OUTPUT_PATH ${PROJECT_BINARY_DIR}/bin)
set(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} --std=c++11 -O2")
# NOTE: Don't include the path in MACOSX_BUNDLE_ICON_FILE -- this is
# the property added to Info.plist
# <key>CFBundleIconFile</key>
# <string>logo</string>
set(MACOSX_BUNDLE_ICON_FILE Logo_ZY)
# And this part tells CMake where to find and install the file itself
set(myApp_ICON ${PROJECT_SOURCE_DIR}/assert/logo/Logo_ZY.icns)
# set_source_files_properties would create Resources folder to store icns
# The folder is at same path where Info.plist located
set_source_files_properties(${myApp_ICON} PROPERTIES MACOSX_PACKAGE_LOCATION "Resources")
#add_executable(gm3d_gui WIN32 MACOSX_BUNDLE ${DIR_SRC})
add_executable(gm3d MACOSX_BUNDLE ${DIR_SRC} ${myApp_ICON})
# 添加openmp的编译命令 设置编译选项
find_package(OpenMP REQUIRED)
if (OpenMP_CXX_FOUND)
message(STATUS "OpenMP Found.")
set(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} ${OpenMP_CXX_FLAGS}")
set(CMAKE_EXE_LINKER_FLAGS "${CMAKE_EXE_LINKER_FLAGS} ${OpenMP_EXE_LINKER_FLAGS}")
set(CMAKE_SHARED_LINKER_FLAGS "${CMAKE_SHARED_LINKER_FLAGS} ${OpenMP_SHARED_LINKER_FLAGS}")
#target_link_libraries(gm3d PUBLIC OpenMP::OpenMP_CXX)
endif()
#find_library(FLTK_LIBRARY libfltk.a /usr/local/Cellar/fltk/1.3.5/lib)
# 添加头文件地址
include_directories(/opt/homebrew/include)
include_directories(/usr/local/include)
# 添加动态库地址
find_library(FLTK_LIBRARY fltk HINTS /opt/homebrew/lib)
target_link_libraries(gm3d ${FLTK_LIBRARY} OpenMP::OpenMP_CXX)

8
src/add_interface_block.cpp
View File

@ -1,5 +1,5 @@
#include "gm3d.h"
//我们读入一个界面数据 插值计算每个块体中心位置的值 然后按情况赋值
//我们读入一个界面数据 插值计算每个块体中心位置的值 然后按情况赋值
int GM3D::AddInterfaceBlock(modelist para_list){
int m, n, xnum, ynum;
double xs,xe,xmin,xmax,dx;
@ -18,7 +18,7 @@ int GM3D::AddInterfaceBlock(modelist para_list){
if (open_infile(infile,filename)) return -1;
while (getline(infile,temp_str)){
//#range必须出现在数据之前
//#range必须出现在数据之前
if (*(temp_str.begin()) == '#'){
if (6 == sscanf(temp_str.c_str(),"# range=%lf/%lf/%lf/%lf/%lf/%lf",&xs,&dx,&xe,&ys,&dy,&ye)){
xmin = MIN(xs,xe); xmax = MAX(xs,xe);
@ -75,7 +75,7 @@ int GM3D::AddInterfaceBlock(modelist para_list){
if (model_cube_[i].cen.z >= temp_topo.z){
if (model_block_val_[i] == BDL_MAX)
model_block_val_[i] = para_list.mod_value; //注意重复赋值的块体会覆盖
model_block_val_[i] = para_list.mod_value; //注意重复赋值的块体会覆盖
else
model_block_val_[i] += para_list.mod_value;
model_added = true;
@ -92,7 +92,7 @@ int GM3D::AddInterfaceBlock(modelist para_list){
if (model_cube_[i].cen.z <= temp_topo.z){
if (model_block_val_[i] == BDL_MAX)
model_block_val_[i] = para_list.mod_value; //注意重复赋值的块体会覆盖
model_block_val_[i] = para_list.mod_value; //注意重复赋值的块体会覆盖
else
model_block_val_[i] += para_list.mod_value;
model_added = true;

8
archive/add_model_gui.cpp → src/add_model_gui.cpp
View File

@ -18,7 +18,7 @@ Fl_Return_Button *can_add_btn=(Fl_Return_Button *)0;
Fl_Button *sig_add_btn=(Fl_Button *)0;
Fl_Input *mod_val_input=(Fl_Input *)0;
// 声明三个字符串用于保存新建的模型类型和赋值类型
// 声明三个字符串用于保存新建的模型类型和赋值类型
char add_mod_type[1024];
char add_val_type[1024];
char agn_val_type[1024];
@ -81,8 +81,8 @@ void cb_bot_val_rbtn(Fl_Button*, void*){
return;
}
// 声明一个字符串用于保存新建的模型参数并添加至列表中
// 模型参数的排列为<模型类型> <赋值类型> <物理参数字符串> <几何参数字符串>
// 声明一个字符串用于保存新建的模型参数并添加至列表中
// 模型参数的排列为<模型类型> <赋值类型> <物理参数字符串> <几何参数字符串>
char add_mod_para[1024];
void cb_sig_add_btn(Fl_Widget*, void*){
@ -101,7 +101,7 @@ void cb_sig_add_btn(Fl_Widget*, void*){
mod_para_brw->add(add_mod_para);
return;
}
//取消添加模型
//取消添加模型
void cb_can_add_btn(Fl_Widget*, void*){
add_mod_win->hide();
return;

2
src/add_models.cpp
View File

@ -10,7 +10,7 @@ int GM3D::AddModels(char* filename){
while(getline(infile,temp_str)){
if (*(temp_str.begin()) == '#') continue;
else{
//按每行5个数据解析 初始化为含观测值与不确定度的观测点
//按每行5个数据解析 初始化为含观测值与不确定度的观测点
if (4 == sscanf(temp_str.c_str(),"%s %s %lf %s",
temp_list.mod_type,temp_list.val_type,&temp_list.mod_value,temp_list.mod_para)){
model_list_.push_back(temp_list);

6
src/add_regular_block.cpp
View File

@ -16,7 +16,7 @@ int GM3D::AddRegularBlock(modelist para_list){
if (model_cube_[i].cen.x >= xmin && model_cube_[i].cen.x <= xmax &&
model_cube_[i].cen.y >= ymin && model_cube_[i].cen.y <= ymax &&
model_cube_[i].cen.z >= zmin && model_cube_[i].cen.z <= zmax){
model_block_val_[i] = para_list.mod_value; //注意重复赋值的块体会覆盖
model_block_val_[i] = para_list.mod_value; //注意重复赋值的块体会覆盖
model_added = true;
}
}
@ -27,7 +27,7 @@ int GM3D::AddRegularBlock(modelist para_list){
model_cube_[i].cen.y >= ymin && model_cube_[i].cen.y <= ymax &&
model_cube_[i].cen.z >= zmin && model_cube_[i].cen.z <= zmax){
if (model_block_val_[i] == BDL_MAX)
model_block_val_[i] = para_list.mod_value; //注意重复赋值的块体会覆盖
model_block_val_[i] = para_list.mod_value; //注意重复赋值的块体会覆盖
else
model_block_val_[i] += para_list.mod_value;
model_added = true;
@ -39,7 +39,7 @@ int GM3D::AddRegularBlock(modelist para_list){
if (model_cube_[i].cen.x >= xmin && model_cube_[i].cen.x <= xmax &&
model_cube_[i].cen.y >= ymin && model_cube_[i].cen.y <= ymax &&
model_cube_[i].cen.z >= zmin && model_cube_[i].cen.z <= zmax){
model_block_val_[i] = BDL_MAX; //注意重复赋值的块体会覆盖
model_block_val_[i] = BDL_MAX; //注意重复赋值的块体会覆盖
model_added = true;
}
}

6
src/add_sphere_block.cpp
View File

@ -20,7 +20,7 @@ int GM3D::AddSphereBlock(modelist para_list){
rad_limit = rad_x*rad_y*rad_z/sqrt(pow(rad_y*rad_z*sin(theta)*cos(phi),2) + pow(rad_x*rad_z*sin(theta)*sin(phi),2) + pow(rad_x*rad_y*cos(theta),2));
if (dist <= rad_limit){
model_block_val_[i] = para_list.mod_value; //注意重复赋值的块体会覆盖
model_block_val_[i] = para_list.mod_value; //注意重复赋值的块体会覆盖
model_added = true;
}
}
@ -39,7 +39,7 @@ int GM3D::AddSphereBlock(modelist para_list){
if (dist <= rad_limit){
if (model_block_val_[i] == BDL_MAX)
model_block_val_[i] = para_list.mod_value; //注意重复赋值的块体会覆盖
model_block_val_[i] = para_list.mod_value; //注意重复赋值的块体会覆盖
else
model_block_val_[i] += para_list.mod_value;
model_added = true;
@ -59,7 +59,7 @@ int GM3D::AddSphereBlock(modelist para_list){
rad_limit = rad_x*rad_y*rad_z/sqrt(pow(rad_y*rad_z*sin(theta)*cos(phi),2) + pow(rad_x*rad_z*sin(theta)*sin(phi),2) + pow(rad_x*rad_y*cos(theta),2));
if (dist <= rad_limit){
model_block_val_[i] = BDL_MAX; //注意重复赋值的块体会覆盖
model_block_val_[i] = BDL_MAX; //注意重复赋值的块体会覆盖
model_added = true;
}
}

12
src/add_tilted_block.cpp
View File

@ -19,7 +19,7 @@ int GM3D::AddTiltedBlock(modelist para_list){
if (!strcmp(para_list.val_type,"replace")){
for (int i = 0; i < model_num_; i++){
//计算当前层的x与y范围
//计算当前层的x与y范围
layer_xmin = (model_cube_[i].cen.z - zmin)*(xmin_2 - xmin_1)/(zmax - zmin) + xmin_1;
layer_xmax = (model_cube_[i].cen.z - zmin)*(xmax_2 - xmax_1)/(zmax - zmin) + xmax_1;
layer_ymin = (model_cube_[i].cen.z - zmin)*(ymin_2 - ymin_1)/(zmax - zmin) + ymin_1;
@ -28,14 +28,14 @@ int GM3D::AddTiltedBlock(modelist para_list){
if (model_cube_[i].cen.x >= layer_xmin && model_cube_[i].cen.x <= layer_xmax &&
model_cube_[i].cen.y >= layer_ymin && model_cube_[i].cen.y <= layer_ymax &&
model_cube_[i].cen.z >= zmin && model_cube_[i].cen.z <= zmax){
model_block_val_[i] = para_list.mod_value; //注意重复赋值的块体会覆盖
model_block_val_[i] = para_list.mod_value; //注意重复赋值的块体会覆盖
model_added = true;
}
}
}
else if (!strcmp(para_list.val_type,"add")){
for (int i = 0; i < model_num_; i++){
//计算当前层的x与y范围
//计算当前层的x与y范围
layer_xmin = (model_cube_[i].cen.z - zmin)*(xmin_2 - xmin_1)/(zmax - zmin) + xmin_1;
layer_xmax = (model_cube_[i].cen.z - zmin)*(xmax_2 - xmax_1)/(zmax - zmin) + xmax_1;
layer_ymin = (model_cube_[i].cen.z - zmin)*(ymin_2 - ymin_1)/(zmax - zmin) + ymin_1;
@ -45,7 +45,7 @@ int GM3D::AddTiltedBlock(modelist para_list){
model_cube_[i].cen.y >= layer_ymin && model_cube_[i].cen.y <= layer_ymax &&
model_cube_[i].cen.z >= zmin && model_cube_[i].cen.z <= zmax){
if (model_block_val_[i] == BDL_MAX)
model_block_val_[i] = para_list.mod_value; //注意重复赋值的块体会覆盖
model_block_val_[i] = para_list.mod_value; //注意重复赋值的块体会覆盖
else
model_block_val_[i] += para_list.mod_value;
model_added = true;
@ -54,7 +54,7 @@ int GM3D::AddTiltedBlock(modelist para_list){
}
else if (!strcmp(para_list.val_type,"erase")){
for (int i = 0; i < model_num_; i++){
//计算当前层的x与y范围
//计算当前层的x与y范围
layer_xmin = (model_cube_[i].cen.z - zmin)*(xmin_2 - xmin_1)/(zmax - zmin) + xmin_1;
layer_xmax = (model_cube_[i].cen.z - zmin)*(xmax_2 - xmax_1)/(zmax - zmin) + xmax_1;
layer_ymin = (model_cube_[i].cen.z - zmin)*(ymin_2 - ymin_1)/(zmax - zmin) + ymin_1;
@ -63,7 +63,7 @@ int GM3D::AddTiltedBlock(modelist para_list){
if (model_cube_[i].cen.x >= layer_xmin && model_cube_[i].cen.x <= layer_xmax &&
model_cube_[i].cen.y >= layer_ymin && model_cube_[i].cen.y <= layer_ymax &&
model_cube_[i].cen.z >= zmin && model_cube_[i].cen.z <= zmax){
model_block_val_[i] = BDL_MAX; //注意重复赋值的块体会覆盖
model_block_val_[i] = BDL_MAX; //注意重复赋值的块体会覆盖
model_added = true;
}
}

20
src/build_regular_grid.cpp
View File

@ -46,26 +46,26 @@ int GM3D::BuildRegularGrid(char* space_para){
while(x >= xmin && x <= xmax){
z = zs;
while(z >= zmin && z <= zmax){
//添加denMod
//添加denMod
temp_cu.cen.id = model_cube_.size();
temp_cu.cen.x = x; temp_cu.cen.y = y; temp_cu.cen.z = z;
//添加mshVert
//添加mshVert
for (int i = 0; i < 8; i++){
temp_cp.id = model_vert_.size(); //添加msh的顶点索引为mshVert的大小
temp_cp.x = temp_cu.cen.x - sign[i][0]*temp_cu.dx; //左下底角
temp_cp.id = model_vert_.size(); //添加msh的顶点索引为mshVert的大小
temp_cp.x = temp_cu.cen.x - sign[i][0]*temp_cu.dx; //左下底角
temp_cp.y = temp_cu.cen.y - sign[i][1]*temp_cu.dy;
temp_cp.z = temp_cu.cen.z - sign[i][2]*temp_cu.dz;
temp_id_str = cpoint_id(temp_cp);
imsp = map_str_point.find(temp_id_str);
//利用map_vert查到当前顶点是否存在,这里需要注意，如果顶点已经存在则只需要将顶点索引置为已存在顶点的索引，不增加顶点计数
//利用map_vert查到当前顶点是否存在,这里需要注意，如果顶点已经存在则只需要将顶点索引置为已存在顶点的索引，不增加顶点计数
if(imsp!=map_str_point.end()){
temp_cu.ids[i] = imsp->second.id;
}
//若为新的顶点则将其增加到两个映射和一个链表中
//若为新的顶点则将其增加到两个映射和一个链表中
else{
temp_cu.ids[i] = temp_cp.id;//新的顶点索引等于顶点集的数量
model_vert_.push_back(temp_cp);//将新产生的顶点保存到顶点链表中
map_str_point[temp_id_str] = temp_cp;//将新产生的顶点保存到顶点位置映射中
temp_cu.ids[i] = temp_cp.id;//新的顶点索引等于顶点集的数量
model_vert_.push_back(temp_cp);//将新产生的顶点保存到顶点链表中
map_str_point[temp_id_str] = temp_cp;//将新产生的顶点保存到顶点位置映射中
}
}
model_cube_.push_back(temp_cu);
@ -82,7 +82,7 @@ int GM3D::BuildRegularGrid(char* space_para){
else{
vert_num_ = model_vert_.size();
model_num_ = model_cube_.size();
model_block_val_.resize(model_num_,BDL_MAX); //初始化模型块体值为BDL_MAX
model_block_val_.resize(model_num_,BDL_MAX); //初始化模型块体值为BDL_MAX
}
return 0;
}

8
src/forward_delta_t.cpp
View File

@ -10,18 +10,18 @@ int GM3D::ForwardDeltaT(char* noise_level,char* mag_para){
double R222,R122,R212,R112,R221,R121,R211,R111;
double G222,G122,G212,G112,G221,G121,G211,G111;
//初始化正演值和不确定度
//初始化正演值和不确定度
for (int i = 0; i < obs_num_; i++){
obs_p_[i].val = obs_p_[i].dev = 0.0;
}
//确定噪声水平
//确定噪声水平
double noise_mean, noise_dev, temp_noise;
if (2 != sscanf(noise_level,"%lf/%lf",&noise_mean,&noise_dev)){
noise_mean = noise_dev = 0.0;
}
//确定磁化参数
//确定磁化参数
if (4 != sscanf(mag_para,"%lf/%lf/%lf/%lf",&I0,&A0,&I,&A)){
I0 = I = 90; A0 = A = 0;
}
@ -39,7 +39,7 @@ int GM3D::ForwardDeltaT(char* noise_level,char* mag_para){
k6=-sin(I0)*sin(I);
}
//添加高斯噪声值
//添加高斯噪声值
default_random_engine generator;
normal_distribution<double> dist(noise_mean, noise_dev);

8
src/forward_delta_tx.cpp
View File

@ -10,18 +10,18 @@ int GM3D::ForwardDeltaTx(char* noise_level,char* mag_para){
double R222,R122,R212,R112,R221,R121,R211,R111;
double G222,G122,G212,G112,G221,G121,G211,G111;
//初始化正演值和不确定度
//初始化正演值和不确定度
for (int i = 0; i < obs_num_; i++){
obs_p_[i].val = obs_p_[i].dev = 0.0;
}
//确定噪声水平
//确定噪声水平
double noise_mean, noise_dev, temp_noise;
if (2 != sscanf(noise_level,"%lf/%lf",&noise_mean,&noise_dev)){
noise_mean = noise_dev = 0.0;
}
//确定磁化参数
//确定磁化参数
if (4 != sscanf(mag_para,"%lf/%lf/%lf/%lf",&I0,&A0,&I,&A)){
I0 = I = 90; A0 = A = 0;
}
@ -39,7 +39,7 @@ int GM3D::ForwardDeltaTx(char* noise_level,char* mag_para){
k6=-sin(I0)*sin(I);
}
//添加高斯噪声值
//添加高斯噪声值
default_random_engine generator;
normal_distribution<double> dist(noise_mean, noise_dev);

8
src/forward_delta_ty.cpp
View File

@ -10,18 +10,18 @@ int GM3D::ForwardDeltaTy(char* noise_level,char* mag_para){
double R222,R122,R212,R112,R221,R121,R211,R111;
double G222,G122,G212,G112,G221,G121,G211,G111;
//初始化正演值和不确定度
//初始化正演值和不确定度
for (int i = 0; i < obs_num_; i++){
obs_p_[i].val = obs_p_[i].dev = 0.0;
}
//确定噪声水平
//确定噪声水平
double noise_mean, noise_dev, temp_noise;
if (2 != sscanf(noise_level,"%lf/%lf",&noise_mean,&noise_dev)){
noise_mean = noise_dev = 0.0;
}
//确定磁化参数
//确定磁化参数
if (4 != sscanf(mag_para,"%lf/%lf/%lf/%lf",&I0,&A0,&I,&A)){
I0 = I = 90; A0 = A = 0;
}
@ -39,7 +39,7 @@ int GM3D::ForwardDeltaTy(char* noise_level,char* mag_para){
k6=-sin(I0)*sin(I);
}
//添加高斯噪声值
//添加高斯噪声值
default_random_engine generator;
normal_distribution<double> dist(noise_mean, noise_dev);

8
src/forward_delta_tz.cpp
View File

@ -11,18 +11,18 @@ int GM3D::ForwardDeltaTz(char* noise_level,char* mag_para){
double R222,R122,R212,R112,R221,R121,R211,R111;
double G222,G122,G212,G112,G221,G121,G211,G111;
//初始化正演值和不确定度
//初始化正演值和不确定度
for (int i = 0; i < obs_num_; i++){
obs_p_[i].val = obs_p_[i].dev = 0.0;
}
//确定噪声水平
//确定噪声水平
double noise_mean, noise_dev, temp_noise;
if (2 != sscanf(noise_level,"%lf/%lf",&noise_mean,&noise_dev)){
noise_mean = noise_dev = 0.0;
}
//确定磁化参数
//确定磁化参数
if (4 != sscanf(mag_para,"%lf/%lf/%lf/%lf",&I0,&A0,&I,&A)){
I0 = I = 90; A0 = A = 0;
}
@ -40,7 +40,7 @@ int GM3D::ForwardDeltaTz(char* noise_level,char* mag_para){
k6=-sin(I0)*sin(I);
}
//添加高斯噪声值
//添加高斯噪声值
default_random_engine generator;
normal_distribution<double> dist(noise_mean, noise_dev);

8
src/forward_hax.cpp
View File

@ -10,18 +10,18 @@ int GM3D::ForwardHax(char* noise_level,char* mag_para){
double R222,R122,R212,R112,R221,R121,R211,R111;
double G222,G122,G212,G112,G221,G121,G211,G111;
//初始化正演值和不确定度
//初始化正演值和不确定度
for (int i = 0; i < obs_num_; i++){
obs_p_[i].val = obs_p_[i].dev = 0.0;
}
//确定噪声水平
//确定噪声水平
double noise_mean, noise_dev, temp_noise;
if (2 != sscanf(noise_level,"%lf/%lf",&noise_mean,&noise_dev)){
noise_mean = noise_dev = 0.0;
}
//确定磁化参数
//确定磁化参数
if (4 != sscanf(mag_para,"%lf/%lf/%lf/%lf",&I0,&A0,&I,&A)){
I = 90; A = 0;
}
@ -34,7 +34,7 @@ int GM3D::ForwardHax(char* noise_level,char* mag_para){
Gamma=sin(I);
}
//添加高斯噪声值
//添加高斯噪声值
default_random_engine generator;
normal_distribution<double> dist(noise_mean, noise_dev);

8
src/forward_hay.cpp
View File

@ -10,18 +10,18 @@ int GM3D::ForwardHay(char* noise_level,char* mag_para){
double R222,R122,R212,R112,R221,R121,R211,R111;
double G222,G122,G212,G112,G221,G121,G211,G111;
//初始化正演值和不确定度
//初始化正演值和不确定度
for (int i = 0; i < obs_num_; i++){
obs_p_[i].val = obs_p_[i].dev = 0.0;
}
//确定噪声水平
//确定噪声水平
double noise_mean, noise_dev, temp_noise;
if (2 != sscanf(noise_level,"%lf/%lf",&noise_mean,&noise_dev)){
noise_mean = noise_dev = 0.0;
}
//确定磁化参数
//确定磁化参数
if (4 != sscanf(mag_para,"%lf/%lf/%lf/%lf",&I0,&A0,&I,&A)){
I = 90; A = 0;
}
@ -34,7 +34,7 @@ int GM3D::ForwardHay(char* noise_level,char* mag_para){
Gamma=sin(I);
}
//添加高斯噪声值
//添加高斯噪声值
default_random_engine generator;
normal_distribution<double> dist(noise_mean, noise_dev);

6
src/forward_vz.cpp
View File

@ -6,18 +6,18 @@ int GM3D::ForwardVz(char* noise_level){
double R222,R122,R212,R112,R221,R121,R211,R111;
double G222,G122,G212,G112,G221,G121,G211,G111;
//初始化正演值和不确定度
//初始化正演值和不确定度
for (int i = 0; i < obs_num_; i++){
obs_p_[i].val = obs_p_[i].dev = 0.0;
}
//确定噪声水平
//确定噪声水平
double noise_mean, noise_dev, temp_noise;
if (2 != sscanf(noise_level,"%lf/%lf",&noise_mean,&noise_dev)){
noise_mean = noise_dev = 0.0;
}
//添加高斯噪声值
//添加高斯噪声值
default_random_engine generator;
normal_distribution<double> dist(noise_mean, noise_dev);

4
src/forward_vzx.cpp
View File

@ -6,7 +6,7 @@ int GM3D::ForwardVzx(char* noise_level){
double R222,R122,R212,R112,R221,R121,R211,R111;
double G222,G122,G212,G112,G221,G121,G211,G111;
//初始化正演值和不确定度
//初始化正演值和不确定度
for (int i = 0; i < obs_num_; i++){
obs_p_[i].val = obs_p_[i].dev = 0.0;
}
@ -16,7 +16,7 @@ int GM3D::ForwardVzx(char* noise_level){
noise_mean = noise_dev = 0.0;
}
//添加高斯噪声值
//添加高斯噪声值
default_random_engine generator;
normal_distribution<double> dist(noise_mean, noise_dev);

4
src/forward_vzy.cpp
View File

@ -6,7 +6,7 @@ int GM3D::ForwardVzy(char* noise_level){
double R222,R122,R212,R112,R221,R121,R211,R111;
double G222,G122,G212,G112,G221,G121,G211,G111;
//初始化正演值和不确定度
//初始化正演值和不确定度
for (int i = 0; i < obs_num_; i++){
obs_p_[i].val = obs_p_[i].dev = 0.0;
}
@ -16,7 +16,7 @@ int GM3D::ForwardVzy(char* noise_level){
noise_mean = noise_dev = 0.0;
}
//添加高斯噪声值
//添加高斯噪声值
default_random_engine generator;
normal_distribution<double> dist(noise_mean, noise_dev);

4
src/forward_vzz.cpp
View File

@ -6,7 +6,7 @@ int GM3D::ForwardVzz(char* noise_level){
double R222,R122,R212,R112,R221,R121,R211,R111;
double G222,G122,G212,G112,G221,G121,G211,G111;
//初始化正演值和不确定度
//初始化正演值和不确定度
for (int i = 0; i < obs_num_; i++){
obs_p_[i].val = obs_p_[i].dev = 0.0;
}
@ -16,7 +16,7 @@ int GM3D::ForwardVzz(char* noise_level){
noise_mean = noise_dev = 0.0;
}
//添加高斯噪声值
//添加高斯噪声值
default_random_engine generator;
normal_distribution<double> dist(noise_mean, noise_dev);

8
src/forward_za.cpp
View File

@ -10,18 +10,18 @@ int GM3D::ForwardZa(char* noise_level,char* mag_para){
double R222,R122,R212,R112,R221,R121,R211,R111;
double G222,G122,G212,G112,G221,G121,G211,G111;
//初始化正演值和不确定度
//初始化正演值和不确定度
for (int i = 0; i < obs_num_; i++){
obs_p_[i].val = obs_p_[i].dev = 0.0;
}
//确定噪声水平
//确定噪声水平
double noise_mean, noise_dev, temp_noise;
if (2 != sscanf(noise_level,"%lf/%lf",&noise_mean,&noise_dev)){
noise_mean = noise_dev = 0.0;
}
//确定磁化参数
//确定磁化参数
if (4 != sscanf(mag_para,"%lf/%lf/%lf/%lf",&I0,&A0,&I,&A)){
I = 90; A = 0;
}
@ -34,7 +34,7 @@ int GM3D::ForwardZa(char* noise_level,char* mag_para){
Gamma=sin(I);
}
//添加高斯噪声值
//添加高斯噪声值
default_random_engine generator;
normal_distribution<double> dist(noise_mean, noise_dev);

41
src/gm3d.h
View File

@ -7,25 +7,24 @@ class GM3D{
public:
GM3D(){}
~GM3D(){}
int BuildRegularGrid(char*); //初始化反演模型空间
int AddModels(char*); //读取模型块体参数文件
void AddModels_GUI(); //添加模型块体参数 GUI程序调用
int AddRegularBlock(modelist); //添加普通模型块体
int AddTiltedBlock(modelist); //添加倾斜模型块体
int AddSphereBlock(modelist); //添加球体椭球体块体
int AddInterfaceBlock(modelist); //添加密度界面
//模型操作
int BuildRegularGrid(char*); //初始化反演模型空间
int AddModels(char*); //读取模型块体参数文件
void AddModels_GUI(); //添加模型块体参数
int AddRegularBlock(modelist); //添加普通模型块体
int AddTiltedBlock(modelist); //添加倾斜模型块体
int AddSphereBlock(modelist); //添加球体椭球体块体
int AddInterfaceBlock(modelist); //添加密度界面
//模型操作
int ReadModel(char*,char*);
int ReadMeshtoolModel(char*); //读入meshtool 3D格式文件
//输出模型
int RegisteredOuput(bool); //注册输出的块体模型
int OutMshFile(char*,string); //输出模型文件
int OutMeshtoolFile(char*); //输出meshtool 3D格式文件
int OutNeighborFile(char*,char*); //输出模型块体或顶点的相邻关系 暂缓
//观测数据
int ReadModel_mst(char*);
//输出模型
int RegisteredOuput(bool); //注册输出的块体模型
int OutMshFile(char*,string); //输出模型文件
int OutNeighborFile(char*,char*); //输出模型块体或顶点的相邻关系 暂缓
//观测数据
int InitObs(char*);
int OutObs(char*);
//正演函数
//正演函数
int ForwardVz(char*);
int ForwardVzx(char*);
int ForwardVzy(char*);
@ -37,16 +36,16 @@ public:
int ForwardHax(char*,char*);
int ForwardHay(char*,char*);
int ForwardZa(char*,char*);
//gm3d_gui函数
//gm3d_gui函数
void get_model_list(modelistArray);
private:
int obs_num_, model_num_, vert_num_;
//正演数组
//正演数组
obspointArray obs_p_;
_2dArray input_models_;
_1sArray input_model_names_;
_1dArray forward_model_;
//模型数据
//模型数据
cubeArray model_cube_;
cpointArray model_vert_;
_1dArray model_block_val_;
@ -57,7 +56,7 @@ private:
_int2intMap vert_out_map_;
_int2intMap ele_data_out_map_;
_2iArray model_vert_neighbor_; //暂缓
_2iArray model_cube_neighbor_; //暂缓
_2iArray model_vert_neighbor_; //暂缓
_2iArray model_cube_neighbor_; //暂缓
};
#endif

468
src/gm3d_gui.cpp
View File

@ -10,7 +10,8 @@ Fl_Button *mod_para_file_btn=(Fl_Button *)0;
Fl_Input *mod_ele_input_build=(Fl_Input *)0;
Fl_Button *build_mod_btn=(Fl_Button *)0;
Fl_Output *mesh_para_output=(Fl_Output *)0;
Fl_Button *edit_mod_btn=(Fl_Button *)0;
Fl_Button *add_mod_btn=(Fl_Button *)0;
Fl_Button *del_mod_btn=(Fl_Button *)0;
Fl_Check_Button *rm_emp_bok_check=(Fl_Check_Button *)0;
Fl_Browser *mod_para_brw=(Fl_Browser *)0;
Fl_Button *mod_file_out_btn=(Fl_Button *)0;
@ -48,19 +49,10 @@ Fl_Output *res_file_output=(Fl_Output *)0;
char mesh_filename[1024];
char out_msh_filename[1024];
// 声明三个字符串用于保存新建的模型类型和赋值类型
char add_mod_type[1024];
char add_val_type[1024];
char agn_val_type[1024];
// 声明一个字符串用于保存新建的模型参数并添加至列表中
// 模型参数的排列为<模型类型> <赋值类型> <物理参数字符串> <几何参数字符串>
char add_mod_para[1024];
void cb_mesh_para_input(Fl_Input*, void*){
strcpy(mesh_filename,mesh_para_input->value());
mesh_para_output->value(mesh_para_input->value());
return;
strcpy(mesh_filename,mesh_para_input->value());
mesh_para_output->value(mesh_para_input->value());
return;
}
void cb_mesh_file_btn(Fl_Widget*, void*){
@ -72,7 +64,7 @@ void cb_mesh_file_btn(Fl_Widget*, void*){
// Block until user picks something.
while(chooser.shown()) { Fl::wait(); }
//从文件的绝对路径中分离文件名并从输出框输出
//从文件的绝对路径中分离文件名并从输出框输出
if (chooser.value() != NULL)
{
strcpy(mesh_filename,chooser.value());
@ -82,9 +74,8 @@ void cb_mesh_file_btn(Fl_Widget*, void*){
std::stringstream temp_ss;
temp_ss << abs_filename.substr(pos+1);
//const char* disp_filename = temp_ss.str().c_str();
//mesh_para_output->value(disp_filename);
mesh_para_output->value(temp_ss.str().c_str());
const char* disp_filename = temp_ss.str().c_str();
mesh_para_output->value(disp_filename);
}
else
{
@ -101,7 +92,7 @@ void cb_mod_para_file_btn(Fl_Widget*, void*){
// Block until user picks something.
while(chooser.shown()) { Fl::wait(); }
//下面我们读入文件并添加模型参数到列表显示
//下面我们读入文件并添加模型参数到列表显示
std::ifstream modpara_in;
modpara_in.open(chooser.value());
if (!modpara_in)
@ -126,132 +117,17 @@ void cb_mod_para_file_btn(Fl_Widget*, void*){
return;
}
void cb_del_mod_btn(Fl_Widget*, void*){
int index = mod_para_brw->value();
mod_para_brw->remove(index);
del_mod_btn->deactivate();
}
void cb_mod_para_brw(Fl_Browser*, void*){
std::stringstream temp_ss;
modelist temp_list;
int index;
for ( int t=1; t<=mod_para_brw->size(); t++ ) {
if ( mod_para_brw->selected(t) ) {
//printf("%d) '%s'\n", t, mod_para_brw->text(t));
index = mod_para_brw->value();
temp_ss.clear(); temp_ss.str("");
temp_ss << mod_para_brw->text(index);
temp_ss >> temp_list.mod_type >> temp_list.val_type
>> temp_list.mod_value >> temp_list.mod_para;
// 重置模型类型按钮
sscanf(temp_list.mod_type, "%s", add_mod_type);
if (!strcmp(add_mod_type, "regular_block"))
{
reg_bok_rbtn->value(1);
til_bok_rbtn->value(0);
sph_rbtn->value(0);
int_face_rbtn->value(0);
top_val_btn->value(0);
bot_val_btn->value(0);
agn_part_group->deactivate();
}
else if (!strcmp(add_mod_type, "tilted_block"))
{
reg_bok_rbtn->value(0);
til_bok_rbtn->value(1);
sph_rbtn->value(0);
int_face_rbtn->value(0);
top_val_btn->value(0);
bot_val_btn->value(0);
agn_part_group->deactivate();
}
else if (!strcmp(add_mod_type, "sphere"))
{
reg_bok_rbtn->value(0);
til_bok_rbtn->value(0);
sph_rbtn->value(1);
int_face_rbtn->value(0);
top_val_btn->value(0);
bot_val_btn->value(0);
agn_part_group->deactivate();
}
else if (!strcmp(add_mod_type, "interface"))
{
reg_bok_rbtn->value(0);
til_bok_rbtn->value(0);
sph_rbtn->value(0);
int_face_rbtn->value(1);
top_val_btn->value(0);
bot_val_btn->value(0);
agn_part_group->activate();
}
if (2 == sscanf(temp_list.val_type, "%[^/]/%s", add_val_type,
agn_val_type))
{
// 重置模型数据类型按钮
if (!strcmp(add_val_type, "add"))
{
app_val_rbtn->value(1);
rep_val_rbtn->value(0);
era_val_rbtn->value(0);
}
else if (!strcmp(add_val_type, "replace"))
{
app_val_rbtn->value(0);
rep_val_rbtn->value(1);
era_val_rbtn->value(0);
}
else if (!strcmp(add_val_type, "erase"))
{
app_val_rbtn->value(0);
rep_val_rbtn->value(0);
era_val_rbtn->value(1);
}
if (!strcmp(agn_val_type, "top"))
{
top_val_btn->value(1);
bot_val_btn->value(0);
}
else if (!strcmp(agn_val_type, "bot"))
{
top_val_btn->value(0);
bot_val_btn->value(1);
}
}
else
{
sscanf(temp_list.val_type, "%s", add_val_type);
// 重置模型数据类型按钮
if (!strcmp(add_val_type, "add"))
{
app_val_rbtn->value(1);
rep_val_rbtn->value(0);
era_val_rbtn->value(0);
}
else if (!strcmp(add_val_type, "replace"))
{
app_val_rbtn->value(0);
rep_val_rbtn->value(1);
era_val_rbtn->value(0);
}
else if (!strcmp(add_val_type, "erase"))
{
app_val_rbtn->value(0);
rep_val_rbtn->value(0);
era_val_rbtn->value(1);
}
}
char tmp_char[1024];
temp_ss.clear(); temp_ss.str("");
temp_ss << temp_list.mod_value;
temp_ss >> tmp_char;
mod_val_input->value(tmp_char);
sig_mod_para_input->value(temp_list.mod_para);
sig_delete_btn->activate();
sig_replace_btn->activate();
del_mod_btn->activate();
}
}
return;
@ -273,7 +149,7 @@ void cb_mod_file_out_btn(Fl_Widget*, void*){
// Block until user picks something.
while(chooser.shown()) { Fl::wait(); }
//从文件的绝对路径中分离文件名并从输出框输出
//从文件的绝对路径中分离文件名并从输出框输出
if (chooser.value() != NULL)
{
strcpy(out_msh_filename,chooser.value());
@ -283,9 +159,8 @@ void cb_mod_file_out_btn(Fl_Widget*, void*){
std::stringstream temp_ss;
temp_ss << abs_filename.substr(pos+1);
//const char* disp_filename = temp_ss.str().c_str();
//mod_out_file_output->value(disp_filename);
mod_out_file_output->value(temp_ss.str().c_str());
const char* disp_filename = temp_ss.str().c_str();
mod_out_file_output->value(disp_filename);
}
else
{
@ -295,17 +170,17 @@ void cb_mod_file_out_btn(Fl_Widget*, void*){
}
void cb_build_mod_btn(Fl_Widget*, void*){
//首先构建一个GM3D实例
//首先构建一个GM3D实例
GM3D gm3d_instance;
//设置初始化变量
//设置初始化变量
char dimension[1024] = "10/20/990/10/20/990/10/20/490";
char out_mshname[1024] = "Untitled.msh";
char elename[1024] = "Untitled";
//将mesh_para_output拷贝到dimension
//将mesh_para_output拷贝到dimension
strcpy(dimension, mesh_filename);
//std::cout << dimension << std::endl;
//将mod_para_brw中的模型参数列表拷贝到一个临时参数列表
//将mod_para_brw中的模型参数列表拷贝到一个临时参数列表
std::stringstream temp_ss;
modelist temp_list;
modelistArray brw_model_list;
@ -321,21 +196,21 @@ void cb_build_mod_btn(Fl_Widget*, void*){
}
}
//将brw_model_list拷贝到GM3D中
//将brw_model_list拷贝到GM3D中
gm3d_instance.get_model_list(brw_model_list);
//将mod_out_file_input拷贝到out_mshname
//将mod_out_file_input拷贝到out_mshname
strcpy(out_mshname, out_msh_filename);
if (!strcmp(out_mshname,"")){
fl_message("Output file's name can't be empty !");
return;
}
//将mod_ele_input_build拷贝到elename
//将mod_ele_input_build拷贝到elename
strcpy(elename, mod_ele_input_build->value());
if (!strcmp(elename,"")){
fl_message("Element data's name can't be empty !");
return;
}
//构建模型网络
//构建模型网络
if(gm3d_instance.BuildRegularGrid(dimension)){
fl_message("Mesh Parameters Load Error !");
return;
@ -344,32 +219,32 @@ void cb_build_mod_btn(Fl_Widget*, void*){
gm3d_instance.RegisteredOuput(remove_null);
if (gm3d_instance.OutMshFile(out_mshname,elename)) return;
//这里需要一个输出信息窗口
//这里需要一个输出信息窗口
fl_message("Model Construction Completed !");
return;
}
/***********************************以下是正演计算GUI函数******************************/
/***********************************以下是正演计算GUI函数******************************/
char in_msh_filename[1024];
char in_obs_filename[1024];
char out_res_filename[1024];
void cb_mod_file_input(Fl_Input*, void*){
strcpy(in_msh_filename,mod_file_input->value());
mod_file_output->value(mod_file_input->value());
return;
strcpy(in_msh_filename,mod_file_input->value());
mod_file_output->value(mod_file_input->value());
return;
}
void cb_mod_file_btn(Fl_Widget*, void*){
// Create the file chooser, and show it
// directory, filter, chooser type, title
Fl_File_Chooser chooser(".", "*.msh", Fl_File_Chooser::SINGLE, "Create a model file.");
Fl_File_Chooser chooser(".", "*.msh, *.mst", Fl_File_Chooser::SINGLE, "Create a model file.");
chooser.show();
// Block until user picks something.
while(chooser.shown()) { Fl::wait(); }
//从文件的绝对路径中分离文件名并从输出框输出
//从文件的绝对路径中分离文件名并从输出框输出
if (chooser.value() != NULL)
{
strcpy(in_msh_filename,chooser.value());
@ -379,9 +254,8 @@ void cb_mod_file_btn(Fl_Widget*, void*){
std::stringstream temp_ss;
temp_ss << abs_filename.substr(pos+1);
//const char* disp_filename = temp_ss.str().c_str();
//mod_file_output->value(disp_filename);
mod_file_output->value(temp_ss.str().c_str());
const char* disp_filename = temp_ss.str().c_str();
mod_file_output->value(disp_filename);
}
else
{
@ -391,9 +265,9 @@ void cb_mod_file_btn(Fl_Widget*, void*){
}
void cb_obs_file_input(Fl_Input*, void*){
strcpy(in_obs_filename,obs_file_input->value());
obs_file_output->value(obs_file_input->value());
return;
strcpy(in_obs_filename,obs_file_input->value());
obs_file_output->value(obs_file_input->value());
return;
}
void cb_obs_file_btn(Fl_Widget*, void*){
@ -405,7 +279,7 @@ void cb_obs_file_btn(Fl_Widget*, void*){
// Block until user picks something.
while(chooser.shown()) { Fl::wait(); }
//从文件的绝对路径中分离文件名并从输出框输出
//从文件的绝对路径中分离文件名并从输出框输出
if (chooser.value() != NULL)
{
strcpy(in_obs_filename,chooser.value());
@ -415,9 +289,8 @@ void cb_obs_file_btn(Fl_Widget*, void*){
std::stringstream temp_ss;
temp_ss << abs_filename.substr(pos+1);
//const char* disp_filename = temp_ss.str().c_str();
//obs_file_output->value(disp_filename);
obs_file_output->value(temp_ss.str().c_str());
const char* disp_filename = temp_ss.str().c_str();
obs_file_output->value(disp_filename);
}
else
{
@ -427,9 +300,9 @@ void cb_obs_file_btn(Fl_Widget*, void*){
}
void cb_res_file_input(Fl_Input*, void*){
strcpy(out_res_filename,res_file_input->value());
res_file_output->value(res_file_input->value());
return;
strcpy(out_res_filename,res_file_input->value());
res_file_output->value(res_file_input->value());
return;
}
void cb_res_file_btn(Fl_Widget*, void*){
@ -441,7 +314,7 @@ void cb_res_file_btn(Fl_Widget*, void*){
// Block until user picks something.
while(chooser.shown()) { Fl::wait(); }
//从文件的绝对路径中分离文件名并从输出框输出
//从文件的绝对路径中分离文件名并从输出框输出
if (chooser.value() != NULL)
{
strcpy(out_res_filename,chooser.value());
@ -451,9 +324,8 @@ void cb_res_file_btn(Fl_Widget*, void*){
std::stringstream temp_ss;
temp_ss << abs_filename.substr(pos+1);
//const char* disp_filename = temp_ss.str().c_str();
//res_file_output->value(disp_filename);
res_file_output->value(temp_ss.str().c_str());
const char* disp_filename = temp_ss.str().c_str();
res_file_output->value(disp_filename);
}
else
{
@ -497,9 +369,9 @@ void cb_mag_data_check(Fl_Check_Button*, void*){
}
void cb_cal_btn(Fl_Button*, void*){
//首先构建一个GM3D实例
//首先构建一个GM3D实例
GM3D gm3d_instance;
//设置初始化变量
//设置初始化变量
char in_mshname[1024] = "Untitled.msh";
char in_obspara[1024] = "Untitled.txt";
char res_outfile[1024] = "Untitled";
@ -514,11 +386,18 @@ void cb_cal_btn(Fl_Button*, void*){
strcpy(ele_name,mod_ele_input->value());
strcpy(noise_para,noise_para_input->value());
strcpy(mag_para,mag_para_input->value());
if (gm3d_instance.ReadModel(in_mshname,ele_name))
{
fl_message("Model Load Error !");
return;
if(in_mshname[strlen(in_mshname)-1] == 'h'){
if (gm3d_instance.ReadModel(in_mshname,ele_name))
{
fl_message("Model Load Error !");
return;
}
}else if(in_mshname[strlen(in_mshname)-1] == 't'){
if (gm3d_instance.ReadModel_mst(in_mshname))
{
fl_message("Model Load Error !");
return;
}
}
if (gm3d_instance.InitObs(in_obspara))
@ -661,227 +540,8 @@ void cb_cal_btn(Fl_Button*, void*){
}
}
}
fl_message("Forward calculation completed !");
return;
}
Fl_Double_Window *edit_mod_win=(Fl_Double_Window *)0;
Fl_Group *mod_type_group=(Fl_Group *)0;
Fl_Radio_Round_Button *reg_bok_rbtn=(Fl_Radio_Round_Button *)0;
Fl_Radio_Round_Button *til_bok_rbtn=(Fl_Radio_Round_Button *)0;
Fl_Radio_Round_Button *sph_rbtn=(Fl_Radio_Round_Button *)0;
Fl_Radio_Round_Button *int_face_rbtn=(Fl_Radio_Round_Button *)0;
Fl_Group *val_type_group=(Fl_Group *)0;
Fl_Radio_Round_Button *app_val_rbtn=(Fl_Radio_Round_Button *)0;
Fl_Radio_Round_Button *rep_val_rbtn=(Fl_Radio_Round_Button *)0;
Fl_Radio_Round_Button *era_val_rbtn=(Fl_Radio_Round_Button *)0;
Fl_Group *agn_part_group=(Fl_Group *)0;
Fl_Radio_Round_Button *top_val_btn=(Fl_Radio_Round_Button *)0;
Fl_Radio_Round_Button *bot_val_btn=(Fl_Radio_Round_Button *)0;
Fl_Input *sig_mod_para_input=(Fl_Input *)0;
Fl_Input *mod_val_input=(Fl_Input *)0;
Fl_Button *sig_add_btn=(Fl_Button *)0;
Fl_Button *sig_replace_btn=(Fl_Button *)0;
Fl_Button *sig_delete_btn=(Fl_Button *)0;
Fl_Return_Button *can_add_btn=(Fl_Return_Button *)0;
void cb_reg_bok_rbtn(Fl_Button*, void*){
strcpy(add_mod_type,"regular_block");
top_val_btn->value(0);
bot_val_btn->value(0);
agn_part_group->deactivate();
//printf("%s\n", "Regular block button selected!");
return;
}
void cb_til_bok_rbtn(Fl_Button*, void*){
strcpy(add_mod_type,"tilted_block");
top_val_btn->value(0);
bot_val_btn->value(0);
agn_part_group->deactivate();
//printf("%s\n", "Tilted block button selected!");
return;
}
void cb_sph_rbtn(Fl_Button*, void*){
strcpy(add_mod_type,"sphere");
top_val_btn->value(0);
bot_val_btn->value(0);
agn_part_group->deactivate();
//printf("%s\n", "Sphere button selected!");
return;
}
void cb_int_face_rbtn(Fl_Button*, void*){
strcpy(add_mod_type,"interface");
agn_part_group->activate();
//printf("%s\n", "Interface button selected!");
return;
}
void cb_app_val_rbtn(Fl_Button*, void*){
strcpy(add_val_type,"add");
//printf("%s\n", "append button selected!");
return;
}
void cb_rep_val_rbtn(Fl_Button*, void*){
strcpy(add_val_type,"replace");
//printf("%s\n", "replace button selected!");
return;
}
void cb_era_val_rbtn(Fl_Button*, void*){
strcpy(add_val_type,"erase");
//printf("%s\n", "Erase button selected!");
return;
}
void cb_top_val_rbtn(Fl_Button*, void*){
strcpy(agn_val_type,"top");
//printf("%s\n", "Erase button selected!");
return;
}
void cb_bot_val_rbtn(Fl_Button*, void*){
strcpy(agn_val_type,"bot");
//printf("%s\n", "Erase button selected!");
return;
}
void cb_sig_add_btn(Fl_Widget*, void*){
strcpy(add_mod_para,add_mod_type);
strcat(add_mod_para," ");
strcat(add_mod_para,add_val_type);
if (!strcmp(add_mod_type,"interface"))
{
strcat(add_mod_para,"/");
strcat(add_mod_para,agn_val_type);
}
strcat(add_mod_para," ");
strcat(add_mod_para,mod_val_input->value());
strcat(add_mod_para," ");
strcat(add_mod_para,sig_mod_para_input->value());
mod_para_brw->add(add_mod_para);
return;
}
void cb_sig_replace_btn(Fl_Widget*, void*)
{
int index = mod_para_brw->value();
strcpy(add_mod_para,add_mod_type);
strcat(add_mod_para," ");
strcat(add_mod_para,add_val_type);
if (!strcmp(add_mod_type,"interface"))
{
strcat(add_mod_para,"/");
strcat(add_mod_para,agn_val_type);
}
strcat(add_mod_para," ");
strcat(add_mod_para,mod_val_input->value());
strcat(add_mod_para," ");
strcat(add_mod_para,sig_mod_para_input->value());
mod_para_brw->replace(index, add_mod_para);
sig_delete_btn->deactivate();
sig_replace_btn->deactivate();
return;
}
void cb_sig_del_btn(Fl_Widget*, void*)
{
int index = mod_para_brw->value();
mod_para_brw->remove(index);
sig_delete_btn->deactivate();
sig_replace_btn->deactivate();
return;
}
//取消添加模型
void cb_can_add_btn(Fl_Widget*, void*){
edit_mod_win->hide();
return;
}
void cb_edit_mod_btn(Fl_Widget*, void*)
{
{ edit_mod_win = new Fl_Double_Window(314, 320, "Model editor (gm3d)");
{ mod_type_group = new Fl_Group(20, 25, 250, 53, "Model Type :");
mod_type_group->align(Fl_Align(FL_ALIGN_TOP_LEFT));
{ reg_bok_rbtn = new Fl_Radio_Round_Button(20, 25, 110, 28, "Regular Block");
reg_bok_rbtn->down_box(FL_ROUND_DOWN_BOX);
reg_bok_rbtn->callback((Fl_Callback*)cb_reg_bok_rbtn);
reg_bok_rbtn->tooltip("model parameter:\n<xmin>/<xmax>/<ymin>/<ymax>/<zmin>/<zmax>");
} // Fl_Radio_Round_Button* reg_bok_rbtn
{ til_bok_rbtn = new Fl_Radio_Round_Button(160, 25, 100, 28, "Tilted Block");
til_bok_rbtn->down_box(FL_ROUND_DOWN_BOX);
til_bok_rbtn->callback((Fl_Callback*)cb_til_bok_rbtn);
til_bok_rbtn->tooltip("model parameter:\n<xmin_z>/<xmax_z>/<ymin_z>/<ymax_z>/<zmin>/<xmin_Z>/<xmax_Z>/<ymin_Z>/<ymax_Z>/<zmax>");
} // Fl_Radio_Round_Button* til_bok_rbtn
{ sph_rbtn = new Fl_Radio_Round_Button(20, 50, 70, 28, "Sphere");
sph_rbtn->down_box(FL_ROUND_DOWN_BOX);
sph_rbtn->callback((Fl_Callback*)cb_sph_rbtn);
sph_rbtn->tooltip("model parameter:\n<x_c>/<y_c>/<z_c>/<x_radius>/<y_radius>/<z_radius>");
} // Fl_Radio_Round_Button* sph_rbtn
{ int_face_rbtn = new Fl_Radio_Round_Button(160, 50, 80, 28, "Interface");
int_face_rbtn->down_box(FL_ROUND_DOWN_BOX);
int_face_rbtn->callback((Fl_Callback*)cb_int_face_rbtn);
int_face_rbtn->tooltip("model parameter:\n<grid-filename>");
} // Fl_Radio_Round_Button* int_face_rbtn
mod_type_group->end();
} // Fl_Group* mod_type_group
{ val_type_group = new Fl_Group(20, 100, 250, 30, "Value Type :");
val_type_group->align(Fl_Align(FL_ALIGN_TOP_LEFT));
{ app_val_rbtn = new Fl_Radio_Round_Button(20, 100, 80, 28, "Append");
app_val_rbtn->down_box(FL_ROUND_DOWN_BOX);
app_val_rbtn->callback((Fl_Callback*)cb_app_val_rbtn);
} // Fl_Radio_Round_Button* app_val_rbtn
{ rep_val_rbtn = new Fl_Radio_Round_Button(110, 100, 80, 28, "Replace");
rep_val_rbtn->down_box(FL_ROUND_DOWN_BOX);
rep_val_rbtn->callback((Fl_Callback*)cb_rep_val_rbtn);
} // Fl_Radio_Round_Button* rep_val_rbtn
{ era_val_rbtn = new Fl_Radio_Round_Button(200, 100, 60, 28, "Erase");
era_val_rbtn->down_box(FL_ROUND_DOWN_BOX);
era_val_rbtn->callback((Fl_Callback*)cb_era_val_rbtn);
} // Fl_Radio_Round_Button* era_val_rbtn
val_type_group->end();
} // Fl_Group* val_type_group
{ agn_part_group = new Fl_Group(20, 150, 160, 28, "Assign Part :");
agn_part_group->align(Fl_Align(FL_ALIGN_TOP_LEFT));
{ top_val_btn = new Fl_Radio_Round_Button(20, 150, 50, 28, "Top");
top_val_btn->down_box(FL_ROUND_DOWN_BOX);
top_val_btn->callback((Fl_Callback*)cb_top_val_rbtn);
} // Fl_Radio_Round_Button* top_val_btn
{ bot_val_btn = new Fl_Radio_Round_Button(90, 150, 70, 28, "Bottom");
bot_val_btn->down_box(FL_ROUND_DOWN_BOX);
bot_val_btn->callback((Fl_Callback*)cb_bot_val_rbtn);
} // Fl_Radio_Round_Button* bot_val_btn
agn_part_group->deactivate();
agn_part_group->end();
} // Fl_Group* agn_part_group
{ sig_mod_para_input = new Fl_Input(20, 200, 275, 28, "Model Parameter :");
sig_mod_para_input->align(Fl_Align(FL_ALIGN_TOP_LEFT));
} // Fl_Input* sig_mod_para_input
{ mod_val_input = new Fl_Input(180, 150, 115, 28, "Model Value :");
mod_val_input->align(Fl_Align(FL_ALIGN_TOP_LEFT));
} // Fl_Input* mod_val_input
{ sig_add_btn = new Fl_Button(19, 240, 80, 28, "Add");
sig_add_btn->callback((Fl_Callback*)cb_sig_add_btn);
} // Fl_Button* sig_add_btn
{ sig_replace_btn = new Fl_Button(118, 240, 80, 28, "Replace");
sig_replace_btn->callback((Fl_Callback*)cb_sig_replace_btn);
sig_replace_btn->deactivate();
} // Fl_Button* sig_replace_btn
{ sig_delete_btn = new Fl_Button(215, 240, 80, 28, "Remove");
sig_delete_btn->callback((Fl_Callback*)cb_sig_del_btn);
sig_delete_btn->deactivate();
} // Fl_Button* sig_delete_btn
{ can_add_btn = new Fl_Return_Button(215, 280, 80, 28, "Cancel");
can_add_btn->callback((Fl_Callback*)cb_can_add_btn);
} // Fl_Return_Button* can_add_btn
edit_mod_win->end();
edit_mod_win->show();
} // Fl_Double_Window* edit_mod_win
}
}

61
src/gm3d_gui.h
View File

@ -10,10 +10,9 @@
#include <FL/Fl_Button.H>
#include <FL/Fl_Output.H>
#include <FL/Fl_Check_Button.H>
#include <FL/Fl_Radio_Round_Button.H>
#include <FL/Fl_Browser.H>
#include <FL/Fl_File_Chooser.H>
#include <FL/Fl_Radio_Round_Button.H>
#include <FL/Fl_Return_Button.H>
#include <iostream>
#include <string.h>
@ -22,24 +21,34 @@
extern void cb_mesh_para_input(Fl_Input*, void*);
extern void cb_mesh_file_btn(Fl_Widget*, void*);
extern void cb_add_mod_btn(Fl_Widget*, void*);
extern void cb_mod_para_file_btn(Fl_Widget*, void*);
extern void cb_build_mod_btn(Fl_Widget*, void*);
extern void cb_edit_mod_btn(Fl_Widget*, void*);
extern void cb_rm_emp_bok_check(Fl_Check_Button*, void*);
extern void cb_del_mod_btn(Fl_Widget*, void*);
extern void cb_mod_para_brw(Fl_Browser*, void*);
extern void cb_rm_emp_bok_check(Fl_Check_Button*, void*);
extern void cb_mod_file_out_btn(Fl_Widget*, void*);
extern void cb_mod_file_btn(Fl_Widget*, void*);
extern void cb_obs_file_btn(Fl_Widget*, void*);
extern void cb_build_mod_btn(Fl_Widget*, void*);
extern void cb_mod_file_input(Fl_Input*, void*);
extern void cb_mod_file_btn(Fl_Widget*, void*);
extern void cb_obs_file_input(Fl_Input*, void*);
extern void cb_obs_file_btn(Fl_Widget*, void*);
extern void cb_res_file_input(Fl_Input*, void*);
extern void cb_res_file_btn(Fl_Widget*, void*);
extern void cb_noise_check(Fl_Check_Button*, void*);
extern void cb_mag_data_check(Fl_Check_Button*, void*);
extern void cb_cal_btn(Fl_Button*, void*);
extern void cb_noise_check(Fl_Check_Button*, void*);
extern void cb_res_file_input(Fl_Input*, void*);
extern void cb_res_file_btn(Fl_Widget*, void*);
extern void cb_reg_bok_rbtn(Fl_Button*, void*);
extern void cb_til_bok_rbtn(Fl_Button*, void*);
extern void cb_sph_rbtn(Fl_Button*, void*);
extern void cb_int_face_rbtn(Fl_Button*, void*);
extern void cb_app_val_rbtn(Fl_Button*, void*);
extern void cb_rep_val_rbtn(Fl_Button*, void*);
extern void cb_era_val_rbtn(Fl_Button*, void*);
extern void cb_top_val_rbtn(Fl_Button*, void*);
extern void cb_bot_val_rbtn(Fl_Button*, void*);
extern void cb_sig_add_btn(Fl_Widget*, void*);
extern void cb_can_add_btn(Fl_Widget*, void*);
extern Fl_Double_Window *main_window;
extern Fl_Tabs *main_tabs;
@ -50,7 +59,8 @@ extern Fl_Button *mod_para_file_btn;
extern Fl_Input *mod_ele_input_build;
extern Fl_Button *build_mod_btn;
extern Fl_Output *mesh_para_output;
extern Fl_Button *edit_mod_btn;
extern Fl_Button *add_mod_btn;
extern Fl_Button *del_mod_btn;
extern Fl_Check_Button *rm_emp_bok_check;
extern Fl_Browser *mod_para_brw;
extern Fl_Button *mod_file_out_btn;
@ -85,21 +95,7 @@ extern Fl_Input *res_file_input;
extern Fl_Button *res_file_btn;
extern Fl_Output *res_file_output;
extern void cb_reg_bok_rbtn(Fl_Button*, void*);
extern void cb_til_bok_rbtn(Fl_Button*, void*);
extern void cb_sph_rbtn(Fl_Button*, void*);
extern void cb_int_face_rbtn(Fl_Button*, void*);
extern void cb_app_val_rbtn(Fl_Button*, void*);
extern void cb_rep_val_rbtn(Fl_Button*, void*);
extern void cb_era_val_rbtn(Fl_Button*, void*);
extern void cb_top_val_rbtn(Fl_Button*, void*);
extern void cb_bot_val_rbtn(Fl_Button*, void*);
extern void cb_sig_add_btn(Fl_Widget*, void*);
extern void cb_sig_replace_btn(Fl_Widget*, void*);
extern void cb_sig_del_btn(Fl_Widget*, void*);
extern void cb_can_add_btn(Fl_Widget*, void*);
extern Fl_Double_Window *edit_mod_win;
extern Fl_Double_Window *add_mod_win;
extern Fl_Group *mod_type_group;
extern Fl_Radio_Round_Button *reg_bok_rbtn;
extern Fl_Radio_Round_Button *til_bok_rbtn;
@ -113,10 +109,7 @@ extern Fl_Group *agn_part_group;
extern Fl_Radio_Round_Button *top_val_btn;
extern Fl_Radio_Round_Button *bot_val_btn;
extern Fl_Input *sig_mod_para_input;
extern Fl_Input *mod_val_input;
extern Fl_Button *sig_add_btn;
extern Fl_Button *sig_replace_btn;
extern Fl_Button *sig_delete_btn;
extern Fl_Return_Button *can_add_btn;
#endif //_GM3D_GUI_H
extern Fl_Button *sig_add_btn;
extern Fl_Input *mod_val_input;
#endif

30
src/head_func.cpp
View File

@ -1,6 +1,6 @@
#include "head_func.h"
/*************************数据结构函数********************************/
//cpoint减法
/*************************数据结构函数********************************/
//cpoint减法
cpoint operator -(cpoint a, cpoint b){
cpoint m;
m.x=a.x-b.x;
@ -8,25 +8,25 @@ cpoint operator -(cpoint a, cpoint b){
m.z=a.z-b.z;
return m;
}
//cpoint模长
//cpoint模长
double modCpoint(cpoint v){
return sqrt(v.x*v.x+v.y*v.y+v.z*v.z);
}
/*************************全局函数********************************/
//正负分离的atan函数 正数返回atan 负数返回atan+pi
/*************************全局函数********************************/
//正负分离的atan函数 正数返回atan 负数返回atan+pi
double arctg(double v){
double ang;
if(v>=0) ang=atan(v);
else if(v<0) ang=atan(v)+Pi;
return ang;
}
//将string转换为stringstream
//将string转换为stringstream
stringstream str2ss(string s){
stringstream sstr;
sstr.str(""); sstr.clear(); sstr.str(s);
return sstr;
}
//返回一个cpoint的位置id字符串
//返回一个cpoint的位置id字符串
string cpoint_id(cpoint c){
string vert_id, mid_id;
stringstream sstemp;
@ -40,7 +40,7 @@ string cpoint_id(cpoint c){
vert_id = vert_id + " " + mid_id;
return vert_id;
}
//测试打开输入文件 如果成功则返回0并输出信息 否则返回1
//测试打开输入文件 如果成功则返回0并输出信息 否则返回1
int open_infile(ifstream &infile,char* filename){
infile.open(filename);
if (!infile){
@ -49,7 +49,7 @@ int open_infile(ifstream &infile,char* filename){
}
return 0;
}
//测试打开输出文件 如果成功则返回0并输出信息 否则返回1
//测试打开输出文件 如果成功则返回0并输出信息 否则返回1
int open_outfile(ofstream &outfile,char* filename){
outfile.open(filename);
if (!outfile){
@ -59,8 +59,8 @@ int open_outfile(ofstream &outfile,char* filename){
return 0;
}
//规则网络插值 长方形内数据插值 距离平方反比
/*长方体示意图*/
//规则网络插值 长方形内数据插值 距离平方反比
/*长方体示意图*/
// y
// |
// |
@ -69,10 +69,10 @@ int open_outfile(ofstream &outfile,char* filename){
// | |
// | |
// 0------------1--->x
// 左下角坐标x0 y0
// 块体尺寸dx dy
// 插值点坐标x y
// 四个角点值
// 左下角坐标x0 y0
// 块体尺寸dx dy
// 插值点坐标x y
// 四个角点值
double grid_interpolate(double x0,double y0,double dx,double dy,double x,double y,
double d0,double d1,double d2,double d3)
{

96
src/head_func.h
View File

@ -18,42 +18,42 @@
using namespace std;
//数学常量
#define BDL_MAX 1e+30 ///< 定义变量的最大值
#define BDL_MIN -1e+30 ///< 定义变量的最小值
#define PRECISION 16 ///< 定义小数点后需要使用的位数
#define ZERO 1e-16 ///< 定义零值
//数学常量
#define BDL_MAX 1e+30 ///< 定义变量的最大值
#define BDL_MIN -1e+30 ///< 定义变量的最小值
#define PRECISION 16 ///< 定义小数点后需要使用的位数
#define ZERO 1e-16 ///< 定义零值
//物理常量
#define Pi (4.0*atan(1.0)) ///< 圆周率
#define G0 6.67408e-3 ///< 万有引力常数。注意这里的量级本来应该是e-11，考虑到单位转换，取维度单位为m，密度单位为g/cm^3，乘以G0则重力单位即为mGal
#define T0 5.0e+4 ///< 地磁场平均强度
//宏函数
#define MAX(a,b) (a>b?a:b) ///< 返回a与b的最大值
#define MIN(a,b) (a<b?a:b) ///< 返回a与b的最小值
#define SetToBox(a,b,in) (MAX(a,MIN(b,in))) ///< 返回a与b之间的一个值，若in在a与b之间则直接返回，否则返回较近的边界值
//终端显示控制符
#define BOLDRED "\033[1m\033[31m" ///< 设置后续字符字体为红色加粗
#define BOLDGREEN "\033[1m\033[32m" ///< 设置后续字符字体为绿色加粗
#define BOLDYELLOW "\033[1m\033[33m" ///< 设置后续字符字体为黄色加粗
#define BOLDBLUE "\033[1m\033[34m" ///< 设置后续字符字体为蓝色加粗
#define UNDERLINE "\033[1m\033[4m" ///< 设置后续字符为添加下划线
#define RESET "\033[0m" ///< 重置字符设置
#define MOVEUP(x) printf("\033[%dA", (x)) ///< 将光标向上挪x行
#define MOVEDOWN(x) printf("\033[%dB", (x)) ///< 将光标向下娜x行
#define MOVELEFT(x) printf("\033[%dD", (x)) ///< 将光标向左娜x字符
#define MOVERIGHT(x) printf("\033[%dC", (x)) ///< 将光标向右娜x字符
#define MOVETO(y,x) printf("\033[%d;%dH", (y), (x)) ///< 将光标向右娜动y字符,向上挪动x字符
#define CLEARLINE "\033[K" ///< 清除本行
#define CLEARALL "\033[2J" ///< 清除终端满屏
//数据结构
typedef vector<int> _1iArray; ///< 整形一维向量
typedef vector<double> _1dArray; ///< 双精度浮点一维向量
typedef vector<string> _1sArray; ///< 字符串一维向量
typedef vector<vector<int> > _2iArray; ///< 整形浮点二维向量
typedef vector<vector<double> > _2dArray; ///< 双精度浮点二维向量
typedef map<int,int> _int2intMap; ///< 整型到整形的映射
//模型块体参数
//物理常量
#define Pi (4.0*atan(1.0)) ///< 圆周率
#define G0 6.67408e-3 ///< 万有引力常数。注意这里的量级本来应该是e-11，考虑到单位转换，取维度单位为m，密度单位为g/cm^3，乘以G0则重力单位即为mGal
#define T0 5.0e+4 ///< 地磁场平均强度
//宏函数
#define MAX(a,b) (a>b?a:b) ///< 返回a与b的最大值
#define MIN(a,b) (a<b?a:b) ///< 返回a与b的最小值
#define SetToBox(a,b,in) (MAX(a,MIN(b,in))) ///< 返回a与b之间的一个值，若in在a与b之间则直接返回，否则返回较近的边界值
//终端显示控制符
#define BOLDRED "\033[1m\033[31m" ///< 设置后续字符字体为红色加粗
#define BOLDGREEN "\033[1m\033[32m" ///< 设置后续字符字体为绿色加粗
#define BOLDYELLOW "\033[1m\033[33m" ///< 设置后续字符字体为黄色加粗
#define BOLDBLUE "\033[1m\033[34m" ///< 设置后续字符字体为蓝色加粗
#define UNDERLINE "\033[1m\033[4m" ///< 设置后续字符为添加下划线
#define RESET "\033[0m" ///< 重置字符设置
#define MOVEUP(x) printf("\033[%dA", (x)) ///< 将光标向上挪x行
#define MOVEDOWN(x) printf("\033[%dB", (x)) ///< 将光标向下娜x行
#define MOVELEFT(x) printf("\033[%dD", (x)) ///< 将光标向左娜x字符
#define MOVERIGHT(x) printf("\033[%dC", (x)) ///< 将光标向右娜x字符
#define MOVETO(y,x) printf("\033[%d;%dH", (y), (x)) ///< 将光标向右娜动y字符,向上挪动x字符
#define CLEARLINE "\033[K" ///< 清除本行
#define CLEARALL "\033[2J" ///< 清除终端满屏
//数据结构
typedef vector<int> _1iArray; ///< 整形一维向量
typedef vector<double> _1dArray; ///< 双精度浮点一维向量
typedef vector<string> _1sArray; ///< 字符串一维向量
typedef vector<vector<int> > _2iArray; ///< 整形浮点二维向量
typedef vector<vector<double> > _2dArray; ///< 双精度浮点二维向量
typedef map<int,int> _int2intMap; ///< 整型到整形的映射
//模型块体参数
struct modelist{
char mod_type[1024];
char val_type[1024];
@ -61,33 +61,33 @@ struct modelist{
double mod_value;
};
typedef vector<modelist> modelistArray;
//直角坐标系点
//直角坐标系点
struct cpoint{
int id = -1;
double x = BDL_MAX; double y = BDL_MAX; double z = BDL_MAX;
};
typedef vector<cpoint> cpointArray;
typedef map<string,cpoint> _str2pointMap;
//观测点
//观测点
struct obspoint : public cpoint{
double val = BDL_MAX; double dev = BDL_MAX;
};
typedef vector<obspoint> obspointArray;
//块体
//块体
struct cube{
cpoint cen;
int ids[8] = {-1,-1,-1,-1,-1,-1,-1,-1};
double dx = BDL_MAX; double dy = BDL_MAX; double dz = BDL_MAX;
};
typedef vector<cube> cubeArray;
/*************************数据结构函数********************************/
cpoint operator -(cpoint,cpoint); //矢量减法
double modCpoint(cpoint); //矢量模
/*************************全局函数********************************/
double arctg(double); //正负分离的atan函数 正数返回atan 负数返回atan+pi
stringstream str2ss(string); //将string转换为stringstream
string cpoint_id(cpoint); //返回一个cpoint的位置id
int open_infile(ifstream&,char*); //测试打开输入文件 如果成功则返回0并输出信息 否则返回1
int open_outfile(ofstream&,char*); //测试打开输出文件 如果成功则返回0并输出信息 否则返回1
double grid_interpolate(double,double,double,double,double,double,double,double,double,double); //规则网络插值
/*************************数据结构函数********************************/
cpoint operator -(cpoint,cpoint); //矢量减法
double modCpoint(cpoint); //矢量模
/*************************全局函数********************************/
double arctg(double); //正负分离的atan函数 正数返回atan 负数返回atan+pi
stringstream str2ss(string); //将string转换为stringstream
string cpoint_id(cpoint); //返回一个cpoint的位置id
int open_infile(ifstream&,char*); //测试打开输入文件 如果成功则返回0并输出信息 否则返回1
int open_outfile(ofstream&,char*); //测试打开输出文件 如果成功则返回0并输出信息 否则返回1
double grid_interpolate(double,double,double,double,double,double,double,double,double,double); //规则网络插值
#endif

6
src/init_obs.cpp
View File

@ -8,7 +8,7 @@ int GM3D::InitObs(char* obs_para){
double xmin,xmax,ymin,ymax;
double xs,xe,ys,ye,eleva,dx,dy;
//按格式解析参数 初始化观测位置 用于正演计算
//按格式解析参数 初始化观测位置 用于正演计算
if (7 == sscanf(obs_para,"%lf/%lf/%lf/%lf/%lf/%lf/%lf",&xs,&dx,&xe,&ys,&dy,&ye,&eleva)){
xmin = MIN(xs,xe); xmax = MAX(xs,xe);
ymin = MIN(ys,ye); ymax = MAX(ys,ye);
@ -26,7 +26,7 @@ int GM3D::InitObs(char* obs_para){
x += dx;
}
}
//解析失败 按文件读入 用于反演使用或者正演计算
//解析失败 按文件读入 用于反演使用或者正演计算
else{
ifstream infile;
if (open_infile(infile,obs_para)) return -1;
@ -34,7 +34,7 @@ int GM3D::InitObs(char* obs_para){
while(getline(infile,temp_str)){
if (*(temp_str.begin()) == '#') continue;
else{
//按每行3个数据解析 初始化为用于正演的观测点
//按每行3个数据解析 初始化为用于正演的观测点
if (3 == sscanf(temp_str.c_str(),"%lf %lf %lf",&temp_obs.y,&temp_obs.x,&temp_obs.z)){
temp_obs.z *= -1.0;
temp_obs.id = obs_p_.size();

27
src/main.cpp
View File

@ -3,31 +3,25 @@
#include "gm3d_gui.h"
int main(int argc, char **argv) {
// use high resolution mode for "retina display"
// Fl::use_high_res_GL(1);
{ main_window = new Fl_Double_Window(500, 600, "gm3d");
{ main_tabs = new Fl_Tabs(10, 10, 480, 580);
{ model_tab = new Fl_Group(10, 40, 480, 550, "Build Model");
//model_tab->hide();
{ mesh_para_input = new Fl_Input(40, 70, 300, 28, "Input Mesh Parameters :");
mesh_para_input->callback((Fl_Callback*)cb_mesh_para_input);
mesh_para_input->align(Fl_Align(FL_ALIGN_TOP_LEFT));
mesh_para_input->tooltip("mesh parameters:\n<xmin>/<dx>/<xmax>/<ymin>/<dy>/<ymax>/<zmin>/<dz>/<zmax>");
} // Fl_Input* mesh_para_input
{ mesh_file_btn = new Fl_Button(360, 70, 100, 28, "Mesh File");
mesh_file_btn->callback((Fl_Callback*)cb_mesh_file_btn);
mesh_file_btn->tooltip("read mesh parameters from an input file.");
} // Fl_Button* mesh_file_btn
{ mod_para_file_btn = new Fl_Button(170, 160, 160, 28, "Add Model From File");
mod_para_file_btn->callback((Fl_Callback*)cb_mod_para_file_btn);
mod_para_file_btn->tooltip("read model parameters from an input file.");
} // Fl_Button* mod_para_file_btn
{ mod_ele_input_build = new Fl_Input(40, 445, 215, 28, "Input Model Element Data Name:");
mod_ele_input_build->align(Fl_Align(FL_ALIGN_TOP_LEFT));
mod_ele_input_build->tooltip("a string. e.g. my_model");
} // Fl_Input* mod_ele_input_build
{ build_mod_btn = new Fl_Button(40, 545, 420, 28, "Build Model !");
build_mod_btn->callback((Fl_Callback*)cb_build_mod_btn);
build_mod_btn->tooltip("callback a gm3d function to build the model.");
} // Fl_Button* build_mod_btn
{ mesh_para_output = new Fl_Output(40, 120, 300, 28, "Mesh Parameters :");
mesh_para_output->value("Unset");
@ -35,14 +29,16 @@ int main(int argc, char **argv) {
mesh_para_output->color(FL_BACKGROUND_COLOR);
mesh_para_output->align(Fl_Align(FL_ALIGN_TOP_LEFT));
} // Fl_Output* mesh_para_output
{ edit_mod_btn = new Fl_Button(40, 160, 100, 28, "Model Editor");
edit_mod_btn->callback((Fl_Callback*)cb_edit_mod_btn);
edit_mod_btn->tooltip("add or edit model in a gui interface.");
} // Fl_Button* edit_mod_btn
{ add_mod_btn = new Fl_Button(40, 160, 100, 28, "Add Model");
add_mod_btn->callback((Fl_Callback*)cb_add_mod_btn);
} // Fl_Button* add_mod_btn
{ del_mod_btn = new Fl_Button(360, 160, 100, 28, "Delete Model");
del_mod_btn->callback((Fl_Callback*)cb_del_mod_btn);
del_mod_btn->deactivate();
} // Fl_Button* del_mod_btn
{ rm_emp_bok_check = new Fl_Check_Button(280, 447, 180, 28, "Remove Empty Blocks");
rm_emp_bok_check->down_box(FL_DOWN_BOX);
rm_emp_bok_check->callback((Fl_Callback*)cb_rm_emp_bok_check);
rm_emp_bok_check->tooltip("suppress the output of blocks with zero values.");
} // Fl_Check_Button* rm_emp_bok_check
{ mod_para_brw = new Fl_Browser(40, 205, 420, 210);
mod_para_brw->callback((Fl_Callback*)cb_mod_para_brw);
@ -50,7 +46,6 @@ int main(int argc, char **argv) {
} // Fl_Browser* mod_para_brw
{ mod_file_out_btn = new Fl_Button(360, 500, 100, 28, "Model File");
mod_file_out_btn->callback((Fl_Callback*)cb_mod_file_out_btn);
mod_file_out_btn->tooltip("create or choose a filename for outputting.");
} // Fl_Button* mod_file_out_btn
{ mod_out_file_output = new Fl_Output(40, 500, 300, 28, "Output File Name :");
mod_out_file_output->box(FL_FLAT_BOX);
@ -96,6 +91,7 @@ int main(int argc, char **argv) {
{ obs_file_input = new Fl_Input(40, 130, 300, 28, "Input Observe FIle :");
obs_file_input->callback((Fl_Callback*)cb_obs_file_input);
obs_file_input->align(Fl_Align(FL_ALIGN_TOP_LEFT));
obs_file_input->tooltip("<xmin>/<dx>/<xmax>/<ymin>/<dy>/<ymax>/<elevation>|<filename>");
} // Fl_Input* obs_file_input
{ obs_file_output = new Fl_Output(260, 240, 200, 28, "Chosen Observe File :");
obs_file_output->box(FL_FLAT_BOX);
@ -136,7 +132,7 @@ int main(int argc, char **argv) {
Za_check->deactivate();
} // Fl_Check_Button* Za_check
{ mag_para_input = new Fl_Input(220, 440, 240, 28, "Magnetization Parameters : ");
mag_para_input->tooltip("<I0>/<D0>/<I>/<D>");
mag_para_input->tooltip("<geo-inclination>/<geo-declination>/<mag-inclination>/<mag-declination>");
mag_para_input->deactivate();
} // Fl_Input* mag_para_input
mag_group->end();
@ -149,7 +145,7 @@ int main(int argc, char **argv) {
cal_btn->callback((Fl_Callback*)cb_cal_btn);
} // Fl_Button* cal_btn
{ noise_para_input = new Fl_Input(260, 370, 200, 28);
noise_para_input->tooltip("<noise-mean>/<nosie-deviation>");
noise_para_input->tooltip("<noise-mean>/<noise-deviation>");
noise_para_input->align(Fl_Align(FL_ALIGN_TOP_LEFT));
noise_para_input->deactivate();
} // Fl_Input* noise_para_input
@ -176,7 +172,6 @@ int main(int argc, char **argv) {
} // Fl_Tabs* main_tabs
main_window->end();
} // Fl_Double_Window* main_window
// show software interface
main_window->show(argc, argv);
return Fl::run();
}

7
src/out_meshtool_file.cpp
View File

@ -1,7 +0,0 @@
#include "gm3d.h"
int GM3D::OutMeshtoolFile(char* filename)
{
//
return 0;
}

4
src/out_msh_file.cpp
View File

@ -6,9 +6,9 @@ int GM3D::OutMshFile(char* filename,string data_name){
ofstream outfile;
if (open_outfile(outfile,filename)) return -1;
//好啦 我们这里输出的模型类型应该是块体
//好啦 我们这里输出的模型类型应该是块体
outfile<<"$MeshFormat"<<endl<<"2.2 0 8"<<endl<<"$EndMeshFormat"<<endl<<"$Nodes"<<endl<< out_vert_ids_.size() <<endl;
//输出顶点 最后一个不输出
//输出顶点 最后一个不输出
for (int i = 0; i < out_vert_ids_.size(); i++){
outfile << i << " "
<< setprecision(16) << model_vert_[out_vert_ids_[i]].x << " "

18
src/out_neighbor_file.cpp
View File

@ -3,22 +3,22 @@
int GM3D::OutNeighborFile(char* v_name,char* b_name){
/*
if (strcmp(v_name,"NULL") || strcmp(b_name,"NULL")){
//整理块体间的相邻关系 先初始化顶点相邻数组
//整理块体间的相邻关系 先初始化顶点相邻数组
model_vert_neighbor_.resize(vert_num_);
for (int i = 0; i < vert_num_; i++)
model_vert_neighbor_[i].resize(8,-1);
//遍历所有块体整理顶点相邻关系
//遍历所有块体整理顶点相邻关系
for (int i = 0; i < model_num_; i++){
for (int j = 0; j < 8; j++){
model_vert_neighbor_[model_cube_[i].ids[j]][j] = model_cube_[i].cen.id;
}
}
//遍历所有顶点相邻关系 所有共点的块体都被认定为相邻块体
//遍历所有顶点相邻关系 所有共点的块体都被认定为相邻块体
model_cube_neighbor_.resize(model_num_);
//循环顶点相邻列表 都不为-1则相互添加
//循环顶点相邻列表 都不为-1则相互添加
for (int i = 0; i < vert_num_; i++){
for (int n = 0; n < 8; n++){
if (model_vert_neighbor_[i][n] != -1){
@ -32,14 +32,14 @@ int GM3D::OutNeighborFile(char* v_name,char* b_name){
}
}
vector<int>::iterator pos; //整型向量的迭代器
vector<int>::iterator pos; //整型向量的迭代器
for (int i = 0; i < model_num_; i++){
sort(model_cube_neighbor_[i].begin(),model_cube_neighbor_[i].end()); //对顶点序列由小到大排序
pos = unique(model_cube_neighbor_[i].begin(),model_cube_neighbor_[i].end()); //获取重复序列开始的位置
model_cube_neighbor_[i].erase(pos,model_cube_neighbor_[i].end()); //删除重复点
sort(model_cube_neighbor_[i].begin(),model_cube_neighbor_[i].end()); //对顶点序列由小到大排序
pos = unique(model_cube_neighbor_[i].begin(),model_cube_neighbor_[i].end()); //获取重复序列开始的位置
model_cube_neighbor_[i].erase(pos,model_cube_neighbor_[i].end()); //删除重复点
}
//清理数组
//清理数组
for (int i = 0; i < vert_num_; i++){
model_vert_neighbor_[i].clear();
vector <int>().swap(model_vert_neighbor_[i]);

2
src/progress_bar.cpp
View File

@ -47,7 +47,7 @@ int ProgressBar::GetConsoleWidth(){
width = csbi.srWindow.Right - csbi.srWindow.Left;
#else
struct winsize win;
//注意！当我们使用pipe here-doc等通道获取程序参数时无法正确的获取窗口大小 此时我们将使用预定值
//注意！当我们使用pipe here-doc等通道获取程序参数时无法正确的获取窗口大小 此时我们将使用预定值
if (ioctl(0, TIOCGWINSZ, &win) != -1)
width = win.ws_col;
else width = 100;

15
src/progress_bar.h
View File

@ -1,23 +1,12 @@
#ifndef _PROGRESS_BAR_
#ifndef _PROGRESS_BAR_
#define _PROGRESS_BAR_
//#include </*sys/ioctl.h*/>
#include <iostream>
#include <iomanip>
#include <cstring>
#include <thread>
#include <chrono>
#if defined _WINDOWS || __WIN32__
#include <windows.h>
#else
#include <sys/ioctl.h>
#include <unistd.h>
#endif
#define TOTAL_PERCENTAGE 100.0
#define CHARACTER_WIDTH_PERCENTAGE 4

6
src/read_meshtool_model.cpp
View File

@ -1,6 +0,0 @@
#include "gm3d.h"
int GM3D::ReadMeshtoolModel(char* filename)
{
return 0;
}

44
src/read_model.cpp
View File

@ -10,14 +10,14 @@ int GM3D::ReadModel(char* filename,char* input_forward_model_name){
stringstream temp_ss;
ifstream mshin;
if (open_infile(mshin,filename)) return -1; //检查并打开模型文件
if (open_infile(mshin,filename)) return -1; //检查并打开模型文件
while(getline(mshin,temp_str)){
//读入模型空间顶点集 msh文件版本为2.2
//读入模型空间顶点集 msh文件版本为2.2
if (temp_str == "$Nodes"){
getline(mshin,temp_str);
temp_ss = str2ss(temp_str);
temp_ss >> vert_num_; //第一个数为顶点的个数
model_vert_.resize(vert_num_); //开辟空间
temp_ss >> vert_num_; //第一个数为顶点的个数
model_vert_.resize(vert_num_); //开辟空间
for (int i = 0; i < vert_num_; i++){
getline(mshin,temp_str);
temp_ss = str2ss(temp_str);
@ -25,19 +25,19 @@ int GM3D::ReadModel(char* filename,char* input_forward_model_name){
model_vert_[i] = temp_vert;
}
}
//读入模型空间单元体
//读入模型空间单元体
else if (temp_str == "$Elements"){
getline(mshin,temp_str);
temp_ss = str2ss(temp_str);
temp_ss >> model_num_; //第一个数为总元素的个数 包含了所有类型的元素 比如三角形 四边形 块体等
temp_ss >> model_num_; //第一个数为总元素的个数 包含了所有类型的元素 比如三角形 四边形 块体等
model_cube_.resize(model_num_);
for (int i = 0; i < model_num_; i++){
getline(mshin,temp_str);
temp_ss = str2ss(temp_str);
temp_ss >> temp_cu.cen.id >> ele_type;
//只读入块体
//只读入块体
if (ele_type == 5){
temp_ss >> attri_num; //跳过模型单元的几何组与物理组等信息 以后可能会有用
temp_ss >> attri_num; //跳过模型单元的几何组与物理组等信息 以后可能会有用
for (int a = 0; a < attri_num; a++)
temp_ss >> temp_attri;
for (int a = 0; a < 8; a++)
@ -46,44 +46,44 @@ int GM3D::ReadModel(char* filename,char* input_forward_model_name){
}
}
}
else continue; //不能识别的单元都被忽略了
else continue; //不能识别的单元都被忽略了
}
mshin.close();
//第二次读入模型文件 初始化模型单元属性
if (open_infile(mshin,filename)) return -1; //检查并打开模型文件
//第二次读入模型文件 初始化模型单元属性
if (open_infile(mshin,filename)) return -1; //检查并打开模型文件
while(getline(mshin,temp_str)){
//读入模型单元属性 注意因为msh文件中$ElementData并未注明所属元素类型
//所以可能会将其他元素类型的属性值也读入 但因为其在pyIdMap中并未注册 所以属性值会全为0 在后续使用时我们需要通过名称辨别
//读入模型单元属性 注意因为msh文件中$ElementData并未注明所属元素类型
//所以可能会将其他元素类型的属性值也读入 但因为其在pyIdMap中并未注册 所以属性值会全为0 在后续使用时我们需要通过名称辨别
if (temp_str == "$ElementData"){
temp_model.resize(model_num_,0.0); //初始化temp_model 为读入模型单元属性做准备
for (int i = 0; i < 2; i++) //先读入元素块的名称 添加到数组
temp_model.resize(model_num_,0.0); //初始化temp_model 为读入模型单元属性做准备
for (int i = 0; i < 2; i++) //先读入元素块的名称 添加到数组
getline(mshin,temp_str);
input_model_names_.push_back(temp_str);
for (int i = 0; i < 6; i++) //跳过元素属性前面的值 最后一次为当前元素块的个数
for (int i = 0; i < 6; i++) //跳过元素属性前面的值 最后一次为当前元素块的个数
getline(mshin,temp_str);
temp_ss = str2ss(temp_str);
temp_ss >> temp_int;
for (int i = 0; i < temp_int; i++){
getline(mshin,temp_str);
temp_ss = str2ss(temp_str);
temp_ss >> temp_id >> temp_val; //读入单元体索引与属性值
temp_ss >> temp_id >> temp_val; //读入单元体索引与属性值
temp_model[temp_id] = temp_val;
}
input_models_.push_back(temp_model);
temp_model.clear();
}
else continue; //不能识别的单元都被忽略了
else continue; //不能识别的单元都被忽略了
}
mshin.close();
//清理映射
//清理映射
temp_model.clear();
vector <double>().swap(temp_model);
//初始化数组
//初始化数组
forward_model_.resize(model_num_,0.0);
//匹配数据名称
//匹配数据名称
char forward_model_name[1024] = "\"";
strcat(forward_model_name,input_forward_model_name);
strcat(forward_model_name,"\"");
@ -94,7 +94,7 @@ int GM3D::ReadModel(char* filename,char* input_forward_model_name){
}
}
//计算块体的中心位置和尺寸
//计算块体的中心位置和尺寸
cpoint corner[8];
for (int i = 0; i < model_num_; i++){
for (int j = 0; j < 8; j++){

51
src/read_model_mst.cpp Normal file
View File

@ -0,0 +1,51 @@
#include "gm3d.h"
int GM3D::ReadModel_mst(char* filename){
int temp_int,ele_type,attri_num,temp_attri,temp_id;
double temp_val;
_1dArray temp_model;
cpoint temp_vert;
cube temp_cu;
string temp_str;
stringstream temp_ss;
ifstream mstin;
if (open_infile(mstin,filename)) return -1; //¼ì²é²¢´ò¿ªÄ£ÐÍÎļþ
double ynum, xnum, znum;
double ystart, xstart, zstart;
double ylength, xlength, zlength;
char temp;
mstin>>ynum>>xnum>>znum;
mstin>>ystart>>xstart>>zstart;
mstin>>ynum>>temp>>ylength;
mstin>>xnum>>temp>>xlength;
mstin>>znum>>temp>>zlength;
model_num_ = ynum*xnum*znum;
model_cube_.resize(model_num_);
//³õʼ»¯Êý×é
forward_model_.resize(model_num_,0.0);
double f_mdl_tmp;
int i, j, k;
int tmp_num;
for(i = 0; i < xnum; i++){
for(j = 0; j < ynum; j++){
for(k = 0; k < znum; k++){
tmp_num = i*ynum*znum + j*znum + k;
mstin>>f_mdl_tmp;
forward_model_[tmp_num] = f_mdl_tmp;
model_cube_[tmp_num].cen.x = xstart + xlength*i + xlength*0.5;
model_cube_[tmp_num].cen.y = ystart + ylength*j + ylength*0.5;
model_cube_[tmp_num].cen.z = zstart + zlength*k + zlength*0.5;
model_cube_[tmp_num].dx = xlength;
model_cube_[tmp_num].dy = ylength;
model_cube_[tmp_num].dz = zlength;
}
}
}
mstin.close();
return 0;
}

26
src/registered_output.cpp
View File

@ -2,10 +2,10 @@
int GM3D::RegisteredOuput(bool remove_empty_element){
int count;
//统计输出模型单元块体和顶点 以及输出的块体数据列表
//统计输出模型单元块体和顶点 以及输出的块体数据列表
if (remove_empty_element){
count = 0;
//遍历所有块体数据 注册有值的块体
//遍历所有块体数据 注册有值的块体
for (int i = 0; i < model_num_; i++){
if (model_block_val_[i] != BDL_MAX){
out_ele_data_ids_.push_back(i);
@ -15,38 +15,38 @@ int GM3D::RegisteredOuput(bool remove_empty_element){
}
}
//遍历所有注册的块体 添加顶点
//遍历所有注册的块体 添加顶点
for (int i = 0; i < out_ele_ids_.size(); i++){
for (int j = 0; j < 8; j++){
out_vert_ids_.push_back(model_cube_[out_ele_ids_[i]].ids[j]);
}
}
//去除输出顶点中的重复部分
vector<int>::iterator pos; //整型向量的迭代器
sort(out_vert_ids_.begin(),out_vert_ids_.end()); //对顶点序列由小到大排序
pos = unique(out_vert_ids_.begin(),out_vert_ids_.end()); //获取重复序列开始的位置
out_vert_ids_.erase(pos,out_vert_ids_.end()); //删除重复点
//去除输出顶点中的重复部分
vector<int>::iterator pos; //整型向量的迭代器
sort(out_vert_ids_.begin(),out_vert_ids_.end()); //对顶点序列由小到大排序
pos = unique(out_vert_ids_.begin(),out_vert_ids_.end()); //获取重复序列开始的位置
out_vert_ids_.erase(pos,out_vert_ids_.end()); //删除重复点
//将需要输出的模型顶点序号与它们的排序做一个对应 保证在输出文件中顶点索引号始终是从0开始的连续的序列
//将需要输出的模型顶点序号与它们的排序做一个对应 保证在输出文件中顶点索引号始终是从0开始的连续的序列
for (int i = 0; i < out_vert_ids_.size(); i++){
vert_out_map_[out_vert_ids_[i]] = i;
}
}
else{
//输出的模型块体为所有
//输出的模型块体为所有
out_ele_ids_.resize(model_num_);
for (int i = 0; i < model_num_; i++){
out_ele_ids_[i] = i;
}
//输出的模型顶点为所有
//输出的模型顶点为所有
out_vert_ids_.resize(vert_num_);
for (int i = 0; i < vert_num_; i++){
vert_out_map_[i] = out_vert_ids_[i] = i;
}
//注册所有有值的块体数据
//注册所有有值的块体数据
count = 0;
for (int i = 0; i < model_num_; i++){
if (model_block_val_[i] != BDL_MAX){
@ -56,6 +56,6 @@ int GM3D::RegisteredOuput(bool remove_empty_element){
}
}
}
//输出所有模型单元块体和顶点 只统计输出的块体数据列表
//输出所有模型单元块体和顶点 只统计输出的块体数据列表
return 0;
}