Re: [deal.II] Supported matrix/vector to LinearOperator

2022-04-19 Thread hkch...@gmail.com
Thanks, Bruno and Jean-Paul. I will not use the TrilinosWrapper::SolverDirect in inverse_operator, replaced by an iterative solver. But Bruno, this example initializes a precondition by a matrix preconditioner.initialize(system_matrix, data); While I can't get the entity of the matrix, I want

Re: [deal.II] simple application of anti-periodic boundary conditions

2022-04-19 Thread Wolfgang Bangerth
On 4/19/22 13:17, Paras Kumar wrote: One follow-up question: The mesh (which is initially periodic at the boundaries)  is refined adaptively during the simulation and it cannot be guaranteed if both the left and the right boundary faces would be refined at the same time. Will

[deal.II] Re: Python bindings

2022-04-19 Thread Bruno Turcksin
Corbin, On Tuesday, April 19, 2022 at 3:04:39 PM UTC-4 corbin@gmail.com wrote: > I've built deal.ii version 9.3.0 with the > -DDEAL_II_COMPONENT_PYTHON_BINDING=ON configuration, yet when I run the > notebook, the Jupyter kernel is unable to locate the PyDealII module, that > is, the

Re: [deal.II] simple application of anti-periodic boundary conditions

2022-04-19 Thread Paras Kumar
Dear Wolfgang, Thank you for the hints. The second approach seems easy to start with. One follow-up question: The mesh (which is initially periodic at the boundaries) is refined adaptively during the simulation and it cannot be guaranteed if both the left and the right boundary faces would be

[deal.II] Python bindings

2022-04-19 Thread Corbin Foucart
Hello everyone, I'd like to write some custom Python bindings for deal.ii, and as a first step, I've been trying to reproduce the ipython notebook for step-49. I've built deal.ii version 9.3.0 with the -DDEAL_II_COMPONENT_PYTHON_BINDING=ON configuration, yet when I run the notebook, the

Re: [deal.II] Supported matrix/vector to LinearOperator

2022-04-19 Thread Jean-Paul Pelteret
Dear Chen, Bruno pointed you to the exact example that I was going to. It is certainly possible to use linear operators as preconditioners for solvers. It looks like support for TrilinosWrappers::SolverDirect has not yet been implemented:

[deal.II] deal.II Newsletter #208

2022-04-19 Thread 'Rene Gassmoeller' via deal.II User Group
Hello everyone! This is deal.II newsletter #208. 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: #13624: fixing typo of the type 'the the' (proposed by cdev911)

Re: [deal.II] Supported matrix/vector to LinearOperator

2022-04-19 Thread Bruno Turcksin
Chen, I am not sure if TrilinosWrappers::SolverDirect is supported but for your second point, you can use the linear_operator as a preconditioner. That's exactly what's done here for example. The

Re: [deal.II] Issue using component_wise()

2022-04-19 Thread Hermes Sampedro
Dear Wolfgang, That is totally right, thank you very much for the suggestions. Regards, Hermes El martes, 19 de abril de 2022 a las 17:02:31 UTC+2, Wolfgang Bangerth escribió: > On 4/19/22 08:39, Hermes Sampedro wrote: > > > > Thank you, could you please suggest me one in particular? > >

Re: [deal.II] Issue using component_wise()

2022-04-19 Thread Wolfgang Bangerth
On 4/19/22 08:39, Hermes Sampedro wrote: Thank you, could you please suggest me one in particular? Hermes -- they are not actually very difficult to find. Take a look at https://dealii.org/developer/doxygen/deal.II/namespaceDoFTools.html and then just read the section headings for different

Re: [deal.II] Issue using component_wise()

2022-04-19 Thread Hermes Sampedro
Dear Wolfgang, Thank you, could you please suggest me one in particular? Regards H El martes, 19 de abril de 2022 a las 16:28:04 UTC+2, Wolfgang Bangerth escribió: > On 4/19/22 08:05, Hermes Sampedro wrote: > > > > I have not been able to find the issue. I am considering looking for an > >

Re: [deal.II] Supported matrix/vector to LinearOperator

2022-04-19 Thread hkch...@gmail.com
Hi everyone, More questions: 1. where can I find the supported linear solver of Trilinos of linear_operator? TrilinosWrappers::SolverDirect seems not to be supported by the linear_operator 2. Does linear_operator can be a preconditioner? for example, SparseMatrix A, B, C, D, E, F; Matrix

Re: [deal.II] Issue using component_wise()

2022-04-19 Thread Wolfgang Bangerth
On 4/19/22 08:05, Hermes Sampedro wrote: I have not been able to find the issue. I am considering looking for an alternative. I am using 2 components to describe a complex solution (for real and imaginary parts). Assuming that I keep the ordering by default and that I export my solution

Re: [deal.II] Supported matrix/vector to LinearOperator

2022-04-19 Thread hkch...@gmail.com
One more question: where I can find the support linear solver of Trilinos of linear_operator? TrilinosWrappers::SolverDirect seem not be supported by the linear_operator best chen 在2022年4月18日星期一 UTC+8 22:07:36 写道: > Bruno Turcksin > Thank you! it really helpful! > > best > chen >

Re: [deal.II] Issue using component_wise()

2022-04-19 Thread Hermes Sampedro
Dear Wolfgang, I have not been able to find the issue. I am considering looking for an alternative. I am using 2 components to describe a complex solution (for real and imaginary parts). Assuming that I keep the ordering by default and that I export my solution vector (to work in an external

Re: [deal.II] Issues with step-1 and step-2

2022-04-19 Thread Wolfgang Bangerth
On 4/19/22 04:52, ankit tyagi wrote: Thank you so much. I think I got it now. Just one last thing. If the output stops plotting the parent edge, and plots only the child edges, the visual discrepancy can be avoided. Right? Will that be a problem in some other way? Yes, one could only show the

Re: [deal.II] Issues with step-1 and step-2

2022-04-19 Thread ankit tyagi
Thank you so much. I think I got it now. Just one last thing. If the output stops plotting the parent edge, and plots only the child edges, the visual discrepancy can be avoided. Right? Will that be a problem in some other way? On Tuesday, April 19, 2022 at 1:31:25 AM UTC+5:30 Wolfgang