[Nek5000-users] userf/userq in .usr file

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Fri Aug 12 08:18:21 CDT 2011


Thank you for the clarification Aleks.  I wanted to make sure I wasn't
unknowingly integrating my forcing function twice :).

Josh

On Thu, Aug 11, 2011 at 5:07 PM,  <nek5000-users at lists.mcs.anl.gov> wrote:
> Hi Josh,
>
> Yes, you are right that the variable SOURCE in userq seems to be outdated -- only qvol is used in nekuq().
>
> As far as multiplication of ffx,... by masss matrix is concerned -- it's done in makeuf() (after calling nekuf() which calls userf()), so there is no need to do the multiplication in userf():
>
>
>      subroutine makeuf
> ...
>      CALL NEKUF   (BFX,BFY,BFZ)
>      CALL OPCOLV (BFX,BFY,BFZ,BM1)
>
>
> Best,
> Aleks
>
>
> ----- Original Message -----
> From: nek5000-users at lists.mcs.anl.gov
> To: "nek5000-users" <nek5000-users at lists.mcs.anl.gov>
> Sent: Tuesday, August 9, 2011 4:43:15 PM
> Subject: [Nek5000-users] userf/userq in .usr file
>
> Hello Neks,
>
> I had a couple of questions regarding userf and userq.
>
> 1)  If we add something to, for example, ffx in userf, is that
> "something" subsequently multiplied by the mass matrix later on in the
> source code, or is that on the user to multiply ffx(i,j,k,e) by its
> appropriate bm1 entry?
>
> 2)  In userq, is "source" used for anything?  I tried looking through
> the source code, and I can see where "qvol" is used, but not where
> "source" is used.
>
>
> Thanks,
>
> --
> Josh Camp
>
> "All that is necessary for the triumph of evil is that good men do
> nothing" -- Edmund Burke
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>



-- 
Josh Camp

"All that is necessary for the triumph of evil is that good men do
nothing" -- Edmund Burke



More information about the Nek5000-users mailing list