Question on mesh partitions for post processing.
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.
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<mailto: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<mailto: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 s.sherwin@imperial.ac.uk<mailto:s.sherwin@imperial.ac.uk> +44 (0) 20 759 45052
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
s.sherwin@imperial.ac.uk +44 (0) 20 759 45052
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
s.sherwin@imperial.ac.uk +44 (0) 20 759 45052
participants (2)
-
ashok jallepalli
-
Sherwin, Spencer J