Lines Matching refs:f_eq
186 subroutine ShashiFormFunction(an_r, f_eq) argument
198 PetscScalar an_r(26), k_eq(23), f_eq(26)
279 f_eq(1) = atom_h_init &
285 f_eq(2) = atom_o_init &
292 f_eq(3) = an_r(2) - 1.0d-150
294 f_eq(4) = atom_c_init &
309 f_eq(5) = k_eq(11)*an_r(1)*an_r(3)**a_io2 &
313 f_eq(6) = atom_n_init &
318 f_eq(7) = part_p(11) &
320 f_eq(8) = part_p(8) &
323 f_eq(9) = part_p(7) &
326 f_eq(10) = part_p(10) &
330 f_eq(11) = part_p(9) &
333 f_eq(12) = part_p(5) &
337 f_eq(13) = part_p(4) &
341 f_eq(14) = part_p(15) &
344 f_eq(15) = part_p(16) &
348 f_eq(16) = part_p(12) &
352 f_eq(17) = part_p(14)*part_p(18)**2 &
355 f_eq(18) = (part_p(13)**2) &
357 print *, f_eq(18), part_p(3), part_p(18), part_p(13), k_eq(16)
359 f_eq(19) = part_p(19)*part_p(3) - k_eq(17)*part_p(13)*part_p(10)
361 f_eq(20) = part_p(21)*part_p(20) - k_eq(18)*part_p(19)*part_p(8)
363 f_eq(21) = part_p(21)*part_p(23) - k_eq(19)*part_p(7)*part_p(8)
365 f_eq(22) = part_p(5)*part_p(11) - k_eq(20)*part_p(21)*part_p(22)
367 f_eq(23) = part_p(24) - k_eq(21)*part_p(21)*part_p(3)
369 f_eq(24) = part_p(3)*part_p(25) - k_eq(22)*part_p(24)*part_p(8)
371 f_eq(25) = part_p(26) - k_eq(23)*part_p(21)*part_p(10)
373 f_eq(26) = -(an_r(20) + an_r(22) + an_r(23)) &
377 write (44, *) i, f_eq(i)