Hi Professor Spencer,

       I have added the debugger output. Please let me know if you need any information to be able to comment.


#0  0x00007ffff1a670c7 in raise () from /lib64/libc.so.6

#1  0x00007ffff1a68478 in abort () from /lib64/libc.so.6

#2  0x00007ffff2388b7d in __gnu_cxx::__verbose_terminate_handler() () from /usr/lib64/libstdc++.so.6

#3  0x00007ffff2386c06 in ?? () from /usr/lib64/libstdc++.so.6

#4  0x00007ffff2386c51 in std::terminate() () from /usr/lib64/libstdc++.so.6

#5  0x00007ffff2386e68 in __cxa_throw () from /usr/lib64/libstdc++.so.6

#6  0x00007ffff23ae7af in std::__throw_out_of_range(char const*) () from /usr/lib64/libstdc++.so.6

#7  0x00007ffff6d7a641 in std::vector<boost::shared_ptr<Nektar::SpatialDomains::Geometry>, std::allocator<boost::shared_ptr<Nektar::SpatialDomains::Geometry> > >::_M_range_check (this=0x7fff61f2a530, __n=0) at /usr/include/c++/4.8/bits/stl_vector.h:794

#8  0x00007ffff6d78553 in std::vector<boost::shared_ptr<Nektar::SpatialDomains::Geometry>, std::allocator<boost::shared_ptr<Nektar::SpatialDomains::Geometry> > >::at (this=0x7fff61f2a530, __n=0) at /usr/include/c++/4.8/bits/stl_vector.h:812

#9  0x00007ffff5cb7226 in Nektar::SpatialDomains::MeshGraph::WriteGeometry (this=0xb048d0, doc=...)

    at /home/sci/ashokj/gitlab/nektar/library/SpatialDomains/MeshGraph.cpp:1928

#10 0x00007ffff5cb7ef0 in Nektar::SpatialDomains::MeshGraph::WriteGeometry (this=0xb048d0, outfilename=...)

    at /home/sci/ashokj/gitlab/nektar/library/SpatialDomains/MeshGraph.cpp:1979

#11 0x00007ffff796b528 in Nektar::FieldUtils::OutputXml::Process (this=0xb044d0, vm=...)

    at /home/sci/ashokj/gitlab/nektar/library/FieldUtils/OutputModules/OutputXml.cpp:76

#12 0x0000000000474db4 in main (argc=5, argv=0x7fffffffe0c8) at /home/sci/ashokj/gitlab/nektar/utilities/FieldConvert/FieldConvert.cpp:386


Ashok.


On Sun, Nov 27, 2016 at 1:41 PM, ashok jallepalli <ash.nani@gmail.com> wrote:
Thank you, Professor Spencer.
          It works on 2D files but given the following error when applied on 3D mesh files.

FieldConvert -r 2,2.5,0.05,0.25,0.5,0.75 naca.xml naca_part.xml

Error:

what():  vector::_M_range_check

Aborted


Thank you,

Ashok.




On Sun, Nov 27, 2016 at 10:57 AM, Sherwin, Spencer J <s.sherwin@imperial.ac.uk> wrote:
Hi Ashok,

You might be able to achieve this using the range option in FieldConvert so possibly try

FieldConvert -r a,b,c,d  infield.xml outfield.xml

Cheers,
Spencer.

On 27 Nov 2016, at 12:36, ashok jallepalli <ash.nani@gmail.com> wrote:

Hi Everyone,
     Does any know a smart way to partition to mesh(Nektar XML) file to contain all element in a bounding box?

Example: If i have 10x10 element structure mesh on domain of [0,1]x[0,1]. I need to partition the mesh such that, the new partition file has all elements in the region (a,b) to (c,d). Assuming a,b,c,d are in [0,1].
     - If there is no command line option for the above scenario, How can I accomplish it using the API?

I have tried using the part-only and part-only-overlapping. I am unable to accomplish what I want.

Thanks in advance,
Ashok.
_______________________________________________
Nektar-users mailing list
Nektar-users@imperial.ac.uk
https://mailman.ic.ac.uk/mailman/listinfo/nektar-users

Spencer  Sherwin
McLaren Racing/Royal Academy of Engineering Research Chair, 
Professor of Computational Fluid Mechanics,
Department of Aeronautics,
Imperial College London
South Kensington Campus
London SW7 2AZ

+44 (0) 20 759 45052