[deal.II] Re: Using constraints for already assembled RHS

2016-09-12 Thread RAJAT ARORA
Hi Daniel, Thanks for the reply. This indeed is one way of doing it. I was wondering if there is any other way in which this can be done. Thanks again. On Monday, September 12, 2016 at 5:04:47 AM UTC-4, Daniel Arndt wrote: > > Rajat, > > The first equation is the usual equilibrium equation.

[deal.II] More power outages

2016-09-12 Thread Guido Kanschat
Dear all, We are expecting more power outages on September 26 to 28. They are trying to fix the emergency supply such that the recent failures don't repeat. Let's keep our fingers crossed. The dealii web site will be down during these outages, which may last a few hours. Please be patient. I apo

Re: [deal.II] Dummy questions on parallelizing DG without meshworker

2016-09-12 Thread yoczhang
Hi, Weixiong Zheng I am wondering where is the DG code written by Praveen C, using the slope limiter? > One real quick question: when trying to see how to do slope limiter, I > visited your DG code, which is awesome, though it's a bit hard to follow > for me (newbie on DG and slope limiting)

Re: [deal.II] Dummy questions on parallelizing DG without meshworker

2016-09-12 Thread yoczhang
Hi, Weixiong Zheng I am wondering where is the DG code written by Praveen C, using the slope limiter? Thanks. Best, Yongchao Zhang 在 2015年11月4日星期三 UTC+8下午1:44:22,Weixiong Zheng写道: > > Thanks for the reply, > > > > 在 2015年11月3日星期二 UTC-6下午11:15:14,Praveen C写道: >> >> Hi >> >> I am using uniform

Re: [deal.II] Help with Vector template instantiations

2016-09-12 Thread dragnikolic
Hello Wolfgang, Here is the patch. I added one more change in functions Vector::reinit: if (omit_zeroing_entries == false) *this = static_cast(0); The line above fails if the Number type is not primitive. If it is a class it can't cast. So it could be: if (omit_zeroing_entries == false

[deal.II] Re: Using constraints for already assembled RHS

2016-09-12 Thread Daniel Arndt
Rajat, The first equation is the usual equilibrium equation. It is solved to give > displacements. The stiffness matrix is then multiplied by this global > displacement vector to give the global traction vector which forms the rhs > for the second equation. > The easiest solution to deal with c