706 |
|
|
707 |
/* set nodes: */ |
/* set nodes: */ |
708 |
/* INTERNAL/BOUNDARY NODES */ |
/* INTERNAL/BOUNDARY NODES */ |
709 |
|
k=0; |
710 |
#pragma omp parallel for private(i0,i1,i2,k) |
#pragma omp parallel for private(i0,i1,i2,k) |
711 |
for (k=0,i2=0;i2<N2;i2++) { |
for (i2=0;i2<N2;i2++) { |
712 |
for (i1=0;i1<N1;i1++) { |
for (i1=0;i1<N1;i1++) { |
713 |
for (i0=0;i0<numNodesLocal-domLeft*periodic[0];i0++,k++) { |
for (i0=0;i0<numNodesLocal-domLeft*periodic[0];i0++,k++) { |
714 |
out->Nodes->Coordinates[INDEX2(0,k,3)]=DBLE(i0+firstNode)/DBLE(N0-1)*Length[0]; |
out->Nodes->Coordinates[INDEX2(0,k,3)]=DBLE(i0+firstNode)/DBLE(N0-1)*Length[0]; |