Lines Matching refs:ratio
69 static PetscErrorCode TSMPRKGenerateTableau2(PetscInt ratio, PetscInt s, const PetscReal Abase[], c… in TSMPRKGenerateTableau2() argument
74 for (k = 0; k < ratio; k++) { in TSMPRKGenerateTableau2()
78 A1[(k * s + i) * ratio * s + k * s + j] = Abase[i * s + j]; in TSMPRKGenerateTableau2()
79 A2[(k * s + i) * ratio * s + k * s + j] = Abase[i * s + j] / ratio; in TSMPRKGenerateTableau2()
84 for (j = 0; j < s; j++) A2[(k * s + i) * ratio * s + l * s + j] = bbase[j] / ratio; in TSMPRKGenerateTableau2()
86 b1[k * s + j] = bbase[j] / ratio; in TSMPRKGenerateTableau2()
87 b2[k * s + j] = bbase[j] / ratio; in TSMPRKGenerateTableau2()
93 static PetscErrorCode TSMPRKGenerateTableau3(PetscInt ratio, PetscInt s, const PetscReal Abase[], c… in TSMPRKGenerateTableau3() argument
98 for (k = 0; k < ratio; k++) { /* diagonal blocks of size ratio*s by ratio*s */ in TSMPRKGenerateTableau3()
99 for (l = 0; l < ratio; l++) /* diagonal sub-blocks of size s by s */ in TSMPRKGenerateTableau3()
102 … A1[((k * ratio + l) * s + i) * ratio * ratio * s + (k * ratio + l) * s + j] = Abase[i * s + j]; in TSMPRKGenerateTableau3()
103 …A2[((k * ratio + l) * s + i) * ratio * ratio * s + (k * ratio + l) * s + j] = Abase[i * s + j] / r… in TSMPRKGenerateTableau3()
104 …A3[((k * ratio + l) * s + i) * ratio * ratio * s + (k * ratio + l) * s + j] = Abase[i * s + j] / r… in TSMPRKGenerateTableau3()
107 for (m = 0; m < ratio; m++) in TSMPRKGenerateTableau3()
108 for (n = 0; n < ratio; n++) in TSMPRKGenerateTableau3()
111 …A2[((k * ratio + m) * s + i) * ratio * ratio * s + (l * ratio + n) * s + j] = bbase[j] / ratio / r… in TSMPRKGenerateTableau3()
112 …A3[((k * ratio + m) * s + i) * ratio * ratio * s + (l * ratio + n) * s + j] = bbase[j] / ratio / r… in TSMPRKGenerateTableau3()
114 for (m = 0; m < ratio; m++) in TSMPRKGenerateTableau3()
117 … (j = 0; j < s; j++) A3[((k * ratio + m) * s + i) * ratio * ratio * s + (k * ratio + n) * s + j] =… in TSMPRKGenerateTableau3()
118 for (n = 0; n < ratio; n++) in TSMPRKGenerateTableau3()
120 b1[(k * ratio + n) * s + j] = bbase[j] / ratio / ratio; in TSMPRKGenerateTableau3()
121 b2[(k * ratio + n) * s + j] = bbase[j] / ratio / ratio; in TSMPRKGenerateTableau3()
122 b3[(k * ratio + n) * s + j] = bbase[j] / ratio / ratio; in TSMPRKGenerateTableau3()