Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Fix pyramid element ref coords #566

Merged
merged 1 commit into from
May 31, 2023
Merged

Conversation

sandro-elsweijer
Copy link
Collaborator

@sandro-elsweijer sandro-elsweijer commented May 24, 2023

Describe your changes here:
The implementation of pyramid element reference coordinates was faulty:
image (1)
This can only be observed, if the pyramids are written out as curved and if the nonlinear subdivision level in ParaView is >0.
After the fix, it looks like this:
grafik

Here is a quick code for generating pyramids and write them as curved elements:

#include <t8.h>
#include <t8_cmesh/t8_cmesh_examples.h>
#include <t8_forest/t8_forest_general.h>
#include <t8_forest/t8_forest_geometrical.h>
#include <t8_forest/t8_forest_profiling.h>
#include <t8_forest/t8_forest_io.h>
#include <t8_forest/t8_forest_vtk.h>
#include <t8_schemes/t8_default/t8_default_cxx.hxx>

int main(int argc, char **argv)
{ 
  int mpiret = sc_MPI_Init (&argc, &argv);
  SC_CHECK_MPI (mpiret);

  sc_init (sc_MPI_COMM_WORLD, 1, 1, NULL, SC_LP_ESSENTIAL);
  t8_init (SC_LP_DEFAULT);
  t8_cmesh_t cmesh = t8_cmesh_new_hypercube(T8_ECLASS_PYRAMID, sc_MPI_COMM_WORLD, 0, 0, 0);
  t8_forest_t forest =  t8_forest_new_uniform (cmesh, t8_scheme_new_default_cxx(), 1, 0, sc_MPI_COMM_WORLD);
  t8_forest_write_vtk_ext(forest, "testfile", 1, 1, 1, 1, 0, 1, 0, 0, NULL);
  t8_forest_unref(&forest);
  sc_finalize ();

  mpiret = sc_MPI_Finalize ();
  SC_CHECK_MPI (mpiret);
}

All these boxes must be checked by the reviewers before merging the pull request:

As a reviewer please read through all the code lines and make sure that the code is fully understood, bug free, well-documented and well-structured.

General

  • The reviewer executed the new code features at least once and checked the results manually

  • The code follows the t8code coding guidelines

  • New source/header files are properly added to the Makefiles

  • The code is well documented

  • All function declarations, structs/classes and their members have a proper doxygen documentation

  • All new algorithms and data structures are sufficiently optimal in terms of memory and runtime (If this should be merged, but there is still potential for optimization, create a new issue)

Tests

  • The code is covered in an existing or new test case using Google Test

Github action

  • The code compiles without warning in debugging and release mode, with and without MPI (this should be executed automatically in a github action)

  • All tests pass (in various configurations, this should be executed automatically in a github action)

    If the Pull request introduces code that is not covered by the github action (for example coupling with a new library):

    • Should this use case be added to the github action?
    • If not, does the specific use case compile and all tests pass (check manually)

Scripts and Wiki

  • If a new directory with source-files is added, it must be covered by the script/find_all_source_files.scp to check the indentation of these files.
  • New Datatypes are added to t8indent_custom_datatypes.txt
  • If this PR introduces a new feature, it must be covered in an example/tutorial and a Wiki article.

Licence

  • The author added a BSD statement to doc/ (or already has one)

@sandro-elsweijer sandro-elsweijer added Bug For a bug in the Code shouldn't take long Can be resolved in under 30 mins labels May 24, 2023
@Davknapp Davknapp self-requested a review May 31, 2023 08:58
@Davknapp Davknapp merged commit c2bab86 into main May 31, 2023
3 checks passed
@Davknapp Davknapp deleted the fix-pyramid_element_ref_coords branch May 31, 2023 08:59
jfussbro pushed a commit to jfussbro/t8code that referenced this pull request Jun 27, 2023
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Bug For a bug in the Code shouldn't take long Can be resolved in under 30 mins
Projects
None yet
Development

Successfully merging this pull request may close these issues.

None yet

2 participants