Lines Matching refs:C
1 C Copyright (c) 2017-2026, Lawrence Livermore National Security, LLC and other CEED contributors
2 C All Rights Reserved. See the top-level COPYRIGHT and NOTICE files for details.
3 C
4 C SPDX-License-Identifier: (BSD-2-Clause)
5 C
6 C This file is part of CEED: http://github.com/ceed
8 C> @file
9 C> Mass and diffusion operators examples using Nek5000
13 C-----------------------------------------------------------------------
17 C Set up mass operator
18 C Input: u1,u2,u3,q Output: v1,v2,ierr
30 C Quadrature Point Loop
58 C Rho
61 C RHS
70 C-----------------------------------------------------------------------
74 C Apply mass operator
75 C Input: u1,u2,q Output: v1,ierr
82 C Quadrature Point Loop
89 C-----------------------------------------------------------------------
93 C Set up diffusion operator
94 C Input: u1,u2,u3,q Output: v1,v2,ierr
107 C Quadrature Point Loop
146 C Geometric factors
147 C Stored in Voigt convention
148 C 0 5 4
149 C 5 1 3
150 C 4 3 2
158 C RHS
169 C-----------------------------------------------------------------------
173 C Apply diffusion operator
174 C Input: u1,u2,q Output: v1,ierr
181 C Quadrature Point Loop
193 C-----------------------------------------------------------------------
195 C Set h2 as rhoJac
196 C Input: bmq,nxq Output: h2
217 C-----------------------------------------------------------------------
219 C Set distance initial condition for BP1
220 C Input: Output: e
240 C-----------------------------------------------------------------------
242 C Set sine initial condition for BP3
243 C Input: Output: e
275 C-----------------------------------------------------------------------
281 C e = gllel(eg)
288 C-----------------------------------------------------------------------
290 C
291 C Note: this is an acceleration term, NOT a force!
292 C Thus, ffx will subsequently be multiplied by rho(x,t).
293 C
298 C e = gllel(eg)
306 C-----------------------------------------------------------------------
318 C-----------------------------------------------------------------------
320 C NOTE ::: This subroutine MAY NOT be called by every process
333 C-----------------------------------------------------------------------
347 C-----------------------------------------------------------------------
354 C-----------------------------------------------------------------------
369 C-----------------------------------------------------------------------
376 C-----------------------------------------------------------------------
378 C Apply zero Dirichlet boundary conditions
379 C Input: h2,nel Output: r1
391 C-----------------------------------------------------------------------
393 C Compute Linfty norm of (x-y)
394 C Input: x,y Output: n
420 C-----------------------------------------------------------------------
422 C 3D transpose of local gradient
423 C Input: u,N,D,Dt Output: ur,us,ut
441 C 3D transpose of local gradient
442 C Input: ur,us,ut,N,D,Dt Output: u
462 C-----------------------------------------------------------------------
464 C Routine to generate elemental geometric matrices on mesh 1
465 C (Gauss-Legendre Lobatto mesh).
554 C-----------------------------------------------------------------------
556 C Generate diagonal preconditioner for Helmholtz operator
557 C Input: h1,h2 Output: d
595 C
596 C Add cross terms if element is deformed.
597 C
631 C If axisymmetric, add a diagonal term in the radial direction (ISD=2)
662 C-----------------------------------------------------------------------
664 C Generate dummy diagonal preconditioner for Helmholtz operator
665 C Input: h1,h2 Output: d
676 C-----------------------------------------------------------------------
695 C-----------------------------------------------------------------------
697 C Solution to BP1 using libCEED
738 C Set up coordinates
752 C Init ceed library
758 C Set up Nek geometry data
762 C Set up true soln
766 C Set up solver parameters
773 C Create ceed basis for mesh and computation
783 C Create ceed element restrictions for mesh and computation
795 C Create ceed vectors
806 C Create ceed qfunctions for masssetupf and massf
831 C Create ceed operators
854 C Compute setup data
860 C Set up true RHS
863 C Set up algebraic RHS with libCEED
872 C Set up algebraic RHS with Nek5000
878 C Solve true RHS
885 C Output
912 C Solve libCEED algebraic RHS
920 C Output
947 C Solve Nek5000 algebraic RHS
955 C Output
985 C Destroy ceed handles
1004 C-----------------------------------------------------------------------
1006 C Solution to BP3 using libCEED
1048 C Set up coordinates and mask
1069 C Init ceed library
1075 C Set up Nek geometry data
1078 C Set up true soln
1083 C Set up solver parameters
1090 C Create ceed basis for mesh and computation
1100 C Create ceed element restrictions for mesh and computation
1114 C Create ceed vectors
1125 C Create ceed qfunctions for diffsetupf and diffusionf
1150 C Create ceed operators
1173 C Compute setup data
1179 C Set up true RHS
1183 C Set up algebraic RHS with libCEED
1193 C Set up algebraic RHS with Nek5000
1200 C Solve true RHS
1207 C Output
1234 C Solve libCEED algebraic RHS
1242 C Output
1269 C Solve Nek5000 algebraic RHS
1277 C Output
1307 C Destroy ceed handles
1326 C-----------------------------------------------------------------------
1329 C Scalar conjugate gradient iteration for solution of uncoupled
1330 C Helmholtz equations
1331 C Input: r1,h1,h2,rmult,binv,tin,ceed,ceed_op,vec_p1,vec_ap1,bpname
1332 C Output: u1,maxit
1339 C INPUT: rhs1 - rhs
1340 C h1 - exact solution
1392 C tolerance check here
1401 C if (nio.eq.0) write(6,1) ifield,istep,iter,nx,(wv(k),k=1,1)
1410 C if (nio.eq.0) write(6,2) iter,enorm,alph(1),pap(1),'alpha'
1435 C-----------------------------------------------------------------------
1438 C Vector conjugate gradient matvec for solution of uncoupled
1439 C Helmholtz equations
1440 C Input: pap,p1,h1,h2,bpname,ceed,ceed_op,vec_ap1,vec_p1
1441 C Output: ap1
1480 C-----------------------------------------------------------------------
1482 C Local matrix-vector for solution of BP3 (stiffness matrix)
1483 C Input: u,g,h1,h2,b,ju,us,ut Output: w
1503 C-----------------------------------------------------------------------
1505 C Vector conjugate gradient matvec for solution of BP1 (mass matrix)
1506 C Input: pap,p1,h1,h2 Output: ap1
1542 C-----------------------------------------------------------------------
1544 C Local matrix-vector for solution of BP3 (stiffness matrix)
1545 C Input: u,g,ur,us,ut,wk Output: w
1576 C-----------------------------------------------------------------------
1578 C Vector conjugate gradient matvec for solution of BP3 (stiffness matrix)
1579 C Input: pap,p1,h1,h2 Output: ap1
1610 C-----------------------------------------------------------------------
1612 C Vector conjugate gradient matvec for solution of uncoupled
1613 C Helmholtz equations
1614 C Input: pap,p1,h1,h2,bpname Output: ap1
1639 C-----------------------------------------------------------------------
1641 C Get BP to run
1642 C Input: Output: bp
1658 C-----------------------------------------------------------------------
1660 C Get CEED backend specification
1661 C Input: Output: spec
1672 C-----------------------------------------------------------------------
1674 C Get test mode flag
1675 C Input: Output: test
1689 C-----------------------------------------------------------------------