[MOAB-dev] problem with Tqdcfr.cpp

Tim Tautges (ANL) tautges at mcs.anl.gov
Wed Feb 19 16:31:10 CST 2014


Before writing the file, try the cubit command 'compress ids all', which renumbers geometric entities so they're 
contiguous starting from 1.

Iulian, for reverse engineering what's in the cub file, see the util/CubitFile* in CGM, those even have information on 
how mesh-related stuff is written.

- tim


On 02/19/2014 03:57 PM, Grindeanu, Iulian R. wrote:
> OK, while looking at this, there is another problem with cubit 14 models, that I cannot explain yet.
> So the current code, with cubit 14 created model, will not read the global id for geometric sets. :(
> I will log a ticket.
>
> The old code (<=cubit 13) will not read new models, created with cubit 14. This is a bug I was trying to fix. It looks
> like I did not fix it well, because
> the global ids are not read, for geometric sets. (actually, this is another bug)  We need to reverse engineer more to
> guess the format of cubit file :(
>
> the current code will still read global ids for old models, created with cubit 13 or less.
>
> This  is related to your problem, Shengyong, but it is different. In your case, the sat file does not have the global ids.
>
> Thoughts?
>
>
> Iulian
>
> ------------------------------------------------------------------------------------------------------------------------
> *From:* Shengyong Cai [shengyongcai at gmail.com]
> *Sent:* Wednesday, February 19, 2014 12:56 PM
> *To:* Grindeanu, Iulian R.; Tautges, Timothy J.
> *Subject:* Re: problem with Tqdcfr.cpp
>
> It seems that user is not allowed to reassign the Global IDs in Cubit.
>
> Kind regards.
> Shengyong Cai
>
>
> On Wed, Feb 19, 2014 at 12:46 PM, Grindeanu, Iulian R. <iulian at anl.gov <mailto:iulian at anl.gov>> wrote:
>
>     So I have verified that in your scenario:
>     import the sat file
>     delete some volumes (leave vol 6)
>     mesh one face
>
>     save the cub file
>
>     moab will read geometric  global id as 0
>
>     With a different sat file (one from sharp, unic) the same scenario works fine
>
>     So it could be something with the sat file.
>     Maybe there is a way to give global ids in cubit, if they are missing
>     These are read as ACIS attributes. Still looking ...
>
>     Iulian
>     ------------------------------------------------------------------------------------------------------------------------
>     *From:* Shengyong Cai [shengyongcai at gmail.com <mailto:shengyongcai at gmail.com>]
>     *Sent:* Wednesday, February 19, 2014 12:27 PM
>     *To:* Grindeanu, Iulian R.; Tautges, Timothy J.
>
>     *Subject:* Re: problem with Tqdcfr.cpp
>
>     Yes! I import the sat file into cubit directly.
>     The attached file is the whole model.
>
>     Kind regards.
>     Shengyong Cai
>
>
>     On Wed, Feb 19, 2014 at 12:23 PM, Grindeanu, Iulian R. <iulian at anl.gov <mailto:iulian at anl.gov>> wrote:
>
>         my guess is that he is importing in cubit directly.
>         Shengyong, can you send me the sat file? Does it have global ids or not for geometry?
>
>         Iulian
>         ________________________________________
>         From: Tautges, Timothy J.
>         Sent: Wednesday, February 19, 2014 12:22 PM
>         To: Shengyong Cai; Grindeanu, Iulian R.
>         Subject: Re: problem with Tqdcfr.cpp
>
>         Wait, you're reading the .sat file into MOAB through ReadCGM in MOAB, or from .cub file through Tqdcfr?  That's
>         important, since there have been changes to both readers recently (by different sets of people)...
>
>         - tim
>
>         On 02/19/2014 12:17 PM, Shengyong Cai wrote:
>          > Hello Iulian,
>          > I import from sat file which is from Tim.
>          >
>          > Kind regards.
>          > Shengyong Cai
>          >
>          >
>          > On Wed, Feb 19, 2014 at 12:04 PM, Grindeanu, Iulian R. <iulian at anl.gov <mailto:iulian at anl.gov>
>         <mailto:iulian at anl.gov <mailto:iulian at anl.gov>>> wrote:
>          >
>          >     So entity sets global IDs are 0; vertices, edges and quads have global id, different from 0, in your file.
>          >     I am not sure how is this geometry created. It could be that global ids were not obtained.
>          >
>          >     I thought we give some global ids to entity sets, based on geometry ids.
>          >     I need to look at the code.
>          >
>          >     But how was this geometry created? Was it imported or something?
>          >
>          >     Tim, do you know how to set some global ids to geometry ents? I think that if the geometry is created in
>         cubit, it
>          >     will get a global id by default.
>          >
>          >     Iulian
>          >
>          >     ________________________________________
>          >     From: Shengyong Cai [shengyongcai at gmail.com <mailto:shengyongcai at gmail.com>
>         <mailto:shengyongcai at gmail.com <mailto:shengyongcai at gmail.com>>]
>          >     Sent: Wednesday, February 19, 2014 11:46 AM
>          >     To: Grindeanu, Iulian R.
>          >     Subject: problem with Tqdcfr.cpp
>          >
>          >     Hello Iulian,
>          >     I saw there are some commits from you for Tqdcfr.cpp. Right now I have
>          >     some problems from Tqdcfr.cpp. When MOAB reads a cub file with mesh, all
>          >     the global IDs are zero. Check the attached file for details.
>          >
>          >     FYI, cub file is from cubit 13.1.
>          >
>          >     Could you help me check that? Thanks!
>          >
>          >     --
>          >     Kind regards.
>          >     Shengyong Cai
>          >     ---------------------------------------------------------------
>          >     Room 422, Engineering Research Building
>          >     1500 Engineering Drive
>          >     Madison, WI, 53706
>          >     U.S.A.
>          >     Phone: 001-608-265-8086
>          >     Email: scai2 at wisc.edu <mailto:scai2 at wisc.edu> <mailto:scai2 at wisc.edu <mailto:scai2 at wisc.edu>> or
>         shengyongcai at gmail.com <mailto:shengyongcai at gmail.com> <mailto:shengyongcai at gmail.com
>         <mailto:shengyongcai at gmail.com>>
>          >     ---------------------------------------------------------------
>          >
>          >
>
>         --
>         ================================================================
>         "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 <mailto:tautges at mcs.anl.gov>)      (telecommuting from UW-Madison)
>            phone (gvoice): (608) 354-1459 <tel:%28608%29%20354-1459>      1500 Engineering Dr.
>                       fax: (608) 263-4499 <tel:%28608%29%20263-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 (gvoice): (608) 354-1459      1500 Engineering Dr.
             fax: (608) 263-4499      Madison, WI 53706



More information about the moab-dev mailing list