/[escript]/trunk/paso/src/SchurComplement.c
ViewVC logotype

Diff of /trunk/paso/src/SchurComplement.c

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 4285 by jfenwick, Tue Jan 22 09:30:23 2013 UTC revision 4286 by caltinay, Thu Mar 7 04:28:11 2013 UTC
# Line 85  void Paso_Solver_updateIncompleteSchurCo Line 85  void Paso_Solver_updateIncompleteSchurCo
85                  where_p=(index_t*)bsearch(&col_FC,index_CC,index_CC_len,sizeof(index_t),Paso_comparIndex);                  where_p=(index_t*)bsearch(&col_FC,index_CC,index_CC_len,sizeof(index_t),Paso_comparIndex);
86                  if (where_p!=NULL) {                  if (where_p!=NULL) {
87                      iPtr_CC_2=iPtr_CC+(index_t)(where_p-index_CC);                      iPtr_CC_2=iPtr_CC+(index_t)(where_p-index_CC);
88                      /* this calculutes A_CF*invA_FF(i,col_CF) */                      /* this calculates A_CF*invA_FF(i,col_CF) */
89                      if (set_A) {                      if (set_A) {
90                         A_CF_11=A_CF->val[4*iPtr_CF  ];                         A_CF_11=A_CF->val[4*iPtr_CF  ];
91                         A_CF_21=A_CF->val[4*iPtr_CF+1];                         A_CF_21=A_CF->val[4*iPtr_CF+1];
# Line 135  void Paso_Solver_updateIncompleteSchurCo Line 135  void Paso_Solver_updateIncompleteSchurCo
135                  where_p=(index_t*)bsearch(&col_FC,index_CC,index_CC_len,sizeof(index_t),Paso_comparIndex);                  where_p=(index_t*)bsearch(&col_FC,index_CC,index_CC_len,sizeof(index_t),Paso_comparIndex);
136                  if (where_p!=NULL) {                  if (where_p!=NULL) {
137                      iPtr_CC_2=iPtr_CC+(index_t)(where_p-index_CC);                      iPtr_CC_2=iPtr_CC+(index_t)(where_p-index_CC);
138                      /* this calculutes A_CF*invA_FF(i,col_CF) */                      /* this calculates A_CF*invA_FF(i,col_CF) */
139                      if (set_A) {                      if (set_A) {
140                         A_CF_11=A_CF->val[9*iPtr_CF  ];                         A_CF_11=A_CF->val[9*iPtr_CF  ];
141                         A_CF_21=A_CF->val[9*iPtr_CF+1];                         A_CF_21=A_CF->val[9*iPtr_CF+1];
# Line 202  void Paso_Solver_updateIncompleteSchurCo Line 202  void Paso_Solver_updateIncompleteSchurCo
202   * Merge of development branch dev-02 back to main trunk on 2005-09-15   * Merge of development branch dev-02 back to main trunk on 2005-09-15
203   *   *
204   * Revision 1.1.2.1  2005/09/05 06:29:50  gross   * Revision 1.1.2.1  2005/09/05 06:29:50  gross
205   * These files have been extracted from finley to define a stand alone libray for iterative   * These files have been extracted from finley to define a stand alone library for iterative
206   * linear solvers on the ALTIX. main entry through Paso_solve. this version compiles but   * linear solvers on the ALTIX. main entry through Paso_solve. this version compiles but
207   * has not been tested yet.   * has not been tested yet.
208   *   *

Legend:
Removed from v.4285  
changed lines
  Added in v.4286

  ViewVC Help
Powered by ViewVC 1.1.26