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

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

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

revision 2784 by artak, Thu Nov 26 05:09:14 2009 UTC revision 2802 by artak, Thu Dec 3 01:51:55 2009 UTC
# Line 58  Paso_SparseMatrix* Paso_SparseMatrix_get Line 58  Paso_SparseMatrix* Paso_SparseMatrix_get
58      for (i=0;i<n;++i) {      for (i=0;i<n;++i) {
59        if (mis_marker[i]) {        if (mis_marker[i]) {
60            for (iptr=W->pattern->ptr[k];iptr<W->pattern->ptr[k+1]; ++iptr) {            for (iptr=W->pattern->ptr[k];iptr<W->pattern->ptr[k+1]; ++iptr) {
61            j=W->pattern->index[iptr];              j=W->pattern->index[iptr];
62            Paso_IndexList_insertIndex(&(index_list[i]),j);              Paso_IndexList_insertIndex(&(index_list[i]),j);
63          }            }
64          k++;            k++;
65        }        }
66        else {        else {
67            Paso_IndexList_insertIndex(&(index_list[i]),i-k);            Paso_IndexList_insertIndex(&(index_list[i]),i-k);
# Line 73  Paso_SparseMatrix* Paso_SparseMatrix_get Line 73  Paso_SparseMatrix* Paso_SparseMatrix_get
73            
74      k=0;      k=0;
75    
76      for (i=0;i<out->numRows;++i) {      for (i=0;i<n;++i) {
77        if (mis_marker[i]) {        if (mis_marker[i]) {
78          wptr=W->pattern->ptr[k];          wptr=W->pattern->ptr[k];
79          for (iptr=out->pattern->ptr[i];iptr<out->pattern->ptr[i+1]; ++iptr) {          for (iptr=out->pattern->ptr[i];iptr<out->pattern->ptr[i+1]; ++iptr) {

Legend:
Removed from v.2784  
changed lines
  Added in v.2802

  ViewVC Help
Powered by ViewVC 1.1.26