Re: [deal.II] Matrix-free DG with h-adaptive mesh refinement

2022-08-16 Thread Marc Fehling
Hello! You are calling prepare_for_coarsening_and_refinement() twice which is not allowed. Thanks for reporting this on the mailing list -- we will improve the error message that you encountered. I spotted a few issues in the order of events of your program. Tutorials step-32, step-42,

[deal.II] get_generalized_support_points() returns only a vector of size 12 ?

2022-08-16 Thread LY XXXiao
Dear deal.ii group users, I am debugging a program using the function *'get_generalized_support_points()' *( where has_support_points()=0, while has_generalized_support_points()=1*)*. My FE system is defined as *'FESystem<3>fe(FE_Nedelec<3>(0), 2);'*, therefore, each active cell has

Re: [deal.II] interface between dealii and matlab for data transfer (using mex files?)

2022-08-16 Thread 'Andreas Hegendörfer' via deal . II User Group
Hi Simon, maybe you can solve your problem with Matlab's C++ interface ( https://de.mathworks.com/help/matlab/matlab_external/call-matlab-functions-from-c-1.html ). The interface allows you to access Matlab's command line via C++ code, and, thus, you can use the full capabilities of Matlab within

Re: [deal.II] Refine per direction

2022-08-16 Thread Wolfgang Bangerth
On 8/16/22 07:23, Uclus Heis wrote: GridGenerator::hyper_rectangle(triangulation,Point(0,0,0),Point(2,1,1),false ); triangulation.refine_global(2); This code generates. rectangle with 4 cells per direction. How can I perform a refinement so that I get more elements at the long dimension of

[deal.II] Refine per direction

2022-08-16 Thread Uclus Heis
Good afternoon, I am using the following function to create a mesh of a rectangle 2x1x1m GridGenerator::hyper_rectangle(triangulation,Point(0,0,0),Point(2,1,1),false ); triangulation.refine_global(2); This code generates. rectangle with 4 cells per direction. How can I perform a refinement

[deal.II] deal.II Newsletter #223

2022-08-16 Thread 'Rene Gassmoeller' via deal.II User Group
Hello everyone! This is deal.II newsletter #223. It automatically reports recently merged features and discussions about the deal.II finite element library. ## Below you find a list of recently proposed or merged features: #14190: FEEvaluationData: add virtual destructor (proposed by