commit d22e9eeb93c6a6092d409633eb510adc86f56ad3 Author: rowanc1 Date: Tue May 13 09:25:25 2014 -0700 Initial Commit. diff --git a/.travis.yml b/.travis.yml new file mode 100644 index 00000000..08e8ee46 --- /dev/null +++ b/.travis.yml @@ -0,0 +1,26 @@ +language: python +python: + - "2.7" +virtualenv: + system_site_packages: true +before_install: + - sudo apt-get install -qq gcc gfortran libblas-dev liblapack-dev python-numpy python-scipy python-matplotlib python-pip + - sudo pip install scipy --upgrade + - sudo pip install numpy --upgrade + - cd ../ + - git clone https://github.com/simpeg/simpeg.git + - cd simpeg/SimPEG/ + - python setup.py + - cd ../../ + - echo export PYTHONPATH=$PYTHONPATH:/home/travis/build/simpeg/simpeg >> .bashrc + - source .bashrc + - cd simpegmt +# command to install dependencies +install: "pip install -r requirements.txt --use-mirrors" +# command to run tests +script: nosetests -v + +notifications: + email: + - rowanc1@gmail.com + - gkrosen@gmail.com diff --git a/LICENSE b/LICENSE new file mode 100644 index 00000000..f94a23fd --- /dev/null +++ b/LICENSE @@ -0,0 +1,20 @@ +The MIT License (MIT) + +Copyright (c) 2013-2014 SimPEG Developers + +Permission is hereby granted, free of charge, to any person obtaining a copy of +this software and associated documentation files (the "Software"), to deal in +the Software without restriction, including without limitation the rights to +use, copy, modify, merge, publish, distribute, sublicense, and/or sell copies of +the Software, and to permit persons to whom the Software is furnished to do so, +subject to the following conditions: + +The above copyright notice and this permission notice shall be included in all +copies or substantial portions of the Software. + +THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR +IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS +FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR +COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER +IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN +CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE. diff --git a/README.md b/README.md new file mode 100644 index 00000000..1f8dc810 --- /dev/null +++ b/README.md @@ -0,0 +1,21 @@ +simpegem +======== + +A electromagnetic forward modelling and inversion package for SimPEG. + + + +Documentation: +[http://simpegem.readthedocs.org/en/latest/](http://simpegem.readthedocs.org/en/latest/) + +Code: +[https://github.com/simpeg/simpegem](https://github.com/simpeg/simpegem) + +Tests: +[https://travis-ci.org/simpeg/simpegem](https://travis-ci.org/simpeg/simpegem) + +Build Status: +[![Build Status](https://travis-ci.org/simpeg/simpegem.svg?branch=master)](https://travis-ci.org/simpeg/simpegem) + +Bugs & Issues: +[https://github.com/simpeg/simpegem/issues](https://github.com/simpeg/simpegem/issues) diff --git a/docs/Makefile b/docs/Makefile new file mode 100644 index 00000000..bee80244 --- /dev/null +++ b/docs/Makefile @@ -0,0 +1,153 @@ +# Makefile for Sphinx documentation +# + +# You can set these variables from the command line. +SPHINXOPTS = +SPHINXBUILD = sphinx-build +PAPER = +BUILDDIR = _build + +# Internal variables. +PAPEROPT_a4 = -D latex_paper_size=a4 +PAPEROPT_letter = -D latex_paper_size=letter +ALLSPHINXOPTS = -d $(BUILDDIR)/doctrees $(PAPEROPT_$(PAPER)) $(SPHINXOPTS) . +# the i18n builder cannot share the environment and doctrees with the others +I18NSPHINXOPTS = $(PAPEROPT_$(PAPER)) $(SPHINXOPTS) . + +.PHONY: help clean html dirhtml singlehtml pickle json htmlhelp qthelp devhelp epub latex latexpdf text man changes linkcheck doctest gettext + +help: + @echo "Please use \`make ' where is one of" + @echo " html to make standalone HTML files" + @echo " dirhtml to make HTML files named index.html in directories" + @echo " singlehtml to make a single large HTML file" + @echo " pickle to make pickle files" + @echo " json to make JSON files" + @echo " htmlhelp to make HTML files and a HTML help project" + @echo " qthelp to make HTML files and a qthelp project" + @echo " devhelp to make HTML files and a Devhelp project" + @echo " epub to make an epub" + @echo " latex to make LaTeX files, you can set PAPER=a4 or PAPER=letter" + @echo " latexpdf to make LaTeX files and run them through pdflatex" + @echo " text to make text files" + @echo " man to make manual pages" + @echo " texinfo to make Texinfo files" + @echo " info to make Texinfo files and run them through makeinfo" + @echo " gettext to make PO message catalogs" + @echo " changes to make an overview of all changed/added/deprecated items" + @echo " linkcheck to check all external links for integrity" + @echo " doctest to run all doctests embedded in the documentation (if enabled)" + +clean: + -rm -rf $(BUILDDIR)/* + +html: + $(SPHINXBUILD) -b html $(ALLSPHINXOPTS) $(BUILDDIR)/html + @echo + @echo "Build finished. The HTML pages are in $(BUILDDIR)/html." + +dirhtml: + $(SPHINXBUILD) -b dirhtml $(ALLSPHINXOPTS) $(BUILDDIR)/dirhtml + @echo + @echo "Build finished. The HTML pages are in $(BUILDDIR)/dirhtml." + +singlehtml: + $(SPHINXBUILD) -b singlehtml $(ALLSPHINXOPTS) $(BUILDDIR)/singlehtml + @echo + @echo "Build finished. The HTML page is in $(BUILDDIR)/singlehtml." + +pickle: + $(SPHINXBUILD) -b pickle $(ALLSPHINXOPTS) $(BUILDDIR)/pickle + @echo + @echo "Build finished; now you can process the pickle files." + +json: + $(SPHINXBUILD) -b json $(ALLSPHINXOPTS) $(BUILDDIR)/json + @echo + @echo "Build finished; now you can process the JSON files." + +htmlhelp: + $(SPHINXBUILD) -b htmlhelp $(ALLSPHINXOPTS) $(BUILDDIR)/htmlhelp + @echo + @echo "Build finished; now you can run HTML Help Workshop with the" \ + ".hhp project file in $(BUILDDIR)/htmlhelp." + +qthelp: + $(SPHINXBUILD) -b qthelp $(ALLSPHINXOPTS) $(BUILDDIR)/qthelp + @echo + @echo "Build finished; now you can run "qcollectiongenerator" with the" \ + ".qhcp project file in $(BUILDDIR)/qthelp, like this:" + @echo "# qcollectiongenerator $(BUILDDIR)/qthelp/SimPEG.qhcp" + @echo "To view the help file:" + @echo "# assistant -collectionFile $(BUILDDIR)/qthelp/SimPEG.qhc" + +devhelp: + $(SPHINXBUILD) -b devhelp $(ALLSPHINXOPTS) $(BUILDDIR)/devhelp + @echo + @echo "Build finished." + @echo "To view the help file:" + @echo "# mkdir -p $$HOME/.local/share/devhelp/SimPEG" + @echo "# ln -s $(BUILDDIR)/devhelp $$HOME/.local/share/devhelp/SimPEG" + @echo "# devhelp" + +epub: + $(SPHINXBUILD) -b epub $(ALLSPHINXOPTS) $(BUILDDIR)/epub + @echo + @echo "Build finished. The epub file is in $(BUILDDIR)/epub." + +latex: + $(SPHINXBUILD) -b latex $(ALLSPHINXOPTS) $(BUILDDIR)/latex + @echo + @echo "Build finished; the LaTeX files are in $(BUILDDIR)/latex." + @echo "Run \`make' in that directory to run these through (pdf)latex" \ + "(use \`make latexpdf' here to do that automatically)." + +latexpdf: + $(SPHINXBUILD) -b latex $(ALLSPHINXOPTS) $(BUILDDIR)/latex + @echo "Running LaTeX files through pdflatex..." + $(MAKE) -C $(BUILDDIR)/latex all-pdf + @echo "pdflatex finished; the PDF files are in $(BUILDDIR)/latex." + +text: + $(SPHINXBUILD) -b text $(ALLSPHINXOPTS) $(BUILDDIR)/text + @echo + @echo "Build finished. The text files are in $(BUILDDIR)/text." + +man: + $(SPHINXBUILD) -b man $(ALLSPHINXOPTS) $(BUILDDIR)/man + @echo + @echo "Build finished. The manual pages are in $(BUILDDIR)/man." + +texinfo: + $(SPHINXBUILD) -b texinfo $(ALLSPHINXOPTS) $(BUILDDIR)/texinfo + @echo + @echo "Build finished. The Texinfo files are in $(BUILDDIR)/texinfo." + @echo "Run \`make' in that directory to run these through makeinfo" \ + "(use \`make info' here to do that automatically)." + +info: + $(SPHINXBUILD) -b texinfo $(ALLSPHINXOPTS) $(BUILDDIR)/texinfo + @echo "Running Texinfo files through makeinfo..." + make -C $(BUILDDIR)/texinfo info + @echo "makeinfo finished; the Info files are in $(BUILDDIR)/texinfo." + +gettext: + $(SPHINXBUILD) -b gettext $(I18NSPHINXOPTS) $(BUILDDIR)/locale + @echo + @echo "Build finished. The message catalogs are in $(BUILDDIR)/locale." + +changes: + $(SPHINXBUILD) -b changes $(ALLSPHINXOPTS) $(BUILDDIR)/changes + @echo + @echo "The overview file is in $(BUILDDIR)/changes." + +linkcheck: + $(SPHINXBUILD) -b linkcheck $(ALLSPHINXOPTS) $(BUILDDIR)/linkcheck + @echo + @echo "Link check complete; look for any errors in the above output " \ + "or in $(BUILDDIR)/linkcheck/output.txt." + +doctest: + $(SPHINXBUILD) -b doctest $(ALLSPHINXOPTS) $(BUILDDIR)/doctest + @echo "Testing of doctests in the sources finished, look at the " \ + "results in $(BUILDDIR)/doctest/output.txt." diff --git a/docs/conf.py b/docs/conf.py new file mode 100644 index 00000000..86aacce4 --- /dev/null +++ b/docs/conf.py @@ -0,0 +1,244 @@ +# -*- coding: utf-8 -*- +# +# SimPEG documentation build configuration file, created by +# sphinx-quickstart on Fri Aug 30 18:42:44 2013. +# +# This file is execfile()d with the current directory set to its containing dir. +# +# Note that not all possible configuration values are present in this +# autogenerated file. +# +# All configuration values have a default; values that are commented out +# serve to show the default. + +import sys, os + +# If extensions (or modules to document with autodoc) are in another directory, +# add these directories to sys.path here. If the directory is relative to the +# documentation root, use os.path.abspath to make it absolute, like shown here. +#sys.path.insert(0, os.path.abspath('.')) + +sys.path.append('../') + +# -- General configuration ----------------------------------------------------- + +# If your documentation needs a minimal Sphinx version, state it here. +#needs_sphinx = '1.0' + +# Add any Sphinx extension module names here, as strings. They can be extensions +# coming with Sphinx (named 'sphinx.ext.*') or your custom ones. +extensions = ['sphinx.ext.todo', 'sphinx.ext.mathjax', 'sphinx.ext.viewcode', 'sphinx.ext.autodoc', 'matplotlib.sphinxext.plot_directive'] + +# Add any paths that contain templates here, relative to this directory. +templates_path = ['_templates'] + +# The suffix of source filenames. +source_suffix = '.rst' + +# The encoding of source files. +#source_encoding = 'utf-8-sig' + +# The master toctree document. +master_doc = 'index' + +# General information about the project. +project = u'SimPEG' +copyright = u'2013, SimPEG Developers' + +# The version info for the project you're documenting, acts as replacement for +# |version| and |release|, also used in various other places throughout the +# built documents. +# +# The short X.Y version. +version = '0.0.1' +# The full version, including alpha/beta/rc tags. +release = '0.0.1' + +# The language for content autogenerated by Sphinx. Refer to documentation +# for a list of supported languages. +#language = None + +# There are two options for replacing |today|: either, you set today to some +# non-false value, then it is used: +#today = '' +# Else, today_fmt is used as the format for a strftime call. +#today_fmt = '%B %d, %Y' + +# List of patterns, relative to source directory, that match files and +# directories to ignore when looking for source files. +exclude_patterns = ['_build'] + +# The reST default role (used for this markup: `text`) to use for all documents. +#default_role = None + +# If true, '()' will be appended to :func: etc. cross-reference text. +#add_function_parentheses = True + +# If true, the current module name will be prepended to all description +# unit titles (such as .. function::). +#add_module_names = True + +# If true, sectionauthor and moduleauthor directives will be shown in the +# output. They are ignored by default. +#show_authors = False + +# The name of the Pygments (syntax highlighting) style to use. +pygments_style = 'sphinx' + +# A list of ignored prefixes for module index sorting. +#modindex_common_prefix = [] + + +# -- Options for HTML output --------------------------------------------------- + +# The theme to use for HTML and HTML Help pages. See the documentation for +# a list of builtin themes. +html_theme = 'default' + +# Theme options are theme-specific and customize the look and feel of a theme +# further. For a list of options available for each theme, see the +# documentation. +#html_theme_options = {} + +# Add any paths that contain custom themes here, relative to this directory. +#html_theme_path = [] + +# The name for this set of Sphinx documents. If None, it defaults to +# " v documentation". +#html_title = None + +# A shorter title for the navigation bar. Default is the same as html_title. +#html_short_title = None + +# The name of an image file (relative to this directory) to place at the top +# of the sidebar. +#html_logo = None + +# The name of an image file (within the static path) to use as favicon of the +# docs. This file should be a Windows icon file (.ico) being 16x16 or 32x32 +# pixels large. +#html_favicon = None + +# Add any paths that contain custom static files (such as style sheets) here, +# relative to this directory. They are copied after the builtin static files, +# so a file named "default.css" will overwrite the builtin "default.css". +html_static_path = ['_static'] + +# If not '', a 'Last updated on:' timestamp is inserted at every page bottom, +# using the given strftime format. +#html_last_updated_fmt = '%b %d, %Y' + +# If true, SmartyPants will be used to convert quotes and dashes to +# typographically correct entities. +#html_use_smartypants = True + +# Custom sidebar templates, maps document names to template names. +#html_sidebars = {} + +# Additional templates that should be rendered to pages, maps page names to +# template names. +#html_additional_pages = {} + +# If false, no module index is generated. +#html_domain_indices = True + +# If false, no index is generated. +#html_use_index = True + +# If true, the index is split into individual pages for each letter. +#html_split_index = False + +# If true, links to the reST sources are added to the pages. +#html_show_sourcelink = True + +# If true, "Created using Sphinx" is shown in the HTML footer. Default is True. +#html_show_sphinx = True + +# If true, "(C) Copyright ..." is shown in the HTML footer. Default is True. +#html_show_copyright = True + +# If true, an OpenSearch description file will be output, and all pages will +# contain a tag referring to it. The value of this option must be the +# base URL from which the finished HTML is served. +#html_use_opensearch = '' + +# This is the file name suffix for HTML files (e.g. ".xhtml"). +#html_file_suffix = None + +# Output file base name for HTML help builder. +htmlhelp_basename = 'SimPEGdoc' + + +# -- Options for LaTeX output -------------------------------------------------- + +latex_elements = { +# The paper size ('letterpaper' or 'a4paper'). +#'papersize': 'letterpaper', + +# The font size ('10pt', '11pt' or '12pt'). +#'pointsize': '10pt', + +# Additional stuff for the LaTeX preamble. +#'preamble': '', +} + +# Grouping the document tree into LaTeX files. List of tuples +# (source start file, target name, title, author, documentclass [howto/manual]). +latex_documents = [ + ('index', 'SimPEG.tex', u'SimPEG Documentation', + u'Rowan Cockett', 'manual'), +] + +# The name of an image file (relative to this directory) to place at the top of +# the title page. +#latex_logo = None + +# For "manual" documents, if this is true, then toplevel headings are parts, +# not chapters. +#latex_use_parts = False + +# If true, show page references after internal links. +#latex_show_pagerefs = False + +# If true, show URL addresses after external links. +#latex_show_urls = False + +# Documents to append as an appendix to all manuals. +#latex_appendices = [] + +# If false, no module index is generated. +#latex_domain_indices = True + + +# -- Options for manual page output -------------------------------------------- + +# One entry per manual page. List of tuples +# (source start file, name, description, authors, manual section). +man_pages = [ + ('index', 'simpeg', u'SimPEG Documentation', + [u'Rowan Cockett'], 1) +] + +# If true, show URL addresses after external links. +#man_show_urls = False + + +# -- Options for Texinfo output ------------------------------------------------ + +# Grouping the document tree into Texinfo files. List of tuples +# (source start file, target name, title, author, +# dir menu entry, description, category) +texinfo_documents = [ + ('index', 'SimPEG', u'SimPEG Documentation', + u'Rowan Cockett', 'SimPEG', 'One line description of project.', + 'Miscellaneous'), +] + +# Documents to append as an appendix to all manuals. +#texinfo_appendices = [] + +# If false, no module index is generated. +#texinfo_domain_indices = True + +# How to display URL addresses: 'footnote', 'no', or 'inline'. +#texinfo_show_urls = 'footnote' diff --git a/docs/index.rst b/docs/index.rst new file mode 100644 index 00000000..69dcd3a0 --- /dev/null +++ b/docs/index.rst @@ -0,0 +1,34 @@ +.. image:: simpeg-logo.png + :width: 300 px + :alt: SimPEG + :align: center + + +SimPEG for Magnetotellurics +=========================== + +SimPEG (Simulation and Parameter Estimation in Geophysics) is a python +package for simulation and gradient based parameter estimation in the +context of geoscience applications. + +simpegMT uses SimPEG as the framework for the forward and inverse +magnetotellurics geophysical problems. + + +Testing simpegMT +================ + +* Master Branch + .. image:: https://travis-ci.org/simpeg/simpegmt.svg?branch=master + :target: https://travis-ci.org/simpeg/simpegmt + :alt: Master Branch + :align: center + + +Project Index & Search +====================== + +* :ref:`genindex` +* :ref:`modindex` +* :ref:`search` + diff --git a/docs/make.bat b/docs/make.bat new file mode 100644 index 00000000..2ac3df69 --- /dev/null +++ b/docs/make.bat @@ -0,0 +1,190 @@ +@ECHO OFF + +REM Command file for Sphinx documentation + +if "%SPHINXBUILD%" == "" ( + set SPHINXBUILD=sphinx-build +) +set BUILDDIR=_build +set ALLSPHINXOPTS=-d %BUILDDIR%/doctrees %SPHINXOPTS% . +set I18NSPHINXOPTS=%SPHINXOPTS% . +if NOT "%PAPER%" == "" ( + set ALLSPHINXOPTS=-D latex_paper_size=%PAPER% %ALLSPHINXOPTS% + set I18NSPHINXOPTS=-D latex_paper_size=%PAPER% %I18NSPHINXOPTS% +) + +if "%1" == "" goto help + +if "%1" == "help" ( + :help + echo.Please use `make ^` where ^ is one of + echo. html to make standalone HTML files + echo. dirhtml to make HTML files named index.html in directories + echo. singlehtml to make a single large HTML file + echo. pickle to make pickle files + echo. json to make JSON files + echo. htmlhelp to make HTML files and a HTML help project + echo. qthelp to make HTML files and a qthelp project + echo. devhelp to make HTML files and a Devhelp project + echo. epub to make an epub + echo. latex to make LaTeX files, you can set PAPER=a4 or PAPER=letter + echo. text to make text files + echo. man to make manual pages + echo. texinfo to make Texinfo files + echo. gettext to make PO message catalogs + echo. changes to make an overview over all changed/added/deprecated items + echo. linkcheck to check all external links for integrity + echo. doctest to run all doctests embedded in the documentation if enabled + goto end +) + +if "%1" == "clean" ( + for /d %%i in (%BUILDDIR%\*) do rmdir /q /s %%i + del /q /s %BUILDDIR%\* + goto end +) + +if "%1" == "html" ( + %SPHINXBUILD% -b html %ALLSPHINXOPTS% %BUILDDIR%/html + if errorlevel 1 exit /b 1 + echo. + echo.Build finished. The HTML pages are in %BUILDDIR%/html. + goto end +) + +if "%1" == "dirhtml" ( + %SPHINXBUILD% -b dirhtml %ALLSPHINXOPTS% %BUILDDIR%/dirhtml + if errorlevel 1 exit /b 1 + echo. + echo.Build finished. The HTML pages are in %BUILDDIR%/dirhtml. + goto end +) + +if "%1" == "singlehtml" ( + %SPHINXBUILD% -b singlehtml %ALLSPHINXOPTS% %BUILDDIR%/singlehtml + if errorlevel 1 exit /b 1 + echo. + echo.Build finished. The HTML pages are in %BUILDDIR%/singlehtml. + goto end +) + +if "%1" == "pickle" ( + %SPHINXBUILD% -b pickle %ALLSPHINXOPTS% %BUILDDIR%/pickle + if errorlevel 1 exit /b 1 + echo. + echo.Build finished; now you can process the pickle files. + goto end +) + +if "%1" == "json" ( + %SPHINXBUILD% -b json %ALLSPHINXOPTS% %BUILDDIR%/json + if errorlevel 1 exit /b 1 + echo. + echo.Build finished; now you can process the JSON files. + goto end +) + +if "%1" == "htmlhelp" ( + %SPHINXBUILD% -b htmlhelp %ALLSPHINXOPTS% %BUILDDIR%/htmlhelp + if errorlevel 1 exit /b 1 + echo. + echo.Build finished; now you can run HTML Help Workshop with the ^ +.hhp project file in %BUILDDIR%/htmlhelp. + goto end +) + +if "%1" == "qthelp" ( + %SPHINXBUILD% -b qthelp %ALLSPHINXOPTS% %BUILDDIR%/qthelp + if errorlevel 1 exit /b 1 + echo. + echo.Build finished; now you can run "qcollectiongenerator" with the ^ +.qhcp project file in %BUILDDIR%/qthelp, like this: + echo.^> qcollectiongenerator %BUILDDIR%\qthelp\SimPEG.qhcp + echo.To view the help file: + echo.^> assistant -collectionFile %BUILDDIR%\qthelp\SimPEG.ghc + goto end +) + +if "%1" == "devhelp" ( + %SPHINXBUILD% -b devhelp %ALLSPHINXOPTS% %BUILDDIR%/devhelp + if errorlevel 1 exit /b 1 + echo. + echo.Build finished. + goto end +) + +if "%1" == "epub" ( + %SPHINXBUILD% -b epub %ALLSPHINXOPTS% %BUILDDIR%/epub + if errorlevel 1 exit /b 1 + echo. + echo.Build finished. The epub file is in %BUILDDIR%/epub. + goto end +) + +if "%1" == "latex" ( + %SPHINXBUILD% -b latex %ALLSPHINXOPTS% %BUILDDIR%/latex + if errorlevel 1 exit /b 1 + echo. + echo.Build finished; the LaTeX files are in %BUILDDIR%/latex. + goto end +) + +if "%1" == "text" ( + %SPHINXBUILD% -b text %ALLSPHINXOPTS% %BUILDDIR%/text + if errorlevel 1 exit /b 1 + echo. + echo.Build finished. The text files are in %BUILDDIR%/text. + goto end +) + +if "%1" == "man" ( + %SPHINXBUILD% -b man %ALLSPHINXOPTS% %BUILDDIR%/man + if errorlevel 1 exit /b 1 + echo. + echo.Build finished. The manual pages are in %BUILDDIR%/man. + goto end +) + +if "%1" == "texinfo" ( + %SPHINXBUILD% -b texinfo %ALLSPHINXOPTS% %BUILDDIR%/texinfo + if errorlevel 1 exit /b 1 + echo. + echo.Build finished. The Texinfo files are in %BUILDDIR%/texinfo. + goto end +) + +if "%1" == "gettext" ( + %SPHINXBUILD% -b gettext %I18NSPHINXOPTS% %BUILDDIR%/locale + if errorlevel 1 exit /b 1 + echo. + echo.Build finished. The message catalogs are in %BUILDDIR%/locale. + goto end +) + +if "%1" == "changes" ( + %SPHINXBUILD% -b changes %ALLSPHINXOPTS% %BUILDDIR%/changes + if errorlevel 1 exit /b 1 + echo. + echo.The overview file is in %BUILDDIR%/changes. + goto end +) + +if "%1" == "linkcheck" ( + %SPHINXBUILD% -b linkcheck %ALLSPHINXOPTS% %BUILDDIR%/linkcheck + if errorlevel 1 exit /b 1 + echo. + echo.Link check complete; look for any errors in the above output ^ +or in %BUILDDIR%/linkcheck/output.txt. + goto end +) + +if "%1" == "doctest" ( + %SPHINXBUILD% -b doctest %ALLSPHINXOPTS% %BUILDDIR%/doctest + if errorlevel 1 exit /b 1 + echo. + echo.Testing of doctests in the sources finished, look at the ^ +results in %BUILDDIR%/doctest/output.txt. + goto end +) + +:end diff --git a/docs/simpeg-logo.png b/docs/simpeg-logo.png new file mode 100644 index 00000000..a0ec740f Binary files /dev/null and b/docs/simpeg-logo.png differ diff --git a/requirements.txt b/requirements.txt new file mode 100644 index 00000000..6b7df820 --- /dev/null +++ b/requirements.txt @@ -0,0 +1,4 @@ +numpy +scipy +ipython +matplotlib diff --git a/setup.py b/setup.py new file mode 100644 index 00000000..95fb3b72 --- /dev/null +++ b/setup.py @@ -0,0 +1,8 @@ +import os +try: + import SimPEG +except ImportError, e: + os.system('git clone https://github.com/simpeg/simpeg.git') + os.system('mv simpeg/SimPEG temp') + os.system('rm -rf simpeg') + os.system('mv temp SimPEG') diff --git a/simpegMT/Base.py b/simpegMT/Base.py new file mode 100644 index 00000000..da7c97fd --- /dev/null +++ b/simpegMT/Base.py @@ -0,0 +1,68 @@ +from SimPEG import Survey, Problem, Utils, np, sp, Solver as SimpegSolver +from scipy.constants import mu_0 + +class BaseEMProblem(Problem.BaseProblem): + + def __init__(self, mesh, **kwargs): + Problem.BaseProblem.__init__(self, mesh, **kwargs) + + solType = None + storeTheseFields = ['e', 'b'] + + surveyPair = Survey.BaseSurvey + dataPair = Survey.Data + + Solver = SimpegSolver + solverOpts = {} + + #################################################### + # Mass Matrices + #################################################### + + @property + def MfMui(self): + #TODO: assuming constant mu + if getattr(self, '_MfMui', None) is None: + self._MfMui = self.mesh.getFaceInnerProduct(1/mu_0) + return self._MfMui + + @property + def Me(self): + if getattr(self, '_Me', None) is None: + self._Me = self.mesh.getEdgeInnerProduct() + return self._Me + + @property + def MeSigma(self): + #TODO: hardcoded to sigma as the model + if getattr(self, '_MeSigma', None) is None: + sigma = self.curTModel + self._MeSigma = self.mesh.getEdgeInnerProduct(sigma) + return self._MeSigma + + @property + def MeSigmaI(self): + #TODO: hardcoded to sigma as the model + if getattr(self, '_MeSigmaI', None) is None: + sigma = self.curTModel + self._MeSigmaI = self.mesh.getEdgeInnerProduct(sigma, invMat=True) + return self._MeSigmaI + + curModel = Utils.dependentProperty('_curModel', None, ['_MeSigma', '_MeSigmaI', '_curTModel', '_curTModelDeriv'], 'Sets the current model, and removes dependent mass matrices.') + + @property + def curTModel(self): + if getattr(self, '_curTModel', None) is None: + self._curTModel = self.mapping.transform(self.curModel) + return self._curTModel + + @property + def curTModelDeriv(self): + if getattr(self, '_curTModelDeriv', None) is None: + self._curTModelDeriv = self.mapping.transformDeriv(self.curModel) + return self._curTModelDeriv + + def fields(self, m): + self.curModel = m + F = self.forward(m, self.getRHS, self.calcFields) + return F diff --git a/simpegMT/FDEM/FDEM.py b/simpegMT/FDEM/FDEM.py new file mode 100644 index 00000000..16571fc3 --- /dev/null +++ b/simpegMT/FDEM/FDEM.py @@ -0,0 +1,296 @@ +from SimPEG import Survey, Problem, Utils, np, sp, Solver as SimpegSolver +from scipy.constants import mu_0 +from SurveyFDEM import SurveyFDEM, FieldsFDEM +from simpegEM.Utils import Sources +from simpegEM.Base import BaseEMProblem + +def omega(freq): + """Change frequency to angular frequency, omega""" + return 2.*np.pi*freq + +class BaseFDEMProblem(BaseEMProblem): + """ + We start by looking at Maxwell's equations in the electric field \\(\\vec{E}\\) and the magnetic flux density \\(\\vec{B}\\): + + .. math:: + + \\nabla \\times \\vec{E} + i \\omega \\vec{B} = 0 \\\\ + \\nabla \\times \\mu^{-1} \\vec{B} - \\sigma \\vec{E} = \\vec{J_s} + + """ + + surveyPair = SurveyFDEM + + def forward(self, m, RHS, CalcFields): + + F = FieldsFDEM(self.mesh, self.survey) + + for freq in self.survey.freqs: + A = self.getA(freq) + rhs = RHS(freq) + solver = self.Solver(A, **self.solverOpts) + sol = solver.solve(rhs) + for fieldType in self.storeTheseFields: + Txs = self.survey.getTransmitters(freq) + F[Txs, fieldType] = CalcFields(sol, freq, fieldType) + + return F + + def Jvec(self, m, v, u=None): + if u is None: + u = self.fields(m) + + self.curModel = m + + Jv = self.dataPair(self.survey) + + for freq in self.survey.freqs: + A = self.getA(freq) + solver = self.Solver(A, **self.solverOpts) + + for tx in self.survey.getTransmitters(freq): + u_tx = u[tx, self.solType] + w = self.getADeriv(freq, u_tx, v) + Ainvw = solver.solve(w) + for rx in tx.rxList: + fAinvw = self.calcFields(Ainvw, freq, rx.projField) + P = lambda v: rx.projectFieldsDeriv(tx, self.mesh, u, v) + + df_dm = self.calcFieldsDeriv(u_tx, freq, rx.projField, v) + if df_dm is None: + Jv[tx, rx] = - P(fAinvw) + else: + Jv[tx, rx] = - P(fAinvw) + P(df_dm) + + return Utils.mkvc(Jv) + + def Jtvec(self, m, v, u=None): + if u is None: + u = self.fields(m) + + self.curModel = m + + # Ensure v is a data object. + if not isinstance(v, self.dataPair): + v = self.dataPair(self.survey, v) + + Jtv = np.zeros(self.mapping.nP) + + for freq in self.survey.freqs: + AT = self.getA(freq).T + solver = self.Solver(AT, **self.solverOpts) + + for tx in self.survey.getTransmitters(freq): + u_tx = u[tx, self.solType] + + for rx in tx.rxList: + PTv = rx.projectFieldsDeriv(tx, self.mesh, u, v[tx, rx], adjoint=True) + fPTv = self.calcFields(PTv, freq, rx.projField, adjoint=True) + + w = solver.solve( fPTv ) + Jtv_rx = - self.getADeriv(freq, u_tx, w, adjoint=True) + + df_dm = self.calcFieldsDeriv(u_tx, freq, rx.projField, PTv, adjoint=True) + + if df_dm is not None: + Jtv_rx += df_dm + + real_or_imag = rx.projComp + if real_or_imag == 'real': + Jtv += Jtv_rx.real + elif real_or_imag == 'imag': + Jtv += - Jtv_rx.real + else: + raise Exception('Must be real or imag') + + return Jtv + + +class ProblemFDEM_e(BaseFDEMProblem): + """ + By eliminating the magnetic flux density using + + .. math:: + + \\vec{B} = \\frac{-1}{i\\omega}\\nabla\\times\\vec{E}, + + we can write Maxwell's equations as a second order system in \\ \\vec{E} \\ only: + + .. math:: + + \\nabla \\times \\mu^{-1} \\nabla \\times \\vec{E} + i \\omega \\sigma \\vec{E} = \\vec{J_s} + + This is the definition of the Forward Problem using the E-formulation of Maxwell's equations. + + + """ + solType = 'e' + + def __init__(self, model, **kwargs): + BaseFDEMProblem.__init__(self, model, **kwargs) + + def getA(self, freq): + """ + :param float freq: Frequency + :rtype: scipy.sparse.csr_matrix + :return: A + """ + mui = self.MfMui + sig = self.MeSigma + C = self.mesh.edgeCurl + + return C.T*mui*C + 1j*omega(freq)*sig + + def getADeriv(self, freq, u, v, adjoint=False): + sig = self.curTModel + dsig_dm = self.curTModelDeriv + dMe_dsig = self.mesh.getEdgeInnerProductDeriv(sig, v=u) + + if adjoint: + return 1j * omega(freq) * ( dsig_dm.T * ( dMe_dsig.T * v ) ) + + return 1j * omega(freq) * ( dMe_dsig * ( dsig_dm * v ) ) + + def getRHS(self, freq): + """ + :param float freq: Frequency + :rtype: numpy.ndarray (nE, nTx) + :return: RHS + """ + Txs = self.survey.getTransmitters(freq) + rhs = range(len(Txs)) + for i, tx in enumerate(Txs): + if tx.txType == 'VMD': + src = Sources.MagneticDipoleVectorPotential + else: + raise NotImplemented('%s txType is not implemented' % tx.txType) + SRCx = src(tx.loc, self.mesh.gridEx, 'x') + SRCy = src(tx.loc, self.mesh.gridEy, 'y') + SRCz = src(tx.loc, self.mesh.gridEz, 'z') + rhs[i] = np.concatenate((SRCx, SRCy, SRCz)) + + a = np.concatenate(rhs).reshape((self.mesh.nE, len(Txs)), order='F') + mui = self.MfMui + C = self.mesh.edgeCurl + + j_s = C.T*mui*C*a + return -1j*omega(freq)*j_s + + def calcFields(self, sol, freq, fieldType, adjoint=False): + e = sol + if fieldType == 'e': + return e + elif fieldType == 'b': + if not adjoint: + b = -(1./(1j*omega(freq))) * ( self.mesh.edgeCurl * e ) + else: + b = -(1./(1j*omega(freq))) * ( self.mesh.edgeCurl.T * e ) + return b + raise NotImplementedError('fieldType "%s" is not implemented.' % fieldType) + + def calcFieldsDeriv(self, sol, freq, fieldType, v, adjoint=False): + e = sol + if fieldType == 'e': + return None + elif fieldType == 'b': + return None + raise NotImplementedError('fieldType "%s" is not implemented.' % fieldType) + + +class ProblemFDEM_b(BaseFDEMProblem): + """ + Solving for b! + """ + solType = 'b' + + def __init__(self, model, **kwargs): + BaseFDEMProblem.__init__(self, model, **kwargs) + + def getA(self, freq): + """ + :param float freq: Frequency + :rtype: scipy.sparse.csr_matrix + :return: A + """ + mui = self.MfMui + sigI = self.MeSigmaI + C = self.mesh.edgeCurl + + return mui*C*sigI*C.T*mui + 1j*omega(freq)*mui + + def getADeriv(self, freq, u, v, adjoint=False): + + mui = self.MfMui + C = self.mesh.edgeCurl + sig = self.curTModel + dsig_dm = self.curTModelDeriv + #TODO: This only works if diagonal (no tensors)... + dMeSigmaI_dI = - self.MeSigmaI**2 + + vec = (C.T*(mui*u)) + dMe_dsig = self.mesh.getEdgeInnerProductDeriv(sig, v=vec) + + if adjoint: + return dsig_dm.T * ( dMe_dsig.T * ( dMeSigmaI_dI.T * ( C.T * ( mui.T * v ) ) ) ) + + return mui * ( C * ( dMeSigmaI_dI * ( dMe_dsig * ( dsig_dm * v ) ) ) ) + + def getRHS(self, freq): + """ + :param float freq: Frequency + :rtype: numpy.ndarray (nE, nTx) + :return: RHS + """ + Txs = self.survey.getTransmitters(freq) + rhs = range(len(Txs)) + for i, tx in enumerate(Txs): + if tx.txType == 'VMD': + src = Sources.MagneticDipoleVectorPotential + else: + raise NotImplemented('%s txType is not implemented' % tx.txType) + SRCx = src(tx.loc, self.mesh.gridEx, 'x') + SRCy = src(tx.loc, self.mesh.gridEy, 'y') + SRCz = src(tx.loc, self.mesh.gridEz, 'z') + rhs[i] = np.concatenate((SRCx, SRCy, SRCz)) + + a = np.concatenate(rhs).reshape((self.mesh.nE, len(Txs)), order='F') + mui = self.MfMui + C = self.mesh.edgeCurl + + b_0 = C*a + return -1j*omega(freq)*mui*b_0 + + def calcFields(self, sol, freq, fieldType, adjoint=False): + b = sol + if fieldType == 'e': + if not adjoint: + e = self.MeSigmaI * ( self.mesh.edgeCurl.T * ( self.MfMui * b ) ) + else: + e = self.MfMui.T * ( self.mesh.edgeCurl * ( self.MeSigmaI.T * b ) ) + return e + elif fieldType == 'b': + return b + raise NotImplementedError('fieldType "%s" is not implemented.' % fieldType) + + def calcFieldsDeriv(self, sol, freq, fieldType, v, adjoint=False): + b = sol + if fieldType == 'e': + sig = self.curTModel + dsig_dm = self.curTModelDeriv + + C = self.mesh.edgeCurl + mui = self.MfMui + + #TODO: This only works if diagonal (no tensors)... + dMeSigmaI_dI = - self.MeSigmaI**2 + + vec = C.T * ( mui * b ) + dMe_dsig = self.mesh.getEdgeInnerProductDeriv(sig, v=vec) + if not adjoint: + return dMeSigmaI_dI * ( dMe_dsig * ( dsig_dm * v ) ) + else: + return dsig_dm.T * ( dMe_dsig.T * ( dMeSigmaI_dI.T * v ) ) + elif fieldType == 'b': + return None + raise NotImplementedError('fieldType "%s" is not implemented.' % fieldType) + diff --git a/simpegMT/FDEM/SurveyFDEM.py b/simpegMT/FDEM/SurveyFDEM.py new file mode 100644 index 00000000..67db6674 --- /dev/null +++ b/simpegMT/FDEM/SurveyFDEM.py @@ -0,0 +1,140 @@ +from SimPEG import Survey, Utils, np, sp + +class RxFDEM(Survey.BaseRx): + + knownRxTypes = { + 'exr':['e', 'Ex', 'real'], + 'eyr':['e', 'Ey', 'real'], + 'ezr':['e', 'Ez', 'real'], + 'exi':['e', 'Ex', 'imag'], + 'eyi':['e', 'Ey', 'imag'], + 'ezi':['e', 'Ez', 'imag'], + + 'bxr':['b', 'Fx', 'real'], + 'byr':['b', 'Fy', 'real'], + 'bzr':['b', 'Fz', 'real'], + 'bxi':['b', 'Fx', 'imag'], + 'byi':['b', 'Fy', 'imag'], + 'bzi':['b', 'Fz', 'imag'], + } + + def __init__(self, locs, rxType): + Survey.BaseRx.__init__(self, locs, rxType) + + @property + def projField(self): + """Field Type projection (e.g. e b ...)""" + return self.knownRxTypes[self.rxType][0] + + @property + def projGLoc(self): + """Grid Location projection (e.g. Ex Fy ...)""" + return self.knownRxTypes[self.rxType][1] + + @property + def projComp(self): + """Component projection (real/imag)""" + return self.knownRxTypes[self.rxType][2] + + def projectFields(self, tx, mesh, u): + P = self.getP(mesh) + u_part_complex = u[tx, self.projField] + # get the real or imag component + real_or_imag = self.projComp + u_part = getattr(u_part_complex, real_or_imag) + return P*u_part + + def projectFieldsDeriv(self, tx, mesh, u, v, adjoint=False): + P = self.getP(mesh) + + if not adjoint: + Pv_complex = P * v + real_or_imag = self.projComp + Pv = getattr(Pv_complex, real_or_imag) + elif adjoint: + Pv_real = P.T * v + + real_or_imag = self.projComp + if real_or_imag == 'imag': + Pv = 1j*Pv_real + elif real_or_imag == 'real': + Pv = Pv_real.astype(complex) + else: + raise NotImplementedError('must be real or imag') + + return Pv + + +class TxFDEM(Survey.BaseTx): + + freq = None #: Frequency (float) + + rxPair = RxFDEM + + knownTxTypes = ['VMD'] + + def __init__(self, loc, txType, freq, rxList): + self.freq = float(freq) + Survey.BaseTx.__init__(self, loc, txType, rxList) + + + +class FieldsFDEM(Survey.Fields): + """Fancy Field Storage for a FDEM survey.""" + knownFields = {'b': 'F', 'e': 'E'} + dtype = complex + + +class SurveyFDEM(Survey.BaseSurvey): + """ + docstring for SurveyFDEM + """ + + txPair = TxFDEM + + def __init__(self, txList, **kwargs): + # Sort these by frequency + self.txList = txList + Survey.BaseSurvey.__init__(self, **kwargs) + + _freqDict = {} + for tx in txList: + if tx.freq not in _freqDict: + _freqDict[tx.freq] = [] + _freqDict[tx.freq] += [tx] + + self._freqDict = _freqDict + self._freqs = sorted([f for f in self._freqDict]) + + @property + def freqs(self): + """Frequencies""" + return self._freqs + + @property + def nFreq(self): + """Number of frequencies""" + return len(self._freqDict) + + @property + def nTxByFreq(self): + if getattr(self, '_nTxByFreq', None) is None: + self._nTxByFreq = {} + for freq in self.freqs: + self._nTxByFreq[freq] = len(self.getTransmitters(freq)) + return self._nTxByFreq + + def getTransmitters(self, freq): + """Returns the transmitters associated with a specific frequency.""" + assert freq in self._freqDict, "The requested frequency is not in this survey." + return self._freqDict[freq] + + def projectFields(self, u): + data = Survey.Data(self) + for tx in self.txList: + for rx in tx.rxList: + data[tx, rx] = rx.projectFields(tx, self.mesh, u) + return data + + def projectFieldsDeriv(self, u): + raise Exception('Use Transmitters to project fields deriv.') diff --git a/simpegMT/FDEM/__init__.py b/simpegMT/FDEM/__init__.py new file mode 100644 index 00000000..70124686 --- /dev/null +++ b/simpegMT/FDEM/__init__.py @@ -0,0 +1,2 @@ +from SurveyFDEM import * +from FDEM import ProblemFDEM_e, ProblemFDEM_b diff --git a/simpegMT/Tests/__init__.py b/simpegMT/Tests/__init__.py new file mode 100644 index 00000000..420388ef --- /dev/null +++ b/simpegMT/Tests/__init__.py @@ -0,0 +1,12 @@ +import os +import glob +import unittest + +if __name__ == '__main__': + test_file_strings = glob.glob('test_*.py') + module_strings = [str[0:len(str)-3] for str in test_file_strings] + suites = [unittest.defaultTestLoader.loadTestsFromName(str) for str + in module_strings] + testSuite = unittest.TestSuite(suites) + + unittest.TextTestRunner(verbosity=2).run(testSuite) diff --git a/simpegMT/Tests/test_FieldsObject.py b/simpegMT/Tests/test_FieldsObject.py new file mode 100644 index 00000000..eac83bad --- /dev/null +++ b/simpegMT/Tests/test_FieldsObject.py @@ -0,0 +1,103 @@ +import unittest +from SimPEG import * +import simpegMT as MT + +class FieldsTest(unittest.TestCase): + + def setUp(self): + mesh = Mesh.TensorMesh([np.ones(n)*5 for n in [10,11,12]],[0,0,-30]) + x = np.linspace(5,10,3) + XYZ = Utils.ndgrid(x,x,np.r_[0.]) + txLoc = np.r_[0,0,0.] + rxList0 = MT.FDEM.RxFDEM(XYZ, 'exi') + Tx0 = MT.FDEM.TxFDEM(txLoc, 'VMD', 3., [rxList0]) + rxList1 = MT.FDEM.RxFDEM(XYZ, 'bxi') + Tx1 = MT.FDEM.TxFDEM(txLoc, 'VMD', 3., [rxList1]) + rxList2 = MT.FDEM.RxFDEM(XYZ, 'bxi') + Tx2 = MT.FDEM.TxFDEM(txLoc, 'VMD', 2., [rxList2]) + rxList3 = MT.FDEM.RxFDEM(XYZ, 'bxi') + Tx3 = MT.FDEM.TxFDEM(txLoc, 'VMD', 2., [rxList3]) + Tx4 = MT.FDEM.TxFDEM(txLoc, 'VMD', 1., [rxList0, rxList1, rxList2, rxList3]) + txList = [Tx0,Tx1,Tx2,Tx3,Tx4] + survey = MT.FDEM.SurveyFDEM(txList) + self.F = MT.FDEM.FieldsFDEM(mesh, survey) + self.Tx0 = Tx0 + self.Tx1 = Tx1 + self.mesh = mesh + self.XYZ = XYZ + + def test_SetGet(self): + F = self.F + for freq in F.survey.freqs: + nFreq = F.survey.nTxByFreq[freq] + Txs = F.survey.getTransmitters(freq) + e = np.random.rand(F.mesh.nE, nFreq) + F[Txs, 'e'] = e + b = np.random.rand(F.mesh.nF, nFreq) + F[Txs, 'b'] = b + if nFreq == 1: + F[Txs, 'b'] = Utils.mkvc(b) + if e.shape[1] == 1: + e, b = Utils.mkvc(e), Utils.mkvc(b) + self.assertTrue(np.all(F[Txs, 'e'] == e)) + self.assertTrue(np.all(F[Txs, 'b'] == b)) + F[Txs] = {'b':b,'e':e} + self.assertTrue(np.all(F[Txs, 'e'] == e)) + self.assertTrue(np.all(F[Txs, 'b'] == b)) + + lastFreq = F[Txs] + self.assertTrue(type(lastFreq) is dict) + self.assertTrue(sorted([k for k in lastFreq]) == ['b','e']) + self.assertTrue(np.all(lastFreq['b'] == b)) + self.assertTrue(np.all(lastFreq['e'] == e)) + + Tx_f3 = F.survey.getTransmitters(3.) + self.assertTrue(F[Tx_f3,'b'].shape == (F.mesh.nF, 2)) + + b = np.random.rand(F.mesh.nF, 2) + Tx_f0 = F.survey.getTransmitters(self.Tx0.freq) + F[Tx_f0,'b'] = b + self.assertTrue(F[self.Tx0]['b'].shape == (F.mesh.nF,)) + self.assertTrue(F[self.Tx0,'b'].shape == (F.mesh.nF,)) + self.assertTrue(np.all(F[self.Tx0,'b'] == b[:,0])) + self.assertTrue(np.all(F[self.Tx1,'b'] == b[:,1])) + + def test_assertions(self): + freq = self.F.survey.freqs[0] + Txs = self.F.survey.getTransmitters(freq) + bWrongSize = np.random.rand(self.F.mesh.nE, self.F.survey.nTxByFreq[freq]) + def fun(): self.F[Txs, 'b'] = bWrongSize + self.assertRaises(ValueError, fun) + def fun(): self.F[-999.] + self.assertRaises(KeyError, fun) + def fun(): self.F['notRight'] + self.assertRaises(KeyError, fun) + def fun(): self.F[Txs,'notThere'] + self.assertRaises(KeyError, fun) + + def test_FieldProjections(self): + F = self.F + for freq in F.survey.freqs: + nFreq = F.survey.nTxByFreq[freq] + Txs = F.survey.getTransmitters(freq) + e = np.random.rand(F.mesh.nE, nFreq) + b = np.random.rand(F.mesh.nF, nFreq) + F[Txs] = {'b':b,'e':e} + + Txs = F.survey.getTransmitters(freq) + for ii, tx in enumerate(Txs): + for jj, rx in enumerate(tx.rxList): + dat = rx.projectFields(tx, self.mesh, F) + self.assertTrue(dat.dtype == float) + fieldType = rx.projField + u = {'b':b[:,ii], 'e': e[:,ii]}[fieldType] + real_or_imag = rx.projComp + u = getattr(u, real_or_imag) + gloc = rx.projGLoc + d = self.mesh.getInterpolationMat(self.XYZ, gloc)*u + self.assertTrue(np.all(dat == d)) + + + +if __name__ == '__main__': + unittest.main() diff --git a/simpegMT/Utils/__init__.py b/simpegMT/Utils/__init__.py new file mode 100644 index 00000000..e69de29b diff --git a/simpegMT/__init__.py b/simpegMT/__init__.py new file mode 100644 index 00000000..a18e6230 --- /dev/null +++ b/simpegMT/__init__.py @@ -0,0 +1,4 @@ +# from EM import * +import Utils +import FDEM +import Base