| /petsc/src/ts/event/tests/ |
| H A D | ex2.c | 8 "y_dot = -0.2*x\n" 9 … "Using two event functions = piecewise-polynomials, zeros = 1 (rank-0), 9 (last rank)\n" 11 "-dir d : zero-crossing direction for events\n" 12 "-flg : additional output in Postevent\n" 13 … "-errtol e : error tolerance, for printing 'pass/fail' for located events (1e-5 by default)\n" 14 "-restart : flag for TSRestartStep() in PostEvent\n" 15 … "-dtpost x : if x > 0, then on even PostEvent calls 1st-post-event-step = x is set,\n" 16 … " on odd PostEvent calls 1st-post-event-step = PETSC_DECIDE is set,\n" 26 PetscReal evres[MAX_NEV]; // times of found zero-crossings 27 PetscReal ref[MAX_NEV]; // reference times of zero-crossings, for checking [all …]
|
| H A D | ex3.c | 8 "y_dot = -0.2*x\n" 11 … "- two polynomial event functions on rank-0 and last-rank (with zeros: 1.05, 9.05[terminating])\n" 12 … "- one event function on rank = '1%size', equal to V*sin(pi*t), zeros = 1,...,10\n" 16 "-dir d : zero-crossing direction for events: 0, 1, -1\n" 17 "-flg : additional output in Postevent\n" 18 … "-errtol e : error tolerance, for printing 'pass/fail' for located events (1e-5 by default)\n" 19 "-restart : flag for TSRestartStep() in PostEvent\n" 20 … "-dtpost x : if x > 0, then on even PostEvent calls 1st-post-event-step = x is set,\n" 21 … " on odd PostEvent calls 1st-post-event-step = PETSC_DECIDE is set,\n" 23 …"-v {float}: scaling of the sin() event function; for small v this event is triggered by the funct… [all …]
|
| H A D | ex1sin.c | 8 "y_dot = -0.2*x\n" 11 "-dir d : zero-crossing direction for events\n" 12 "-flg : additional output in Postevent\n" 13 … "-errtol e : error tolerance, for printing 'pass/fail' for located events (1e-5 by default)\n" 14 "-restart : flag for TSRestartStep() in PostEvent\n" 15 … "-dtpost x : if x > 0, then on even PostEvent calls 1st-post-event-step = x is set,\n" 16 … " on odd PostEvent calls 1st-post-event-step = PETSC_DECIDE is set,\n" 26 PetscReal evres[MAX_NEV]; // times of found zero-crossings 27 PetscReal ref[MAX_NEV]; // reference times of zero-crossings, for checking 31 …errtol; // error tolerance, for printing 'pass/fail' for located events (1e-5 by default) [all …]
|
| H A D | ex4.c | 8 "y_dot = -0.2*x\n" 9 "Using one or several event functions (on rank-0)\n" 10 …"This program is mostly intended to test the Anderson-Bjorck iteration with challenging event-func… 12 "-dir d : zero-crossing direction for events\n" 13 "-flg : additional output in Postevent\n" 14 … "-errtol e : error tolerance, for printing 'pass/fail' for located events (1e-5 by default)\n" 15 "-restart : flag for TSRestartStep() in PostEvent\n" 16 … "-dtpost x : if x > 0, then on even PostEvent calls 1st-post-event-step = x is set,\n" 17 … " on odd PostEvent calls 1st-post-event-step = PETSC_DECIDE is set,\n" 19 …"-func F : selects the event function [0, ..., 11], if F == -1 (default) is set, all event funct… [all …]
|
| H A D | ex5.c | 8 "y_dot = -0.2*x\n" 10 … "7 polynomials (dir=+1) with zeros: 1+2^i, i=-3,...3, on ranks=(i+3)%size\n" 11 … "7 polynomials (dir=-1) with zeros: 1+(8-2^i), i=-3,...3, on ranks=(i+3)%size\n" 12 "(t-5)^2 * sin(pi*t), with zeros = 1,2,...10, on rank-0\n" 13 " 0.5 * cos(pi*t), with zeros = 0.5,1.5,...9.5, on last-rank\n" 15 "-dir d : zero-crossing direction for events\n" 16 "-flg : additional output in Postevent\n" 17 … "-errtol e : error tolerance, for printing 'pass/fail' for located events (1e-5 by default)\n" 18 "-restart : flag for TSRestartStep() in PostEvent\n" 19 … "-dtpost x : if x > 0, then on even PostEvent calls 1st-post-event-step = x is set,\n" [all …]
|
| /petsc/src/vec/is/sf/tutorials/ |
| H A D | ex1.c | 1 static const char help[] = "Test star forest communication (PetscSF)\n\n"; 6 …ommon communication patterns can be expressed as updates of rootdata using leafdata and vice-versa. 31 …tscCall(PetscViewerASCIISynchronizedPrintf(viewer, "[%d] %" PetscInt_FMT " <- (%d,%" PetscInt_FMT … in PetscSFViewCustomLocals_Private() 54 PetscOptionsBegin(PETSC_COMM_WORLD, "", "PetscSF Test Options", "none"); in main() 56 …PetscCall(PetscOptionsBool("-test_all", "Test all SF communications", "", test_all, &test_all, NUL… in main() 58 PetscCall(PetscOptionsBool("-test_bcast", "Test broadcast", "", test_bcast, &test_bcast, NULL)); in main() 60 …PetscCall(PetscOptionsBool("-test_bcastop", "Test broadcast and reduce", "", test_bcastop, &test_b… in main() 62 …PetscCall(PetscOptionsBool("-test_reduce", "Test reduction", "", test_reduce, &test_reduce, NULL)); in main() 64 …PetscCall(PetscOptionsBool("-test_char", "Test signed char, unsigned char, and char", "", test_cha… in main() 67 …PetscCall(PetscOptionsString("-test_op", "Designate which MPI_Op to use", "", opstring, opstring, … in main() [all …]
|
| H A D | ex3.c | 1 static const char help[] = "Test freeing of MPI types in PetscSF\n\n"; 30 …PetscCall(PetscOptionsBool("-test_dupped_type", "Test dupped input type", "", test_dupped_type, &t… in main() 88 /*TEST 90 test: 91 suffix: basic 92 args: -sf_type basic 94 test: 96 args: -test_dupped_type -sf_type basic 98 test: 100 filter: grep -v "type" | grep -v "sort" [all …]
|
| H A D | ex2.c | 1 static const char help[] = "Test overlapped communication on a single star forest (PetscSF)\n\n"; 95 /*TEST 97 test: 98 suffix: basic 100 filter: grep -v "type" | grep -v "sort" 101 args: -sf_type basic 103 test: 106 filter: grep -v "type" | grep -v "sort" 108 …args: -sf_type window -sf_window_sync {{fence active lock}} -sf_window_flavor {{create dynamic all… 111 …# The nightly test suite with MPICH uses ch3:sock, which is broken when winsize == 0 in some of th… [all …]
|
| /petsc/src/vec/is/ao/tests/ex3d/ |
| H A D | ex3.cxx | 1 static char help[] = "AO test contributed by Sebastian Steiger <steiger@purdue.edu>, March 2011\n\n… 5 mpiexec -n 12 ./ex3 6 mpiexec -n 30 ./ex3 -ao_type basic 30 PetscCall(PetscOptionsGetString(NULL, NULL, "-datafiles", datafiles, sizeof(datafiles), &flg)); in main() 31 PetscCheck(flg, PETSC_COMM_WORLD, PETSC_ERR_USER, "Must specify -datafiles ${DATAFILESPATH}/ao"); in main() 38 int tmp = -1; in main() 40 tmp = -1; in main() 42 if (tmp == -1) break; in main() 69 /*TEST 74 test: [all …]
|
| /petsc/src/snes/tutorials/ |
| H A D | ex12.c | 10 …-dm_adapt_view hdf5:$PWD/adapt.h5 -sol_adapt_view hdf5:$PWD/adapt.h5::append -dm_adapt_pre_view hd… 14 -info :~sys,vec,is,mat,ksp,snes,ts 56 /* Reproducing tests from SISC 40(3), pp. A1473-A1493, 2018 */ 86 -\Delta u + f = -4 + 4 = 0 90 -\nabla u \cdot -\hat y |_{y=0} = (2y)|_{y=0} = 0 (bottom) 91 -\nabla u \cdot \hat y |_{y=1} = -(2y)|_{y=1} = -2 (top) 92 -\nabla u \cdot -\hat x |_{x=0} = (2x)|_{x=0} = 0 (left) 93 -\nabla u \cdot \hat x |_{x=1} = -(2x)|_{x=1} = -2 (right) 118 const PetscReal r2 = PetscSqr(x[0] - 0.5) + PetscSqr(x[1] - 0.5); in ball_u_2d() 119 const PetscReal xi = alpha * (radius2 - r2); in ball_u_2d() [all …]
|
| /petsc/src/ts/tutorials/ |
| H A D | ex40.c | 1 static char help[] = "Serial bouncing ball example to test TS event feature.\n"; 6 u2_t = -9.8 43 u[1] = -0.9 * u[1]; in PostEventFunction() 45 app->nbounces++; in PostEventFunction() 46 if (app->nbounces >= app->maxbounces) { in PostEventFunction() 47 … PetscCall(PetscPrintf(PETSC_COMM_SELF, "Ball bounced %" PetscInt_FMT " times\n", app->nbounces)); in PostEventFunction() 67 f[1] = -9.8; in RHSFunction() 117 f[0] = udot[0] - u[1]; in IFunction() 140 J[0][1] = -1.0; in IJacobian() 170 /* - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - in main() [all …]
|
| H A D | ex44.c | 1 static char help[] = "Parallel bouncing ball example formulated as a second-order system to test TS… 6 u_tt = -9.8 - 1/2 Cd (u_t)^2 sign(u_t) 54 v[0] = -app->Cr * v[0]; in PostEvent() 57 app->bounces++; in PostEvent() 59 if (app->bounces >= app->maxbounces) { // 'app->bounces' may be different on different processes in PostEvent() 60 …tf(PETSC_COMM_SELF, "Processor [%d]: Ball bounced %" PetscInt_FMT " times\n", rank, app->bounces)); in PostEvent() 78 Res = a[0] + 9.8 + 0.5 * app->Cd * v[0] * v[0] * PetscSignReal(PetscRealPart(v[0])); in I2Function() 100 Jac = shiftA + shiftV * app->Cd * v[0]; in I2Jacobian() 138 PetscCall(PetscOptionsReal("-Cd", "Drag coefficient", "", app.Cd, &app.Cd, NULL)); in main() 139 PetscCall(PetscOptionsReal("-Cr", "Restitution coefficient", "", app.Cr, &app.Cr, NULL)); in main() [all …]
|
| /petsc/src/ts/tutorials/output/ |
| H A D | ex40_n_ts_max_step_rejections-2.out | 1 …TSAdapt basic alpha 0: step 0 rejected t=0 + 2.500e-01 dt=2.500e-02 wlte= 538 wltea=… 2 …TSAdapt basic alpha 0: step 0 rejected t=0 + 2.500e-02 dt=2.500e-03 wlte= 20.5 wltea=… 3 …TSAdapt basic alpha 0: step 0 accepted t=0 + 2.500e-03 dt=4.166e-03 wlte=0.292 wltea=… 20 type: basic 27 maximum solution absolute value to be ignored -1. 31 tolerances: relative=1e-08, absolute=1e-50, solution=1e-08 38 alpha=1.000000e-04 39 maxlambda=1.000000e+00, minlambda=1.000000e-12 40 tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08 44 …restart=30, using classical (unmodified) Gram-Schmidt orthogonalization with no iterative refineme… [all …]
|
| H A D | ex40_n_ts_max_step_rejections-1.out | 1 …TSAdapt basic alpha 0: step 0 rejected t=0 + 2.500e-01 dt=2.500e-02 wlte= 538 wltea=… 2 …TSAdapt basic alpha 0: step 0 rejected t=0 + 2.500e-02 dt=2.500e-03 wlte= 20.5 wltea=… 19 type: basic 26 maximum solution absolute value to be ignored -1. 30 tolerances: relative=1e-08, absolute=1e-50, solution=1e-08 37 alpha=1.000000e-04 38 maxlambda=1.000000e+00, minlambda=1.000000e-12 39 tolerances: relative=1.000000e-08, absolute=1.000000e-15, lambda=1.000000e-08 43 …restart=30, using classical (unmodified) Gram-Schmidt orthogonalization with no iterative refineme… 44 happy breakdown tolerance=1e-30 [all …]
|
| /petsc/doc/developers/ |
| H A D | testing.md | 5 The PETSc test system consists of 7 - Formatted comments at the bottom of the tutorials and test source files that describes the tests … 8 - The *test generator* (`config/gmakegentest.py`) that parses the tutorial and test source files an… 10 - The *PETSc test harness* that consists of makefile and shell scripts that runs the executables wi… 14 ## PETSc Test Description Language 27 /*TEST 35 args: -debug -fields v1,v2,v3 36 test: 37 test: 38 args: -foo bar [all …]
|
| /petsc/src/ts/tests/ |
| H A D | ex17.c | 7 return p ? t * poly(p - 1, t) : 1.0; in poly() 12 return p > 0 ? (PetscReal)p * poly(p - 1, t) : 0.0; in dpoly() 41 PetscCall(MatShift(A, (PetscReal)-1)); in CreateMat() 115 …0 here, but arkimex errors with 1.e-14 discrepancy when using TSRResize without restart on some ma… in main() 120 PetscCall(PetscOptionsGetInt(NULL, NULL, "-order", &order, NULL)); in main() 121 PetscCall(PetscOptionsGetReal(NULL, NULL, "-tol", &tol, NULL)); in main() 158 err = PetscAbsScalar(results[0][i] - results[1][i]); in main() 163 err = PetscAbsScalar(results[0][i] - results[2][i]); in main() 178 /*TEST 181 args: -vec_mdot_use_gemv 0 -vec_maxpy_use_gemv 0 [all …]
|
| /petsc/src/snes/tutorials/output/ |
| H A D | ex12_fas_newton_1.out | 3 1 SNES Function norm 6.223e-10 4 0 SNES Function norm 4.482e-10 5 1 SNES Function norm < 1.e-11 6 0 SNES Function norm < 1.e-11 7 1 SNES Function norm < 1.e-11 8 0 SNES Function norm < 1.e-11 9 1 SNES Function norm < 1.e-11 10 0 SNES Function norm 5.875e-10 11 1 SNES Function norm < 1.e-11 12 1 SNES Function norm < 1.e-11 [all …]
|
| H A D | ex12_p4est_full_q2_conformal_serial.out | 3 1 SNES Function norm 8.08715e-05 4 0 SNES Function norm 6.89596e-05 5 1 SNES Function norm 4.693e-10 6 0 SNES Function norm 4.180e-10 7 1 SNES Function norm < 1.e-11 8 0 SNES Function norm 4.663e-10 9 1 SNES Function norm < 1.e-11 10 0 SNES Function norm 7.21962e-05 11 1 SNES Function norm 5.178e-10 12 1 SNES Function norm 5.178e-10 [all …]
|
| H A D | ex12_p4est_full_q2_conformal_serial_plex.out | 3 1 SNES Function norm 8.08715e-05 4 0 SNES Function norm 6.89596e-05 5 1 SNES Function norm 4.693e-10 6 0 SNES Function norm 4.180e-10 7 1 SNES Function norm < 1.e-11 8 0 SNES Function norm 4.663e-10 9 1 SNES Function norm < 1.e-11 10 0 SNES Function norm 7.21962e-05 11 1 SNES Function norm 5.178e-10 12 1 SNES Function norm 5.178e-10 [all …]
|
| /petsc/src/dm/tests/ |
| H A D | ex49.c | 1 static char help[] = "Test basic DMProduct operations.\n\n"; 14 PetscCall(PetscOptionsGetInt(NULL, NULL, "-dim", &dim, NULL)); in main() 25 /*TEST 27 test: 31 TEST*/
|
| /petsc/src/vec/is/ao/tests/output/ |
| H A D | ex1_1.out | 2 type: basic 4 PETSc->App App->PETSc 10 [0] 2,1,9,7 PetscToApplication 2 1 -1 -1 13 Test AOCreateMemoryScalable: 17 PETSc->App App->PETSc 25 Test AOCreateMemoryScalable with ispetsc=NULL: 29 PETSc->App App->PETSc 37 Test general API: 38 [0] 2,1,9,7 PetscToApplication 2 1 -1 -1
|
| H A D | ex1_2.out | 2 type: basic 4 PETSc->App App->PETSc 20 Test AOCreateMemoryScalable: 24 PETSc->App App->PETSc 38 Test AOCreateMemoryScalable with ispetsc=NULL: 42 PETSc->App App->PETSc 56 Test general API:
|
| /petsc/src/sys/objects/device/tests/ |
| H A D | ex2.c | 12 /* basic creation and destruction */ in main() 17 /* double free is no-op */ in main() 21 /* test global context returns a valid context */ in main() 25 /* test locally setting to null doesn't clobber the global */ in main() 30 /* test duplicate */ in main() 46 /*TEST 52 args: -device_enable {{lazy eager}} 53 test: 56 test: 58 args: -root_device_context_device_type host [all …]
|
| /petsc/ |
| H A D | .gitlab-ci.yml | 2 # stage-1 take only a few minutes; they do not run the full test suite or external packages. 4 # stage-2 runs on MCS systems and may take 10 to 15 minutes. They run the full test suite but with … 6 # stage-3 runs on MCS systems and may take an hour or more. They run the full test suite and heavil… 8 # The stage-(n) tests are only started if all of the stage-(n-1) tests run without error 13 - stage-1 14 - stage-2 15 - stage-3 16 - stage-4 19 GIT_CLEAN_FLAGS: -ffdxq 21 BASE_EXTRA_OPTIONS: -nox -nox_warning -malloc_dump [all …]
|
| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex62.c | 3 code indicates the procedure for setting user-defined subdomains.\n\ 4 See section 'ex62' below for command-line options.\n\ 5 Without -user_set_subdomains, the general PCGASM options are meaningful:\n\ 6 -pc_gasm_total_subdomains\n\ 7 -pc_gasm_print_subdomains\n\ 13 and ex2.c for more detailed comments on the basic usage of KSP 16 The GASM preconditioner is fully parallel. The user-space routine 29 petscsys.h - base PETSc routines petscvec.h - vectors 30 petscmat.h - matrices 31 petscis.h - index sets petscksp.h - Krylov subspace methods [all …]
|