Re: [deal.II] Why is this asserted to be impossible?

On 09/18/2016 03:45 AM, Alex Zimmerman wrote: I question all of the asserts in this file: https://github.com/dealii/dealii/blob/master/source/numerics/vector_tools_rhs.cc My particular case is the one on line 36, for dim = 1 and spacedim = 1. Thanks for looking into this. This looks like we

Re: [deal.II] Re: Can we say that the higher order method, the more accurate?

On 09/18/2016 05:34 PM, JAEKWANG KIM wrote: I also wonder whether their is a limitation in higher order method that deal offer. If no, then we can use Q10 method, just putting degree value as 10? Yes. It will work, but you will find practical restrictions such as the fact that matrices

Re: [deal.II] Re: Can we say that the higher order method, the more accurate?

On 09/18/2016 03:43 PM, JAEKWANG KIM wrote: To, Daniel Arndt. Also, I didn't understand what you meant by _"Are you using a Mapping of the same order as your velocity ansatz space is?"_ You will want to look at the step-10 and step-11 tutorial programs to see what mappings are. You can't

Re: [deal.II] A*diag(V) with mmult?

On 09/18/2016 10:01 AM, 'Franck Kalala' via deal.II User Group wrote: How to perform the matrix multiplication A*diag(V) where diag(V) is a diagonal matrix with the vector V in the main diagonal. is there any such function in dealii? how to create the SparseMatrix diag(V) for a given vector

[deal.II] Re: Can we say that the higher order method, the more accurate?

I also wonder whether their is a limitation in higher order method that deal offer. If no, then we can use Q10 method, just putting degree value as 10? -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see

[deal.II] Re: Can we say that the higher order method, the more accurate?

This is how my mesh looks like. I am using half domain of 2D sphere assuming that it would symmetrical along phi domain. Tho my domain is 2d, I am

Re: [deal.II] A*diag(V) with mmult?

Franck, > How to perform the matrix multiplication > > A*diag(V) There is no function inside the library to do this, but you can write it yourself by scaling the entry a_ij by v_i. Code should be something like this (not tested): for (SparseMatrix::iterator it=A.begin();it!=A.end(); ++it)

[deal.II] Re: A*diag(V) with mmult?

Franck, One possibility is to use A.mmult(C,B,V) where B is an IdentityMatrix. Another way is to create an IdentityMatrix B, modify its entries suitably and use A.mmult(C,B). A third approach would be to just write the the scaling yourself. Best, Daniel Am Sonntag, 18. September 2016 18:01:17

[deal.II] Re: Can we say that the higher order method, the more accurate?

Jaekwank, You would indeed expect significantly better results with higher polynomial degree on the same mesh, if the solution is sufficiently regular. What kind of mesh are you using? Are you resolving all the boundary layers suitably? Are you using a Mapping of the same order as your velocity

[deal.II] A*diag(V) with mmult?

I am having a SparseMatrix A which has a certain sparsity pattern and a vector V. How to perform the matrix multiplication A*diag(V) where diag(V) is a diagonal matrix with the vector V in the main diagonal. is there any such function in dealii? how to create the SparseMatrix diag(V) for a

Re: [deal.II] Why is this asserted to be impossible?

I question all of the asserts in this file: https://github.com/dealii/dealii/blob/master/source/numerics/vector_tools_rhs.cc My particular case is the one on line 36, for dim = 1 and spacedim = 1. Thanks for looking into this. On Sunday, September 18, 2016 at 2:42:01 AM UTC+2, Wolfgang