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

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

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

revision 414 by jgs, Wed Nov 9 02:02:19 2005 UTC revision 415 by gross, Wed Jan 4 05:37:33 2006 UTC
# Line 19  Line 19 
19  void Paso_SystemMatrix_copy(Paso_SystemMatrix* in,double* array) {  void Paso_SystemMatrix_copy(Paso_SystemMatrix* in,double* array) {
20    dim_t i,j;    dim_t i,j;
21    index_t iptr;    index_t iptr;
22      index_t index_offset=(in->type & MATRIX_FORMAT_OFFSET1 ? 1:0);
23    #pragma omp parallel for private(i,iptr,j) schedule(static)    #pragma omp parallel for private(i,iptr,j) schedule(static)
24    for (i=0;i< in->pattern->n_ptr;++i) {    for (i=0;i< in->pattern->n_ptr;++i) {
25       for (iptr=(in->pattern->ptr[i])-PTR_OFFSET;iptr<(in->pattern->ptr[i+1])-PTR_OFFSET; ++iptr) {       for (iptr=(in->pattern->ptr[i])-index_offset;iptr<(in->pattern->ptr[i+1])-index_offset; ++iptr) {
26           for (j=0;j<in->block_size;j++) array[iptr*(in->block_size)+j]=in->val[iptr*(in->block_size)+j];           for (j=0;j<in->block_size;j++) array[iptr*(in->block_size)+j]=in->val[iptr*(in->block_size)+j];
27       }       }
28    }    }

Legend:
Removed from v.414  
changed lines
  Added in v.415

  ViewVC Help
Powered by ViewVC 1.1.26