Skip to content

Conversation

@Algiane
Copy link
Member

@Algiane Algiane commented Dec 5, 2023

This PR adds a suitable handling of edge length histogram. Now we differentiate the cases where:

  • we have no metric at all - in this case, the histo cannot be printed because we cannot compute edge lengths;
  • we fail due to a memory issue;
  • we have no metric on some MPI processes but a metrics on others (empty group) : we can print the histo;
  • we have a metric everywhere : we can print the histo.

It also fixes an error in the computation of edge lengths for an anisotropic metric with ridge storage.

  - Calls the suitable function in PMMG_computePrilen when an aniso metric
 is provided with a specific metric storage at ridges (metRidTyp==1)
  - Check if metric is allocated before computing the edge lengths (to avoid memory errors)
 but inside the `MMG3D|PMMG_computePrilen` functions to be able to differentiate the cases where:
    - computePrilen fails due to memory error
    - We have a metric on at least one MPI proc so we can compute the edge lengths
    - We don't have any metric (and we cannot compute the edge length histo). In this case and
      with a centralized mesh on 1 proc and -v 10 verbosity, we were previously printing infinit edge lengths.

It makes the PMMG_prilen function callable everywhere in the software.
@Algiane Algiane added kind: bug Something isn't working part: parmmg specific to parmmg code part labels Dec 5, 2023
@Algiane Algiane self-assigned this Dec 5, 2023
@Algiane Algiane merged commit eaebcbe into develop Dec 8, 2023
@Algiane Algiane deleted the feature/fix-edgeLength branch August 21, 2024 18:45
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

kind: bug Something isn't working part: parmmg specific to parmmg code part

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants