[MOAB-dev] [Fwd: Fw: user: resolve_shared_ents]

[Lukasz Kaczmarczyk]

Hi

This is different problem. In attachment you can find a test file with  
generating a bug.

   MBInterface gMB = new MBCore();
   MBParallelComm *pcMB = new MBParallelComm(gMB);

   gMB->load_mesh(argv[1]);	

   MBRange tets;
   result = gMB->get_entities_by_type( 0, MBTET, tets );
   result = pcMB->resolve_shared_ents( 0, tets, 3, 0 );

If I find a bug I send solution.

Regards
Lukasz


-------------- next part --------------
A non-text attachment was scrubbed...
Name: aaa.h5m
Type: application/octet-stream
Size: 56664 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/moab-dev/attachments/20090619/00ae5014/attachment-0001.obj>
-------------- next part --------------





> -------- Original Message --------
> Subject: Fw: [MOAB-dev] user: resolve_shared_ents
> Date: Thu, 18 Jun 2009 16:59:21 -0700 (PDT)
> From: Brandon Smith <bsmith82 at yahoo.com>
> To: Tim Tautges <tautges at mcs.anl.gov>
> CC: kraftche at cae.wisc.edu
>
>
> Hi Tim,
>
> I found a similar problem (the same problem?) a while ago and fixed  
> it by adding this to MBSkinner.cpp at line 312.
>
> dum_elems = dum_elems.intersect(source_entities);
>
> This limits the dum_elems to entities that were input to the Skinner.
>
> Brandon
>
>
>
> ----- Forwarded Message ----
> From: Lukasz Kaczmarczyk <likask at civil.gla.ac.uk>
> To: moab-dev at mcs.anl.gov
> Sent: Thursday, June 18, 2009 4:43:02 PM
> Subject: [MOAB-dev] user: resolve_shared_ents
>
> Hello
> I am new user of moab.  I use Zoltan to partition a mesh. After that  
> I like to resolve shared ents. Function
> MBErrorCode resolve_shared_ents(MBEntityHandle this_set,
>                                  MBRange &proc_ents,
>                                  int resolve_dim = -1,
>                                  int shared_dim = -1);
> works ok for some meshes. However for tetrahedral mesh imported form  
> Gmsh I have following error:
> Resolving shared entities.
> Assertion failed: (0 < dum_elems.size() && 3 > dum_elems.size()),  
> function find_skin, file MBSkinner.cpp, line 315.
> For my case dum_elems.size()==4.
>
> Can You give me a tip how to solve this problem?
>
> Regards,
> Lukasz
>
>
>
>
>
>
> -- 
> ================================================================
> "You will keep in perfect peace him whose mind is
>  steadfast, because he trusts in you."               Isaiah 26:3
>
>             Tim Tautges            Argonne National Laboratory
>         (tautges at mcs.anl.gov)      (telecommuting from UW-Madison)
>         phone: (608) 263-8485      1500 Engineering Dr.
>           fax: (608) 263-4499      Madison, WI 53706
>
>