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

[Tim Tautges]

FYI, I haven't answered in detail because I'm working on a re-design of 
the whole parallel shared entity and ghosting stuff.  Hopefully your 
problem will go away once it's done.  I'll keep this example as a test.

- tim

Lukasz Kaczmarczyk wrote:
> 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
> 
> 
> 
> 
> 
> 
> 
>> -------- 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
>>
>>
> 

-- 
================================================================
"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