Re: [deal.II] Setting boundary id's for points instead of faces

2022-10-11 Thread Wasim Niyaz Munshi ce21d400
I was using a constraint object to apply the boundary conditions and component mask to deal with hanging nodes and apply x and y displacements separately. Now, I want to solve an elasticity problem on a unit square wherein I want to apply a displacement in y direction on top face and have a pin and

Re: [deal.II] Re: Trying to increase precision of floating point values in step3 of tutorial

2022-10-11 Thread Bruno Turcksin
You want to do something like this: DataOut<2> data_out; data_out.attach_dof_handler(dof_handler); data_out.add_data_vector(solution, "solution"); data_out.build_patches(); std::ofstream output("solution.vtk"); output.set_precisition(12); data_out.write_vtk(output); Best, Bruno L

Re: [deal.II] Re: Trying to increase precision of floating point values in step3 of tutorial

2022-10-11 Thread Wasim Niyaz Munshi ce21d400
Did you mean output.set_precision(12); It gives the following error error: ‘std::ofstream’ {aka ‘class std::basic_ofstream’} has no member named ‘set_precision’; did you mean ‘std::streamsize std::ios_base::_M_precision’? (not accessible from this context) I also tried output.setprecision(12)

Re: [deal.II] Re: Trying to increase precision of floating point values in step3 of tutorial

2022-10-11 Thread Bruno Turcksin
Sorry, it should be output.precision(12) Bruno Le mar. 11 oct. 2022 à 09:03, Wasim Niyaz Munshi ce21d400 < ce21d...@smail.iitm.ac.in> a écrit : > Did you mean > output.set_precision(12); > > It gives the following error > error: ‘std::ofstream’ {aka ‘class std::basic_ofstream’} has no > member

Re: [deal.II] Setting boundary id's for points instead of faces

2022-10-11 Thread Wolfgang Bangerth
Wasim, I imagine you are calling your setup_boundary_values() function before you call dof_handler.distribute_dofs()? Best W. On 10/11/22 05:16, Wasim Niyaz Munshi ce21d400 wrote: *** Caution: EXTERNAL Sender *** I was using a constraint object to apply the boundary conditions and componen

Re: [deal.II] Regarding boundary condition

2022-10-11 Thread Deepika Kushwah
Thank you Dr. Wolfgang Bangerth, I have successfully applied boundary conditions at the extreme node of 1-D bar but not able to impose it on intermediate points. I have following doubts- 1. How to check boundary ids and assign boundary values corresponding to any intermediate node? 2. In the solut

[deal.II] Re: Segmentation fault in Trilinos direct solver

2022-10-11 Thread Marc Fehling
Hello, I can't identify the problem right away, but here are some tips in debugging it: - Does the problem persist if you disable both h- and p-refinement? If so, than an issue might be that your data structures haven't been reinitialized after you updated your discretization. - Do

[deal.II] Re: Segmentation fault in Trilinos direct solver

2022-10-11 Thread 'yy.wayne' via deal.II User Group
Thank for your reply, Actually I just kill this error by writing the code again from scratch. I think the problem might be "system_matrix.clear()" is missed befor system_matrix.reinit(). 在2022年10月12日星期三 UTC+8 04:05:32 写道: > Hello, > > I can't identify the problem right away, but here are some

[deal.II] Re: Segmentation fault in Trilinos direct solver

2022-10-11 Thread 'yy.wayne' via deal.II User Group
... sadly after some modification the errors show again(I adapt step-29 into this framework(involves FESystem), with matrix-based p-multigrid). I think the error is not on direct solver, since the coarse matrix is initialized successfully in Klu direct solver. However the worst thing is *no err

Re: [deal.II] Re: Segmentation fault in Trilinos direct solver

2022-10-11 Thread Wolfgang Bangerth
On 10/11/22 23:19, 'yy.wayne' via deal.II User Group wrote: ** Specifically, when a Trilinos::SolverDirect object calls solve(), the line         solver_control.check(0,0) is not executed. I've debuged with both QT and gdb but error information is similiar, both only returns a segmentation f