| /petsc/src/tao/leastsquares/tests/ |
| H A D | chwirut2.c | 33 PetscReal t[NOBSERVATIONS]; /* array of independent variables of observation */ member
130 PetscReal *y = user->y, *f, *t = user->t; in EvaluateFunction() local
137 …for (i = 0; i < NOBSERVATIONS; i++) f[i] = y[i] - PetscExpScalar(-x[0] * t[i]) / (x[1] + x[2] * t[… in EvaluateFunction()
150 PetscReal *t = user->t; in EvaluateJacobian() local
157 base = PetscExpScalar(-x[0] * t[i]) / (x[1] + x[2] * t[i]); in EvaluateJacobian()
159 user->j[i][0] = t[i] * base; in EvaluateJacobian()
160 user->j[i][1] = base / (x[1] + x[2] * t[i]); in EvaluateJacobian()
161 user->j[i][2] = base * t[i] / (x[1] + x[2] * t[i]); in EvaluateJacobian()
191 PetscReal *t = user->t, *y = user->y; in InitializeData() local
196 t[i++] = 0.5000; in InitializeData()
[all …]
|
| /petsc/src/tao/leastsquares/tutorials/ |
| H A D | chwirut1.c | 32 PetscReal t[NOBSERVATIONS]; /* array of independent variables of observation */ member
108 PetscReal *y = user->y, *f, *t = user->t; in EvaluateFunction() local
114 …for (i = 0; i < NOBSERVATIONS; i++) f[i] = y[i] - PetscExpScalar(-x[0] * t[i]) / (x[1] + x[2] * t[… in EvaluateFunction()
128 PetscReal *t = user->t; in EvaluateJacobian() local
134 base = PetscExpScalar(-x[0] * t[i]) / (x[1] + x[2] * t[i]); in EvaluateJacobian()
136 user->j[i][0] = t[i] * base; in EvaluateJacobian()
137 user->j[i][1] = base / (x[1] + x[2] * t[i]); in EvaluateJacobian()
138 user->j[i][2] = base * t[i] / (x[1] + x[2] * t[i]); in EvaluateJacobian()
168 PetscReal *t = user->t, *y = user->y; in InitializeData() local
173 t[i++] = 0.5000; in InitializeData()
[all …]
|
| H A D | chwirut2.c | 35 PetscReal t[NOBSERVATIONS]; /* array of independent variables of observation */ member
166 PetscReal *t = user->t, *y = user->y; in InitializeData() local
171 t[i++] = 0.5000; in InitializeData()
173 t[i++] = 0.6250; in InitializeData()
175 t[i++] = 0.7500; in InitializeData()
177 t[i++] = 0.8750; in InitializeData()
179 t[i++] = 1.0000; in InitializeData()
181 t[i++] = 1.2500; in InitializeData()
183 t[i++] = 1.7500; in InitializeData()
185 t[i++] = 2.2500; in InitializeData()
[all …]
|
| H A D | chwirut1f.F90 | 19 PetscReal t(0:213)
53 f_v(i + 1) = y(i) - exp(-x_v(1)*t(i))/(x_v(2) + x_v(3)*t(i))
77 y(i) = 92.9000; t(i) = 0.5000; i = i + 1
78 y(i) = 78.7000; t(i) = 0.6250; i = i + 1
79 y(i) = 64.2000; t(i) = 0.7500; i = i + 1
80 y(i) = 64.9000; t(i) = 0.8750; i = i + 1
81 y(i) = 57.1000; t(i) = 1.0000; i = i + 1
82 y(i) = 43.3000; t(i) = 1.2500; i = i + 1
83 y(i) = 31.1000; t(i) = 1.7500; i = i + 1
84 y(i) = 23.6000; t(i) = 2.2500; i = i + 1
[all …]
|
| H A D | chwirut2f.F90 | 14 PetscReal t(0:213)
28 y(i) = 92.9000; t(i) = 0.5000; i = i + 1
29 y(i) = 78.7000; t(i) = 0.6250; i = i + 1
30 y(i) = 64.2000; t(i) = 0.7500; i = i + 1
31 y(i) = 64.9000; t(i) = 0.8750; i = i + 1
32 y(i) = 57.1000; t(i) = 1.0000; i = i + 1
33 y(i) = 43.3000; t(i) = 1.2500; i = i + 1
34 y(i) = 31.1000; t(i) = 1.7500; i = i + 1
35 y(i) = 23.6000; t(i) = 2.2500; i = i + 1
36 y(i) = 31.0500; t(i) = 1.7500; i = i + 1
[all …]
|
| /petsc/src/ts/tests/output/ |
| H A D | ex2_euler.out | 1 At t = 0.000000e+00 u = 1.000000e+00 1.000000e+00 1.000000e+00
2 At t = 0.000000e+00 errors = 0.000000e+00 0.000000e+00 0.000000e+00
3 At t = 1.000000e-03 u = 1.003000e+00 1.004000e+00 1.003000e+00
4 At t = 1.000000e-03 errors = -5.005672e-06 -7.008007e-06 -5.005672e-06
5 At t = 2.000000e-03 u = 1.006010e+00 1.008014e+00 1.006010e+00
6 At t = 2.000000e-03 errors = -1.004541e-05 -1.406411e-05 -1.004541e-05
7 At t = 3.000000e-03 u = 1.009030e+00 1.012042e+00 1.009030e+00
8 At t = 3.000000e-03 errors = -1.511939e-05 -2.116855e-05 -1.511939e-05
9 At t = 4.000000e-03 u = 1.012060e+00 1.016084e+00 1.012060e+00
10 At t = 4.000000e-03 errors = -2.022779e-05 -2.832159e-05 -2.022779e-05
[all …]
|
| H A D | ex2_be_adapt.out | 1 At t = 0.000000e+00 u = 1.000000e+00 1.000000e+00 1.000000e+00
2 At t = 0.000000e+00 errors = 0.000000e+00 0.000000e+00 0.000000e+00
3 At t = 1.000000e-03 u = 1.003010e+00 1.004014e+00 1.003010e+00
4 At t = 1.000000e-03 errors = 5.028445e-06 7.040158e-06 5.028445e-06
5 At t = 2.000000e-03 u = 1.006030e+00 1.008042e+00 1.006030e+00
6 At t = 2.000000e-03 errors = 1.009117e-05 1.412871e-05 1.009117e-05
7 At t = 4.000000e-03 u = 1.012111e+00 1.016155e+00 1.012111e+00
8 At t = 4.000000e-03 errors = 3.052675e-05 4.274366e-05 3.052675e-05
9 At t = 8.000000e-03 u = 1.024437e+00 1.032612e+00 1.024437e+00
10 At t = 8.000000e-03 errors = 1.139107e-04 1.595211e-04 1.139107e-04
[all …]
|
| H A D | ex2_beuler.out | 1 At t = 0.000000e+00 u = 1.000000e+00 1.000000e+00 1.000000e+00
2 At t = 0.000000e+00 errors = 0.000000e+00 0.000000e+00 0.000000e+00
3 At t = 1.000000e-03 u = 1.003010e+00 1.004014e+00 1.003010e+00
4 At t = 1.000000e-03 errors = 5.028445e-06 7.040158e-06 5.028445e-06
5 At t = 2.000000e-03 u = 1.006030e+00 1.008042e+00 1.006030e+00
6 At t = 2.000000e-03 errors = 1.009117e-05 1.412871e-05 1.009117e-05
7 At t = 3.000000e-03 u = 1.009060e+00 1.012084e+00 1.009060e+00
8 At t = 3.000000e-03 errors = 1.518836e-05 2.126592e-05 1.518836e-05
9 At t = 4.000000e-03 u = 1.012101e+00 1.016141e+00 1.012101e+00
10 At t = 4.000000e-03 errors = 2.032017e-05 2.845202e-05 2.032017e-05
[all …]
|
| H A D | ex2_rk.out | 1 At t = 0.000000e+00 u = 1.000000e+00 1.000000e+00 1.000000e+00
2 At t = 0.000000e+00 errors = 0.000000e+00 0.000000e+00 0.000000e+00
3 …TSAdapt basic rk 0:3bs step 0 accepted t=0 + 1.000e-03 dt=1.000e-02 wlte=4.09e-06 wlt…
4 At t = 1.000000e-03 u = 1.003005e+00 1.004007e+00 1.003005e+00
5 At t = 1.000000e-03 errors = -4.836798e-12 -6.837864e-12 -4.836798e-12
6 …TSAdapt basic rk 0:3bs step 1 accepted t=0.001 + 1.000e-02 dt=5.597e-02 wlte=0.00416 wlte…
7 At t = 1.100000e-02 u = 1.033613e+00 1.044858e+00 1.033613e+00
8 At t = 1.100000e-02 errors = -4.883664e-08 -6.904508e-08 -4.883664e-08
9 …TSAdapt basic rk 0:3bs step 2 accepted t=0.011 + 5.597e-02 dt=5.465e-02 wlte=0.783 wltea=…
10 At t = 6.696691e-02 u = 1.225076e+00 1.301734e+00 1.225076e+00
[all …]
|
| /petsc/src/ts/interface/ |
| H A D | tscreate.c | 31 TS t; in TSCreate() local
37 PetscCall(PetscHeaderCreate(t, TS_CLASSID, "TS", "Time stepping", "TS", comm, TSDestroy, TSView)); in TSCreate()
40 t->problem_type = TS_NONLINEAR; in TSCreate()
41 t->equation_type = TS_EQ_UNSPECIFIED; in TSCreate()
43 t->ptime = 0.0; in TSCreate()
44 t->time_step = 0.1; in TSCreate()
45 t->initial_time_step = t->time_step; in TSCreate()
46 t->max_time = PETSC_MAX_REAL; in TSCreate()
47 t->exact_final_time = TS_EXACTFINALTIME_UNSPECIFIED; in TSCreate()
48 t->steps = 0; in TSCreate()
[all …]
|
| /petsc/src/ts/tutorials/output/ |
| H A D | ex24_2.out | 1 0 t= 0.0e+00 dt= 1.0e-05 f= 2.0e+06 df= 2.6e+05 it=( 0, 0)
2 1 t= 1.0e-05 dt= 1.1e-05 f= 1.5e+06 df= 2.1e+05 it=( 3, 3)
3 2 t= 2.1e-05 dt= 1.5e-05 f= 1.2e+06 df= 1.7e+05 it=( 3, 3)
4 3 t= 3.6e-05 dt= 2.0e-05 f= 8.5e+05 df= 1.4e+05 it=( 3, 3)
5 4 t= 5.6e-05 dt= 2.9e-05 f= 6.0e+05 df= 1.0e+05 it=( 3, 3)
6 5 t= 8.6e-05 dt= 4.4e-05 f= 4.0e+05 df= 7.7e+04 it=( 3, 3)
7 6 t= 1.3e-04 dt= 6.9e-05 f= 2.4e+05 df= 5.4e+04 it=( 3, 3)
8 7 t= 2.0e-04 dt= 1.2e-04 f= 1.4e+05 df= 3.5e+04 it=( 3, 3)
9 8 t= 3.2e-04 dt= 2.1e-04 f= 6.9e+04 df= 2.1e+04 it=( 3, 3)
10 9 t= 5.3e-04 dt= 4.3e-04 f= 3.0e+04 df= 1.1e+04 it=( 4, 4)
[all …]
|
| H A D | ex41_d.out | 1 Ball hit the ground at t = 4.08 seconds -> Processor[0]
2 Ball hit the ground at t = 4.08 seconds -> Processor[0]
3 Ball hit the ground at t = 4.13 seconds -> Processor[1]
4 Ball hit the ground at t = 4.13 seconds -> Processor[1]
5 Ball hit the ground at t = 7.76 seconds -> Processor[0]
6 Ball hit the ground at t = 7.76 seconds -> Processor[0]
7 Ball hit the ground at t = 7.95 seconds -> Processor[1]
8 Ball hit the ground at t = 7.95 seconds -> Processor[1]
9 Ball hit the ground at t = 11.06 seconds -> Processor[0]
10 Ball hit the ground at t = 11.06 seconds -> Processor[0]
[all …]
|
| H A D | ex41_c.out | 1 Ball hit the ground at t = 4.08 seconds -> Processor[0]
2 Ball hit the ground at t = 4.08 seconds -> Processor[0]
3 Ball hit the ground at t = 4.13 seconds -> Processor[1]
4 Ball hit the ground at t = 4.13 seconds -> Processor[1]
5 Ball hit the ground at t = 7.76 seconds -> Processor[0]
6 Ball hit the ground at t = 7.76 seconds -> Processor[0]
7 Ball hit the ground at t = 7.95 seconds -> Processor[1]
8 Ball hit the ground at t = 7.95 seconds -> Processor[1]
9 Ball hit the ground at t = 11.06 seconds -> Processor[0]
10 Ball hit the ground at t = 11.06 seconds -> Processor[0]
[all …]
|
| H A D | ex41_f.out | 1 Ball hit the ground at t = 4.08 seconds -> Processor[0]
2 Ball hit the ground at t = 4.08 seconds -> Processor[0]
3 Ball hit the ground at t = 4.13 seconds -> Processor[1]
4 Ball hit the ground at t = 4.13 seconds -> Processor[1]
5 Ball hit the ground at t = 7.76 seconds -> Processor[0]
6 Ball hit the ground at t = 7.76 seconds -> Processor[0]
7 Ball hit the ground at t = 7.95 seconds -> Processor[1]
8 Ball hit the ground at t = 7.95 seconds -> Processor[1]
9 Ball hit the ground at t = 11.06 seconds -> Processor[0]
10 Ball hit the ground at t = 11.06 seconds -> Processor[0]
[all …]
|
| H A D | ex41_b.out | 1 Ball hit the ground at t = 4.08 seconds -> Processor[0]
2 Ball hit the ground at t = 4.08 seconds -> Processor[0]
3 Ball hit the ground at t = 4.13 seconds -> Processor[1]
4 Ball hit the ground at t = 4.13 seconds -> Processor[1]
5 Ball hit the ground at t = 7.76 seconds -> Processor[0]
6 Ball hit the ground at t = 7.76 seconds -> Processor[0]
7 Ball hit the ground at t = 7.95 seconds -> Processor[1]
8 Ball hit the ground at t = 7.95 seconds -> Processor[1]
9 Ball hit the ground at t = 11.06 seconds -> Processor[0]
10 Ball hit the ground at t = 11.06 seconds -> Processor[0]
[all …]
|
| H A D | ex41_e.out | 1 Ball hit the ground at t = 4.08 seconds -> Processor[0]
2 Ball hit the ground at t = 4.08 seconds -> Processor[0]
3 Ball hit the ground at t = 4.13 seconds -> Processor[1]
4 Ball hit the ground at t = 4.13 seconds -> Processor[1]
5 Ball hit the ground at t = 7.76 seconds -> Processor[0]
6 Ball hit the ground at t = 7.76 seconds -> Processor[0]
7 Ball hit the ground at t = 7.95 seconds -> Processor[1]
8 Ball hit the ground at t = 7.95 seconds -> Processor[1]
9 Ball hit the ground at t = 11.06 seconds -> Processor[0]
10 Ball hit the ground at t = 11.06 seconds -> Processor[0]
[all …]
|
| H A D | ex41_g.out | 1 Ball hit the ground at t = 4.08 seconds -> Processor[0]
2 Ball hit the ground at t = 4.08 seconds -> Processor[0]
3 Ball hit the ground at t = 4.13 seconds -> Processor[1]
4 Ball hit the ground at t = 4.13 seconds -> Processor[1]
5 Ball hit the ground at t = 7.76 seconds -> Processor[0]
6 Ball hit the ground at t = 7.76 seconds -> Processor[0]
7 Ball hit the ground at t = 7.95 seconds -> Processor[1]
8 Ball hit the ground at t = 7.95 seconds -> Processor[1]
9 Ball hit the ground at t = 11.06 seconds -> Processor[0]
10 Ball hit the ground at t = 11.06 seconds -> Processor[0]
[all …]
|
| H A D | ex41_a.out | 1 Ball hit the ground at t = 4.08 seconds -> Processor[0]
2 Ball hit the ground at t = 4.08 seconds -> Processor[0]
3 Ball hit the ground at t = 4.13 seconds -> Processor[1]
4 Ball hit the ground at t = 4.13 seconds -> Processor[1]
5 Ball hit the ground at t = 7.76 seconds -> Processor[0]
6 Ball hit the ground at t = 7.76 seconds -> Processor[0]
7 Ball hit the ground at t = 7.95 seconds -> Processor[1]
8 Ball hit the ground at t = 7.95 seconds -> Processor[1]
9 Ball hit the ground at t = 11.06 seconds -> Processor[0]
10 Ball hit the ground at t = 11.06 seconds -> Processor[0]
[all …]
|
| H A D | ex31_4.out | 2 …TSAdapt glee glee 0:35 step 0 accepted t=0 + 1.000e-02 dt=1.000e-01 wlte=1.44e-06 wlt…
3 …TSAdapt glee glee 0:35 step 1 accepted t=0.01 + 1.000e-01 dt=2.296e-01 wlte=0.0112 wltea…
4 …TSAdapt glee glee 0:35 step 2 accepted t=0.11 + 2.296e-01 dt=3.604e-01 wlte=0.0451 wltea…
5 …TSAdapt glee glee 0:35 step 3 rejected t=0.33956 + 3.604e-01 dt=1.899e-01 wlte= 2.93 wltea=…
6 …TSAdapt glee glee 0:35 step 3 accepted t=0.33956 + 1.899e-01 dt=2.195e-01 wlte=0.138 wltea=…
7 …TSAdapt glee glee 0:35 step 4 rejected t=0.529423 + 2.195e-01 dt=1.604e-01 wlte=0.773 wltea=…
8 …TSAdapt glee glee 0:35 step 4 accepted t=0.529423 + 1.604e-01 dt=1.473e-01 wlte=0.319 wltea=…
9 …TSAdapt glee glee 0:35 step 5 rejected t=0.689836 + 1.473e-01 dt=1.196e-01 wlte=0.487 wltea=…
10 …TSAdapt glee glee 0:35 step 5 rejected t=0.689836 + 1.196e-01 dt=1.025e-01 wlte=0.397 wltea=…
11 …TSAdapt glee glee 0:35 step 5 rejected t=0.689836 + 1.025e-01 dt=8.995e-02 wlte=0.364 wltea=…
[all …]
|
| H A D | ex31_4_alt_2.out | 2 …TSAdapt glee glee 0:35 step 0 accepted t=0 + 1.000e-02 dt=1.000e-01 wlte=1.44e-06 wlt…
3 …TSAdapt glee glee 0:35 step 1 accepted t=0.01 + 1.000e-01 dt=2.296e-01 wlte=0.0112 wltea…
4 …TSAdapt glee glee 0:35 step 2 accepted t=0.11 + 2.296e-01 dt=3.604e-01 wlte=0.0451 wltea…
5 …TSAdapt glee glee 0:35 step 3 rejected t=0.33956 + 3.604e-01 dt=1.899e-01 wlte= 2.93 wltea=…
6 …TSAdapt glee glee 0:35 step 3 accepted t=0.33956 + 1.899e-01 dt=2.195e-01 wlte=0.138 wltea=…
7 …TSAdapt glee glee 0:35 step 4 rejected t=0.529423 + 2.195e-01 dt=1.604e-01 wlte=0.773 wltea=…
8 …TSAdapt glee glee 0:35 step 4 accepted t=0.529423 + 1.604e-01 dt=1.473e-01 wlte=0.319 wltea=…
9 …TSAdapt glee glee 0:35 step 5 rejected t=0.689836 + 1.473e-01 dt=1.196e-01 wlte=0.487 wltea=…
10 …TSAdapt glee glee 0:35 step 5 rejected t=0.689836 + 1.196e-01 dt=1.025e-01 wlte=0.397 wltea=…
11 …TSAdapt glee glee 0:35 step 5 rejected t=0.689836 + 1.025e-01 dt=8.995e-02 wlte=0.364 wltea=…
[all …]
|
| H A D | ex31_4_alt.out | 2 …TSAdapt glee glee 0:35 step 0 accepted t=0 + 1.000e-02 dt=1.000e-01 wlte=1.44e-06 wlt…
3 …TSAdapt glee glee 0:35 step 1 accepted t=0.01 + 1.000e-01 dt=2.296e-01 wlte=0.0112 wltea…
4 …TSAdapt glee glee 0:35 step 2 accepted t=0.11 + 2.296e-01 dt=3.604e-01 wlte=0.0451 wltea…
5 …TSAdapt glee glee 0:35 step 3 rejected t=0.33956 + 3.604e-01 dt=1.899e-01 wlte= 2.93 wltea=…
6 …TSAdapt glee glee 0:35 step 3 accepted t=0.33956 + 1.899e-01 dt=2.195e-01 wlte=0.138 wltea=…
7 …TSAdapt glee glee 0:35 step 4 rejected t=0.529423 + 2.195e-01 dt=1.604e-01 wlte=0.773 wltea=…
8 …TSAdapt glee glee 0:35 step 4 accepted t=0.529423 + 1.604e-01 dt=1.473e-01 wlte=0.319 wltea=…
9 …TSAdapt glee glee 0:35 step 5 rejected t=0.689836 + 1.473e-01 dt=1.196e-01 wlte=0.487 wltea=…
10 …TSAdapt glee glee 0:35 step 5 rejected t=0.689836 + 1.196e-01 dt=1.025e-01 wlte=0.397 wltea=…
11 …TSAdapt glee glee 0:35 step 5 rejected t=0.689836 + 1.025e-01 dt=8.995e-02 wlte=0.364 wltea=…
[all …]
|
| /petsc/src/ts/event/tests/output/ |
| H A D | ex3span_adaptmonitor.out | 1 TSAdapt none beuler 0: step 0 accepted t=0.01 + 9.900e-02 dt=5.050e-02
2 TSAdapt none beuler 0: step 1 accepted t=0.109 + 5.050e-02 dt=5.050e-02
3 TSAdapt none beuler 0: step 2 accepted t=0.1595 + 5.050e-02 dt=9.900e-02
4 TSAdapt none beuler 0: step 3 accepted t=0.21 + 9.900e-02 dt=9.900e-02
5 TSAdapt none beuler 0: step 4 accepted t=0.309 + 9.900e-02 dt=9.900e-02
6 TSAdapt none beuler 0: step 5 accepted t=0.408 + 9.900e-02 dt=9.900e-02
7 TSAdapt none beuler 0: step 6 accepted t=0.507 + 9.900e-02 dt=9.900e-02
8 TSAdapt none beuler 0: step 7 accepted t=0.606 + 9.900e-02 dt=9.900e-02
9 TSAdapt none beuler 0: step 8 accepted t=0.705 + 9.900e-02 dt=9.900e-02
10 TSAdapt none beuler 0: step 9 accepted t=0.804 + 9.900e-02 dt=5.350e-02
[all …]
|
| /petsc/src/mat/impls/baij/seq/ |
| H A D | baijsolvtran7.c | 12 PetscScalar s1, s2, s3, s4, s5, s6, s7, x1, x2, x3, x4, x5, x6, x7, *x, *t; in MatSolveTranspose_SeqBAIJ_7_inplace() local
18 t = a->solve_work; in MatSolveTranspose_SeqBAIJ_7_inplace()
29 t[ii] = b[ic]; in MatSolveTranspose_SeqBAIJ_7_inplace()
30 t[ii + 1] = b[ic + 1]; in MatSolveTranspose_SeqBAIJ_7_inplace()
31 t[ii + 2] = b[ic + 2]; in MatSolveTranspose_SeqBAIJ_7_inplace()
32 t[ii + 3] = b[ic + 3]; in MatSolveTranspose_SeqBAIJ_7_inplace()
33 t[ii + 4] = b[ic + 4]; in MatSolveTranspose_SeqBAIJ_7_inplace()
34 t[ii + 5] = b[ic + 5]; in MatSolveTranspose_SeqBAIJ_7_inplace()
35 t[ii + 6] = b[ic + 6]; in MatSolveTranspose_SeqBAIJ_7_inplace()
44 x1 = t[idx]; in MatSolveTranspose_SeqBAIJ_7_inplace()
[all …]
|
| H A D | baijsolvtran6.c | 12 PetscScalar s1, s2, s3, s4, s5, s6, x1, x2, x3, x4, x5, x6, *x, *t; in MatSolveTranspose_SeqBAIJ_6_inplace() local
18 t = a->solve_work; in MatSolveTranspose_SeqBAIJ_6_inplace()
29 t[ii] = b[ic]; in MatSolveTranspose_SeqBAIJ_6_inplace()
30 t[ii + 1] = b[ic + 1]; in MatSolveTranspose_SeqBAIJ_6_inplace()
31 t[ii + 2] = b[ic + 2]; in MatSolveTranspose_SeqBAIJ_6_inplace()
32 t[ii + 3] = b[ic + 3]; in MatSolveTranspose_SeqBAIJ_6_inplace()
33 t[ii + 4] = b[ic + 4]; in MatSolveTranspose_SeqBAIJ_6_inplace()
34 t[ii + 5] = b[ic + 5]; in MatSolveTranspose_SeqBAIJ_6_inplace()
43 x1 = t[idx]; in MatSolveTranspose_SeqBAIJ_6_inplace()
44 x2 = t[1 + idx]; in MatSolveTranspose_SeqBAIJ_6_inplace()
[all …]
|
| H A D | baijsolvtran5.c | 12 PetscScalar s1, s2, s3, s4, s5, x1, x2, x3, x4, x5, *x, *t; in MatSolveTranspose_SeqBAIJ_5_inplace() local
18 t = a->solve_work; in MatSolveTranspose_SeqBAIJ_5_inplace()
29 t[ii] = b[ic]; in MatSolveTranspose_SeqBAIJ_5_inplace()
30 t[ii + 1] = b[ic + 1]; in MatSolveTranspose_SeqBAIJ_5_inplace()
31 t[ii + 2] = b[ic + 2]; in MatSolveTranspose_SeqBAIJ_5_inplace()
32 t[ii + 3] = b[ic + 3]; in MatSolveTranspose_SeqBAIJ_5_inplace()
33 t[ii + 4] = b[ic + 4]; in MatSolveTranspose_SeqBAIJ_5_inplace()
42 x1 = t[idx]; in MatSolveTranspose_SeqBAIJ_5_inplace()
43 x2 = t[1 + idx]; in MatSolveTranspose_SeqBAIJ_5_inplace()
44 x3 = t[2 + idx]; in MatSolveTranspose_SeqBAIJ_5_inplace()
[all …]
|