Lines Matching refs:f
73 PetscScalar *f; in RHSFunction_Hull1972A1() local
78 PetscCall(VecGetArray(F, &f)); in RHSFunction_Hull1972A1()
79 f[0] = -y[0]; in RHSFunction_Hull1972A1()
81 PetscCall(VecRestoreArray(F, &f)); in RHSFunction_Hull1972A1()
103 PetscScalar *f; in IFunction_Hull1972A1() local
107 PetscCall(VecGetArray(F, &f)); in IFunction_Hull1972A1()
108 f[0] = -y[0]; in IFunction_Hull1972A1()
110 PetscCall(VecRestoreArray(F, &f)); in IFunction_Hull1972A1()
136 PetscScalar *f; in RHSFunction_Hull1972A2() local
140 PetscCall(VecGetArray(F, &f)); in RHSFunction_Hull1972A2()
141 f[0] = -0.5 * y[0] * y[0] * y[0]; in RHSFunction_Hull1972A2()
143 PetscCall(VecRestoreArray(F, &f)); in RHSFunction_Hull1972A2()
165 PetscScalar *f; in IFunction_Hull1972A2() local
170 PetscCall(VecGetArray(F, &f)); in IFunction_Hull1972A2()
171 f[0] = -0.5 * y[0] * y[0] * y[0]; in IFunction_Hull1972A2()
173 PetscCall(VecRestoreArray(F, &f)); in IFunction_Hull1972A2()
200 PetscScalar *f; in RHSFunction_Hull1972A3() local
204 PetscCall(VecGetArray(F, &f)); in RHSFunction_Hull1972A3()
205 f[0] = y[0] * PetscCosReal(t); in RHSFunction_Hull1972A3()
207 PetscCall(VecRestoreArray(F, &f)); in RHSFunction_Hull1972A3()
228 PetscScalar *f; in IFunction_Hull1972A3() local
233 PetscCall(VecGetArray(F, &f)); in IFunction_Hull1972A3()
234 f[0] = y[0] * PetscCosReal(t); in IFunction_Hull1972A3()
236 PetscCall(VecRestoreArray(F, &f)); in IFunction_Hull1972A3()
261 PetscScalar *f; in RHSFunction_Hull1972A4() local
266 PetscCall(VecGetArray(F, &f)); in RHSFunction_Hull1972A4()
267 f[0] = (0.25 * y[0]) * (1.0 - 0.05 * y[0]); in RHSFunction_Hull1972A4()
269 PetscCall(VecRestoreArray(F, &f)); in RHSFunction_Hull1972A4()
291 PetscScalar *f; in IFunction_Hull1972A4() local
296 PetscCall(VecGetArray(F, &f)); in IFunction_Hull1972A4()
297 f[0] = (0.25 * y[0]) * (1.0 - 0.05 * y[0]); in IFunction_Hull1972A4()
299 PetscCall(VecRestoreArray(F, &f)); in IFunction_Hull1972A4()
325 PetscScalar *f; in RHSFunction_Hull1972A5() local
330 PetscCall(VecGetArray(F, &f)); in RHSFunction_Hull1972A5()
331 f[0] = (y[0] - t) / (y[0] + t); in RHSFunction_Hull1972A5()
333 PetscCall(VecRestoreArray(F, &f)); in RHSFunction_Hull1972A5()
355 PetscScalar *f; in IFunction_Hull1972A5() local
360 PetscCall(VecGetArray(F, &f)); in IFunction_Hull1972A5()
361 f[0] = (y[0] - t) / (y[0] + t); in IFunction_Hull1972A5()
363 PetscCall(VecRestoreArray(F, &f)); in IFunction_Hull1972A5()
389 PetscScalar *f; in RHSFunction_Hull1972B1() local
394 PetscCall(VecGetArray(F, &f)); in RHSFunction_Hull1972B1()
395 f[0] = 2.0 * (y[0] - y[0] * y[1]); in RHSFunction_Hull1972B1()
396 f[1] = -(y[1] - y[0] * y[1]); in RHSFunction_Hull1972B1()
398 PetscCall(VecRestoreArray(F, &f)); in RHSFunction_Hull1972B1()
404 PetscScalar *f; in IFunction_Hull1972B1() local
409 PetscCall(VecGetArray(F, &f)); in IFunction_Hull1972B1()
410 f[0] = 2.0 * (y[0] - y[0] * y[1]); in IFunction_Hull1972B1()
411 f[1] = -(y[1] - y[0] * y[1]); in IFunction_Hull1972B1()
413 PetscCall(VecRestoreArray(F, &f)); in IFunction_Hull1972B1()
442 PetscScalar *f; in RHSFunction_Hull1972B2() local
447 PetscCall(VecGetArray(F, &f)); in RHSFunction_Hull1972B2()
448 f[0] = -y[0] + y[1]; in RHSFunction_Hull1972B2()
449 f[1] = y[0] - 2.0 * y[1] + y[2]; in RHSFunction_Hull1972B2()
450 f[2] = y[1] - y[2]; in RHSFunction_Hull1972B2()
452 PetscCall(VecRestoreArray(F, &f)); in RHSFunction_Hull1972B2()
458 PetscScalar *f; in IFunction_Hull1972B2() local
463 PetscCall(VecGetArray(F, &f)); in IFunction_Hull1972B2()
464 f[0] = -y[0] + y[1]; in IFunction_Hull1972B2()
465 f[1] = y[0] - 2.0 * y[1] + y[2]; in IFunction_Hull1972B2()
466 f[2] = y[1] - y[2]; in IFunction_Hull1972B2()
468 PetscCall(VecRestoreArray(F, &f)); in IFunction_Hull1972B2()
502 PetscScalar *f; in RHSFunction_Hull1972B3() local
507 PetscCall(VecGetArray(F, &f)); in RHSFunction_Hull1972B3()
508 f[0] = -y[0]; in RHSFunction_Hull1972B3()
509 f[1] = y[0] - y[1] * y[1]; in RHSFunction_Hull1972B3()
510 f[2] = y[1] * y[1]; in RHSFunction_Hull1972B3()
512 PetscCall(VecRestoreArray(F, &f)); in RHSFunction_Hull1972B3()
518 PetscScalar *f; in IFunction_Hull1972B3() local
523 PetscCall(VecGetArray(F, &f)); in IFunction_Hull1972B3()
524 f[0] = -y[0]; in IFunction_Hull1972B3()
525 f[1] = y[0] - y[1] * y[1]; in IFunction_Hull1972B3()
526 f[2] = y[1] * y[1]; in IFunction_Hull1972B3()
528 PetscCall(VecRestoreArray(F, &f)); in IFunction_Hull1972B3()
562 PetscScalar *f; in RHSFunction_Hull1972B4() local
567 PetscCall(VecGetArray(F, &f)); in RHSFunction_Hull1972B4()
568 f[0] = -y[1] - y[0] * y[2] / PetscSqrtScalar(y[0] * y[0] + y[1] * y[1]); in RHSFunction_Hull1972B4()
569 f[1] = y[0] - y[1] * y[2] / PetscSqrtScalar(y[0] * y[0] + y[1] * y[1]); in RHSFunction_Hull1972B4()
570 f[2] = y[0] / PetscSqrtScalar(y[0] * y[0] + y[1] * y[1]); in RHSFunction_Hull1972B4()
572 PetscCall(VecRestoreArray(F, &f)); in RHSFunction_Hull1972B4()
578 PetscScalar *f; in IFunction_Hull1972B4() local
583 PetscCall(VecGetArray(F, &f)); in IFunction_Hull1972B4()
584 f[0] = -y[1] - y[0] * y[2] / PetscSqrtScalar(y[0] * y[0] + y[1] * y[1]); in IFunction_Hull1972B4()
585 f[1] = y[0] - y[1] * y[2] / PetscSqrtScalar(y[0] * y[0] + y[1] * y[1]); in IFunction_Hull1972B4()
586 f[2] = y[0] / PetscSqrtScalar(y[0] * y[0] + y[1] * y[1]); in IFunction_Hull1972B4()
588 PetscCall(VecRestoreArray(F, &f)); in IFunction_Hull1972B4()
624 PetscScalar *f; in RHSFunction_Hull1972B5() local
629 PetscCall(VecGetArray(F, &f)); in RHSFunction_Hull1972B5()
630 f[0] = y[1] * y[2]; in RHSFunction_Hull1972B5()
631 f[1] = -y[0] * y[2]; in RHSFunction_Hull1972B5()
632 f[2] = -0.51 * y[0] * y[1]; in RHSFunction_Hull1972B5()
634 PetscCall(VecRestoreArray(F, &f)); in RHSFunction_Hull1972B5()
640 PetscScalar *f; in IFunction_Hull1972B5() local
645 PetscCall(VecGetArray(F, &f)); in IFunction_Hull1972B5()
646 f[0] = y[1] * y[2]; in IFunction_Hull1972B5()
647 f[1] = -y[0] * y[2]; in IFunction_Hull1972B5()
648 f[2] = -0.51 * y[0] * y[1]; in IFunction_Hull1972B5()
650 PetscCall(VecRestoreArray(F, &f)); in IFunction_Hull1972B5()
684 PetscScalar *f; in RHSFunction_Kulikov2013I() local
689 PetscCall(VecGetArray(F, &f)); in RHSFunction_Kulikov2013I()
690 f[0] = 2. * t * PetscPowScalar(y[1], 1. / 5.) * y[3]; in RHSFunction_Kulikov2013I()
691 f[1] = 10. * t * y[3] * PetscExpScalar(5.0 * (y[2] - 1.)); in RHSFunction_Kulikov2013I()
692 f[2] = 2. * t * y[3]; in RHSFunction_Kulikov2013I()
693 f[3] = -2. * t * PetscLogScalar(y[0]); in RHSFunction_Kulikov2013I()
695 PetscCall(VecRestoreArray(F, &f)); in RHSFunction_Kulikov2013I()
737 PetscScalar *f; in IFunction_Kulikov2013I() local
742 PetscCall(VecGetArray(F, &f)); in IFunction_Kulikov2013I()
743 f[0] = 2. * t * PetscPowScalar(y[1], 1. / 5.) * y[3]; in IFunction_Kulikov2013I()
744 f[1] = 10. * t * y[3] * PetscExpScalar(5.0 * (y[2] - 1.)); in IFunction_Kulikov2013I()
745 f[2] = 2. * t * y[3]; in IFunction_Kulikov2013I()
746 f[3] = -2. * t * PetscLogScalar(y[0]); in IFunction_Kulikov2013I()
748 PetscCall(VecRestoreArray(F, &f)); in IFunction_Kulikov2013I()
794 PetscScalar *f; in RHSFunction_Hull1972C1() local
801 PetscCall(VecGetArray(F, &f)); in RHSFunction_Hull1972C1()
802 f[0] = -y[0]; in RHSFunction_Hull1972C1()
803 for (i = 1; i < N - 1; i++) f[i] = y[i - 1] - y[i]; in RHSFunction_Hull1972C1()
804 f[N - 1] = y[N - 2]; in RHSFunction_Hull1972C1()
806 PetscCall(VecRestoreArray(F, &f)); in RHSFunction_Hull1972C1()
812 PetscScalar *f; in IFunction_Hull1972C1() local
819 PetscCall(VecGetArray(F, &f)); in IFunction_Hull1972C1()
820 f[0] = -y[0]; in IFunction_Hull1972C1()
821 for (i = 1; i < N - 1; i++) f[i] = y[i - 1] - y[i]; in IFunction_Hull1972C1()
822 f[N - 1] = y[N - 2]; in IFunction_Hull1972C1()
824 PetscCall(VecRestoreArray(F, &f)); in IFunction_Hull1972C1()
869 PetscScalar *f; in RHSFunction_Hull1972C2() local
875 PetscCall(VecGetArray(F, &f)); in RHSFunction_Hull1972C2()
876 f[0] = -y[0]; in RHSFunction_Hull1972C2()
877 for (i = 1; i < N - 1; i++) f[i] = (PetscReal)i * y[i - 1] - (PetscReal)(i + 1) * y[i]; in RHSFunction_Hull1972C2()
878 f[N - 1] = (PetscReal)(N - 1) * y[N - 2]; in RHSFunction_Hull1972C2()
880 PetscCall(VecRestoreArray(F, &f)); in RHSFunction_Hull1972C2()
886 PetscScalar *f; in IFunction_Hull1972C2() local
893 PetscCall(VecGetArray(F, &f)); in IFunction_Hull1972C2()
894 f[0] = -y[0]; in IFunction_Hull1972C2()
895 for (i = 1; i < N - 1; i++) f[i] = (PetscReal)i * y[i - 1] - (PetscReal)(i + 1) * y[i]; in IFunction_Hull1972C2()
896 f[N - 1] = (PetscReal)(N - 1) * y[N - 2]; in IFunction_Hull1972C2()
898 PetscCall(VecRestoreArray(F, &f)); in IFunction_Hull1972C2()
942 PetscScalar *f; in RHSFunction_Hull1972C34() local
949 PetscCall(VecGetArray(F, &f)); in RHSFunction_Hull1972C34()
950 f[0] = -2.0 * y[0] + y[1]; in RHSFunction_Hull1972C34()
951 for (i = 1; i < N - 1; i++) f[i] = y[i - 1] - 2.0 * y[i] + y[i + 1]; in RHSFunction_Hull1972C34()
952 f[N - 1] = y[N - 2] - 2.0 * y[N - 1]; in RHSFunction_Hull1972C34()
954 PetscCall(VecRestoreArray(F, &f)); in RHSFunction_Hull1972C34()
960 PetscScalar *f; in IFunction_Hull1972C34() local
967 PetscCall(VecGetArray(F, &f)); in IFunction_Hull1972C34()
968 f[0] = -2.0 * y[0] + y[1]; in IFunction_Hull1972C34()
969 for (i = 1; i < N - 1; i++) f[i] = y[i - 1] - 2.0 * y[i] + y[i + 1]; in IFunction_Hull1972C34()
970 f[N - 1] = y[N - 2] - 2.0 * y[N - 1]; in IFunction_Hull1972C34()
972 PetscCall(VecRestoreArray(F, &f)); in IFunction_Hull1972C34()