43   const PetscScalar *u, *udot;  in IFunctionPassive()  local
48   PetscCall(VecGetArrayRead(Udot, &udot));  in IFunctionPassive()
50   f[0] = udot[0] + ctx->k * u[0] * u[1];  in IFunctionPassive()
51   f[1] = udot[1] + ctx->k * u[0] * u[1];  in IFunctionPassive()
52   f[2] = udot[2] - ctx->k * u[0] * u[1];  in IFunctionPassive()
54   PetscCall(VecRestoreArrayRead(Udot, &udot));  in IFunctionPassive()
65   const PetscScalar *u, *udot;  in IFunctionActive1()  local
73   PetscCall(VecGetArrayRead(Udot, &udot));  in IFunctionActive1()
81   f_a[0] = udot[0] + ctx->k * u_a[0] * u_a[1];  in IFunctionActive1()
82   f_a[1] = udot[1] + ctx->k * u_a[0] * u_a[1];  in IFunctionActive1()
83   f_a[2] = udot[2] - ctx->k * u_a[0] * u_a[1];  in IFunctionActive1()
91   PetscCall(VecRestoreArrayRead(Udot, &udot));  in IFunctionActive1()
102   const PetscScalar *u, *udot;  in IFunctionActive2()  local
110   PetscCall(VecGetArrayRead(Udot, &udot));  in IFunctionActive2()
115   udot_a[0] <<= udot[0];  in IFunctionActive2()
116   udot_a[1] <<= udot[1];  in IFunctionActive2()
117   udot_a[2] <<= udot[2]; /* Mark independence */  in IFunctionActive2()
128   PetscCall(VecRestoreArrayRead(Udot, &udot));  in IFunctionActive2()