diff --git a/README.rst b/README.rst index 17c34cc..1ae0637 100644 --- a/README.rst +++ b/README.rst @@ -1,42 +1,42 @@ +============== +Stokes Dynamic +============== +Stokes equation is a linearised form of the Navier-Stokes equations in the limit of small Reynolds number. +This equation, also named creeping flow, is a type of fluid flow where advective inertial forces are small compared with viscous forces. +The dynamic form of the Stokes equation is implemented in this example. -================= -Example Name Here -================= - -Short synopsis of the example. Building the example ==================== -If this not a pure Python example place the instructions here on how to configure and build with CMake:: +The fortran version of the example can be configured and built with CMake:: + + git clone https://github.com/OpenCMISS-Examples/stokes_dynamic.git + mkdir stokes_dynamic-build + cd stokes_dynamic-build + cmake -DOpenCMISSLibs_DIR=~/opencmiss/install/ ../stokes_dynamic/ + make - git clone https://github.com/you/your-example.git # maybe your example is somewhere else like: OpenCMISS-Examples - mkdir build - cmake -DOpenCMISSLibs_DIR=/path/to/opencmisslib/install ../your-example - make # cmake --build . will also work here and is much more platform agnostic. Running the example =================== -Explain how the example is run:: +The example may be executed with the following commands:: - cd build - ./src/fortran/XXXXXXXX + cd ./src/fortran/ + ./stokes_dynamic_fortran -or maybe it is a Python only example:: +The results can be visualised by running `visualiseStokesDynamic.cmgui <./src/fortran/visualiseStokesDynamic.cmgui>`_ with the `Cmgui visualiser `_. - source /path/to/opencmisslibs/install/virtaul_environments/oclibs_venv_pyXY_release/bin/activate - python src/python/XXXXXXXX.py - -where the XY in the path are the Python major and minor versions respectively. Prerequisites ============= -Are there any external sources that are required over and above CMake, OpenCMISS libraries. Sources like meshes which might be stored outside of the example itself. +There are no additional input files required for this example as it is self-contained. + License ======= -A line on the license applicable to this example. \ No newline at end of file +License applicable to this example is described in `LICENSE <./LICENSE>`_. diff --git a/docs/Makefile b/docs/Makefile deleted file mode 100644 index 897a3a5..0000000 --- a/docs/Makefile +++ /dev/null @@ -1,192 +0,0 @@ -# Makefile for Sphinx documentation -# - -# You can set these variables from the command line. -SPHINXOPTS = -SPHINXBUILD = sphinx-build -PAPER = -BUILDDIR = _build - -# User-friendly check for sphinx-build -ifeq ($(shell which $(SPHINXBUILD) >/dev/null 2>&1; echo $$?), 1) -$(error The '$(SPHINXBUILD)' command was not found. Make sure you have Sphinx installed, then set the SPHINXBUILD environment variable to point to the full path of the '$(SPHINXBUILD)' executable. Alternatively you can add the directory with the executable to your PATH. If you don't have Sphinx installed, grab it from http://sphinx-doc.org/) -endif - -# 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 coverage 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 " applehelp to make an Apple Help Book" - @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 " latexpdfja to make LaTeX files and run them through platex/dvipdfmx" - @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 " xml to make Docutils-native XML files" - @echo " pseudoxml to make pseudoxml-XML files for display purposes" - @echo " linkcheck to check all external links for integrity" - @echo " doctest to run all doctests embedded in the documentation (if enabled)" - @echo " coverage to run coverage check of 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/XXXXXXXX.qhcp" - @echo "To view the help file:" - @echo "# assistant -collectionFile $(BUILDDIR)/qthelp/XXXXXXXX.qhc" - -applehelp: - $(SPHINXBUILD) -b applehelp $(ALLSPHINXOPTS) $(BUILDDIR)/applehelp - @echo - @echo "Build finished. The help book is in $(BUILDDIR)/applehelp." - @echo "N.B. You won't be able to view it unless you put it in" \ - "~/Library/Documentation/Help or install it in your application" \ - "bundle." - -devhelp: - $(SPHINXBUILD) -b devhelp $(ALLSPHINXOPTS) $(BUILDDIR)/devhelp - @echo - @echo "Build finished." - @echo "To view the help file:" - @echo "# mkdir -p $$HOME/.local/share/devhelp/XXXXXXXX" - @echo "# ln -s $(BUILDDIR)/devhelp $$HOME/.local/share/devhelp/XXXXXXXX" - @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." - -latexpdfja: - $(SPHINXBUILD) -b latex $(ALLSPHINXOPTS) $(BUILDDIR)/latex - @echo "Running LaTeX files through platex and dvipdfmx..." - $(MAKE) -C $(BUILDDIR)/latex all-pdf-ja - @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." - -coverage: - $(SPHINXBUILD) -b coverage $(ALLSPHINXOPTS) $(BUILDDIR)/coverage - @echo "Testing of coverage in the sources finished, look at the " \ - "results in $(BUILDDIR)/coverage/python.txt." - -xml: - $(SPHINXBUILD) -b xml $(ALLSPHINXOPTS) $(BUILDDIR)/xml - @echo - @echo "Build finished. The XML files are in $(BUILDDIR)/xml." - -pseudoxml: - $(SPHINXBUILD) -b pseudoxml $(ALLSPHINXOPTS) $(BUILDDIR)/pseudoxml - @echo - @echo "Build finished. The pseudo-XML files are in $(BUILDDIR)/pseudoxml." diff --git a/docs/conf.py b/docs/conf.py deleted file mode 100644 index bf40d6d..0000000 --- a/docs/conf.py +++ /dev/null @@ -1,286 +0,0 @@ -# -*- coding: utf-8 -*- -# -# XXXXXXXX documentation build configuration file, created by -# sphinx-quickstart on Tue Sep 26 16:27:22 2017. -# -# 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 -import os -import shlex - -# 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('.')) - -# -- 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.mathjax', -] - -# Add any paths that contain templates here, relative to this directory. -templates_path = ['_templates'] - -# The suffix(es) of source filenames. -# You can specify multiple suffix as a list of string: -# source_suffix = ['.rst', '.md'] -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'XXXXXXXX' -copyright = u'2017, YYYYYYYY' -author = u'YYYYYYYY' - -# 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 = 'X.Y.Z' -# The full version, including alpha/beta/rc tags. -release = 'X.Y.Z' - -# The language for content autogenerated by Sphinx. Refer to documentation -# for a list of supported languages. -# -# This is also used if you do content translation via gettext catalogs. -# Usually you set "language" from the command line for these cases. -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 = [] - -# If true, keep warnings as "system message" paragraphs in the built documents. -#keep_warnings = False - -# If true, `todo` and `todoList` produce output, else they produce nothing. -todo_include_todos = False - - -# -- 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 = 'alabaster' - -# 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'] - -# Add any extra paths that contain custom files (such as robots.txt or -# .htaccess) here, relative to this directory. These files are copied -# directly to the root of the documentation. -#html_extra_path = [] - -# 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 - -# Language to be used for generating the HTML full-text search index. -# Sphinx supports the following languages: -# 'da', 'de', 'en', 'es', 'fi', 'fr', 'hu', 'it', 'ja' -# 'nl', 'no', 'pt', 'ro', 'ru', 'sv', 'tr' -#html_search_language = 'en' - -# A dictionary with options for the search language support, empty by default. -# Now only 'ja' uses this config value -#html_search_options = {'type': 'default'} - -# The name of a javascript file (relative to the configuration directory) that -# implements a search results scorer. If empty, the default will be used. -#html_search_scorer = 'scorer.js' - -# Output file base name for HTML help builder. -htmlhelp_basename = 'XXXXXXXXdoc' - -# -- 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': '', - -# Latex figure (float) alignment -#'figure_align': 'htbp', -} - -# Grouping the document tree into LaTeX files. List of tuples -# (source start file, target name, title, -# author, documentclass [howto, manual, or own class]). -latex_documents = [ - (master_doc, 'XXXXXXXX.tex', u'XXXXXXXX Documentation', - u'YYYYYYYY', '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 = [ - (master_doc, 'xxxxxxxx', u'XXXXXXXX Documentation', - [author], 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 = [ - (master_doc, 'XXXXXXXX', u'XXXXXXXX Documentation', - author, 'XXXXXXXX', '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' - -# If true, do not generate a @detailmenu in the "Top" node's menu. -#texinfo_no_detailmenu = False diff --git a/docs/index.rst b/docs/index.rst deleted file mode 100644 index c6e8d46..0000000 --- a/docs/index.rst +++ /dev/null @@ -1,22 +0,0 @@ -.. XXXXXXXX documentation master file, created by - sphinx-quickstart on Tue Sep 26 16:27:22 2017. - You can adapt this file completely to your liking, but it should at least - contain the root `toctree` directive. - -Welcome to XXXXXXXX's documentation! -==================================== - -Contents: - -.. toctree:: - :maxdepth: 2 - - - -Indices and tables -================== - -* :ref:`genindex` -* :ref:`modindex` -* :ref:`search` - diff --git a/docs/make.bat b/docs/make.bat deleted file mode 100644 index 9744509..0000000 --- a/docs/make.bat +++ /dev/null @@ -1,263 +0,0 @@ -@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. xml to make Docutils-native XML files - echo. pseudoxml to make pseudoxml-XML files for display purposes - echo. linkcheck to check all external links for integrity - echo. doctest to run all doctests embedded in the documentation if enabled - echo. coverage to run coverage check of 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 -) - - -REM Check if sphinx-build is available and fallback to Python version if any -%SPHINXBUILD% 2> nul -if errorlevel 9009 goto sphinx_python -goto sphinx_ok - -:sphinx_python - -set SPHINXBUILD=python -m sphinx.__init__ -%SPHINXBUILD% 2> nul -if errorlevel 9009 ( - echo. - echo.The 'sphinx-build' command was not found. Make sure you have Sphinx - echo.installed, then set the SPHINXBUILD environment variable to point - echo.to the full path of the 'sphinx-build' executable. Alternatively you - echo.may add the Sphinx directory to PATH. - echo. - echo.If you don't have Sphinx installed, grab it from - echo.http://sphinx-doc.org/ - exit /b 1 -) - -:sphinx_ok - - -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\XXXXXXXX.qhcp - echo.To view the help file: - echo.^> assistant -collectionFile %BUILDDIR%\qthelp\XXXXXXXX.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" == "latexpdf" ( - %SPHINXBUILD% -b latex %ALLSPHINXOPTS% %BUILDDIR%/latex - cd %BUILDDIR%/latex - make all-pdf - cd %~dp0 - echo. - echo.Build finished; the PDF files are in %BUILDDIR%/latex. - goto end -) - -if "%1" == "latexpdfja" ( - %SPHINXBUILD% -b latex %ALLSPHINXOPTS% %BUILDDIR%/latex - cd %BUILDDIR%/latex - make all-pdf-ja - cd %~dp0 - echo. - echo.Build finished; the PDF 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 -) - -if "%1" == "coverage" ( - %SPHINXBUILD% -b coverage %ALLSPHINXOPTS% %BUILDDIR%/coverage - if errorlevel 1 exit /b 1 - echo. - echo.Testing of coverage in the sources finished, look at the ^ -results in %BUILDDIR%/coverage/python.txt. - goto end -) - -if "%1" == "xml" ( - %SPHINXBUILD% -b xml %ALLSPHINXOPTS% %BUILDDIR%/xml - if errorlevel 1 exit /b 1 - echo. - echo.Build finished. The XML files are in %BUILDDIR%/xml. - goto end -) - -if "%1" == "pseudoxml" ( - %SPHINXBUILD% -b pseudoxml %ALLSPHINXOPTS% %BUILDDIR%/pseudoxml - if errorlevel 1 exit /b 1 - echo. - echo.Build finished. The pseudo-XML files are in %BUILDDIR%/pseudoxml. - goto end -) - -:end diff --git a/src/fortran/expected_results/output/MainTime_0.part0.exelem b/src/fortran/expected_results/output/MainTime_0.part0.exelem new file mode 100644 index 0000000..6adc3cc --- /dev/null +++ b/src/fortran/expected_results/output/MainTime_0.part0.exelem @@ -0,0 +1,201 @@ + Group name: StokesRegion + Shape. Dimension=2 + #Scale factor sets= 1 + l.Lagrange*l.Lagrange, #Scale factors=4 + #Nodes= 4 + #Fields=4 + 1) Coordinate, coordinate, rectangular cartesian, #Components=2 + x. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + y. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2) U, field, rectangular cartesian, #Components=3 + 1. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3) del U_del n, field, rectangular cartesian, #Components=3 + 1. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4) Materials, field, rectangular cartesian, #Components=2 + 1. constant*constant, no modify, grid based. + #xi1=0, #xi2=0 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + Element: 1 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 1 2 5 6 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 2 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 2 3 6 7 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 3 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 3 4 7 8 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 4 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 5 6 9 10 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 5 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 6 7 10 11 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 6 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 7 8 11 12 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 7 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 9 10 13 14 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 8 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 10 11 14 15 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 9 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 11 12 15 16 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 diff --git a/src/fortran/expected_results/output/MainTime_0.part0.exnode b/src/fortran/expected_results/output/MainTime_0.part0.exnode new file mode 100644 index 0000000..32ccb3a --- /dev/null +++ b/src/fortran/expected_results/output/MainTime_0.part0.exnode @@ -0,0 +1,192 @@ + Group name: StokesRegion + #Fields=4 + 1) Coordinate, coordinate, rectangular cartesian, #Components=2 + x. Value index= 1, #Derivatives= 0 + y. Value index= 2, #Derivatives= 0 + 2) U, field, rectangular cartesian, #Components=3 + 1. Value index= 3, #Derivatives= 0 + 2. Value index= 4, #Derivatives= 0 + 3. Value index= 5, #Derivatives= 0 + 3) del U_del n, field, rectangular cartesian, #Components=3 + 1. Value index= 6, #Derivatives= 0 + 2. Value index= 7, #Derivatives= 0 + 3. Value index= 8, #Derivatives= 0 + 4) Materials, field, rectangular cartesian, #Components=2 + 1. Value index= 9, #Derivatives= 0 + 2. Value index= 10, #Derivatives= 0 + Node: 1 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 2 + 3.3333333333333331E-01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 3 + 6.6666666666666663E-01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 4 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 5 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 6 + 3.3333333333333331E-01 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 7 + 6.6666666666666663E-01 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 8 + 1.0000000000000000E+00 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 9 + 0.0000000000000000E+00 + 2.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 10 + 3.3333333333333331E-01 + 2.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 11 + 6.6666666666666663E-01 + 2.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 12 + 1.0000000000000000E+00 + 2.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 13 + 0.0000000000000000E+00 + 3.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 14 + 3.3333333333333331E-01 + 3.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 15 + 6.6666666666666663E-01 + 3.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 16 + 1.0000000000000000E+00 + 3.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 diff --git a/src/fortran/expected_results/output/MainTime_1.part0.exelem b/src/fortran/expected_results/output/MainTime_1.part0.exelem new file mode 100644 index 0000000..6adc3cc --- /dev/null +++ b/src/fortran/expected_results/output/MainTime_1.part0.exelem @@ -0,0 +1,201 @@ + Group name: StokesRegion + Shape. Dimension=2 + #Scale factor sets= 1 + l.Lagrange*l.Lagrange, #Scale factors=4 + #Nodes= 4 + #Fields=4 + 1) Coordinate, coordinate, rectangular cartesian, #Components=2 + x. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + y. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2) U, field, rectangular cartesian, #Components=3 + 1. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3) del U_del n, field, rectangular cartesian, #Components=3 + 1. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4) Materials, field, rectangular cartesian, #Components=2 + 1. constant*constant, no modify, grid based. + #xi1=0, #xi2=0 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + Element: 1 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 1 2 5 6 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 2 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 2 3 6 7 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 3 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 3 4 7 8 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 4 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 5 6 9 10 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 5 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 6 7 10 11 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 6 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 7 8 11 12 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 7 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 9 10 13 14 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 8 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 10 11 14 15 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 9 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 11 12 15 16 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 diff --git a/src/fortran/expected_results/output/MainTime_1.part0.exnode b/src/fortran/expected_results/output/MainTime_1.part0.exnode new file mode 100644 index 0000000..6a9906a --- /dev/null +++ b/src/fortran/expected_results/output/MainTime_1.part0.exnode @@ -0,0 +1,192 @@ + Group name: StokesRegion + #Fields=4 + 1) Coordinate, coordinate, rectangular cartesian, #Components=2 + x. Value index= 1, #Derivatives= 0 + y. Value index= 2, #Derivatives= 0 + 2) U, field, rectangular cartesian, #Components=3 + 1. Value index= 3, #Derivatives= 0 + 2. Value index= 4, #Derivatives= 0 + 3. Value index= 5, #Derivatives= 0 + 3) del U_del n, field, rectangular cartesian, #Components=3 + 1. Value index= 6, #Derivatives= 0 + 2. Value index= 7, #Derivatives= 0 + 3. Value index= 8, #Derivatives= 0 + 4) Materials, field, rectangular cartesian, #Components=2 + 1. Value index= 9, #Derivatives= 0 + 2. Value index= 10, #Derivatives= 0 + Node: 1 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.1812122013384210E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 2 + 3.3333333333333331E-01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 2.9846130637575932E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 3 + 6.6666666666666663E-01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 2.9846130637587621E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 4 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.1812122013397113E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 5 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 2.6286825726128612E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 6 + 3.3333333333333331E-01 + 1.0000000000000000E+00 + 4.0196780935444903E-11 + 9.9999999999176392E-01 + 2.1200746481173198E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 7 + 6.6666666666666663E-01 + 1.0000000000000000E+00 + 4.1106091599225940E-11 + 1.0000000000059759E+00 + 2.1200746481194802E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 8 + 1.0000000000000000E+00 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 2.6286825726150116E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 9 + 0.0000000000000000E+00 + 2.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 5.5263166816708642E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 10 + 3.3333333333333331E-01 + 2.0000000000000000E+00 + -7.8633580547472011E-11 + 1.0000000000002784E+00 + 7.8366250373396911E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 11 + 6.6666666666666663E-01 + 2.0000000000000000E+00 + -8.0591616372742983E-11 + 1.0000000000008473E+00 + 7.8366250372513431E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 12 + 1.0000000000000000E+00 + 2.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 5.5263166815782201E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 13 + 0.0000000000000000E+00 + 3.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 5.8443881702088829E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 14 + 3.3333333333333331E-01 + 3.0000000000000000E+00 + 4.5552450068894528E-11 + 1.0000000000138978E+00 + 1.6892339923448090E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 15 + 6.6666666666666663E-01 + 3.0000000000000000E+00 + 4.5218765955018175E-11 + 9.9999999998926870E-01 + 1.6892339923617283E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 16 + 1.0000000000000000E+00 + 3.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 5.8443881702222544E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 diff --git a/src/fortran/expected_results/output/MainTime_2.part0.exelem b/src/fortran/expected_results/output/MainTime_2.part0.exelem new file mode 100644 index 0000000..6adc3cc --- /dev/null +++ b/src/fortran/expected_results/output/MainTime_2.part0.exelem @@ -0,0 +1,201 @@ + Group name: StokesRegion + Shape. Dimension=2 + #Scale factor sets= 1 + l.Lagrange*l.Lagrange, #Scale factors=4 + #Nodes= 4 + #Fields=4 + 1) Coordinate, coordinate, rectangular cartesian, #Components=2 + x. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + y. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2) U, field, rectangular cartesian, #Components=3 + 1. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3) del U_del n, field, rectangular cartesian, #Components=3 + 1. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4) Materials, field, rectangular cartesian, #Components=2 + 1. constant*constant, no modify, grid based. + #xi1=0, #xi2=0 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + Element: 1 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 1 2 5 6 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 2 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 2 3 6 7 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 3 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 3 4 7 8 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 4 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 5 6 9 10 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 5 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 6 7 10 11 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 6 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 7 8 11 12 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 7 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 9 10 13 14 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 8 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 10 11 14 15 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 9 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 11 12 15 16 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 diff --git a/src/fortran/expected_results/output/MainTime_2.part0.exnode b/src/fortran/expected_results/output/MainTime_2.part0.exnode new file mode 100644 index 0000000..0103c22 --- /dev/null +++ b/src/fortran/expected_results/output/MainTime_2.part0.exnode @@ -0,0 +1,192 @@ + Group name: StokesRegion + #Fields=4 + 1) Coordinate, coordinate, rectangular cartesian, #Components=2 + x. Value index= 1, #Derivatives= 0 + y. Value index= 2, #Derivatives= 0 + 2) U, field, rectangular cartesian, #Components=3 + 1. Value index= 3, #Derivatives= 0 + 2. Value index= 4, #Derivatives= 0 + 3. Value index= 5, #Derivatives= 0 + 3) del U_del n, field, rectangular cartesian, #Components=3 + 1. Value index= 6, #Derivatives= 0 + 2. Value index= 7, #Derivatives= 0 + 3. Value index= 8, #Derivatives= 0 + 4) Materials, field, rectangular cartesian, #Components=2 + 1. Value index= 9, #Derivatives= 0 + 2. Value index= 10, #Derivatives= 0 + Node: 1 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0901513302457076E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 2 + 3.3333333333333331E-01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 2.7592262631258805E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 3 + 6.6666666666666663E-01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 2.7592262631148039E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 4 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0901513302380593E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 5 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 2.4126556016713060E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 6 + 3.3333333333333331E-01 + 1.0000000000000000E+00 + 1.1061237735346829E-11 + 1.0000000000010958E+00 + 1.9402123505118976E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 7 + 6.6666666666666663E-01 + 1.0000000000000000E+00 + 1.1115054301577958E-11 + 9.9999999999955824E-01 + 1.9402123504877174E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 8 + 1.0000000000000000E+00 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 2.4126556016474076E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 9 + 0.0000000000000000E+00 + 2.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 4.9713204783704992E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 10 + 3.3333333333333331E-01 + 2.0000000000000000E+00 + -4.0673088599273896E-12 + 1.0000000000029283E+00 + 7.1783011959370659E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 11 + 6.6666666666666663E-01 + 2.0000000000000000E+00 + -3.7386164374112150E-12 + 9.9999999999631439E-01 + 7.1783011960619092E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 12 + 1.0000000000000000E+00 + 2.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 4.9713204784162279E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 13 + 0.0000000000000000E+00 + 3.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 5.6091042225424621E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 14 + 3.3333333333333331E-01 + 3.0000000000000000E+00 + -1.0758781928043634E-11 + 9.9999999999433287E-01 + 1.5056138636990473E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 15 + 6.6666666666666663E-01 + 3.0000000000000000E+00 + -1.1788039474953334E-11 + 1.0000000000063791E+00 + 1.5056138639048604E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 16 + 1.0000000000000000E+00 + 3.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 5.6091042227001591E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 diff --git a/src/fortran/expected_results/output/MainTime_3.part0.exelem b/src/fortran/expected_results/output/MainTime_3.part0.exelem new file mode 100644 index 0000000..6adc3cc --- /dev/null +++ b/src/fortran/expected_results/output/MainTime_3.part0.exelem @@ -0,0 +1,201 @@ + Group name: StokesRegion + Shape. Dimension=2 + #Scale factor sets= 1 + l.Lagrange*l.Lagrange, #Scale factors=4 + #Nodes= 4 + #Fields=4 + 1) Coordinate, coordinate, rectangular cartesian, #Components=2 + x. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + y. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2) U, field, rectangular cartesian, #Components=3 + 1. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3) del U_del n, field, rectangular cartesian, #Components=3 + 1. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4) Materials, field, rectangular cartesian, #Components=2 + 1. constant*constant, no modify, grid based. + #xi1=0, #xi2=0 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + Element: 1 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 1 2 5 6 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 2 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 2 3 6 7 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 3 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 3 4 7 8 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 4 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 5 6 9 10 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 5 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 6 7 10 11 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 6 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 7 8 11 12 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 7 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 9 10 13 14 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 8 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 10 11 14 15 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 9 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 11 12 15 16 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 diff --git a/src/fortran/expected_results/output/MainTime_3.part0.exnode b/src/fortran/expected_results/output/MainTime_3.part0.exnode new file mode 100644 index 0000000..0103c22 --- /dev/null +++ b/src/fortran/expected_results/output/MainTime_3.part0.exnode @@ -0,0 +1,192 @@ + Group name: StokesRegion + #Fields=4 + 1) Coordinate, coordinate, rectangular cartesian, #Components=2 + x. Value index= 1, #Derivatives= 0 + y. Value index= 2, #Derivatives= 0 + 2) U, field, rectangular cartesian, #Components=3 + 1. Value index= 3, #Derivatives= 0 + 2. Value index= 4, #Derivatives= 0 + 3. Value index= 5, #Derivatives= 0 + 3) del U_del n, field, rectangular cartesian, #Components=3 + 1. Value index= 6, #Derivatives= 0 + 2. Value index= 7, #Derivatives= 0 + 3. Value index= 8, #Derivatives= 0 + 4) Materials, field, rectangular cartesian, #Components=2 + 1. Value index= 9, #Derivatives= 0 + 2. Value index= 10, #Derivatives= 0 + Node: 1 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0901513302457076E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 2 + 3.3333333333333331E-01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 2.7592262631258805E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 3 + 6.6666666666666663E-01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 2.7592262631148039E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 4 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0901513302380593E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 5 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 2.4126556016713060E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 6 + 3.3333333333333331E-01 + 1.0000000000000000E+00 + 1.1061237735346829E-11 + 1.0000000000010958E+00 + 1.9402123505118976E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 7 + 6.6666666666666663E-01 + 1.0000000000000000E+00 + 1.1115054301577958E-11 + 9.9999999999955824E-01 + 1.9402123504877174E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 8 + 1.0000000000000000E+00 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 2.4126556016474076E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 9 + 0.0000000000000000E+00 + 2.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 4.9713204783704992E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 10 + 3.3333333333333331E-01 + 2.0000000000000000E+00 + -4.0673088599273896E-12 + 1.0000000000029283E+00 + 7.1783011959370659E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 11 + 6.6666666666666663E-01 + 2.0000000000000000E+00 + -3.7386164374112150E-12 + 9.9999999999631439E-01 + 7.1783011960619092E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 12 + 1.0000000000000000E+00 + 2.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 4.9713204784162279E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 13 + 0.0000000000000000E+00 + 3.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 5.6091042225424621E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 14 + 3.3333333333333331E-01 + 3.0000000000000000E+00 + -1.0758781928043634E-11 + 9.9999999999433287E-01 + 1.5056138636990473E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 15 + 6.6666666666666663E-01 + 3.0000000000000000E+00 + -1.1788039474953334E-11 + 1.0000000000063791E+00 + 1.5056138639048604E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 16 + 1.0000000000000000E+00 + 3.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 5.6091042227001591E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 diff --git a/src/fortran/expected_results/output/MainTime_4.part0.exelem b/src/fortran/expected_results/output/MainTime_4.part0.exelem new file mode 100644 index 0000000..6adc3cc --- /dev/null +++ b/src/fortran/expected_results/output/MainTime_4.part0.exelem @@ -0,0 +1,201 @@ + Group name: StokesRegion + Shape. Dimension=2 + #Scale factor sets= 1 + l.Lagrange*l.Lagrange, #Scale factors=4 + #Nodes= 4 + #Fields=4 + 1) Coordinate, coordinate, rectangular cartesian, #Components=2 + x. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + y. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2) U, field, rectangular cartesian, #Components=3 + 1. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3) del U_del n, field, rectangular cartesian, #Components=3 + 1. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4) Materials, field, rectangular cartesian, #Components=2 + 1. constant*constant, no modify, grid based. + #xi1=0, #xi2=0 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 3. #Values=1 + Value indices: 1 + Scale factor indices: -1 + 4. #Values=1 + Value indices: 1 + Scale factor indices: -1 + Element: 1 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 1 2 5 6 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 2 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 2 3 6 7 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 3 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 3 4 7 8 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 4 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 5 6 9 10 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 5 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 6 7 10 11 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 6 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 7 8 11 12 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 7 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 9 10 13 14 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 8 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 10 11 14 15 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 9 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 11 12 15 16 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 diff --git a/src/fortran/expected_results/output/MainTime_4.part0.exnode b/src/fortran/expected_results/output/MainTime_4.part0.exnode new file mode 100644 index 0000000..0103c22 --- /dev/null +++ b/src/fortran/expected_results/output/MainTime_4.part0.exnode @@ -0,0 +1,192 @@ + Group name: StokesRegion + #Fields=4 + 1) Coordinate, coordinate, rectangular cartesian, #Components=2 + x. Value index= 1, #Derivatives= 0 + y. Value index= 2, #Derivatives= 0 + 2) U, field, rectangular cartesian, #Components=3 + 1. Value index= 3, #Derivatives= 0 + 2. Value index= 4, #Derivatives= 0 + 3. Value index= 5, #Derivatives= 0 + 3) del U_del n, field, rectangular cartesian, #Components=3 + 1. Value index= 6, #Derivatives= 0 + 2. Value index= 7, #Derivatives= 0 + 3. Value index= 8, #Derivatives= 0 + 4) Materials, field, rectangular cartesian, #Components=2 + 1. Value index= 9, #Derivatives= 0 + 2. Value index= 10, #Derivatives= 0 + Node: 1 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0901513302457076E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 2 + 3.3333333333333331E-01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 2.7592262631258805E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 3 + 6.6666666666666663E-01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 2.7592262631148039E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 4 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0901513302380593E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 5 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 2.4126556016713060E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 6 + 3.3333333333333331E-01 + 1.0000000000000000E+00 + 1.1061237735346829E-11 + 1.0000000000010958E+00 + 1.9402123505118976E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 7 + 6.6666666666666663E-01 + 1.0000000000000000E+00 + 1.1115054301577958E-11 + 9.9999999999955824E-01 + 1.9402123504877174E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 8 + 1.0000000000000000E+00 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 2.4126556016474076E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 9 + 0.0000000000000000E+00 + 2.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 4.9713204783704992E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 10 + 3.3333333333333331E-01 + 2.0000000000000000E+00 + -4.0673088599273896E-12 + 1.0000000000029283E+00 + 7.1783011959370659E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 11 + 6.6666666666666663E-01 + 2.0000000000000000E+00 + -3.7386164374112150E-12 + 9.9999999999631439E-01 + 7.1783011960619092E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 12 + 1.0000000000000000E+00 + 2.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 4.9713204784162279E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 13 + 0.0000000000000000E+00 + 3.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 5.6091042225424621E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 14 + 3.3333333333333331E-01 + 3.0000000000000000E+00 + -1.0758781928043634E-11 + 9.9999999999433287E-01 + 1.5056138636990473E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 15 + 6.6666666666666663E-01 + 3.0000000000000000E+00 + -1.1788039474953334E-11 + 1.0000000000063791E+00 + 1.5056138639048604E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 16 + 1.0000000000000000E+00 + 3.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 5.6091042227001591E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 diff --git a/src/fortran/expected_results/results.txt b/src/fortran/expected_results/results.txt deleted file mode 100644 index 84c0069..0000000 --- a/src/fortran/expected_results/results.txt +++ /dev/null @@ -1 +0,0 @@ -working. diff --git a/src/fortran/input/CMHEART.inp b/src/fortran/input/CMHEART.inp deleted file mode 100644 index 5e4cd14..0000000 --- a/src/fortran/input/CMHEART.inp +++ /dev/null @@ -1,3 +0,0 @@ -'../../../TOOLS/NavierStokesMeshes/HEX-M2-V2-P1_FE.aux' -'../../../TOOLS/NavierStokesMeshes/HEX-M2-V2-P1_FE.C' -'../../../TOOLS/NavierStokesMeshes/HEX-M2-V2-P1_FE.M' diff --git a/src/fortran/stokes_dynamic.F90 b/src/fortran/stokes_dynamic.F90 index a5bad45..1f7f565 100644 --- a/src/fortran/stokes_dynamic.F90 +++ b/src/fortran/stokes_dynamic.F90 @@ -1,36 +1,24 @@ -PROGRAM STOKESDYNAMICEXAMPLE +PROGRAM stokes_dynamic ! !================================================================================================================================ ! - !PROGRAM LIBRARIES - USE OpenCMISS USE OpenCMISS_Iron #ifndef NOMPIMOD USE MPI #endif - - ! - !================================================================================================================================ - ! - - !PROGRAM VARIABLES AND TYPES - IMPLICIT NONE - #ifdef NOMPIMOD #include "mpif.h" #endif !Test program parameters - REAL(CMISSRP), PARAMETER :: HEIGHT=1.0_CMISSRP + REAL(CMISSRP), PARAMETER :: HEIGHT=3.0_CMISSRP REAL(CMISSRP), PARAMETER :: WIDTH=1.0_CMISSRP REAL(CMISSRP), PARAMETER :: LENGTH=1.0_CMISSRP - !Test program parameters - INTEGER(CMISSIntg), PARAMETER :: CoordinateSystemUserNumber=1 INTEGER(CMISSIntg), PARAMETER :: RegionUserNumber=2 INTEGER(CMISSIntg), PARAMETER :: MeshUserNumber=3 @@ -42,6 +30,7 @@ PROGRAM STOKESDYNAMICEXAMPLE INTEGER(CMISSIntg), PARAMETER :: EquationsSetUserNumberStokes=9 INTEGER(CMISSIntg), PARAMETER :: ProblemUserNumber=10 INTEGER(CMISSIntg), PARAMETER :: GeneratedMeshUserNumber=11 + INTEGER(CMISSIntg), PARAMETER :: AnalyticFieldUserNumber=12 INTEGER(CMISSIntg), PARAMETER :: DomainUserNumber=2 INTEGER(CMISSIntg), PARAMETER :: SolverStokesUserNumber=1 @@ -77,7 +66,7 @@ PROGRAM STOKESDYNAMICEXAMPLE INTEGER(CMISSIntg) :: TOTAL_NUMBER_OF_ELEMENTS INTEGER(CMISSIntg) :: MAXIMUM_ITERATIONS INTEGER(CMISSIntg) :: RESTART_VALUE -! INTEGER(CMISSIntg) :: MPI_IERROR + ! INTEGER(CMISSIntg) :: MPI_IERROR INTEGER(CMISSIntg) :: NUMBER_OF_FIXED_WALL_NODES_STOKES INTEGER(CMISSIntg) :: NUMBER_OF_INLET_WALL_NODES_STOKES @@ -117,55 +106,38 @@ PROGRAM STOKESDYNAMICEXAMPLE !CMISS variables - !Regions TYPE(cmfe_RegionType) :: Region TYPE(cmfe_RegionType) :: WorldRegion - !Coordinate systems TYPE(cmfe_CoordinateSystemType) :: CoordinateSystem TYPE(cmfe_CoordinateSystemType) :: WorldCoordinateSystem - !Basis TYPE(cmfe_BasisType) :: BasisGeometry TYPE(cmfe_BasisType) :: BasisVelocity TYPE(cmfe_BasisType) :: BasisPressure - !Nodes TYPE(cmfe_NodesType) :: Nodes - !Elements TYPE(cmfe_MeshElementsType) :: MeshElementsSpace TYPE(cmfe_MeshElementsType) :: MeshElementsVelocity TYPE(cmfe_MeshElementsType) :: MeshElementsPressure - !Meshes TYPE(cmfe_MeshType) :: Mesh TYPE(cmfe_GeneratedMeshType) :: GeneratedMesh - !Decompositions TYPE(cmfe_DecompositionType) :: Decomposition - !Fields TYPE(cmfe_FieldsType) :: Fields - !Field types - TYPE(cmfe_FieldType) :: GeometricField + TYPE(cmfe_FieldType) :: GeometricField,AnalyticField TYPE(cmfe_FieldType) :: EquationsSetField TYPE(cmfe_FieldType) :: DependentFieldStokes TYPE(cmfe_FieldType) :: MaterialsFieldStokes - !Boundary conditions TYPE(cmfe_BoundaryConditionsType) :: BoundaryConditionsStokes - !Equations sets TYPE(cmfe_EquationsSetType) :: EquationsSetStokes - !Equations TYPE(cmfe_EquationsType) :: EquationsStokes - !Problems TYPE(cmfe_ProblemType) :: Problem - !Control loops TYPE(cmfe_ControlLoopType) :: ControlLoop - !Solvers TYPE(cmfe_SolverType) :: DynamicSolverStokes TYPE(cmfe_SolverType) :: LinearSolverStokes - !Solver equations TYPE(cmfe_SolverEquationsType) :: SolverEquationsStokes !Generic CMISS variables INTEGER(CMISSIntg) :: NumberOfComputationalNodes,ComputationalNodeNumber,BoundaryNodeDomain - INTEGER(CMISSIntg) :: EquationsSetIndex - INTEGER(CMISSIntg) :: Err + INTEGER(CMISSIntg) :: EquationsSetIndex,Err ! !================================================================================================================================ @@ -191,13 +163,13 @@ PROGRAM STOKESDYNAMICEXAMPLE ! !PROBLEM CONTROL PANEL - NUMBER_GLOBAL_X_ELEMENTS=1 + NUMBER_GLOBAL_X_ELEMENTS=3 NUMBER_GLOBAL_Y_ELEMENTS=3 - NUMBER_GLOBAL_Z_ELEMENTS=1 + NUMBER_GLOBAL_Z_ELEMENTS=3 MESH_COMPONENT_NUMBER_SPACE=1 MESH_COMPONENT_NUMBER_VELOCITY=1 MESH_COMPONENT_NUMBER_PRESSURE=1 - NUMBER_OF_DIMENSIONS=3 + NUMBER_OF_DIMENSIONS=2 BASIS_TYPE=CMFE_BASIS_LINEAR_LAGRANGE_INTERPOLATION BASIS_XI_INTERPOLATION_SPACE=CMFE_BASIS_LINEAR_LAGRANGE_INTERPOLATION BASIS_XI_INTERPOLATION_VELOCITY=CMFE_BASIS_LINEAR_LAGRANGE_INTERPOLATION @@ -211,18 +183,18 @@ PROGRAM STOKESDYNAMICEXAMPLE FIXED_WALL_NODES_STOKES_FLAG=.TRUE. INLET_WALL_NODES_STOKES_FLAG=.TRUE. IF(FIXED_WALL_NODES_STOKES_FLAG) THEN - NUMBER_OF_FIXED_WALL_NODES_STOKES=1 - ALLOCATE(FIXED_WALL_NODES_STOKES(NUMBER_OF_FIXED_WALL_NODES_STOKES)) - FIXED_WALL_NODES_STOKES=[1] + NUMBER_OF_FIXED_WALL_NODES_STOKES=8 + ALLOCATE(FIXED_WALL_NODES_STOKES(NUMBER_OF_FIXED_WALL_NODES_STOKES)) + FIXED_WALL_NODES_STOKES=[1,4,5,8,9,12,13,16] ENDIF IF(INLET_WALL_NODES_STOKES_FLAG) THEN - NUMBER_OF_INLET_WALL_NODES_STOKES=1 - ALLOCATE(INLET_WALL_NODES_STOKES(NUMBER_OF_INLET_WALL_NODES_STOKES)) - INLET_WALL_NODES_STOKES=[2] - !Set initial boundary conditions - BOUNDARY_CONDITIONS_STOKES(1)=0.0_CMISSRP - BOUNDARY_CONDITIONS_STOKES(2)=1.0_CMISSRP - BOUNDARY_CONDITIONS_STOKES(3)=0.0_CMISSRP + NUMBER_OF_INLET_WALL_NODES_STOKES=2 + ALLOCATE(INLET_WALL_NODES_STOKES(NUMBER_OF_INLET_WALL_NODES_STOKES)) + INLET_WALL_NODES_STOKES=[2,3] + !Set initial boundary conditions + BOUNDARY_CONDITIONS_STOKES(1)=0.0_CMISSRP + BOUNDARY_CONDITIONS_STOKES(2)=1.0_CMISSRP + BOUNDARY_CONDITIONS_STOKES(3)=0.0_CMISSRP ENDIF !Set material parameters MU_PARAM_STOKES=1.0_CMISSRP @@ -276,6 +248,7 @@ PROGRAM STOKESDYNAMICEXAMPLE !Start the creation of a new region CALL cmfe_Region_Initialise(Region,Err) CALL cmfe_Region_CreateStart(RegionUserNumber,WorldRegion,Region,Err) + CALL cmfe_Region_LabelSet(Region,"StokesRegion",Err) !Set the regions coordinate system as defined above CALL cmfe_Region_CoordinateSystemSet(Region,CoordinateSystem,Err) !Finish the creation of the region @@ -297,68 +270,71 @@ PROGRAM STOKESDYNAMICEXAMPLE CALL cmfe_Basis_NumberOfXiSet(BasisGeometry,NUMBER_OF_DIMENSIONS,Err) !Set the basis xi interpolation and number of Gauss points IF(NUMBER_OF_DIMENSIONS==2) THEN - CALL cmfe_Basis_InterpolationXiSet(BasisGeometry,[BASIS_XI_INTERPOLATION_SPACE,BASIS_XI_INTERPOLATION_SPACE],Err) - CALL cmfe_Basis_QuadratureNumberOfGaussXiSet(BasisGeometry,[BASIS_XI_GAUSS_SPACE,BASIS_XI_GAUSS_SPACE],Err) - ELSE IF(NUMBER_OF_DIMENSIONS==3) THEN - CALL cmfe_Basis_InterpolationXiSet(BasisGeometry,[BASIS_XI_INTERPOLATION_SPACE,BASIS_XI_INTERPOLATION_SPACE, & - & BASIS_XI_INTERPOLATION_SPACE],Err) - CALL cmfe_Basis_QuadratureNumberOfGaussXiSet(BasisGeometry,[BASIS_XI_GAUSS_SPACE,BASIS_XI_GAUSS_SPACE,BASIS_XI_GAUSS_SPACE],Err) + CALL cmfe_Basis_InterpolationXiSet(BasisGeometry,[BASIS_XI_INTERPOLATION_SPACE,BASIS_XI_INTERPOLATION_SPACE],Err) + CALL cmfe_Basis_QuadratureNumberOfGaussXiSet(BasisGeometry,[BASIS_XI_GAUSS_SPACE,BASIS_XI_GAUSS_SPACE],Err) + ELSE + CALL cmfe_Basis_InterpolationXiSet(BasisGeometry,[BASIS_XI_INTERPOLATION_SPACE,BASIS_XI_INTERPOLATION_SPACE, & + & BASIS_XI_INTERPOLATION_SPACE],Err) + CALL cmfe_Basis_QuadratureNumberOfGaussXiSet(BasisGeometry,[BASIS_XI_GAUSS_SPACE,BASIS_XI_GAUSS_SPACE,BASIS_XI_GAUSS_SPACE], & + & Err) ENDIF !Finish the creation of the basis CALL cmfe_Basis_CreateFinish(BasisGeometry,Err) + !Start the creation of another basis IF(BASIS_XI_INTERPOLATION_VELOCITY==BASIS_XI_INTERPOLATION_SPACE) THEN - BasisVelocity=BasisGeometry + BasisVelocity=BasisGeometry ELSE - MESH_NUMBER_OF_COMPONENTS=MESH_NUMBER_OF_COMPONENTS+1 - !Initialise a new velocity basis - CALL cmfe_Basis_Initialise(BasisVelocity,Err) - !Start the creation of a basis - CALL cmfe_Basis_CreateStart(BASIS_NUMBER_VELOCITY,BasisVelocity,Err) - !Set the basis type (Lagrange/Simplex) - CALL cmfe_Basis_TypeSet(BasisVelocity,BASIS_TYPE,Err) - !Set the basis xi number - CALL cmfe_Basis_NumberOfXiSet(BasisVelocity,NUMBER_OF_DIMENSIONS,Err) - !Set the basis xi interpolation and number of Gauss points - IF(NUMBER_OF_DIMENSIONS==2) THEN - CALL cmfe_Basis_InterpolationXiSet(BasisVelocity,[BASIS_XI_INTERPOLATION_VELOCITY,BASIS_XI_INTERPOLATION_VELOCITY],Err) - CALL cmfe_Basis_QuadratureNumberOfGaussXiSet(BasisVelocity,[BASIS_XI_GAUSS_VELOCITY,BASIS_XI_GAUSS_VELOCITY],Err) - ELSE IF(NUMBER_OF_DIMENSIONS==3) THEN - CALL cmfe_Basis_InterpolationXiSet(BasisVelocity,[BASIS_XI_INTERPOLATION_VELOCITY,BASIS_XI_INTERPOLATION_VELOCITY, & - & BASIS_XI_INTERPOLATION_VELOCITY],Err) - CALL cmfe_Basis_QuadratureNumberOfGaussXiSet(BasisVelocity,[BASIS_XI_GAUSS_VELOCITY,BASIS_XI_GAUSS_VELOCITY, & - & BASIS_XI_GAUSS_VELOCITY],Err) - ENDIF - !Finish the creation of the basis - CALL cmfe_Basis_CreateFinish(BasisVelocity,Err) + MESH_NUMBER_OF_COMPONENTS=MESH_NUMBER_OF_COMPONENTS+1 + !Initialise a new velocity basis + CALL cmfe_Basis_Initialise(BasisVelocity,Err) + !Start the creation of a basis + CALL cmfe_Basis_CreateStart(BASIS_NUMBER_VELOCITY,BasisVelocity,Err) + !Set the basis type (Lagrange/Simplex) + CALL cmfe_Basis_TypeSet(BasisVelocity,BASIS_TYPE,Err) + !Set the basis xi number + CALL cmfe_Basis_NumberOfXiSet(BasisVelocity,NUMBER_OF_DIMENSIONS,Err) + !Set the basis xi interpolation and number of Gauss points + IF(NUMBER_OF_DIMENSIONS==2) THEN + CALL cmfe_Basis_InterpolationXiSet(BasisVelocity,[BASIS_XI_INTERPOLATION_VELOCITY,BASIS_XI_INTERPOLATION_VELOCITY],Err) + CALL cmfe_Basis_QuadratureNumberOfGaussXiSet(BasisVelocity,[BASIS_XI_GAUSS_VELOCITY,BASIS_XI_GAUSS_VELOCITY],Err) + ELSE + CALL cmfe_Basis_InterpolationXiSet(BasisVelocity,[BASIS_XI_INTERPOLATION_VELOCITY,BASIS_XI_INTERPOLATION_VELOCITY, & + & BASIS_XI_INTERPOLATION_VELOCITY],Err) + CALL cmfe_Basis_QuadratureNumberOfGaussXiSet(BasisVelocity,[BASIS_XI_GAUSS_VELOCITY,BASIS_XI_GAUSS_VELOCITY, & + & BASIS_XI_GAUSS_VELOCITY],Err) + ENDIF + !Finish the creation of the basis + CALL cmfe_Basis_CreateFinish(BasisVelocity,Err) ENDIF + !Start the creation of another basis IF(BASIS_XI_INTERPOLATION_PRESSURE==BASIS_XI_INTERPOLATION_SPACE) THEN - BasisPressure=BasisGeometry + BasisPressure=BasisGeometry ELSE IF(BASIS_XI_INTERPOLATION_PRESSURE==BASIS_XI_INTERPOLATION_VELOCITY) THEN - BasisPressure=BasisVelocity + BasisPressure=BasisVelocity ELSE - MESH_NUMBER_OF_COMPONENTS=MESH_NUMBER_OF_COMPONENTS+1 - !Initialise a new pressure basis - CALL cmfe_Basis_Initialise(BasisPressure,Err) - !Start the creation of a basis - CALL cmfe_Basis_CreateStart(BASIS_NUMBER_PRESSURE,BasisPressure,Err) - !Set the basis type (Lagrange/Simplex) - CALL cmfe_Basis_TypeSet(BasisPressure,BASIS_TYPE,Err) - !Set the basis xi number - CALL cmfe_Basis_NumberOfXiSet(BasisPressure,NUMBER_OF_DIMENSIONS,Err) - !Set the basis xi interpolation and number of Gauss points - IF(NUMBER_OF_DIMENSIONS==2) THEN - CALL cmfe_Basis_InterpolationXiSet(BasisPressure,[BASIS_XI_INTERPOLATION_PRESSURE,BASIS_XI_INTERPOLATION_PRESSURE],Err) - CALL cmfe_Basis_QuadratureNumberOfGaussXiSet(BasisPressure,[BASIS_XI_GAUSS_PRESSURE,BASIS_XI_GAUSS_PRESSURE],Err) - ELSE IF(NUMBER_OF_DIMENSIONS==3) THEN - CALL cmfe_Basis_InterpolationXiSet(BasisPressure,[BASIS_XI_INTERPOLATION_PRESSURE,BASIS_XI_INTERPOLATION_PRESSURE, & - & BASIS_XI_INTERPOLATION_PRESSURE],Err) - CALL cmfe_Basis_QuadratureNumberOfGaussXiSet(BasisPressure,[BASIS_XI_GAUSS_PRESSURE,BASIS_XI_GAUSS_PRESSURE, & - & BASIS_XI_GAUSS_PRESSURE],Err) - ENDIF - !Finish the creation of the basis - CALL cmfe_Basis_CreateFinish(BasisPressure,Err) + MESH_NUMBER_OF_COMPONENTS=MESH_NUMBER_OF_COMPONENTS+1 + !Initialise a new pressure basis + CALL cmfe_Basis_Initialise(BasisPressure,Err) + !Start the creation of a basis + CALL cmfe_Basis_CreateStart(BASIS_NUMBER_PRESSURE,BasisPressure,Err) + !Set the basis type (Lagrange/Simplex) + CALL cmfe_Basis_TypeSet(BasisPressure,BASIS_TYPE,Err) + !Set the basis xi number + CALL cmfe_Basis_NumberOfXiSet(BasisPressure,NUMBER_OF_DIMENSIONS,Err) + !Set the basis xi interpolation and number of Gauss points + IF(NUMBER_OF_DIMENSIONS==2) THEN + CALL cmfe_Basis_InterpolationXiSet(BasisPressure,[BASIS_XI_INTERPOLATION_PRESSURE,BASIS_XI_INTERPOLATION_PRESSURE],Err) + CALL cmfe_Basis_QuadratureNumberOfGaussXiSet(BasisPressure,[BASIS_XI_GAUSS_PRESSURE,BASIS_XI_GAUSS_PRESSURE],Err) + ELSE + CALL cmfe_Basis_InterpolationXiSet(BasisPressure,[BASIS_XI_INTERPOLATION_PRESSURE,BASIS_XI_INTERPOLATION_PRESSURE, & + & BASIS_XI_INTERPOLATION_PRESSURE],Err) + CALL cmfe_Basis_QuadratureNumberOfGaussXiSet(BasisPressure,[BASIS_XI_GAUSS_PRESSURE,BASIS_XI_GAUSS_PRESSURE, & + & BASIS_XI_GAUSS_PRESSURE],Err) + ENDIF + !Finish the creation of the basis + CALL cmfe_Basis_CreateFinish(BasisPressure,Err) ENDIF ! @@ -375,9 +351,14 @@ PROGRAM STOKESDYNAMICEXAMPLE !Set the default basis CALL cmfe_GeneratedMesh_BasisSet(GeneratedMesh,BasisGeometry,Err) !Define the mesh on the region - CALL cmfe_GeneratedMesh_ExtentSet(GeneratedMesh,[WIDTH,HEIGHT,LENGTH],Err) - CALL cmfe_GeneratedMesh_NumberOfElementsSet(GeneratedMesh,[NUMBER_GLOBAL_X_ELEMENTS,NUMBER_GLOBAL_Y_ELEMENTS, & - & NUMBER_GLOBAL_Z_ELEMENTS],Err) + IF(NUMBER_OF_DIMENSIONS==2) THEN + CALL cmfe_GeneratedMesh_ExtentSet(GeneratedMesh,[WIDTH,HEIGHT],Err) + CALL cmfe_GeneratedMesh_NumberOfElementsSet(GeneratedMesh,[NUMBER_GLOBAL_X_ELEMENTS,NUMBER_GLOBAL_Y_ELEMENTS],Err) + ELSE + CALL cmfe_GeneratedMesh_ExtentSet(GeneratedMesh,[WIDTH,HEIGHT,LENGTH],Err) + CALL cmfe_GeneratedMesh_NumberOfElementsSet(GeneratedMesh,[NUMBER_GLOBAL_X_ELEMENTS,NUMBER_GLOBAL_Y_ELEMENTS, & + & NUMBER_GLOBAL_Z_ELEMENTS],Err) + ENDIF !Finish the creation of a generated mesh in the region CALL cmfe_Mesh_Initialise(Mesh,Err) CALL cmfe_GeneratedMesh_CreateFinish(GeneratedMesh,MeshUserNumber,Mesh,Err) @@ -386,7 +367,7 @@ PROGRAM STOKESDYNAMICEXAMPLE !================================================================================================================================ ! - !GEOMETRIC FIELD + !MESH DECOMPOSITION !Create a decomposition CALL cmfe_Decomposition_Initialise(Decomposition,Err) @@ -397,6 +378,12 @@ PROGRAM STOKESDYNAMICEXAMPLE !Finish the decomposition CALL cmfe_Decomposition_CreateFinish(Decomposition,Err) + ! + !================================================================================================================================ + ! + + !GEOMETRIC FIELD + !Start to create a default (geometric) field on the region CALL cmfe_Field_Initialise(GeometricField,Err) CALL cmfe_Field_CreateStart(GeometricFieldUserNumber,Region,GeometricField,Err) @@ -407,18 +394,14 @@ PROGRAM STOKESDYNAMICEXAMPLE !Set the scaling to use CALL cmfe_Field_ScalingTypeSet(GeometricField,CMFE_FIELD_NO_SCALING,Err) !Set the mesh component to be used by the field components. - !DO COMPONENT_NUMBER=1,NUMBER_OF_DIMENSIONS - ! CALL cmfe_Field_ComponentMeshComponentSet(GeometricField,CMFE_FIELD_U_VARIABLE_TYPE,COMPONENT_NUMBER, & - ! & MESH_COMPONENT_NUMBER_SPACE,Err) - !ENDDO - CALL cmfe_Field_ComponentMeshComponentSet(GeometricField,CMFE_FIELD_U_VARIABLE_TYPE,1,1,Err) - CALL cmfe_Field_ComponentMeshComponentSet(GeometricField,CMFE_FIELD_U_VARIABLE_TYPE,2,1,Err) - CALL cmfe_Field_ComponentMeshComponentSet(GeometricField,CMFE_FIELD_U_VARIABLE_TYPE,3,1,Err) + DO COMPONENT_NUMBER=1,NUMBER_OF_DIMENSIONS + CALL cmfe_Field_ComponentMeshComponentSet(GeometricField,CMFE_FIELD_U_VARIABLE_TYPE,COMPONENT_NUMBER, & + & MESH_COMPONENT_NUMBER_SPACE,Err) + ENDDO !Finish creating the field CALL cmfe_Field_CreateFinish(GeometricField,Err) - - CALL cmfe_Field_ParameterSetUpdateStart(GeometricField,CMFE_FIELD_U_VARIABLE_TYPE,CMFE_FIELD_VALUES_SET_TYPE,Err) - CALL cmfe_Field_ParameterSetUpdateFinish(GeometricField,CMFE_FIELD_U_VARIABLE_TYPE,CMFE_FIELD_VALUES_SET_TYPE,Err) + !Update the geometric field parameters + CALL cmfe_GeneratedMesh_GeometricParametersCalculate(GeneratedMesh,GeometricField,Err) ! !================================================================================================================================ @@ -429,11 +412,10 @@ PROGRAM STOKESDYNAMICEXAMPLE !Create the equations set for dynamic Stokes CALL cmfe_EquationsSet_Initialise(EquationsSetStokes,Err) CALL cmfe_Field_Initialise(EquationsSetField,Err) - CALL cmfe_EquationsSet_CreateStart(EquationsSetUserNumberStokes,Region,GeometricField,[CMFE_EQUATIONS_SET_FLUID_MECHANICS_CLASS, & - & CMFE_EQUATIONS_SET_STOKES_EQUATION_TYPE,CMFE_EQUATIONS_SET_TRANSIENT_STOKES_SUBTYPE],EquationsSetFieldUserNumber, & - & EquationsSetField,EquationsSetStokes,Err) !Set the equations set to be a dynamic Stokes problem - + CALL cmfe_EquationsSet_CreateStart(EquationsSetUserNumberStokes,Region,GeometricField,[CMFE_EQUATIONS_SET_FLUID_MECHANICS_CLASS, & + & CMFE_EQUATIONS_SET_STOKES_EQUATION_TYPE,CMFE_EQUATIONS_SET_TRANSIENT_STOKES_SUBTYPE],EquationsSetFieldUserNumber, & + & EquationsSetField,EquationsSetStokes,Err) !Finish creating the equations set CALL cmfe_EquationsSet_CreateFinish(EquationsSetStokes,Err) @@ -445,27 +427,27 @@ PROGRAM STOKESDYNAMICEXAMPLE !Create the equations set dependent field variables for dynamic Stokes CALL cmfe_Field_Initialise(DependentFieldStokes,Err) - CALL cmfe_EquationsSet_DependentCreateStart(EquationsSetStokes,DependentFieldUserNumberStokes, & - & DependentFieldStokes,Err) + CALL cmfe_EquationsSet_DependentCreateStart(EquationsSetStokes,DependentFieldUserNumberStokes,DependentFieldStokes,Err) + CALL cmfe_Field_VariableLabelSet(DependentFieldStokes,CMFE_FIELD_U_VARIABLE_TYPE,"U",Err) !Set the mesh component to be used by the field components. DO COMPONENT_NUMBER=1,NUMBER_OF_DIMENSIONS - CALL cmfe_Field_ComponentMeshComponentSet(DependentFieldStokes,CMFE_FIELD_U_VARIABLE_TYPE,COMPONENT_NUMBER, & - & MESH_COMPONENT_NUMBER_VELOCITY,Err) - CALL cmfe_Field_ComponentMeshComponentSet(DependentFieldStokes,CMFE_FIELD_DELUDELN_VARIABLE_TYPE,COMPONENT_NUMBER, & - & MESH_COMPONENT_NUMBER_VELOCITY,Err) + CALL cmfe_Field_ComponentMeshComponentSet(DependentFieldStokes,CMFE_FIELD_U_VARIABLE_TYPE,COMPONENT_NUMBER, & + & MESH_COMPONENT_NUMBER_VELOCITY,Err) + CALL cmfe_Field_ComponentMeshComponentSet(DependentFieldStokes,CMFE_FIELD_DELUDELN_VARIABLE_TYPE,COMPONENT_NUMBER, & + & MESH_COMPONENT_NUMBER_VELOCITY,Err) ENDDO COMPONENT_NUMBER=NUMBER_OF_DIMENSIONS+1 - CALL cmfe_Field_ComponentMeshComponentSet(DependentFieldStokes,CMFE_FIELD_U_VARIABLE_TYPE,COMPONENT_NUMBER, & - & MESH_COMPONENT_NUMBER_PRESSURE,Err) - CALL cmfe_Field_ComponentMeshComponentSet(DependentFieldStokes,CMFE_FIELD_DELUDELN_VARIABLE_TYPE,COMPONENT_NUMBER, & - & MESH_COMPONENT_NUMBER_PRESSURE,Err) + CALL cmfe_Field_ComponentMeshComponentSet(DependentFieldStokes,CMFE_FIELD_U_VARIABLE_TYPE,COMPONENT_NUMBER, & + & MESH_COMPONENT_NUMBER_PRESSURE,Err) + CALL cmfe_Field_ComponentMeshComponentSet(DependentFieldStokes,CMFE_FIELD_DELUDELN_VARIABLE_TYPE,COMPONENT_NUMBER, & + & MESH_COMPONENT_NUMBER_PRESSURE,Err) !Finish the equations set dependent field variables CALL cmfe_EquationsSet_DependentCreateFinish(EquationsSetStokes,Err) !Initialise dependent field DO COMPONENT_NUMBER=1,NUMBER_OF_DIMENSIONS - CALL cmfe_Field_ComponentValuesInitialise(DependentFieldStokes,CMFE_FIELD_U_VARIABLE_TYPE,CMFE_FIELD_VALUES_SET_TYPE, & - & COMPONENT_NUMBER,INITIAL_FIELD_STOKES(COMPONENT_NUMBER),Err) + CALL cmfe_Field_ComponentValuesInitialise(DependentFieldStokes,CMFE_FIELD_U_VARIABLE_TYPE,CMFE_FIELD_VALUES_SET_TYPE, & + & COMPONENT_NUMBER,INITIAL_FIELD_STOKES(COMPONENT_NUMBER),Err) ENDDO @@ -477,14 +459,27 @@ PROGRAM STOKESDYNAMICEXAMPLE !Create the equations set materials field variables for dynamic Stokes CALL cmfe_Field_Initialise(MaterialsFieldStokes,Err) - CALL cmfe_EquationsSet_MaterialsCreateStart(EquationsSetStokes,MaterialsFieldUserNumberStokes, & - & MaterialsFieldStokes,Err) + CALL cmfe_EquationsSet_MaterialsCreateStart(EquationsSetStokes,MaterialsFieldUserNumberStokes,MaterialsFieldStokes,Err) + CALL cmfe_Field_VariableLabelSet(MaterialsFieldStokes,CMFE_FIELD_U_VARIABLE_TYPE,"Materials",Err) !Finish the equations set materials field variables CALL cmfe_EquationsSet_MaterialsCreateFinish(EquationsSetStokes,Err) CALL cmfe_Field_ComponentValuesInitialise(MaterialsFieldStokes,CMFE_FIELD_U_VARIABLE_TYPE,CMFE_FIELD_VALUES_SET_TYPE, & - & MaterialsFieldUserNumberStokesMu,MU_PARAM_STOKES,Err) + & MaterialsFieldUserNumberStokesMu,MU_PARAM_STOKES,Err) CALL cmfe_Field_ComponentValuesInitialise(MaterialsFieldStokes,CMFE_FIELD_U_VARIABLE_TYPE,CMFE_FIELD_VALUES_SET_TYPE, & - & MaterialsFieldUserNumberStokesRho,RHO_PARAM_STOKES,Err) + & MaterialsFieldUserNumberStokesRho,RHO_PARAM_STOKES,Err) + + ! + !================================================================================================================================ + ! + + !ANALYTIC FIELD + + !Create the equations set analytic field variables + !CALL cmfe_Field_Initialise(AnalyticField,Err) + !CALL cmfe_EquationsSet_AnalyticCreateStart(EquationsSetStokes,CMFE_EQUATIONS_SET_NAVIER_STOKES_EQUATION_TWO_DIM_4, & + ! & AnalyticFieldUserNumber,AnalyticField,Err) + !Finish the equations set analytic field variables + !CALL cmfe_EquationsSet_AnalyticCreateFinish(EquationsSetStokes,Err) ! !================================================================================================================================ @@ -514,16 +509,23 @@ PROGRAM STOKESDYNAMICEXAMPLE CALL cmfe_Problem_Initialise(Problem,Err) CALL cmfe_ControlLoop_Initialise(ControlLoop,Err) CALL cmfe_Problem_CreateStart(ProblemUserNumber,[CMFE_PROBLEM_FLUID_MECHANICS_CLASS,CMFE_PROBLEM_STOKES_EQUATION_TYPE, & - & CMFE_PROBLEM_TRANSIENT_STOKES_SUBTYPE],Problem,Err) + & CMFE_PROBLEM_TRANSIENT_STOKES_SUBTYPE],Problem,Err) !Finish the creation of a problem. CALL cmfe_Problem_CreateFinish(Problem,Err) + + ! + !================================================================================================================================ + ! + + !CONTROL LOOP + !Start the creation of the problem control loop CALL cmfe_Problem_ControlLoopCreateStart(Problem,Err) !Get the control loop CALL cmfe_Problem_ControlLoopGet(Problem,CMFE_CONTROL_LOOP_NODE,ControlLoop,Err) !Set the times CALL cmfe_ControlLoop_TimesSet(ControlLoop,DYNAMIC_SOLVER_STOKES_START_TIME,DYNAMIC_SOLVER_STOKES_STOP_TIME, & - & DYNAMIC_SOLVER_STOKES_TIME_INCREMENT,Err) + & DYNAMIC_SOLVER_STOKES_TIME_INCREMENT,Err) !Set the output timing CALL cmfe_ControlLoop_TimeOutputSet(ControlLoop,DYNAMIC_SOLVER_STOKES_OUTPUT_FREQUENCY,Err) !Finish creating the problem control loop @@ -551,15 +553,15 @@ PROGRAM STOKESDYNAMICEXAMPLE CALL cmfe_Solver_OutputTypeSet(LinearSolverStokes,LINEAR_SOLVER_STOKES_OUTPUT_TYPE,Err) !Set the solver settings IF(LINEAR_SOLVER_STOKES_DIRECT_FLAG) THEN - CALL cmfe_Solver_LinearTypeSet(LinearSolverStokes,CMFE_SOLVER_LINEAR_DIRECT_SOLVE_TYPE,Err) - CALL cmfe_Solver_LibraryTypeSet(LinearSolverStokes,CMFE_SOLVER_MUMPS_LIBRARY,Err) + CALL cmfe_Solver_LinearTypeSet(LinearSolverStokes,CMFE_SOLVER_LINEAR_DIRECT_SOLVE_TYPE,Err) + CALL cmfe_Solver_LibraryTypeSet(LinearSolverStokes,CMFE_SOLVER_MUMPS_LIBRARY,Err) ELSE - CALL cmfe_Solver_LinearTypeSet(LinearSolverStokes,CMFE_SOLVER_LINEAR_ITERATIVE_SOLVE_TYPE,Err) - CALL cmfe_Solver_LinearIterativeMaximumIterationsSet(LinearSolverStokes,MAXIMUM_ITERATIONS,Err) - CALL cmfe_Solver_LinearIterativeDivergenceToleranceSet(LinearSolverStokes,DIVERGENCE_TOLERANCE,Err) - CALL cmfe_Solver_LinearIterativeRelativeToleranceSet(LinearSolverStokes,RELATIVE_TOLERANCE,Err) - CALL cmfe_Solver_LinearIterativeAbsoluteToleranceSet(LinearSolverStokes,ABSOLUTE_TOLERANCE,Err) - CALL cmfe_Solver_LinearIterativeGMRESRestartSet(LinearSolverStokes,RESTART_VALUE,Err) + CALL cmfe_Solver_LinearTypeSet(LinearSolverStokes,CMFE_SOLVER_LINEAR_ITERATIVE_SOLVE_TYPE,Err) + CALL cmfe_Solver_LinearIterativeMaximumIterationsSet(LinearSolverStokes,MAXIMUM_ITERATIONS,Err) + CALL cmfe_Solver_LinearIterativeDivergenceToleranceSet(LinearSolverStokes,DIVERGENCE_TOLERANCE,Err) + CALL cmfe_Solver_LinearIterativeRelativeToleranceSet(LinearSolverStokes,RELATIVE_TOLERANCE,Err) + CALL cmfe_Solver_LinearIterativeAbsoluteToleranceSet(LinearSolverStokes,ABSOLUTE_TOLERANCE,Err) + CALL cmfe_Solver_LinearIterativeGMRESRestartSet(LinearSolverStokes,RESTART_VALUE,Err) ENDIF !Finish the creation of the problem solver CALL cmfe_Problem_SolversCreateFinish(Problem,Err) @@ -596,36 +598,37 @@ PROGRAM STOKESDYNAMICEXAMPLE CALL cmfe_SolverEquations_BoundaryConditionsCreateStart(SolverEquationsStokes,BoundaryConditionsStokes,Err) !Set fixed wall nodes IF(FIXED_WALL_NODES_STOKES_FLAG) THEN - DO NODE_COUNTER=1,NUMBER_OF_FIXED_WALL_NODES_STOKES - NODE_NUMBER=FIXED_WALL_NODES_STOKES(NODE_COUNTER) - CONDITION=CMFE_BOUNDARY_CONDITION_FIXED_WALL - CALL cmfe_Decomposition_NodeDomainGet(Decomposition,NODE_NUMBER,1,BoundaryNodeDomain,Err) - IF(BoundaryNodeDomain==ComputationalNodeNumber) THEN - DO COMPONENT_NUMBER=1,NUMBER_OF_DIMENSIONS - VALUE=0.0_CMISSRP - CALL cmfe_BoundaryConditions_SetNode(BoundaryConditionsStokes,DependentFieldStokes,CMFE_FIELD_U_VARIABLE_TYPE,1, & - & CMFE_NO_GLOBAL_DERIV, & - & NODE_NUMBER,COMPONENT_NUMBER,CONDITION,VALUE,Err) - ENDDO - ENDIF - ENDDO + DO NODE_COUNTER=1,NUMBER_OF_FIXED_WALL_NODES_STOKES + NODE_NUMBER=FIXED_WALL_NODES_STOKES(NODE_COUNTER) + CONDITION=CMFE_BOUNDARY_CONDITION_FIXED_WALL + CALL cmfe_Decomposition_NodeDomainGet(Decomposition,NODE_NUMBER,1,BoundaryNodeDomain,Err) + IF(BoundaryNodeDomain==ComputationalNodeNumber) THEN + DO COMPONENT_NUMBER=1,NUMBER_OF_DIMENSIONS + VALUE=0.0_CMISSRP + CALL cmfe_BoundaryConditions_SetNode(BoundaryConditionsStokes,DependentFieldStokes,CMFE_FIELD_U_VARIABLE_TYPE,1, & + & CMFE_NO_GLOBAL_DERIV, & + & NODE_NUMBER,COMPONENT_NUMBER,CONDITION,VALUE,Err) + ENDDO + ENDIF + ENDDO ENDIF !Set velocity boundary conditions IF(INLET_WALL_NODES_STOKES_FLAG) THEN - DO NODE_COUNTER=1,NUMBER_OF_INLET_WALL_NODES_STOKES - NODE_NUMBER=INLET_WALL_NODES_STOKES(NODE_COUNTER) - CONDITION=CMFE_BOUNDARY_CONDITION_FIXED_INLET - CALL cmfe_Decomposition_NodeDomainGet(Decomposition,NODE_NUMBER,1,BoundaryNodeDomain,Err) - IF(BoundaryNodeDomain==ComputationalNodeNumber) THEN - DO COMPONENT_NUMBER=1,NUMBER_OF_DIMENSIONS - VALUE=BOUNDARY_CONDITIONS_STOKES(COMPONENT_NUMBER) - CALL cmfe_BoundaryConditions_SetNode(BoundaryConditionsStokes,DependentFieldStokes,CMFE_FIELD_U_VARIABLE_TYPE,1, & - & CMFE_NO_GLOBAL_DERIV, & - & NODE_NUMBER,COMPONENT_NUMBER,CONDITION,VALUE,Err) - ENDDO - ENDIF - ENDDO + DO NODE_COUNTER=1,NUMBER_OF_INLET_WALL_NODES_STOKES + NODE_NUMBER=INLET_WALL_NODES_STOKES(NODE_COUNTER) + CONDITION=CMFE_BOUNDARY_CONDITION_FIXED_INLET + CALL cmfe_Decomposition_NodeDomainGet(Decomposition,NODE_NUMBER,1,BoundaryNodeDomain,Err) + IF(BoundaryNodeDomain==ComputationalNodeNumber) THEN + DO COMPONENT_NUMBER=1,NUMBER_OF_DIMENSIONS + VALUE=BOUNDARY_CONDITIONS_STOKES(COMPONENT_NUMBER) + CALL cmfe_BoundaryConditions_SetNode(BoundaryConditionsStokes,DependentFieldStokes,CMFE_FIELD_U_VARIABLE_TYPE,1, & + & CMFE_NO_GLOBAL_DERIV, & + & NODE_NUMBER,COMPONENT_NUMBER,CONDITION,VALUE,Err) + ENDDO + ENDIF + ENDDO ENDIF + !CALL cmfe_SolverEquations_BoundaryConditionsAnalytic(SolverEquationsStokes,Err) !Finish the creation of the equations set boundary conditions CALL cmfe_SolverEquations_BoundaryConditionsCreateFinish(SolverEquationsStokes,Err) @@ -633,10 +636,12 @@ PROGRAM STOKESDYNAMICEXAMPLE !================================================================================================================================ ! - !RUN SOLVERS + !Output Analytic analysis + !Call cmfe_AnalyticAnalysis_Output(DependentFieldStokes,"StokesAnalytic",Err) - !Turn of PETSc error handling - !CALL PETSC_ERRORHANDLING_SET_ON(ERR,ERROR,*999) + ! + !================================================================================================================================ + ! !Solve the problem WRITE(*,'(A)') "Solving problem..." @@ -649,15 +654,15 @@ PROGRAM STOKESDYNAMICEXAMPLE !OUTPUT - EXPORT_FIELD_IO=.TRUE. + EXPORT_FIELD_IO=.FALSE. IF(EXPORT_FIELD_IO) THEN - WRITE(*,'(A)') "Exporting fields..." - CALL cmfe_Fields_Initialise(Fields,Err) - CALL cmfe_Fields_Create(Region,Fields,Err) - CALL cmfe_Fields_NodesExport(Fields,"DynamicStokes","FORTRAN",Err) - CALL cmfe_Fields_ElementsExport(Fields,"DynamicStokes","FORTRAN",Err) - CALL cmfe_Fields_Finalise(Fields,Err) - WRITE(*,'(A)') "Field exported!" + WRITE(*,'(A)') "Exporting fields..." + CALL cmfe_Fields_Initialise(Fields,Err) + CALL cmfe_Fields_Create(Region,Fields,Err) + CALL cmfe_Fields_NodesExport(Fields,"stokes_dynamic","FORTRAN",Err) + CALL cmfe_Fields_ElementsExport(Fields,"stokes_dynamic","FORTRAN",Err) + CALL cmfe_Fields_Finalise(Fields,Err) + WRITE(*,'(A)') "Field exported!" ENDIF !Finialise CMISS @@ -665,4 +670,4 @@ PROGRAM STOKESDYNAMICEXAMPLE WRITE(*,'(A)') "Program successfully completed." STOP -END PROGRAM STOKESDYNAMICEXAMPLE +END PROGRAM stokes_dynamic diff --git a/src/fortran/visualiseStokesDynamic.cmgui b/src/fortran/visualiseStokesDynamic.cmgui new file mode 100755 index 0000000..62a98b5 --- /dev/null +++ b/src/fortran/visualiseStokesDynamic.cmgui @@ -0,0 +1,47 @@ +$dt=1; # frequency +$steps=40000; # number of time steps +$cores=1; # number of processors + +$w=500; # width of the graphical window +$h=1000; # height of the graphical window + +#Read in the sequence of nodal positions. +for ($i=0;$i<$steps;$i=$i+$dt) + { + $time = $i/($dt*100) + for ($j=0;$j<$cores;$j=$j+1) + { + $filename = sprintf("./output/MainTime_%01d.part%01d.exnode", $i, $j); + print "Reading $filename time $time\n"; + gfx read node "$filename" time $i; + } + } + +#Read in the element description +for ($k=0;$k<$cores;$k=$k+1) + { + $filename = sprintf("./output/MainTime_0.part%01d.exelem", $k); + gfx read element "$filename"; + } + +gfx define faces egroup StokesRegion + +gfx define field Coordinate.x component Coordinate.x +gfx define field Coordinate.y component Coordinate.y +gfx define field Coordinate.z component Coordinate.z + +gfx def field x_velocity component U.1 +gfx def field y_velocity component U.2 +gfx def field z_velocity component U.3 +gfx def field pressure component U.4 + +gfx define field vector_field coord rectangular_cartesian component U.1 U.2 U.3 U.4 +gfx modify g_element StokesRegion node_points data y_velocity + +gfx edit scene +gfx create window 1; + +#Set the timekeeper playing +gfx timekeeper default set 1.0; +gfx timekeeper default speed 1; +gfx create time_editor diff --git a/src/python/expected_results/output/MainTime_0.part0.exelem b/src/python/expected_results/output/MainTime_0.part0.exelem new file mode 100644 index 0000000..1eab4e3 --- /dev/null +++ b/src/python/expected_results/output/MainTime_0.part0.exelem @@ -0,0 +1,201 @@ + Group name: StokesRegion + Shape. Dimension=2 + #Scale factor sets= 1 + l.Lagrange*l.Lagrange, #Scale factors=4 + #Nodes= 4 + #Fields=4 + 1) Coordinate, coordinate, rectangular cartesian, #Components=2 + x. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + y. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 2) U, field, rectangular cartesian, #Components=3 + 1. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 3. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 3) del U_del n, field, rectangular cartesian, #Components=3 + 1. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 3. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 4) Material, field, rectangular cartesian, #Components=2 + 1. constant*constant, no modify, grid based. + #xi1=0, #xi2=0 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + Element: 1 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 1 2 5 6 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 2 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 2 3 6 7 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 3 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 3 4 7 8 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 4 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 5 6 9 10 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 5 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 6 7 10 11 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 6 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 7 8 11 12 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 7 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 9 10 13 14 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 8 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 10 11 14 15 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 9 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 11 12 15 16 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 diff --git a/src/python/expected_results/output/MainTime_0.part0.exnode b/src/python/expected_results/output/MainTime_0.part0.exnode new file mode 100644 index 0000000..a1321da --- /dev/null +++ b/src/python/expected_results/output/MainTime_0.part0.exnode @@ -0,0 +1,192 @@ + Group name: StokesRegion + #Fields=4 + 1) Coordinate, coordinate, rectangular cartesian, #Components=2 + x. Value index= 1, #Derivatives= 0 + y. Value index= 2, #Derivatives= 0 + 2) U, field, rectangular cartesian, #Components=3 + 1. Value index= 3, #Derivatives= 0 + 2. Value index= 4, #Derivatives= 0 + 3. Value index= 5, #Derivatives= 0 + 3) del U_del n, field, rectangular cartesian, #Components=3 + 1. Value index= 6, #Derivatives= 0 + 2. Value index= 7, #Derivatives= 0 + 3. Value index= 8, #Derivatives= 0 + 4) Material, field, rectangular cartesian, #Components=2 + 1. Value index= 9, #Derivatives= 0 + 2. Value index= 10, #Derivatives= 0 + Node: 1 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 2 + 3.3333333333333331E-01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 3 + 6.6666666666666663E-01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 4 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 5 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 6 + 3.3333333333333331E-01 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 7 + 6.6666666666666663E-01 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 8 + 1.0000000000000000E+00 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 9 + 0.0000000000000000E+00 + 2.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 10 + 3.3333333333333331E-01 + 2.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 11 + 6.6666666666666663E-01 + 2.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 12 + 1.0000000000000000E+00 + 2.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 13 + 0.0000000000000000E+00 + 3.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 14 + 3.3333333333333331E-01 + 3.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 15 + 6.6666666666666663E-01 + 3.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 16 + 1.0000000000000000E+00 + 3.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 diff --git a/src/python/expected_results/output/MainTime_1.part0.exelem b/src/python/expected_results/output/MainTime_1.part0.exelem new file mode 100644 index 0000000..1eab4e3 --- /dev/null +++ b/src/python/expected_results/output/MainTime_1.part0.exelem @@ -0,0 +1,201 @@ + Group name: StokesRegion + Shape. Dimension=2 + #Scale factor sets= 1 + l.Lagrange*l.Lagrange, #Scale factors=4 + #Nodes= 4 + #Fields=4 + 1) Coordinate, coordinate, rectangular cartesian, #Components=2 + x. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + y. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 2) U, field, rectangular cartesian, #Components=3 + 1. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 3. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 3) del U_del n, field, rectangular cartesian, #Components=3 + 1. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 3. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 4) Material, field, rectangular cartesian, #Components=2 + 1. constant*constant, no modify, grid based. + #xi1=0, #xi2=0 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + Element: 1 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 1 2 5 6 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 2 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 2 3 6 7 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 3 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 3 4 7 8 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 4 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 5 6 9 10 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 5 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 6 7 10 11 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 6 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 7 8 11 12 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 7 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 9 10 13 14 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 8 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 10 11 14 15 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 9 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 11 12 15 16 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 diff --git a/src/python/expected_results/output/MainTime_1.part0.exnode b/src/python/expected_results/output/MainTime_1.part0.exnode new file mode 100644 index 0000000..297c81c --- /dev/null +++ b/src/python/expected_results/output/MainTime_1.part0.exnode @@ -0,0 +1,192 @@ + Group name: StokesRegion + #Fields=4 + 1) Coordinate, coordinate, rectangular cartesian, #Components=2 + x. Value index= 1, #Derivatives= 0 + y. Value index= 2, #Derivatives= 0 + 2) U, field, rectangular cartesian, #Components=3 + 1. Value index= 3, #Derivatives= 0 + 2. Value index= 4, #Derivatives= 0 + 3. Value index= 5, #Derivatives= 0 + 3) del U_del n, field, rectangular cartesian, #Components=3 + 1. Value index= 6, #Derivatives= 0 + 2. Value index= 7, #Derivatives= 0 + 3. Value index= 8, #Derivatives= 0 + 4) Material, field, rectangular cartesian, #Components=2 + 1. Value index= 9, #Derivatives= 0 + 2. Value index= 10, #Derivatives= 0 + Node: 1 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.1812122013393381E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 2 + 3.3333333333333331E-01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 2.9846130637595071E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 3 + 6.6666666666666663E-01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 2.9846130637590665E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 4 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.1812122013388382E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 5 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 2.6286825726135184E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 6 + 3.3333333333333331E-01 + 1.0000000000000000E+00 + -9.7908019770077493E-12 + 9.9999999999979972E-01 + 2.1200746481192425E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 7 + 6.6666666666666663E-01 + 1.0000000000000000E+00 + -8.0700587257565004E-12 + 9.9999999999738842E-01 + 2.1200746481177053E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 8 + 1.0000000000000000E+00 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 2.6286825726120352E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 9 + 0.0000000000000000E+00 + 2.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 5.5263166816422791E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 10 + 3.3333333333333331E-01 + 2.0000000000000000E+00 + 1.8382265636428370E-11 + 1.0000000000002103E+00 + 7.8366250373066046E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 11 + 6.6666666666666663E-01 + 2.0000000000000000E+00 + 1.5851453116772888E-11 + 9.9999999999993827E-01 + 7.8366250373119266E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 12 + 1.0000000000000000E+00 + 2.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 5.5263166816489182E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 13 + 0.0000000000000000E+00 + 3.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 5.8443881702858924E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 14 + 3.3333333333333331E-01 + 3.0000000000000000E+00 + -1.1511640515858268E-11 + 1.0000000000005844E+00 + 1.6892339923515256E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 15 + 6.6666666666666663E-01 + 3.0000000000000000E+00 + -1.5510357975943329E-11 + 1.0000000000060576E+00 + 1.6892339923398016E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 16 + 1.0000000000000000E+00 + 3.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 5.8443881702755824E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 diff --git a/src/python/expected_results/output/MainTime_2.part0.exelem b/src/python/expected_results/output/MainTime_2.part0.exelem new file mode 100644 index 0000000..1eab4e3 --- /dev/null +++ b/src/python/expected_results/output/MainTime_2.part0.exelem @@ -0,0 +1,201 @@ + Group name: StokesRegion + Shape. Dimension=2 + #Scale factor sets= 1 + l.Lagrange*l.Lagrange, #Scale factors=4 + #Nodes= 4 + #Fields=4 + 1) Coordinate, coordinate, rectangular cartesian, #Components=2 + x. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + y. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 2) U, field, rectangular cartesian, #Components=3 + 1. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 3. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 3) del U_del n, field, rectangular cartesian, #Components=3 + 1. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 3. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 4) Material, field, rectangular cartesian, #Components=2 + 1. constant*constant, no modify, grid based. + #xi1=0, #xi2=0 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + Element: 1 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 1 2 5 6 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 2 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 2 3 6 7 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 3 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 3 4 7 8 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 4 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 5 6 9 10 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 5 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 6 7 10 11 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 6 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 7 8 11 12 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 7 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 9 10 13 14 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 8 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 10 11 14 15 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 9 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 11 12 15 16 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 diff --git a/src/python/expected_results/output/MainTime_2.part0.exnode b/src/python/expected_results/output/MainTime_2.part0.exnode new file mode 100644 index 0000000..b41806b --- /dev/null +++ b/src/python/expected_results/output/MainTime_2.part0.exnode @@ -0,0 +1,192 @@ + Group name: StokesRegion + #Fields=4 + 1) Coordinate, coordinate, rectangular cartesian, #Components=2 + x. Value index= 1, #Derivatives= 0 + y. Value index= 2, #Derivatives= 0 + 2) U, field, rectangular cartesian, #Components=3 + 1. Value index= 3, #Derivatives= 0 + 2. Value index= 4, #Derivatives= 0 + 3. Value index= 5, #Derivatives= 0 + 3) del U_del n, field, rectangular cartesian, #Components=3 + 1. Value index= 6, #Derivatives= 0 + 2. Value index= 7, #Derivatives= 0 + 3. Value index= 8, #Derivatives= 0 + 4) Material, field, rectangular cartesian, #Components=2 + 1. Value index= 9, #Derivatives= 0 + 2. Value index= 10, #Derivatives= 0 + Node: 1 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0901513302412457E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 2 + 3.3333333333333331E-01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 2.7592262631192661E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 3 + 6.6666666666666663E-01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 2.7592262631213170E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 4 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0901513302425441E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 5 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 2.4126556016574458E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 6 + 3.3333333333333331E-01 + 1.0000000000000000E+00 + -1.8430868678120357E-12 + 1.0000000000002169E+00 + 1.9402123504978896E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 7 + 6.6666666666666663E-01 + 1.0000000000000000E+00 + -1.8401991541106519E-12 + 1.0000000000004894E+00 + 1.9402123505016444E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 8 + 1.0000000000000000E+00 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 2.4126556016613819E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 9 + 0.0000000000000000E+00 + 2.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 4.9713204783941665E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 10 + 3.3333333333333331E-01 + 2.0000000000000000E+00 + 4.9623973825973567E-13 + 9.9999999999908995E-01 + 7.1783011960081957E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 11 + 6.6666666666666663E-01 + 2.0000000000000000E+00 + 8.8221704527868547E-13 + 1.0000000000001938E+00 + 7.1783011959883742E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 12 + 1.0000000000000000E+00 + 2.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 4.9713204783881153E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 13 + 0.0000000000000000E+00 + 3.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 5.6091042226330501E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 14 + 3.3333333333333331E-01 + 3.0000000000000000E+00 + 2.2973527166178013E-12 + 1.0000000000012448E+00 + 1.5056138638196117E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 15 + 6.6666666666666663E-01 + 3.0000000000000000E+00 + 1.3845717500201821E-12 + 9.9999999999923483E-01 + 1.5056138637839029E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 16 + 1.0000000000000000E+00 + 3.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 5.6091042226032224E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 diff --git a/src/python/expected_results/output/MainTime_3.part0.exelem b/src/python/expected_results/output/MainTime_3.part0.exelem new file mode 100644 index 0000000..1eab4e3 --- /dev/null +++ b/src/python/expected_results/output/MainTime_3.part0.exelem @@ -0,0 +1,201 @@ + Group name: StokesRegion + Shape. Dimension=2 + #Scale factor sets= 1 + l.Lagrange*l.Lagrange, #Scale factors=4 + #Nodes= 4 + #Fields=4 + 1) Coordinate, coordinate, rectangular cartesian, #Components=2 + x. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + y. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 2) U, field, rectangular cartesian, #Components=3 + 1. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 3. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 3) del U_del n, field, rectangular cartesian, #Components=3 + 1. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 3. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 4) Material, field, rectangular cartesian, #Components=2 + 1. constant*constant, no modify, grid based. + #xi1=0, #xi2=0 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + Element: 1 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 1 2 5 6 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 2 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 2 3 6 7 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 3 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 3 4 7 8 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 4 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 5 6 9 10 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 5 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 6 7 10 11 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 6 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 7 8 11 12 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 7 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 9 10 13 14 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 8 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 10 11 14 15 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 9 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 11 12 15 16 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 diff --git a/src/python/expected_results/output/MainTime_3.part0.exnode b/src/python/expected_results/output/MainTime_3.part0.exnode new file mode 100644 index 0000000..b41806b --- /dev/null +++ b/src/python/expected_results/output/MainTime_3.part0.exnode @@ -0,0 +1,192 @@ + Group name: StokesRegion + #Fields=4 + 1) Coordinate, coordinate, rectangular cartesian, #Components=2 + x. Value index= 1, #Derivatives= 0 + y. Value index= 2, #Derivatives= 0 + 2) U, field, rectangular cartesian, #Components=3 + 1. Value index= 3, #Derivatives= 0 + 2. Value index= 4, #Derivatives= 0 + 3. Value index= 5, #Derivatives= 0 + 3) del U_del n, field, rectangular cartesian, #Components=3 + 1. Value index= 6, #Derivatives= 0 + 2. Value index= 7, #Derivatives= 0 + 3. Value index= 8, #Derivatives= 0 + 4) Material, field, rectangular cartesian, #Components=2 + 1. Value index= 9, #Derivatives= 0 + 2. Value index= 10, #Derivatives= 0 + Node: 1 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0901513302412457E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 2 + 3.3333333333333331E-01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 2.7592262631192661E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 3 + 6.6666666666666663E-01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 2.7592262631213170E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 4 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0901513302425441E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 5 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 2.4126556016574458E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 6 + 3.3333333333333331E-01 + 1.0000000000000000E+00 + -1.8430868678120357E-12 + 1.0000000000002169E+00 + 1.9402123504978896E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 7 + 6.6666666666666663E-01 + 1.0000000000000000E+00 + -1.8401991541106519E-12 + 1.0000000000004894E+00 + 1.9402123505016444E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 8 + 1.0000000000000000E+00 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 2.4126556016613819E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 9 + 0.0000000000000000E+00 + 2.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 4.9713204783941665E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 10 + 3.3333333333333331E-01 + 2.0000000000000000E+00 + 4.9623973825973567E-13 + 9.9999999999908995E-01 + 7.1783011960081957E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 11 + 6.6666666666666663E-01 + 2.0000000000000000E+00 + 8.8221704527868547E-13 + 1.0000000000001938E+00 + 7.1783011959883742E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 12 + 1.0000000000000000E+00 + 2.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 4.9713204783881153E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 13 + 0.0000000000000000E+00 + 3.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 5.6091042226330501E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 14 + 3.3333333333333331E-01 + 3.0000000000000000E+00 + 2.2973527166178013E-12 + 1.0000000000012448E+00 + 1.5056138638196117E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 15 + 6.6666666666666663E-01 + 3.0000000000000000E+00 + 1.3845717500201821E-12 + 9.9999999999923483E-01 + 1.5056138637839029E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 16 + 1.0000000000000000E+00 + 3.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 5.6091042226032224E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 diff --git a/src/python/expected_results/output/MainTime_4.part0.exelem b/src/python/expected_results/output/MainTime_4.part0.exelem new file mode 100644 index 0000000..1eab4e3 --- /dev/null +++ b/src/python/expected_results/output/MainTime_4.part0.exelem @@ -0,0 +1,201 @@ + Group name: StokesRegion + Shape. Dimension=2 + #Scale factor sets= 1 + l.Lagrange*l.Lagrange, #Scale factors=4 + #Nodes= 4 + #Fields=4 + 1) Coordinate, coordinate, rectangular cartesian, #Components=2 + x. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + y. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 2) U, field, rectangular cartesian, #Components=3 + 1. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 3. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 3) del U_del n, field, rectangular cartesian, #Components=3 + 1. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 3. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 4) Material, field, rectangular cartesian, #Components=2 + 1. constant*constant, no modify, grid based. + #xi1=0, #xi2=0 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + Element: 1 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 1 2 5 6 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 2 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 2 3 6 7 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 3 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 3 4 7 8 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 4 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 5 6 9 10 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 5 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 6 7 10 11 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 6 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 7 8 11 12 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 7 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 9 10 13 14 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 8 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 10 11 14 15 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 9 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 11 12 15 16 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 diff --git a/src/python/expected_results/output/MainTime_4.part0.exnode b/src/python/expected_results/output/MainTime_4.part0.exnode new file mode 100644 index 0000000..b41806b --- /dev/null +++ b/src/python/expected_results/output/MainTime_4.part0.exnode @@ -0,0 +1,192 @@ + Group name: StokesRegion + #Fields=4 + 1) Coordinate, coordinate, rectangular cartesian, #Components=2 + x. Value index= 1, #Derivatives= 0 + y. Value index= 2, #Derivatives= 0 + 2) U, field, rectangular cartesian, #Components=3 + 1. Value index= 3, #Derivatives= 0 + 2. Value index= 4, #Derivatives= 0 + 3. Value index= 5, #Derivatives= 0 + 3) del U_del n, field, rectangular cartesian, #Components=3 + 1. Value index= 6, #Derivatives= 0 + 2. Value index= 7, #Derivatives= 0 + 3. Value index= 8, #Derivatives= 0 + 4) Material, field, rectangular cartesian, #Components=2 + 1. Value index= 9, #Derivatives= 0 + 2. Value index= 10, #Derivatives= 0 + Node: 1 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0901513302412457E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 2 + 3.3333333333333331E-01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 2.7592262631192661E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 3 + 6.6666666666666663E-01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 2.7592262631213170E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 4 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0901513302425441E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 5 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 2.4126556016574458E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 6 + 3.3333333333333331E-01 + 1.0000000000000000E+00 + -1.8430868678120357E-12 + 1.0000000000002169E+00 + 1.9402123504978896E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 7 + 6.6666666666666663E-01 + 1.0000000000000000E+00 + -1.8401991541106519E-12 + 1.0000000000004894E+00 + 1.9402123505016444E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 8 + 1.0000000000000000E+00 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 2.4126556016613819E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 9 + 0.0000000000000000E+00 + 2.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 4.9713204783941665E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 10 + 3.3333333333333331E-01 + 2.0000000000000000E+00 + 4.9623973825973567E-13 + 9.9999999999908995E-01 + 7.1783011960081957E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 11 + 6.6666666666666663E-01 + 2.0000000000000000E+00 + 8.8221704527868547E-13 + 1.0000000000001938E+00 + 7.1783011959883742E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 12 + 1.0000000000000000E+00 + 2.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 4.9713204783881153E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 13 + 0.0000000000000000E+00 + 3.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 5.6091042226330501E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 14 + 3.3333333333333331E-01 + 3.0000000000000000E+00 + 2.2973527166178013E-12 + 1.0000000000012448E+00 + 1.5056138638196117E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 15 + 6.6666666666666663E-01 + 3.0000000000000000E+00 + 1.3845717500201821E-12 + 9.9999999999923483E-01 + 1.5056138637839029E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 16 + 1.0000000000000000E+00 + 3.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 5.6091042226032224E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 diff --git a/src/python/output/MainTime_0.part0.exelem b/src/python/output/MainTime_0.part0.exelem new file mode 100644 index 0000000..1eab4e3 --- /dev/null +++ b/src/python/output/MainTime_0.part0.exelem @@ -0,0 +1,201 @@ + Group name: StokesRegion + Shape. Dimension=2 + #Scale factor sets= 1 + l.Lagrange*l.Lagrange, #Scale factors=4 + #Nodes= 4 + #Fields=4 + 1) Coordinate, coordinate, rectangular cartesian, #Components=2 + x. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + y. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 2) U, field, rectangular cartesian, #Components=3 + 1. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 3. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 3) del U_del n, field, rectangular cartesian, #Components=3 + 1. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 3. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 4) Material, field, rectangular cartesian, #Components=2 + 1. constant*constant, no modify, grid based. + #xi1=0, #xi2=0 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + Element: 1 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 1 2 5 6 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 2 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 2 3 6 7 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 3 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 3 4 7 8 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 4 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 5 6 9 10 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 5 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 6 7 10 11 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 6 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 7 8 11 12 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 7 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 9 10 13 14 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 8 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 10 11 14 15 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 9 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 11 12 15 16 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 diff --git a/src/python/output/MainTime_0.part0.exnode b/src/python/output/MainTime_0.part0.exnode new file mode 100644 index 0000000..a1321da --- /dev/null +++ b/src/python/output/MainTime_0.part0.exnode @@ -0,0 +1,192 @@ + Group name: StokesRegion + #Fields=4 + 1) Coordinate, coordinate, rectangular cartesian, #Components=2 + x. Value index= 1, #Derivatives= 0 + y. Value index= 2, #Derivatives= 0 + 2) U, field, rectangular cartesian, #Components=3 + 1. Value index= 3, #Derivatives= 0 + 2. Value index= 4, #Derivatives= 0 + 3. Value index= 5, #Derivatives= 0 + 3) del U_del n, field, rectangular cartesian, #Components=3 + 1. Value index= 6, #Derivatives= 0 + 2. Value index= 7, #Derivatives= 0 + 3. Value index= 8, #Derivatives= 0 + 4) Material, field, rectangular cartesian, #Components=2 + 1. Value index= 9, #Derivatives= 0 + 2. Value index= 10, #Derivatives= 0 + Node: 1 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 2 + 3.3333333333333331E-01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 3 + 6.6666666666666663E-01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 4 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 5 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 6 + 3.3333333333333331E-01 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 7 + 6.6666666666666663E-01 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 8 + 1.0000000000000000E+00 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 9 + 0.0000000000000000E+00 + 2.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 10 + 3.3333333333333331E-01 + 2.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 11 + 6.6666666666666663E-01 + 2.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 12 + 1.0000000000000000E+00 + 2.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 13 + 0.0000000000000000E+00 + 3.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 14 + 3.3333333333333331E-01 + 3.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 15 + 6.6666666666666663E-01 + 3.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 16 + 1.0000000000000000E+00 + 3.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 diff --git a/src/python/output/MainTime_1.part0.exelem b/src/python/output/MainTime_1.part0.exelem new file mode 100644 index 0000000..1eab4e3 --- /dev/null +++ b/src/python/output/MainTime_1.part0.exelem @@ -0,0 +1,201 @@ + Group name: StokesRegion + Shape. Dimension=2 + #Scale factor sets= 1 + l.Lagrange*l.Lagrange, #Scale factors=4 + #Nodes= 4 + #Fields=4 + 1) Coordinate, coordinate, rectangular cartesian, #Components=2 + x. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + y. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 2) U, field, rectangular cartesian, #Components=3 + 1. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 3. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 3) del U_del n, field, rectangular cartesian, #Components=3 + 1. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 3. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 4) Material, field, rectangular cartesian, #Components=2 + 1. constant*constant, no modify, grid based. + #xi1=0, #xi2=0 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + Element: 1 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 1 2 5 6 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 2 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 2 3 6 7 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 3 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 3 4 7 8 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 4 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 5 6 9 10 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 5 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 6 7 10 11 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 6 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 7 8 11 12 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 7 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 9 10 13 14 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 8 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 10 11 14 15 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 9 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 11 12 15 16 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 diff --git a/src/python/output/MainTime_1.part0.exnode b/src/python/output/MainTime_1.part0.exnode new file mode 100644 index 0000000..297c81c --- /dev/null +++ b/src/python/output/MainTime_1.part0.exnode @@ -0,0 +1,192 @@ + Group name: StokesRegion + #Fields=4 + 1) Coordinate, coordinate, rectangular cartesian, #Components=2 + x. Value index= 1, #Derivatives= 0 + y. Value index= 2, #Derivatives= 0 + 2) U, field, rectangular cartesian, #Components=3 + 1. Value index= 3, #Derivatives= 0 + 2. Value index= 4, #Derivatives= 0 + 3. Value index= 5, #Derivatives= 0 + 3) del U_del n, field, rectangular cartesian, #Components=3 + 1. Value index= 6, #Derivatives= 0 + 2. Value index= 7, #Derivatives= 0 + 3. Value index= 8, #Derivatives= 0 + 4) Material, field, rectangular cartesian, #Components=2 + 1. Value index= 9, #Derivatives= 0 + 2. Value index= 10, #Derivatives= 0 + Node: 1 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.1812122013393381E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 2 + 3.3333333333333331E-01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 2.9846130637595071E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 3 + 6.6666666666666663E-01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 2.9846130637590665E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 4 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.1812122013388382E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 5 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 2.6286825726135184E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 6 + 3.3333333333333331E-01 + 1.0000000000000000E+00 + -9.7908019770077493E-12 + 9.9999999999979972E-01 + 2.1200746481192425E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 7 + 6.6666666666666663E-01 + 1.0000000000000000E+00 + -8.0700587257565004E-12 + 9.9999999999738842E-01 + 2.1200746481177053E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 8 + 1.0000000000000000E+00 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 2.6286825726120352E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 9 + 0.0000000000000000E+00 + 2.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 5.5263166816422791E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 10 + 3.3333333333333331E-01 + 2.0000000000000000E+00 + 1.8382265636428370E-11 + 1.0000000000002103E+00 + 7.8366250373066046E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 11 + 6.6666666666666663E-01 + 2.0000000000000000E+00 + 1.5851453116772888E-11 + 9.9999999999993827E-01 + 7.8366250373119266E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 12 + 1.0000000000000000E+00 + 2.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 5.5263166816489182E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 13 + 0.0000000000000000E+00 + 3.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 5.8443881702858924E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 14 + 3.3333333333333331E-01 + 3.0000000000000000E+00 + -1.1511640515858268E-11 + 1.0000000000005844E+00 + 1.6892339923515256E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 15 + 6.6666666666666663E-01 + 3.0000000000000000E+00 + -1.5510357975943329E-11 + 1.0000000000060576E+00 + 1.6892339923398016E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 16 + 1.0000000000000000E+00 + 3.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 5.8443881702755824E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 diff --git a/src/python/output/MainTime_2.part0.exelem b/src/python/output/MainTime_2.part0.exelem new file mode 100644 index 0000000..1eab4e3 --- /dev/null +++ b/src/python/output/MainTime_2.part0.exelem @@ -0,0 +1,201 @@ + Group name: StokesRegion + Shape. Dimension=2 + #Scale factor sets= 1 + l.Lagrange*l.Lagrange, #Scale factors=4 + #Nodes= 4 + #Fields=4 + 1) Coordinate, coordinate, rectangular cartesian, #Components=2 + x. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + y. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 2) U, field, rectangular cartesian, #Components=3 + 1. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 3. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 3) del U_del n, field, rectangular cartesian, #Components=3 + 1. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 3. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 4) Material, field, rectangular cartesian, #Components=2 + 1. constant*constant, no modify, grid based. + #xi1=0, #xi2=0 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + Element: 1 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 1 2 5 6 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 2 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 2 3 6 7 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 3 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 3 4 7 8 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 4 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 5 6 9 10 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 5 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 6 7 10 11 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 6 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 7 8 11 12 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 7 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 9 10 13 14 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 8 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 10 11 14 15 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 9 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 11 12 15 16 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 diff --git a/src/python/output/MainTime_2.part0.exnode b/src/python/output/MainTime_2.part0.exnode new file mode 100644 index 0000000..b41806b --- /dev/null +++ b/src/python/output/MainTime_2.part0.exnode @@ -0,0 +1,192 @@ + Group name: StokesRegion + #Fields=4 + 1) Coordinate, coordinate, rectangular cartesian, #Components=2 + x. Value index= 1, #Derivatives= 0 + y. Value index= 2, #Derivatives= 0 + 2) U, field, rectangular cartesian, #Components=3 + 1. Value index= 3, #Derivatives= 0 + 2. Value index= 4, #Derivatives= 0 + 3. Value index= 5, #Derivatives= 0 + 3) del U_del n, field, rectangular cartesian, #Components=3 + 1. Value index= 6, #Derivatives= 0 + 2. Value index= 7, #Derivatives= 0 + 3. Value index= 8, #Derivatives= 0 + 4) Material, field, rectangular cartesian, #Components=2 + 1. Value index= 9, #Derivatives= 0 + 2. Value index= 10, #Derivatives= 0 + Node: 1 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0901513302412457E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 2 + 3.3333333333333331E-01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 2.7592262631192661E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 3 + 6.6666666666666663E-01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 2.7592262631213170E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 4 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0901513302425441E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 5 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 2.4126556016574458E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 6 + 3.3333333333333331E-01 + 1.0000000000000000E+00 + -1.8430868678120357E-12 + 1.0000000000002169E+00 + 1.9402123504978896E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 7 + 6.6666666666666663E-01 + 1.0000000000000000E+00 + -1.8401991541106519E-12 + 1.0000000000004894E+00 + 1.9402123505016444E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 8 + 1.0000000000000000E+00 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 2.4126556016613819E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 9 + 0.0000000000000000E+00 + 2.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 4.9713204783941665E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 10 + 3.3333333333333331E-01 + 2.0000000000000000E+00 + 4.9623973825973567E-13 + 9.9999999999908995E-01 + 7.1783011960081957E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 11 + 6.6666666666666663E-01 + 2.0000000000000000E+00 + 8.8221704527868547E-13 + 1.0000000000001938E+00 + 7.1783011959883742E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 12 + 1.0000000000000000E+00 + 2.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 4.9713204783881153E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 13 + 0.0000000000000000E+00 + 3.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 5.6091042226330501E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 14 + 3.3333333333333331E-01 + 3.0000000000000000E+00 + 2.2973527166178013E-12 + 1.0000000000012448E+00 + 1.5056138638196117E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 15 + 6.6666666666666663E-01 + 3.0000000000000000E+00 + 1.3845717500201821E-12 + 9.9999999999923483E-01 + 1.5056138637839029E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 16 + 1.0000000000000000E+00 + 3.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 5.6091042226032224E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 diff --git a/src/python/output/MainTime_3.part0.exelem b/src/python/output/MainTime_3.part0.exelem new file mode 100644 index 0000000..1eab4e3 --- /dev/null +++ b/src/python/output/MainTime_3.part0.exelem @@ -0,0 +1,201 @@ + Group name: StokesRegion + Shape. Dimension=2 + #Scale factor sets= 1 + l.Lagrange*l.Lagrange, #Scale factors=4 + #Nodes= 4 + #Fields=4 + 1) Coordinate, coordinate, rectangular cartesian, #Components=2 + x. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + y. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 2) U, field, rectangular cartesian, #Components=3 + 1. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 3. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 3) del U_del n, field, rectangular cartesian, #Components=3 + 1. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 3. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 4) Material, field, rectangular cartesian, #Components=2 + 1. constant*constant, no modify, grid based. + #xi1=0, #xi2=0 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + Element: 1 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 1 2 5 6 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 2 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 2 3 6 7 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 3 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 3 4 7 8 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 4 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 5 6 9 10 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 5 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 6 7 10 11 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 6 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 7 8 11 12 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 7 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 9 10 13 14 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 8 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 10 11 14 15 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 9 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 11 12 15 16 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 diff --git a/src/python/output/MainTime_3.part0.exnode b/src/python/output/MainTime_3.part0.exnode new file mode 100644 index 0000000..b41806b --- /dev/null +++ b/src/python/output/MainTime_3.part0.exnode @@ -0,0 +1,192 @@ + Group name: StokesRegion + #Fields=4 + 1) Coordinate, coordinate, rectangular cartesian, #Components=2 + x. Value index= 1, #Derivatives= 0 + y. Value index= 2, #Derivatives= 0 + 2) U, field, rectangular cartesian, #Components=3 + 1. Value index= 3, #Derivatives= 0 + 2. Value index= 4, #Derivatives= 0 + 3. Value index= 5, #Derivatives= 0 + 3) del U_del n, field, rectangular cartesian, #Components=3 + 1. Value index= 6, #Derivatives= 0 + 2. Value index= 7, #Derivatives= 0 + 3. Value index= 8, #Derivatives= 0 + 4) Material, field, rectangular cartesian, #Components=2 + 1. Value index= 9, #Derivatives= 0 + 2. Value index= 10, #Derivatives= 0 + Node: 1 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0901513302412457E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 2 + 3.3333333333333331E-01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 2.7592262631192661E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 3 + 6.6666666666666663E-01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 2.7592262631213170E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 4 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0901513302425441E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 5 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 2.4126556016574458E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 6 + 3.3333333333333331E-01 + 1.0000000000000000E+00 + -1.8430868678120357E-12 + 1.0000000000002169E+00 + 1.9402123504978896E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 7 + 6.6666666666666663E-01 + 1.0000000000000000E+00 + -1.8401991541106519E-12 + 1.0000000000004894E+00 + 1.9402123505016444E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 8 + 1.0000000000000000E+00 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 2.4126556016613819E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 9 + 0.0000000000000000E+00 + 2.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 4.9713204783941665E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 10 + 3.3333333333333331E-01 + 2.0000000000000000E+00 + 4.9623973825973567E-13 + 9.9999999999908995E-01 + 7.1783011960081957E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 11 + 6.6666666666666663E-01 + 2.0000000000000000E+00 + 8.8221704527868547E-13 + 1.0000000000001938E+00 + 7.1783011959883742E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 12 + 1.0000000000000000E+00 + 2.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 4.9713204783881153E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 13 + 0.0000000000000000E+00 + 3.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 5.6091042226330501E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 14 + 3.3333333333333331E-01 + 3.0000000000000000E+00 + 2.2973527166178013E-12 + 1.0000000000012448E+00 + 1.5056138638196117E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 15 + 6.6666666666666663E-01 + 3.0000000000000000E+00 + 1.3845717500201821E-12 + 9.9999999999923483E-01 + 1.5056138637839029E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 16 + 1.0000000000000000E+00 + 3.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 5.6091042226032224E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 diff --git a/src/python/output/MainTime_4.part0.exelem b/src/python/output/MainTime_4.part0.exelem new file mode 100644 index 0000000..1eab4e3 --- /dev/null +++ b/src/python/output/MainTime_4.part0.exelem @@ -0,0 +1,201 @@ + Group name: StokesRegion + Shape. Dimension=2 + #Scale factor sets= 1 + l.Lagrange*l.Lagrange, #Scale factors=4 + #Nodes= 4 + #Fields=4 + 1) Coordinate, coordinate, rectangular cartesian, #Components=2 + x. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + y. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 2) U, field, rectangular cartesian, #Components=3 + 1. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 3. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 3) del U_del n, field, rectangular cartesian, #Components=3 + 1. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 3. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + 4) Material, field, rectangular cartesian, #Components=2 + 1. constant*constant, no modify, grid based. + #xi1=0, #xi2=0 + 2. l.Lagrange*l.Lagrange, no modify, standard node based. + #Nodes= 4 + 1. #Values=1 + Value indices: 1 + Scale factor indices: 1 + 2. #Values=1 + Value indices: 1 + Scale factor indices: 2 + 3. #Values=1 + Value indices: 1 + Scale factor indices: 3 + 4. #Values=1 + Value indices: 1 + Scale factor indices: 4 + Element: 1 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 1 2 5 6 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 2 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 2 3 6 7 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 3 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 3 4 7 8 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 4 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 5 6 9 10 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 5 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 6 7 10 11 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 6 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 7 8 11 12 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 7 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 9 10 13 14 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 8 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 10 11 14 15 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 + Element: 9 0 0 + Values: + 1.0000000000000000E+00 + Nodes: + 11 12 15 16 + Scale factors: + 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 1.0000000000000000E+00 diff --git a/src/python/output/MainTime_4.part0.exnode b/src/python/output/MainTime_4.part0.exnode new file mode 100644 index 0000000..b41806b --- /dev/null +++ b/src/python/output/MainTime_4.part0.exnode @@ -0,0 +1,192 @@ + Group name: StokesRegion + #Fields=4 + 1) Coordinate, coordinate, rectangular cartesian, #Components=2 + x. Value index= 1, #Derivatives= 0 + y. Value index= 2, #Derivatives= 0 + 2) U, field, rectangular cartesian, #Components=3 + 1. Value index= 3, #Derivatives= 0 + 2. Value index= 4, #Derivatives= 0 + 3. Value index= 5, #Derivatives= 0 + 3) del U_del n, field, rectangular cartesian, #Components=3 + 1. Value index= 6, #Derivatives= 0 + 2. Value index= 7, #Derivatives= 0 + 3. Value index= 8, #Derivatives= 0 + 4) Material, field, rectangular cartesian, #Components=2 + 1. Value index= 9, #Derivatives= 0 + 2. Value index= 10, #Derivatives= 0 + Node: 1 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0901513302412457E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 2 + 3.3333333333333331E-01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 2.7592262631192661E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 3 + 6.6666666666666663E-01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 2.7592262631213170E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 4 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.0901513302425441E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 5 + 0.0000000000000000E+00 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 2.4126556016574458E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 6 + 3.3333333333333331E-01 + 1.0000000000000000E+00 + -1.8430868678120357E-12 + 1.0000000000002169E+00 + 1.9402123504978896E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 7 + 6.6666666666666663E-01 + 1.0000000000000000E+00 + -1.8401991541106519E-12 + 1.0000000000004894E+00 + 1.9402123505016444E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 8 + 1.0000000000000000E+00 + 1.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 2.4126556016613819E+01 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 9 + 0.0000000000000000E+00 + 2.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 4.9713204783941665E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 10 + 3.3333333333333331E-01 + 2.0000000000000000E+00 + 4.9623973825973567E-13 + 9.9999999999908995E-01 + 7.1783011960081957E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 11 + 6.6666666666666663E-01 + 2.0000000000000000E+00 + 8.8221704527868547E-13 + 1.0000000000001938E+00 + 7.1783011959883742E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 12 + 1.0000000000000000E+00 + 2.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 4.9713204783881153E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 13 + 0.0000000000000000E+00 + 3.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 5.6091042226330501E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 14 + 3.3333333333333331E-01 + 3.0000000000000000E+00 + 2.2973527166178013E-12 + 1.0000000000012448E+00 + 1.5056138638196117E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 15 + 6.6666666666666663E-01 + 3.0000000000000000E+00 + 1.3845717500201821E-12 + 9.9999999999923483E-01 + 1.5056138637839029E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 + Node: 16 + 1.0000000000000000E+00 + 3.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 5.6091042226032224E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 0.0000000000000000E+00 + 1.2345678806304932E+00 + 1.0000000000000000E+00 diff --git a/src/python/stokes_dynamic.py b/src/python/stokes_dynamic.py new file mode 100644 index 0000000..240f31f --- /dev/null +++ b/src/python/stokes_dynamic.py @@ -0,0 +1,353 @@ +######################################################### +# 2D Dynamic Stokes FLow +# 20180529 +######################################################### + +# In[1]: + +# Import the libraries (OpenCMISS,python,numpy,scipy) +import numpy,math,cmath,csv,time,sys,os,pdb +from scipy.fftpack import fft,ifft +from scipy.sparse import linalg +from scipy.linalg import inv,eig +from scipy.special import jn +from opencmiss.iron import iron + +#------------------------------ +#Set problem parameters +#------------------------------ +height = 3.0 +width = 1.0 +length = 1.0 + +mu=1.0 +rho=1.0 +initial_conc=0.0 +screen_output_freq=2 + +CoordinateSystemUserNumber=1 +RegionUserNumber=2 +MeshUserNumber=3 +DecompositionUserNumber=4 +GeometricFieldUserNumber=5 +EquationsSetFieldUserNumber=6 +DependentFieldUserNumberStokes=7 +MaterialsFieldUserNumberStokes=8 +IndependentFieldUserNumberStokes=9 +EquationsSetUserNumberStokes=10 +ProblemUserNumber=11 +GeneratedMeshUserNumber=12 +AnalyticFieldUserNumber=13 +BasisUserNumber=14 + +DomainUserNumber=1 +SolverStokesUserNumber=1 +MaterialsFieldUserNumberStokesMu=1 +MaterialsFieldUserNumberStokesRho=2 + +numberGlobalXElements = 3 +numberGlobalYElements = 3 +numberGlobalZElements = 3 + +#boundary condition nodes +InletWallNodes=[2,3] +FixedWallNodes=[1,4,5,8,9,12,13,16] + +DYNAMIC_SOLVER_STOKES_OUTPUT_FREQUENCY = 1 + +timeIncrement = 1.0 +startTime = 0.0 +stopTime = 5.0 +dynamicSolverStokesTheta = [1.0] + +#--------------------------------- +# Get the computational nodes info +#--------------------------------- +numberOfComputationalNodes = iron.ComputationalNumberOfNodesGet() +computationalNodeNumber = iron.ComputationalNodeNumberGet() + +#------------------ +#Coordinate system +#------------------ +coordinateSystem = iron.CoordinateSystem() +coordinateSystem.CreateStart(CoordinateSystemUserNumber) +coordinateSystem.dimension = 2 +coordinateSystem.CreateFinish() + +#------------------ +#Regin +#------------------ +# Create a region +region = iron.Region() +region.CreateStart(RegionUserNumber,iron.WorldRegion) +region.label = "StokesRegion" +region.coordinateSystem = coordinateSystem +region.CreateFinish() + + +# In[2]: + +#-------- +#Basis +#-------- +# Create a tri-linear lagrange basis +basis = iron.Basis() +basis.CreateStart(BasisUserNumber) +basis.type = iron.BasisTypes.LAGRANGE_HERMITE_TP +basis.numberOfXi = 2 +basis.interpolationXi = [iron.BasisInterpolationSpecifications.LINEAR_LAGRANGE]*2 +basis.quadratureNumberOfGaussXi = [2]*2 +#basis.quadratureLocalFaceGaussEvaluate = True +basis.CreateFinish() + + +# In[3]: + +#---------- +#Mesh +#---------- +# Create a generated mesh +generatedMesh = iron.GeneratedMesh() +generatedMesh.CreateStart(GeneratedMeshUserNumber,region) +generatedMesh.type = iron.GeneratedMeshTypes.REGULAR +generatedMesh.basis = [basis] +generatedMesh.extent = [width,height] +generatedMesh.numberOfElements = [numberGlobalXElements,numberGlobalYElements] + +mesh = iron.Mesh() +generatedMesh.CreateFinish(MeshUserNumber,mesh) + +numberOfElements = mesh.numberOfElements +print("number of elements: " + str(numberOfElements)) + + +# In[4]: + +#------------------ +#Mesh decomposition +#------------------ +decomposition = iron.Decomposition() +decomposition.CreateStart(DecompositionUserNumber,mesh) +decomposition.type = iron.DecompositionTypes.CALCULATED +decomposition.numberOfDomains = numberOfComputationalNodes +decomposition.CreateFinish() + + +# In[5]: + +#-------------- +#Geometry field +#-------------- + +#Create a field for the geometry +geometricField = iron.Field() +geometricField.CreateStart(GeometricFieldUserNumber,region) +geometricField.meshDecomposition = decomposition +geometricField.ComponentMeshComponentSet(iron.FieldVariableTypes.U,1,1) +geometricField.ComponentMeshComponentSet(iron.FieldVariableTypes.U,2,1) +#geometricField.ComponentMeshComponentSet(iron.FieldVariableTypes.U,3,1) +geometricField.CreateFinish() + +# Set geometry from the generated mesh +generatedMesh.GeometricParametersCalculate(geometricField) + + +# In[6]: + +#------------- +#Equation sets +#------------- + +#Create Stokes equations set +equationsSetField = iron.Field() +equationsSet = iron.EquationsSet() +equationsSetSpecification = [iron.EquationsSetClasses.FLUID_MECHANICS, + iron.EquationsSetTypes.STOKES_EQUATION, + iron.EquationsSetSubtypes.TRANSIENT_STOKES] +equationsSet.CreateStart(EquationsSetUserNumberStokes,region,geometricField, + equationsSetSpecification,EquationsSetFieldUserNumber,equationsSetField) +equationsSet.CreateFinish() + + +# In[7]: + +#--------------- +#Dependent field +#--------------- +# Create dependent field +dependentField = iron.Field() +equationsSet.DependentCreateStart(DependentFieldUserNumberStokes,dependentField) +dependentField.VariableLabelSet(iron.FieldVariableTypes.U,'U') +dependentField.ComponentMeshComponentSet(iron.FieldVariableTypes.U,1,1) +dependentField.ComponentMeshComponentSet(iron.FieldVariableTypes.U,2,1) +dependentField.ComponentMeshComponentSet(iron.FieldVariableTypes.U,3,1) +equationsSet.DependentCreateFinish() + + +# In[8]: + +#-------------- +#Material field +#-------------- +# Create material field +materialField = iron.Field() +equationsSet.MaterialsCreateStart(MaterialsFieldUserNumberStokes,materialField) +materialField.VariableLabelSet(iron.FieldVariableTypes.U, "Material") +equationsSet.MaterialsCreateFinish() + +# Mu and Rho? +materialField.ComponentValuesInitialiseDP(iron.FieldVariableTypes.U,iron.FieldParameterSetTypes.VALUES,1,mu) +materialField.ComponentValuesInitialiseDP(iron.FieldVariableTypes.U,iron.FieldParameterSetTypes.VALUES,2,rho) + +# Initialise dependent field +dependentField.ComponentValuesInitialiseDP(iron.FieldVariableTypes.U,iron.FieldParameterSetTypes.VALUES,1,initial_conc) + + +# In[9]: + +#----------- +# EQUATIONS +#----------- +# Create equations +equations = iron.Equations() +equationsSet.EquationsCreateStart(equations) +equations.sparsityType = iron.EquationsSparsityTypes.SPARSE +equations.lumpingType = iron.EquationsLumpingTypes.UNLUMPED +equations.outputType = iron.EquationsOutputTypes.NONE +equationsSet.EquationsCreateFinish() + + +# In[10]: + +#----------------------------------------------------------------------------------------------------------- +#PROBLEM +#----------------------------------------------------------------------------------------------------------- + +# Create Diffusion equation problem +problem = iron.Problem() +problemSpecification = [iron.ProblemClasses.FLUID_MECHANICS, + iron.ProblemTypes.STOKES_EQUATION, + iron.ProblemSubtypes.TRANSIENT_STOKES] +problem.CreateStart(ProblemUserNumber, problemSpecification) +problem.CreateFinish() + +# Create control loops +#problem.ControlLoopCreateStart() +#problem.ControlLoopCreateFinish() +#=============== +# Control Loops +#=============== +TimeLoop = iron.ControlLoop() +problem.ControlLoopCreateStart() +problem.ControlLoopGet([iron.ControlLoopIdentifiers.NODE],TimeLoop) +TimeLoop.LabelSet('Time Loop') +TimeLoop.TimesSet(startTime,stopTime,timeIncrement) +TimeLoop.TimeOutputSet(DYNAMIC_SOLVER_STOKES_OUTPUT_FREQUENCY) +TimeLoop.OutputTypeSet(iron.ControlLoopOutputTypes.NONE) +problem.ControlLoopCreateFinish() + +# In[11]: + +#------- +#SOLVER +#------- +DYNAMIC_SOLVER_STOKES_OUTPUT_TYPE = iron.SolverOutputTypes.NONE +LINEAR_SOLVER_STOKES_OUTPUT_TYPE = iron.SolverOutputTypes.NONE +# Create problem solver +solverLinear = iron.Solver() +solverDynamic = iron.Solver() +problem.SolversCreateStart() +problem.SolverGet([iron.ControlLoopIdentifiers.NODE],SolverStokesUserNumber,solverDynamic) +solverDynamic.OutputTypeSet(DYNAMIC_SOLVER_STOKES_OUTPUT_TYPE) +solverDynamic.DynamicThetaSet(dynamicSolverStokesTheta) +solverDynamic.DynamicLinearSolverGet(solverLinear) +solverLinear.outputType = iron.SolverOutputTypes.SOLVER +solverLinear.linearType = iron.LinearSolverTypes.ITERATIVE +solverLinear.linearIterativeAbsoluteTolerance = 1.0E-10 +solverLinear.linearIterativeRelativeTolerance = 1.0E-10 +problem.SolversCreateFinish() + + +# In[12]: + +#----------------------------------------------------------------------------------------------------------- +#SOLVER EQUATIONS +#----------------------------------------------------------------------------------------------------------- +## Create solver equations and add equations set to solver equations +solver = iron.Solver() +solverEquations = iron.SolverEquations() +problem.SolverEquationsCreateStart() +problem.SolverGet([iron.ControlLoopIdentifiers.NODE],1,solver) +solver.SolverEquationsGet(solverEquations) +solverEquations.sparsityType = iron.SolverEquationsSparsityTypes.SPARSE +equationsSetIndex = solverEquations.EquationsSetAdd(equationsSet) +problem.SolverEquationsCreateFinish() + + +# In[13]: + +#------------------- +#Boundary conditions +#------------------- +# Create boundary conditions and set first and last nodes to 0.0 and 1.0 +boundaryConditions = iron.BoundaryConditions() +solverEquations.BoundaryConditionsCreateStart(boundaryConditions) + +""" +firstNodeNumber=1 +nodes = iron.Nodes() +region.NodesGet(nodes) +lastNodeNumber = nodes.numberOfNodes +firstNodeDomain = decomposition.NodeDomainGet(firstNodeNumber,1) +lastNodeDomain = decomposition.NodeDomainGet(lastNodeNumber,1) +""" +#Set velocity boundary conditions +for i in InletWallNodes: + print i + boundaryConditions.SetNode(dependentField,iron.FieldVariableTypes.U,1,1,i,1,iron.BoundaryConditionsTypes.FIXED,0.0) + boundaryConditions.SetNode(dependentField,iron.FieldVariableTypes.U,1,1,i,2,iron.BoundaryConditionsTypes.FIXED,1.0) + #boundaryConditions.SetNode(dependentField,iron.FieldVariableTypes.U,1,1,i,3,iron.BoundaryConditionsTypes.FIXED,0.0) + +#Set fixed wall nodes +for i in FixedWallNodes: + print i + boundaryConditions.SetNode(dependentField,iron.FieldVariableTypes.U,1,1,i,1,iron.BoundaryConditionsTypes.FIXED_WALL,0.0) + boundaryConditions.SetNode(dependentField,iron.FieldVariableTypes.U,1,1,i,2,iron.BoundaryConditionsTypes.FIXED_WALL,0.0) + +solverEquations.BoundaryConditionsCreateFinish() + + +# In[ ]: +#----------------------------------------------------------------------------------------------------------- +#SOLVE +#----------------------------------------------------------------------------------------------------------- +#problem.Solve() +#++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ +# Solve the problem +print "Solving problem..." +start = time.time() +problem.Solve() +end = time.time() +elapsed = end - start +print "Calculation Time = %3.4f" %elapsed +print "Problem solved!" +print "#" + +#----------------------------------------------------------------------------------------------------------- +#OUTPUT +#----------------------------------------------------------------------------------------------------------- +# Ensure output directories exist +#if not os.path.exists('./output'): +# os.makedirs('./output') + +# Export results +#fields = iron.Fields() +#fields.CreateRegion(region) +#fields.NodesExport("output/DynamicStokes","FORTRAN") +#fields.ElementsExport("output/DynamicStokes","FORTRAN") +#fields.Finalise() + +iron.Finalise() + + +# In[ ]: diff --git a/src/python/visualiseStokes.com b/src/python/visualiseStokes.com new file mode 100755 index 0000000..62a98b5 --- /dev/null +++ b/src/python/visualiseStokes.com @@ -0,0 +1,47 @@ +$dt=1; # frequency +$steps=40000; # number of time steps +$cores=1; # number of processors + +$w=500; # width of the graphical window +$h=1000; # height of the graphical window + +#Read in the sequence of nodal positions. +for ($i=0;$i<$steps;$i=$i+$dt) + { + $time = $i/($dt*100) + for ($j=0;$j<$cores;$j=$j+1) + { + $filename = sprintf("./output/MainTime_%01d.part%01d.exnode", $i, $j); + print "Reading $filename time $time\n"; + gfx read node "$filename" time $i; + } + } + +#Read in the element description +for ($k=0;$k<$cores;$k=$k+1) + { + $filename = sprintf("./output/MainTime_0.part%01d.exelem", $k); + gfx read element "$filename"; + } + +gfx define faces egroup StokesRegion + +gfx define field Coordinate.x component Coordinate.x +gfx define field Coordinate.y component Coordinate.y +gfx define field Coordinate.z component Coordinate.z + +gfx def field x_velocity component U.1 +gfx def field y_velocity component U.2 +gfx def field z_velocity component U.3 +gfx def field pressure component U.4 + +gfx define field vector_field coord rectangular_cartesian component U.1 U.2 U.3 U.4 +gfx modify g_element StokesRegion node_points data y_velocity + +gfx edit scene +gfx create window 1; + +#Set the timekeeper playing +gfx timekeeper default set 1.0; +gfx timekeeper default speed 1; +gfx create time_editor