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

[Lukasz Kaczmarczyk]

FYI, by changing  create_vert_elem_adjs from true to false for function
MBSkinner::find_skin(const MBRange &source_entities,
                                  MBRange &forward_target_entities,
                                  MBRange &reverse_target_entities,
                                  bool create_vert_elem_adjs)
in file MBParallelComm.cpp line 2491 fix the problem. However I am not  
aware of all consequences  of this change.

Lukasz


On 19 Jun 2009, at 20:58, Tim Tautges wrote:

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