Lines Matching +full:- +full:t
4 petsc.h - base PETSc routines petscvec.h - vectors
5 petscsys.h - system routines petscmat.h - matrices
6 petscis.h - index sets petscksp.h - Krylov subspace methods
7 petscviewer.h - viewers petscpc.h - preconditioners
23 static char help[] = "Finds the nonlinear least-squares solution to the model \n\
24 y = exp[-b1*x]/(b2+b3*x) + e \n";
29 /* User-defined application context */
32 PetscReal t[NOBSERVATIONS]; /* array of independent variables of observation */
45 /*--------------------------------------------------------------------*/
54 AppCtx user; /* user-defined work context */
102 /*--------------------------------------------------------------------*/
108 PetscReal *y = user->y, *f, *t = user->t;
114 for (i = 0; i < NOBSERVATIONS; i++) f[i] = y[i] - PetscExpScalar(-x[0] * t[i]) / (x[1] + x[2] * t[i]);
121 /*------------------------------------------------------------*/
128 PetscReal *t = user->t;
134 base = PetscExpScalar(-x[0] * t[i]) / (x[1] + x[2] * t[i]);
136 user->j[i][0] = t[i] * base;
137 user->j[i][1] = base / (x[1] + x[2] * t[i]);
138 user->j[i][2] = base * t[i] / (x[1] + x[2] * t[i]);
142 PetscCall(MatSetValues(J, NOBSERVATIONS, user->idm, NPARAMETERS, user->idn, (PetscReal *)user->j, INSERT_VALUES));
151 /* ------------------------------------------------------------ */
165 /* ---------------------------------------------------------------------- */
168 PetscReal *t = user->t, *y = user->y;
173 t[i++] = 0.5000;
175 t[i++] = 0.6250;
177 t[i++] = 0.7500;
179 t[i++] = 0.8750;
181 t[i++] = 1.0000;
183 t[i++] = 1.2500;
185 t[i++] = 1.7500;
187 t[i++] = 2.2500;
189 t[i++] = 1.7500;
191 t[i++] = 2.2500;
193 t[i++] = 2.7500;
195 t[i++] = 3.2500;
197 t[i++] = 3.7500;
199 t[i++] = 4.2500;
201 t[i++] = 4.7500;
203 t[i++] = 5.2500;
205 t[i++] = 5.7500;
207 t[i++] = 0.5000;
209 t[i++] = 0.6250;
211 t[i++] = 0.7500;
213 t[i++] = 0.8750;
215 t[i++] = 1.0000;
217 t[i++] = 1.2500;
219 t[i++] = 1.7500;
221 t[i++] = 2.2500;
223 t[i++] = 1.7500;
225 t[i++] = 2.2500;
227 t[i++] = 2.7500;
229 t[i++] = 3.2500;
231 t[i++] = 3.7500;
233 t[i++] = 4.2500;
235 t[i++] = 4.7500;
237 t[i++] = 5.2500;
239 t[i++] = 5.7500;
241 t[i++] = 0.5000;
243 t[i++] = 0.6250;
245 t[i++] = 0.7500;
247 t[i++] = 0.8750;
249 t[i++] = 1.0000;
251 t[i++] = 1.2500;
253 t[i++] = 1.7500;
255 t[i++] = 2.2500;
257 t[i++] = 1.7500;
259 t[i++] = 2.2500;
261 t[i++] = 2.7500;
263 t[i++] = 3.2500;
265 t[i++] = 3.7500;
267 t[i++] = 4.2500;
269 t[i++] = 4.7500;
271 t[i++] = 5.2500;
273 t[i++] = 5.7500;
275 t[i++] = 0.5000;
277 t[i++] = 0.6250;
279 t[i++] = 0.7500;
281 t[i++] = 0.8750;
283 t[i++] = 1.0000;
285 t[i++] = 1.2500;
287 t[i++] = 1.7500;
289 t[i++] = 2.2500;
291 t[i++] = 1.7500;
293 t[i++] = 2.2500;
295 t[i++] = 2.7500;
297 t[i++] = 3.2500;
299 t[i++] = 3.7500;
301 t[i++] = 4.2500;
303 t[i++] = 4.7500;
305 t[i++] = 5.2500;
307 t[i++] = 5.7500;
309 t[i++] = .5000;
311 t[i++] = .7500;
313 t[i++] = 1.5000;
315 t[i++] = 3.0000;
317 t[i++] = 3.0000;
319 t[i++] = 3.0000;
321 t[i++] = 6.0000;
323 t[i++] = .5000;
325 t[i++] = .7500;
327 t[i++] = 1.5000;
329 t[i++] = 3.0000;
331 t[i++] = 3.0000;
333 t[i++] = 3.0000;
335 t[i++] = 6.0000;
337 t[i++] = .5000;
339 t[i++] = .7500;
341 t[i++] = 1.5000;
343 t[i++] = 3.0000;
345 t[i++] = 3.0000;
347 t[i++] = 3.0000;
349 t[i++] = 6.0000;
351 t[i++] = .5000;
353 t[i++] = .7500;
355 t[i++] = 1.5000;
357 t[i++] = 3.0000;
359 t[i++] = 6.0000;
361 t[i++] = 3.0000;
363 t[i++] = 3.0000;
365 t[i++] = 6.0000;
367 t[i++] = .5000;
369 t[i++] = .7500;
371 t[i++] = 1.0000;
373 t[i++] = 1.5000;
375 t[i++] = 2.0000;
377 t[i++] = 2.0000;
379 t[i++] = 2.5000;
381 t[i++] = 3.0000;
383 t[i++] = 4.0000;
385 t[i++] = 5.0000;
387 t[i++] = 6.0000;
389 t[i++] = .5000;
391 t[i++] = .7500;
393 t[i++] = 1.0000;
395 t[i++] = 1.5000;
397 t[i++] = 2.0000;
399 t[i++] = 2.0000;
401 t[i++] = 2.5000;
403 t[i++] = 3.0000;
405 t[i++] = 4.0000;
407 t[i++] = 5.0000;
409 t[i++] = 6.0000;
411 t[i++] = .5000;
413 t[i++] = .7500;
415 t[i++] = 1.0000;
417 t[i++] = 1.5000;
419 t[i++] = 2.0000;
421 t[i++] = 2.0000;
423 t[i++] = 2.5000;
425 t[i++] = 3.0000;
427 t[i++] = 4.0000;
429 t[i++] = 5.0000;
431 t[i++] = 6.0000;
433 t[i++] = .5000;
435 t[i++] = .6250;
437 t[i++] = .7500;
439 t[i++] = .8750;
441 t[i++] = 1.0000;
443 t[i++] = 1.2500;
445 t[i++] = 2.2500;
447 t[i++] = 2.2500;
449 t[i++] = 2.7500;
451 t[i++] = 3.2500;
453 t[i++] = 3.7500;
455 t[i++] = 4.2500;
457 t[i++] = 4.7500;
459 t[i++] = 5.2500;
461 t[i++] = 5.7500;
463 t[i++] = 3.0000;
465 t[i++] = 3.0000;
467 t[i++] = 3.0000;
469 t[i++] = 3.0000;
471 t[i++] = 3.0000;
473 t[i++] = 3.0000;
475 t[i++] = .5000;
477 t[i++] = .7500;
479 t[i++] = 1.0000;
481 t[i++] = 1.5000;
483 t[i++] = 2.0000;
485 t[i++] = 2.5000;
487 t[i++] = 2.0000;
489 t[i++] = 2.5000;
491 t[i++] = 3.0000;
493 t[i++] = 4.0000;
495 t[i++] = 5.0000;
497 t[i++] = 6.0000;
499 t[i++] = .5000;
501 t[i++] = .7500;
503 t[i++] = 1.0000;
505 t[i++] = 1.5000;
507 t[i++] = 2.0000;
509 t[i++] = 2.5000;
511 t[i++] = 2.0000;
513 t[i++] = 2.5000;
515 t[i++] = 3.0000;
517 t[i++] = 4.0000;
519 t[i++] = 5.0000;
521 t[i++] = 6.0000;
523 t[i++] = .5000;
525 t[i++] = .7500;
527 t[i++] = 1.0000;
529 t[i++] = 1.5000;
531 t[i++] = 2.0000;
533 t[i++] = 2.5000;
535 t[i++] = 2.0000;
537 t[i++] = 2.5000;
539 t[i++] = 3.0000;
541 t[i++] = 4.0000;
543 t[i++] = 5.0000;
545 t[i++] = 6.0000;
547 t[i++] = 3.0000;
549 t[i++] = .5000;
551 t[i++] = .7500;
553 t[i++] = 1.5000;
555 t[i++] = 3.0000;
557 t[i++] = 6.0000;
559 t[i++] = 3.0000;
561 t[i++] = 6.0000;
563 t[i++] = 3.0000;
565 t[i++] = 3.0000;
567 t[i++] = 3.0000;
569 t[i++] = 1.7500;
571 t[i++] = 1.7500;
573 t[i++] = .5000;
575 t[i++] = .7500;
577 t[i++] = 1.7500;
579 t[i++] = 1.7500;
581 t[i++] = 2.7500;
583 t[i++] = 3.7500;
585 t[i++] = 1.7500;
587 t[i++] = 1.7500;
589 t[i++] = .5000;
591 t[i++] = .7500;
593 t[i++] = 2.7500;
595 t[i++] = 3.7500;
597 t[i++] = 1.7500;
599 t[i++] = 1.7500;
609 args: -tao_monitor_short -tao_max_it 100 -tao_type pounders -tao_gatol 1.e-5
613 args: -tao_monitor_short -tao_max_it 100 -tao_type brgn -tao_brgn_regularization_type l2prox -tao_brgn_regularizer_weight 1e-4 -tao_gatol 1.e-5 -tao_brgn_subsolver_tao_monitor_short
617 args: -tao_monitor_short -tao_max_it 100 -tao_type brgn -tao_brgn_regularization_type l1dict -tao_brgn_regularizer_weight 1e-4 -tao_brgn_l1_smooth_epsilon 1e-6 -tao_gatol 1.e-5 -tao_brgn_subsolver_tao_monitor_short
621 args: -tao_monitor_short -tao_max_it 100 -tao_type brgn -tao_brgn_regularization_type lm -tao_gatol 1.e-5 -tao_brgn_subsolver_tao_type bnls -tao_brgn_subsolver_tao_monitor_short