Re: [deal.II] Vector Equation

where x component of vector_RHS is -(grad_psi)_y / x and y component is (grad_psi) / x (here psi is the another calculated scalar solution) So, I refer to "class AdvectionField : public TensorFunction<1,dim>" in step-9.cc to form the vector_RHS You're approaching this the wrong way. The

Re: [deal.II] Using Sacado with example 15

On 08/24/2017 09:22 AM, Maxi Miller wrote: In the first case, there is nothing we need to do: we are using a continuous finite element, and face terms do not appear in the bilinear form in this case. The second case usually does not lead to face terms either if we enforce hanging

Re: [deal.II] Using Sacado with example 15

On 08/24/2017 09:11 AM, Praveen C wrote: I did this a while back, see here https://bitbucket.org/cpraveen/deal_ii/src/865f1c963f1089a870054bebf7180c9712e93da8/step-15-sacado/?at=master Ah, most excellent! I suspect that this can be written even more concisely with the current development ver

[deal.II] Re: Step-9, iterator over the cell

oh Thanks, Now I understand what it had meant ! Thank you On Thursday, August 24, 2017 at 4:26:59 PM UTC-5, Bruno Turcksin wrote: > > Jaekwong > > On Thursday, August 24, 2017 at 5:01:24 PM UTC-4, Jaekwang Kim wrote: >> >> >>

[deal.II] Re: Step-9, iterator over the cell

Jaekwong On Thursday, August 24, 2017 at 5:01:24 PM UTC-4, Jaekwang Kim wrote: > > > > > > *case 1: * > when we iterate over on cell-4 in my pic, it has four faces a

[deal.II] Step-9, iterator over the cell

Hi, all. I have a questions on an algorithm in step-9 that does some iteration over the cell with its neighbors. I tried to read explanation carefully, but some

Re: [deal.II] Using deal.II logo in publication

On 08/24/2017 01:38 PM, Timo Heister wrote: What about putting the logo and a download link to a high resolution variant on the homepage ont the Info -> License page [1]? That would be great. I also don't think we have a high-resolution (ideally vector art) version yet. Yes, all that sounds go

Re: [deal.II] Using deal.II logo in publication

> What about putting the logo and a download link to a high resolution > variant on the homepage ont the Info -> License page [1]? That would be great. I also don't think we have a high-resolution (ideally vector art) version yet. -- Timo Heister http://www.math.clemson.edu/~heister/ -- The de

Re: [deal.II] Using deal.II logo in publication

What about putting the logo and a download link to a high resolution variant on the homepage ont the Info -> License page [1]? Stating something like deal.II Logo The deal.II logo is copyrighted by ... and can be used under a [CC BY 4.0 license](link to license). You can fin

[deal.II] Re: Output nodal scalar

Thank you Jean-Paul. I did not hear about L2-projection method before, but it seems that the L2-projection method also evaluates the stress tensor at the quadrature points and then finds the cell-averaged stress. The only difference compared with the heuristic method is that it solves AX = Y wh

[deal.II] Re: Output nodal scalar

There's also the MeshWorker framework which has some integrators that can help. But I don't know how that works so I can't help you further other than to point it out to you. On Thursday, August 24, 2017 at 6:3

[deal.II] Re: Output nodal scalar

I can't think of any off of the top of my head, but I think that MatrixCreator::create_mass_matrix() plus VectorTools::create_right_hand_side()

[deal.II] Re: Output nodal scalar

Jean-Paul Any examples programming-wise? Like a deal.II manual/tutorial... Jie On Thursday, August 24, 2017 at 11:27:40 AM UTC-4, Jean-Paul Pelteret wrote: > > Hi Jie, > > Sure, we used it in this paper > >> @Article{vogel2014a-preprint, >> author= {Vogel, F. and Pelteret, J-P. V. and Kaes

[deal.II] Re: Output nodal scalar

Hi Jie, Sure, we used it in this paper > @Article{vogel2014a-preprint, > author= {Vogel, F. and Pelteret, J-P. V. and Kaessmair, S. and > Steinmann, P.}, > title = {Magnetic force and torque on particles subject to a > magnetic field}, > journal = {European Journal of Mechanics

Re: [deal.II] Using Sacado with example 15

Concerning the constraints: I was refering to this line: In the first case, there is nothing we need to do: we are using a > continuous finite element, and face terms do not appear in the bilinear > form in this case. The second case usually does not lead to face terms > either if we enforce ha

Re: [deal.II] Using Sacado with example 15

I did this a while back, see here https://bitbucket.org/cpraveen/deal_ii/src/865f1c963f1089a870054bebf7180c9712e93da8/step-15-sacado/?at=master Best praveen On Thu, Aug 24, 2017 at 8:26 PM, Wolfgang Bangerth wrote: > > I would like to get a feeling in working with Sacado, and would like to >>

[deal.II] Re: Multiply SparseMatrix and a Vector

In fact, that was the issue.. Works well now. :) Thanks ! On Thursday, August 24, 2017 at 4:21:59 PM UTC+2, Jean-Paul Pelteret wrote: > > Dear Nitish, > > Well, in the documentation to vmult >

[deal.II] Re: Output nodal scalar

Hi Jean-Paul Thanks for your help! The specific add_data_vector() function you shared does solve my problem! For your second suggestion on L2-smoothing, I'd like to give it a try. Can you please point me to some detailed references? Examples or documentations would be great! Jie On Thursday,

Re: [deal.II] Using Sacado with example 15

I would like to get a feeling in working with Sacado, and would like to try it on example 15. According to my understanding I now have to treat the boundary conditions on my own, and not by using the ConstraintMatrix. Is that still true, or was there an update since 2008? I suspect that you

Re: [deal.II] Using deal.II logo in publication

On 08/24/2017 08:02 AM, Matthias Maier wrote: Martin -- good question. I don't think anyone has ever thought about this. I would say go ahead and use it. Certainly, as the person who created that picture, I have no objections. What about we put a small licensing blurb on the homepage? Something

[deal.II] Re: Multiply SparseMatrix and a Vector

Dear Nitish, Well, in the documentation to vmult it states that when > the current object represents a parallel distributed matrix (of type > PETScWrappers::MPI::SparseMatrix

[deal.II] Re: Multiply SparseMatrix and a Vector

Errata: PETSc On Thursday, August 24, 2017 at 4:14:56 PM UTC+2, Nitish Anand wrote: > > Hello, > > I wish to multiply a PETSeWrapper::MPI::SparseMatrix with a > Vector. I would be glad if someone can point to the right material > to do this efficiently. > > Thanks in advance. > > --Nitish > --

[deal.II] Multiply SparseMatrix and a Vector

Hello, I wish to multiply a PETSeWrapper::MPI::SparseMatrix with a Vector. I would be glad if someone can point to the right material to do this efficiently. Thanks in advance. --Nitish -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://g

Re: [deal.II] Using deal.II logo in publication

On Thu, Aug 24, 2017, at 08:22 CDT, Wolfgang Bangerth wrote: > Martin -- good question. I don't think anyone has ever thought about > this. I would say go ahead and use it. Certainly, as the person who > created that picture, I have no objections. What about we put a small licensing blurb on t

Re: [deal.II] Using deal.II logo in publication

On 08/23/2017 06:39 PM, Martin Diehl wrote: I've written a short overview article on continuum scale modeling in which I mention deal.II (including the preferred reference and a link to the homepage). For an attractive paper, I want to include the deal.II-Logo. Are there any specific regulatio

[deal.II] Using Sacado with example 15

I would like to get a feeling in working with Sacado, and would like to try it on example 15. According to my understanding I now have to treat the boundary conditions on my own, and not by using the ConstraintMatrix. Is that still true, or was there an update since 2008? Furthermore, now the f