Lines Matching refs:d_eq
382 subroutine ShashiFormJacobian(an_r, d_eq) argument
405 PetscScalar d_eq(26, 26), H_molar(26)
485 d_eq(ib, jb) = 0.0d0
496 d_eq(1, 1) = -an_h(1)
497 d_eq(1, 2) = -an_h_additive
498 d_eq(1, 5) = -2
499 d_eq(1, 10) = -1
500 d_eq(1, 11) = -1
501 d_eq(1, 14) = -2
502 d_eq(1, 16) = -1
503 d_eq(1, 17) = -2
504 d_eq(1, 19) = -1
505 d_eq(1, 20) = -1
506 d_eq(1, 22) = -3
507 d_eq(1, 26) = -1
509 d_eq(2, 2) = -1*an_o_additive
510 d_eq(2, 3) = -2
511 d_eq(2, 4) = -2
512 d_eq(2, 5) = -1
513 d_eq(2, 8) = -1
514 d_eq(2, 9) = -1
515 d_eq(2, 10) = -1
516 d_eq(2, 12) = -1
517 d_eq(2, 13) = -1
518 d_eq(2, 15) = -2
519 d_eq(2, 16) = -2
520 d_eq(2, 20) = -1
521 d_eq(2, 22) = -1
522 d_eq(2, 23) = -1
523 d_eq(2, 24) = -2
524 d_eq(2, 25) = -1
525 d_eq(2, 26) = -1
527 d_eq(6, 6) = -2
528 d_eq(6, 7) = -1
529 d_eq(6, 9) = -1
530 d_eq(6, 12) = -2
531 d_eq(6, 15) = -1
532 d_eq(6, 23) = -1
534 d_eq(4, 1) = -an_c(1)
535 d_eq(4, 2) = -an_c_additive
536 d_eq(4, 4) = -1
537 d_eq(4, 13) = -1
538 d_eq(4, 17) = -2
539 d_eq(4, 18) = -1
540 d_eq(4, 19) = -1
541 d_eq(4, 20) = -1
575 d_eq(5, 1) = &
580 d_eq(5, jj) = d_eq(5, 1)
583 d_eq(5, 1) = d_eq(5, 1) + k_eq(11)*(an_r(3)**ai_o2)
585 d_eq(5, 3) = d_eq(5, 3) + k_eq(11)*(ai_o2*an_r(3)**(ai_o2 - 1)) &
588 d_eq(5, 4) = d_eq(5, 4) - (i_cc*an_r(4)**(i_cc - 1))* &
590 d_eq(5, 5) = d_eq(5, 5) &
594 d_eq(3, 1) = -(an_r(4)**2)*(an_r(5)**3)*(pt/an_t)*(-1.0/an_t)
596 d_eq(3, jj) = d_eq(3, 1)
599 d_eq(3, 2) = d_eq(3, 2) + k_eq(12)*(an_r(3)**3)
601 d_eq(3, 3) = d_eq(3, 3) + k_eq(12)*(3*an_r(3)**2)*an_r(2)
603 d_eq(3, 4) = d_eq(3, 4) - 2*an_r(4)*(an_r(5)**3)*(pt/an_t)
605 d_eq(3, 5) = d_eq(3, 5) - 3*(an_r(5)**2)*(an_r(4)**2)*(pt/an_t)
609 d_eq(3, ii3) = 0.0d0
611 d_eq(3, 2) = 1.0d0
613 d_eq(7, 1) = pt*an_r(11)*(-1.0)/const2 &
617 d_eq(7, jj) = d_eq(7, 1)
620 d_eq(7, 11) = d_eq(7, 11) + pt/an_t
621 d_eq(7, 14) = d_eq(7, 14) &
624 d_eq(8, 1) = pt*an_r(8)*(-1.0)/const2 &
628 d_eq(8, jj) = d_eq(8, 1)
631 d_eq(8, 3) = d_eq(8, 3) &
633 d_eq(8, 8) = d_eq(8, 8) + pt/an_t
635 d_eq(9, 1) = pt*an_r(7)*(-1.0)/const2 &
639 d_eq(9, jj) = d_eq(9, 1)
642 d_eq(9, 7) = d_eq(9, 7) + pt/an_t
643 d_eq(9, 6) = d_eq(9, 6) &
646 d_eq(10, 1) = pt*an_r(10)*(-1.0)/const2 &
650 d_eq(10, jj) = d_eq(10, 1)
653 d_eq(10, 3) = d_eq(10, 3) &
656 d_eq(10, 10) = d_eq(10, 10) + pt/an_t
657 d_eq(10, 14) = d_eq(10, 14) &
661 d_eq(11, 1) = pt*an_r(9)*(-1.0)/const2 &
666 d_eq(11, jj) = d_eq(11, 1)
669 d_eq(11, 3) = d_eq(11, 3) &
672 d_eq(11, 6) = d_eq(11, 6) &
675 d_eq(11, 9) = d_eq(11, 9) + pt/an_t
677 d_eq(12, 1) = pt*an_r(5)*(-1.0)/const2 &
682 d_eq(12, jj) = d_eq(12, 1)
685 d_eq(12, 3) = d_eq(12, 3) &
689 d_eq(12, 5) = d_eq(12, 5) + pt/an_t
690 d_eq(12, 14) = d_eq(12, 14) &
693 d_eq(13, 1) = pt*an_r(4)*(-1.0)/const2 &
698 d_eq(13, jj) = d_eq(13, 1)
701 d_eq(13, 3) = d_eq(13, 3) &
705 d_eq(13, 4) = d_eq(13, 4) + pt/an_t
706 d_eq(13, 13) = d_eq(13, 13) &
709 d_eq(14, 1) = pt*an_r(15)*(-1.0)/const2 &
714 d_eq(14, jj) = d_eq(14, 1)
717 d_eq(14, 3) = d_eq(14, 3) &
720 d_eq(14, 9) = d_eq(14, 9) &
722 d_eq(14, 15) = d_eq(14, 15) + pt/an_t
724 d_eq(15, 1) = pt*an_r(16)*(-1.0)/const2 &
729 d_eq(15, jj) = d_eq(15, 1)
732 d_eq(15, 3) = d_eq(15, 3) &
734 d_eq(15, 14) = d_eq(15, 14) &
737 d_eq(15, 16) = d_eq(15, 16) + pt/an_t
739 d_eq(16, 1) = pt*an_r(12)*(-1.0)/const2 &
744 d_eq(16, jj) = d_eq(16, 1)
747 d_eq(16, 3) = d_eq(16, 3) &
751 d_eq(16, 6) = d_eq(16, 6) &
753 d_eq(16, 12) = d_eq(16, 12) + pt/an_t
758 d_eq(17, 1) = an_r(14)*an_r(18)*an_r(18)*(pt**3)*(-3/const_four) &
761 d_eq(17, jj) = d_eq(17, 1)
763 d_eq(17, 14) = d_eq(17, 14) + an_r(18)*an_r(18)*(pt**3)/const_cube
764 d_eq(17, 17) = d_eq(17, 17) - k_eq(15)*pt/an_t
765 d_eq(17, 18) = d_eq(17, 18) + 2*an_r(18)*an_r(14) &
768 d_eq(18, 1) = an_r(13)*an_r(13)*(pt**2)*(-2/const_cube) &
772 d_eq(18, jj) = d_eq(18, 1)
774 d_eq(18, 3) = d_eq(18, 3) &
776 d_eq(18, 13) = d_eq(18, 13) &
778 d_eq(18, 18) = d_eq(18, 18) - k_eq(16)*an_r(3) &
783 d_eq(19, 1) = an_r(3)*an_r(19)*(pt**2)*(-2/const_cube) &
786 d_eq(19, jj) = d_eq(19, 1)
788 d_eq(19, 13) = d_eq(19, 13) &
790 d_eq(19, 10) = d_eq(19, 10) &
792 d_eq(19, 3) = d_eq(19, 3) + an_r(19)*pt*pt/const2
793 d_eq(19, 19) = d_eq(19, 19) + an_r(3)*pt*pt/const2
796 d_eq(20, 1) = an_r(21)*an_r(20)*(pt**2)*(-2/const_cube) &
799 d_eq(20, jj) = d_eq(20, 1)
801 d_eq(20, 8) = d_eq(20, 8) &
803 d_eq(20, 19) = d_eq(20, 19) &
805 d_eq(20, 20) = d_eq(20, 20) + an_r(21)*pt*pt/const2
806 d_eq(20, 21) = d_eq(20, 21) + an_r(20)*pt*pt/const2
811 d_eq(21, 1) = an_r(21)*an_r(23)*(pt**2)*(-2/const_cube) &
814 d_eq(21, jj) = d_eq(21, 1)
816 d_eq(21, 7) = d_eq(21, 7) &
818 d_eq(21, 8) = d_eq(21, 8) &
820 d_eq(21, 21) = d_eq(21, 21) + an_r(23)*pt*pt/const2
821 d_eq(21, 23) = d_eq(21, 23) + an_r(21)*pt*pt/const2
825 d_eq(22, 1) = an_r(5)*an_r(11)*(pt**2)*(-2/const_cube) &
828 d_eq(22, jj) = d_eq(22, 1)
830 d_eq(22, 21) = d_eq(22, 21) &
832 d_eq(22, 22) = d_eq(22, 22) &
834 d_eq(22, 11) = d_eq(22, 11) + an_r(5)*pt*pt/(const2)
835 d_eq(22, 5) = d_eq(22, 5) + an_r(11)*pt*pt/(const2)
840 d_eq(23, 1) = an_r(24)*(pt)*(-1/const2) &
843 d_eq(23, jj) = d_eq(23, 1)
845 d_eq(23, 3) = d_eq(23, 3) &
847 d_eq(23, 21) = d_eq(23, 21) &
849 d_eq(23, 24) = d_eq(23, 24) + pt/(an_t)
853 d_eq(24, 1) = an_r(3)*an_r(25)*(pt**2)*(-2/const_cube) &
856 d_eq(24, jj) = d_eq(24, 1)
858 d_eq(24, 8) = d_eq(24, 8) &
860 d_eq(24, 24) = d_eq(24, 24) &
862 d_eq(24, 3) = d_eq(24, 3) + an_r(25)*pt*pt/const2
863 d_eq(24, 25) = d_eq(24, 25) + an_r(3)*pt*pt/const2
868 d_eq(25, 1) = an_r(26)*(pt)*(-1/const2) &
871 d_eq(25, jj) = d_eq(25, 1)
873 d_eq(25, 10) = d_eq(25, 10) &
875 d_eq(25, 21) = d_eq(25, 21) &
877 d_eq(25, 26) = d_eq(25, 26) + pt/(an_t)
881 d_eq(26, 20) = -1
882 d_eq(26, 22) = -1
883 d_eq(26, 23) = -1
884 d_eq(26, 21) = 1
885 d_eq(26, 24) = 1
886 d_eq(26, 25) = 1
887 d_eq(26, 26) = 1
891 write (44, *) i, j, d_eq(i, j)