Mshmet and 2D field

User Help
#1

Hi everyone,

I’m starting to use mmg2d and just finished the tutorials.
I tried to use mshmet ( olderVersionThatWorks) on a simple test 2d case in order to compute a first metric for adaptation. Unfortunately i can’t get any output even if the code seems to rune properly…

The output is :
mshmet toto.mesh -eps 0.00001 -hmin 0.01 -iso
– MSHMET, Release 2.3c (Aug. 25, 2009)
Copyright (c) LJLL, 2007-09
2022-10-17 15:50:06

– INPUT DATA
%% toto.mesh OPENED
– READING DATA FILE toto.mesh
%% toto.sol OPENED
– READING DATA FILE toto.sol
NUMBER OF GIVEN VERTICES 5230
NUMBER OF GIVEN TRIANGLES 9638
NUMBER OF GIVEN DATA 5230
– DATA READING COMPLETED. 0.01 sec.

&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
MODULE MSHMET-LJLL : 2.3c (Aug. 25, 2009)
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&

– PHASE 1 : ANALYSIS
** SETTING ADJACENCIES
– PHASE 1 COMPLETED. 0.00 sec.

– PHASE 2 : METRIC
– PHASE 2 COMPLETED. 0.00 sec.

&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
END OF MODULE MSHMET
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&

– WRITING DATA FILE toto.mesh
%% toto.new.sol OPENED
– WRITING COMPLETED

ELAPSED TIME 0.01 SEC. (0.01)

#2

here are the test files : toto.mesh (299.3 KB)
toto.sol (91.4 KB)

#3

and the output file : toto.new.sol (5.2 KB)

#4

Hello,

mshmet computes the size map from the solution Hessian. As the solution in toto.sol is constant, you obtain a constant metric.

Otherwise, the metric is saved in <filename.new.sol> file but the mesh file is not saved (because it has not changed). Tu visualize your output metric you can create a symbolic link toward your mesh file : ln -s toto.mesh toto.new.mesh

Best Regards,
Algiane

#5

Hi Algiane,

I’ve just found (by reading the soucre code) that in fact Meshmet is looking if the MeshVersionFormatted is 2, while it is 3 in my file.
Simply changing this number leads to a correct (i mean non empty) toto.new.sol.

Sincerly
Julien

#6

Thanks for this feedback: I close this topic then!

Have a nice day,
Algiane

#7