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

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

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

revision 2965 by jfenwick, Wed Mar 3 01:53:41 2010 UTC revision 2966 by jfenwick, Wed Mar 3 03:19:39 2010 UTC
# Line 66  void Paso_SystemMatrix_CalcBorderMIS(Pas Line 66  void Paso_SystemMatrix_CalcBorderMIS(Pas
66          for (j=0;j<bordercount;++j) {          for (j=0;j<bordercount;++j) {
67          index_t bnode=border[j];          index_t bnode=border[j];
68          if (ISAVAILABLE(weights[bnode])) {          if (ISAVAILABLE(weights[bnode])) {
69                Paso_Pattern* p=A->mainBlock->pattern;
70              weights[bnode]=MISIN;              weights[bnode]=MISIN;
71              /* Now walk the neighbours and mark them unavailable */              /* Now walk the neighbours and mark them unavailable */
             Paso_Pattern* p=A->mainBlock->pattern;  
72              for (k=p->ptr[bnode];k<p->ptr[bnode+1];++k) {   /* Walk along the row */              for (k=p->ptr[bnode];k<p->ptr[bnode+1];++k) {   /* Walk along the row */
73              if (p->index[k]!=bnode) {   /* ignore diagonal link to self */              if (p->index[k]!=bnode) {   /* ignore diagonal link to self */
74                  weights[p->index[k]]=MISOUT;    /* node can't be in the set */                  weights[p->index[k]]=MISOUT;    /* node can't be in the set */

Legend:
Removed from v.2965  
changed lines
  Added in v.2966

  ViewVC Help
Powered by ViewVC 1.1.26