Hello,
I have a quick question. Should I be expecting FE_Simplex of cubic order to
work in an FESystem? In case the answer is "yes", I'm attaching a minimal
example that I would expect to work, but it fails some assert in the
library. I tried checking, but could not really understand why this is
happening.
Unfortunately, I wouldn't have any comment on which specific commit might
have caused the issue (if it is an issue) as I just switched to my own
build of deal.II and was previously using 9.5.1. I'm not sure if this new
instalation could be the issue, so I'll describe what I did:
I compiled deal.II with
cmake -DCMAKE_INSTALL_PREFIX=/usr/local -DDEAL_II_WITH_MPI=ON
-DCMAKE_EXPORT_COMPILE_COMMANDS=ON -DDEAL_II_WITH_TRILINOS=ON
-DDEAL_II_WITH_PETSC=ON ..
followed by
sudo make -j 6 install
I did not get any notable errors or warnings.
I'm compiling and running the example with
cmake ..
make example
./example
Thanks for any help you can offer.
Regards,
Kyle
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/b4513a35-7437-4b7c-acc3-7c14c9e0b21cn%40googlegroups.com.
#include <deal.II/fe/fe_system.h>
#include <deal.II/fe/fe_simplex_p.h>
#include <iostream>
using namespace dealii;
int main()
{
unsigned int degree = 2;
FESystem<2> fe_system_1(FE_SimplexP<2>(degree),2,FE_SimplexP<2>(degree-1),1);
std::cout << "Success!" << std::endl;
degree = 3;
// Fails assertion inside deal.II function called from FESystem::initialize.
FESystem<2> fe_system_2(FE_SimplexP<2>(degree),2,FE_SimplexP<2>(degree-1),1);
std::cout << "Success!" << std::endl;
return 0;
}
CMAKE_MINIMUM_REQUIRED(VERSION 3.28.3)
FIND_PACKAGE(deal.II 9.7.0 REQUIRED)
DEAL_II_INITIALIZE_CACHED_VARIABLES()
PROJECT(myproject)
ADD_EXECUTABLE(example example.cpp)
DEAL_II_SETUP_TARGET(example)
ADD_CUSTOM_TARGET(debug
COMMAND ${CMAKE_COMMAND} -DCMAKE_BUILD_TYPE=Debug ${CMAKE_SOURCE_DIR}
COMMAND ${CMAKE_COMMAND} --build ${CMAKE_BINARY_DIR} --target all
COMMENT "Switch CMAKE_BUILD_TYPE to Debug"
)
ADD_CUSTOM_TARGET(release
COMMAND ${CMAKE_COMMAND} -DCMAKE_BUILD_TYPE=Release ${CMAKE_SOURCE_DIR}
COMMAND ${CMAKE_COMMAND} --build ${CMAKE_BINARY_DIR} --target all
COMMENT "Switch CMAKE_BUILD_TYPE to Release"
)
