mrcPixelDataSet(&tmpu[i], x, y, z, 0, mrcPixelRePart);
mrcPixelDataSet(&tmpv[i], x, y, z, 0, mrcPixelRePart);
mrcPixelDataSet(&tmpw[i], x, y, z, 0, mrcPixelRePart);
- mrcPixelDataSet(&outFlow, x, y, z, 0, mrcPixelRePart);
+ //mrcPixelDataSet(&outFlow, x, y, z, 0, mrcPixelRePart);
} else {
count = 0;
mrcPixelDataSet(&(tmpu[i]), x, y, z, u, mrcPixelRePart);
mrcPixelDataSet(&(tmpv[i]), x, y, z, v, mrcPixelRePart);
mrcPixelDataSet(&(tmpw[i]), x, y, z, w, mrcPixelRePart);
- mrcPixelDataSet(&outFlow, x, y, z, size, mrcPixelRePart);
+ //mrcPixelDataSet(&outFlow, x, y, z, size, mrcPixelRePart);
u = 2 * u;
v = 2 * v;
mrcPixelDataSet(&(tmpv[i-1]), 2*x, 2*y, 2*z, v, mrcPixelRePart);
mrcPixelDataSet(&(tmpw[i-1]), 2*x, 2*y, 2*z, w, mrcPixelRePart);
-
+ }
+ }
+ }
+ }
//////////////////////////////////計算した速度を2倍して高い解像度のわたす/////////////////////////////////////////////////////////
- for(z = 1; z < In1Shrink[i-1].HeaderN.z; z++){
- for(y = 1; y < In1Shrink[i-1].HeaderN.y; y++){
- for(x = 1; x < In1Shrink[i-1].HeaderN.x; x++){
+ for(z = 1; z < In1Shrink[i-1].HeaderN.z; z++){
+ for(y = 1; y < In1Shrink[i-1].HeaderN.y; y++){
+ for(x = 1; x < In1Shrink[i-1].HeaderN.x; x++){
- if((int)x%2 != 0 || (int)y%2 != 0 || (int)z%2 != 0){
+ if((int)x%2 != 0 || (int)y%2 != 0 || (int)z%2 != 0){
- mrcPixelDataGet(&(tmpu[i-1]), x - 1, y , z , &data1, mrcPixelRePart, mrcPixelHowNearest);
- mrcPixelDataGet(&(tmpu[i-1]), x + 1, y , z , &data2, mrcPixelRePart, mrcPixelHowNearest);
- mrcPixelDataGet(&(tmpu[i-1]), x , y - 1, z , &data3, mrcPixelRePart, mrcPixelHowNearest);
- mrcPixelDataGet(&(tmpu[i-1]), x , y + 1, z , &data4, mrcPixelRePart, mrcPixelHowNearest);
- mrcPixelDataGet(&(tmpu[i-1]), x , y , z - 1, &data5, mrcPixelRePart, mrcPixelHowNearest);
- mrcPixelDataGet(&(tmpu[i-1]), x , y , z + 1, &data6, mrcPixelRePart, mrcPixelHowNearest);
+ mrcPixelDataGet(&(tmpu[i-1]), x - 1, y , z , &data1, mrcPixelRePart, mrcPixelHowNearest);
+ mrcPixelDataGet(&(tmpu[i-1]), x + 1, y , z , &data2, mrcPixelRePart, mrcPixelHowNearest);
+ mrcPixelDataGet(&(tmpu[i-1]), x , y - 1, z , &data3, mrcPixelRePart, mrcPixelHowNearest);
+ mrcPixelDataGet(&(tmpu[i-1]), x , y + 1, z , &data4, mrcPixelRePart, mrcPixelHowNearest);
+ mrcPixelDataGet(&(tmpu[i-1]), x , y , z - 1, &data5, mrcPixelRePart, mrcPixelHowNearest);
+ mrcPixelDataGet(&(tmpu[i-1]), x , y , z + 1, &data6, mrcPixelRePart, mrcPixelHowNearest);
- size1 = (data1 + data2 + data3 + data4 + data5 + data6) / 6 ;
-
- mrcPixelDataGet(&(tmpv[i-1]), x - 1, y , z , &data1, mrcPixelRePart, mrcPixelHowNearest);
- mrcPixelDataGet(&(tmpv[i-1]), x + 1, y , z , &data2, mrcPixelRePart, mrcPixelHowNearest);
- mrcPixelDataGet(&(tmpv[i-1]), x , y - 1, z , &data3, mrcPixelRePart, mrcPixelHowNearest);
- mrcPixelDataGet(&(tmpv[i-1]), x , y + 1, z , &data4, mrcPixelRePart, mrcPixelHowNearest);
- mrcPixelDataGet(&(tmpv[i-1]), x , y , z - 1, &data5, mrcPixelRePart, mrcPixelHowNearest);
- mrcPixelDataGet(&(tmpv[i-1]), x , y , z + 1, &data6, mrcPixelRePart, mrcPixelHowNearest);
+ size1 = (data1 + data2 + data3 + data4 + data5 + data6) / 6 ;
+
+ mrcPixelDataGet(&(tmpv[i-1]), x - 1, y , z , &data1, mrcPixelRePart, mrcPixelHowNearest);
+ mrcPixelDataGet(&(tmpv[i-1]), x + 1, y , z , &data2, mrcPixelRePart, mrcPixelHowNearest);
+ mrcPixelDataGet(&(tmpv[i-1]), x , y - 1, z , &data3, mrcPixelRePart, mrcPixelHowNearest);
+ mrcPixelDataGet(&(tmpv[i-1]), x , y + 1, z , &data4, mrcPixelRePart, mrcPixelHowNearest);
+ mrcPixelDataGet(&(tmpv[i-1]), x , y , z - 1, &data5, mrcPixelRePart, mrcPixelHowNearest);
+ mrcPixelDataGet(&(tmpv[i-1]), x , y , z + 1, &data6, mrcPixelRePart, mrcPixelHowNearest);
- size2 = (data1 + data2 + data3 + data4 + data5 + data6) / 6 ;
-
- mrcPixelDataGet(&(tmpw[i-1]), x - 1, y , z , &data1, mrcPixelRePart, mrcPixelHowNearest);
- mrcPixelDataGet(&(tmpw[i-1]), x + 1, y , z , &data2, mrcPixelRePart, mrcPixelHowNearest);
- mrcPixelDataGet(&(tmpw[i-1]), x , y - 1, z , &data3, mrcPixelRePart, mrcPixelHowNearest);
- mrcPixelDataGet(&(tmpw[i-1]), x , y + 1, z , &data4, mrcPixelRePart, mrcPixelHowNearest);
- mrcPixelDataGet(&(tmpw[i-1]), x , y , z - 1, &data5, mrcPixelRePart, mrcPixelHowNearest);
- mrcPixelDataGet(&(tmpw[i-1]), x , y , z + 1, &data6, mrcPixelRePart, mrcPixelHowNearest);
+ size2 = (data1 + data2 + data3 + data4 + data5 + data6) / 6 ;
+
+ mrcPixelDataGet(&(tmpw[i-1]), x - 1, y , z , &data1, mrcPixelRePart, mrcPixelHowNearest);
+ mrcPixelDataGet(&(tmpw[i-1]), x + 1, y , z , &data2, mrcPixelRePart, mrcPixelHowNearest);
+ mrcPixelDataGet(&(tmpw[i-1]), x , y - 1, z , &data3, mrcPixelRePart, mrcPixelHowNearest);
+ mrcPixelDataGet(&(tmpw[i-1]), x , y + 1, z , &data4, mrcPixelRePart, mrcPixelHowNearest);
+ mrcPixelDataGet(&(tmpw[i-1]), x , y , z - 1, &data5, mrcPixelRePart, mrcPixelHowNearest);
+ mrcPixelDataGet(&(tmpw[i-1]), x , y , z + 1, &data6, mrcPixelRePart, mrcPixelHowNearest);
- size3 = (data1 + data2 + data3 + data4 + data5 + data6) / 6 ;
-
- mrcPixelDataSet(&(tmpu[i-1]), x, y, z, size1, mrcPixelRePart);
- mrcPixelDataSet(&(tmpv[i-1]), x, y, z, size2, mrcPixelRePart);
- mrcPixelDataSet(&(tmpw[i-1]), x, y, z, size3, mrcPixelRePart);
+ size3 = (data1 + data2 + data3 + data4 + data5 + data6) / 6 ;
- }
- }
- }
- }
+ mrcPixelDataSet(&(tmpu[i-1]), x, y, z, size1, mrcPixelRePart);
+ mrcPixelDataSet(&(tmpv[i-1]), x, y, z, size2, mrcPixelRePart);
+ mrcPixelDataSet(&(tmpw[i-1]), x, y, z, size3, mrcPixelRePart);
-////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
}
}
}
}
+
+////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
}
+ // }
+ // }
+// }
gsl_matrix_free(m);
gsl_matrix_free(b);