Lines Matching defs:coords
149 PetscErrorCode constant(PetscInt dim, PetscReal time, const PetscReal coords[], PetscInt Nc, PetscScalar *u, PetscCtx ctx)
161 PetscErrorCode constantDer(PetscInt dim, PetscReal time, const PetscReal coords[], const PetscReal n[], PetscInt Nc, PetscScalar *u, PetscCtx ctx)
175 PetscErrorCode linear(PetscInt dim, PetscReal time, const PetscReal coords[], PetscInt Nc, PetscScalar *u, PetscCtx ctx)
181 for (d = 0; d < Nc; ++d) u[d] = coords[d];
183 u[0] = coords[d];
187 PetscErrorCode linearDer(PetscInt dim, PetscReal time, const PetscReal coords[], const PetscReal n[], PetscInt Nc, PetscScalar *u, PetscCtx ctx)
205 PetscErrorCode quadratic(PetscInt dim, PetscReal time, const PetscReal coords[], PetscInt Nc, PetscScalar *u, PetscCtx ctx)
212 u[0] = coords[0] * coords[1];
213 u[1] = coords[1] * coords[2];
214 u[2] = coords[2] * coords[0];
216 u[0] = coords[0] * coords[0];
217 u[1] = coords[0] * coords[1];
220 u[0] = coords[d] * coords[d];
224 PetscErrorCode quadraticDer(PetscInt dim, PetscReal time, const PetscReal coords[], const PetscReal n[], PetscInt Nc, PetscScalar *u, PetscCtx ctx)
231 u[0] = coords[1] * n[0] + coords[0] * n[1];
232 u[1] = coords[2] * n[1] + coords[1] * n[2];
233 u[2] = coords[2] * n[0] + coords[0] * n[2];
235 u[0] = 2.0 * coords[0] * n[0];
236 u[1] = coords[1] * n[0] + coords[0] * n[1];
239 u[0] = 2.0 * coords[d] * n[d];
245 PetscErrorCode cubic(PetscInt dim, PetscReal time, const PetscReal coords[], PetscInt Nc, PetscScalar *u, PetscCtx ctx)
252 u[0] = coords[0] * coords[0] * coords[1];
253 u[1] = coords[1] * coords[1] * coords[2];
254 u[2] = coords[2] * coords[2] * coords[0];
256 u[0] = coords[0] * coords[0] * coords[0];
257 u[1] = coords[0] * coords[0] * coords[1];
260 u[0] = coords[d] * coords[d] * coords[d];
264 PetscErrorCode cubicDer(PetscInt dim, PetscReal time, const PetscReal coords[], const PetscReal n[], PetscInt Nc, PetscScalar *u, PetscCtx ctx)
271 u[0] = 2.0 * coords[0] * coords[1] * n[0] + coords[0] * coords[0] * n[1];
272 u[1] = 2.0 * coords[1] * coords[2] * n[1] + coords[1] * coords[1] * n[2];
273 u[2] = 2.0 * coords[2] * coords[0] * n[2] + coords[2] * coords[2] * n[0];
275 u[0] = 3.0 * coords[0] * coords[0] * n[0];
276 u[1] = 2.0 * coords[0] * coords[1] * n[0] + coords[0] * coords[0] * n[1];
279 u[0] = 3.0 * coords[d] * coords[d] * n[d];
285 PetscErrorCode quartic(PetscInt dim, PetscReal time, const PetscReal coords[], PetscInt Nc, PetscScalar *u, PetscCtx ctx)
292 u[0] = coords[0] * coords[0] * coords[1] * coords[1];
293 u[1] = coords[1] * coords[1] * coords[2] * coords[2];
294 u[2] = coords[2] * coords[2] * coords[0] * coords[0];
296 u[0] = coords[0] * coords[0] * coords[0] * coords[0];
297 u[1] = coords[0] * coords[0] * coords[1] * coords[1];
300 u[0] = coords[d] * coords[d] * coords[d] * coords[d];
304 PetscErrorCode quarticDer(PetscInt dim, PetscReal time, const PetscReal coords[], const PetscReal n[], PetscInt Nc, PetscScalar *u, PetscCtx ctx)
311 u[0] = 2.0 * coords[0] * coords[1] * coords[1] * n[0] + 2.0 * coords[0] * coords[0] * coords[1] * n[1];
312 u[1] = 2.0 * coords[1] * coords[2] * coords[2] * n[1] + 2.0 * coords[1] * coords[1] * coords[2] * n[2];
313 u[2] = 2.0 * coords[2] * coords[0] * coords[0] * n[2] + 2.0 * coords[2] * coords[2] * coords[0] * n[0];
315 u[0] = 4.0 * coords[0] * coords[0] * coords[0] * n[0];
316 u[1] = 2.0 * coords[0] * coords[1] * coords[1] * n[0] + 2.0 * coords[0] * coords[0] * coords[1] * n[1];
319 u[0] = 4.0 * coords[d] * coords[d] * coords[d] * n[d];
324 PetscErrorCode mytanh(PetscInt dim, PetscReal time, const PetscReal coords[], PetscInt Nc, PetscScalar *u, PetscCtx ctx)
330 for (d = 0; d < Nc; ++d) u[d] = PetscTanhReal(coords[d] - 0.5);
332 u[0] = PetscTanhReal(coords[d] - 0.5);
336 PetscErrorCode mytanhDer(PetscInt dim, PetscReal time, const PetscReal coords[], const PetscReal n[], PetscInt Nc, PetscScalar *u, PetscCtx ctx)
342 for (d = 0; d < Nc; ++d) u[d] = 1.0 / PetscSqr(PetscCoshReal(coords[d] - 0.5)) * n[d];
344 u[0] = 1.0 / PetscSqr(PetscCoshReal(coords[d] - 0.5)) * n[d];
349 PetscErrorCode trig(PetscInt dim, PetscReal time, const PetscReal coords[], PetscInt Nc, PetscScalar *u, PetscCtx ctx)
355 for (d = 0; d < Nc; ++d) u[d] = PetscSinReal(PETSC_PI * m * coords[d]);
357 u[0] = PetscSinReal(PETSC_PI * m * coords[d]);
361 PetscErrorCode trigDer(PetscInt dim, PetscReal time, const PetscReal coords[], const PetscReal n[], PetscInt Nc, PetscScalar *u, PetscCtx ctx)
367 for (d = 0; d < Nc; ++d) u[d] = PETSC_PI * m * PetscCosReal(PETSC_PI * m * coords[d]) * n[d];
369 u[0] = PETSC_PI * m * PetscCosReal(PETSC_PI * m * coords[d]) * n[d];