添加了两个图形界面

ybw2016v · Nov 29, 2018 · ee90fd1 · ee90fd1
1 parent 9696dd6
commit ee90fd1
Show file tree

Hide file tree

Showing 19 changed files with 207 additions and 9 deletions.
diff --git a/cat.stl b/cat.stl
diff --git a/dogui2.py b/dogui2.py
@@ -100,6 +100,7 @@ def _DogCal_fired(self):
         self.nummax=self.NBdog.max()
         self.DogOutShape=str(self.NBdog.shape)
         mlab.pipeline.volume(fileddog,vmin=0,vmax=self.nummax)
+        # mlab.contour3d(self.NBdog)
         mlab.colorbar()
         self.DogNpy()
         self.DogCwj()
@@ -125,13 +126,13 @@ def DogCwj(self):
     static int y$num[]=$ry1;
     static int z$num[]=$rz1;
     int longdog=$l1;
-    for (int i = 0; i  < longdog; i ++)
+    for (int icd = 0; icd  < longdog; icd ++)
     {
-        //p1[x$num[i]*xar+y$num[i]*yar+z$num[i]*zar]=
-        //vx[x$num[i]*xar+y$num[i]*yar+z$num[i]*zar]=
-        //vy[x$num[i]*xar+y$num[i]*yar+z$num[i]*zar]=
-        //vz[x$num[i]*xar+y$num[i]*yar+z$num[i]*zar]=
-        //z0[x$num[i]*xar+y$num[i]*yar+z$num[i]*zar]=
+        //p1[x$num[icd]*xar+y$num[icd]*yar+z$num[icd]*zar]=
+        //vx[x$num[icd]*xar+y$num[icd]*yar+z$num[icd]*zar]=
+        //vy[x$num[icd]*xar+y$num[icd]*yar+z$num[icd]*zar]=
+        //vz[x$num[icd]*xar+y$num[icd]*yar+z$num[icd]*zar]=
+        //z0[x$num[icd]*xar+y$num[icd]*yar+z$num[icd]*zar]=
     }
     '''
 
@@ -141,10 +142,14 @@ def DogCwj(self):
             # print(chDog)
             if self.nuls[chDog]==True:
                 lup=np.where(self.NBdog==chDog)
-                spx=str(list(lup[0]))
-                spy=str(list(lup[1]))
-                spz=str(list(lup[2]))
+                # llxx=len()
+                spx1=str(list(lup[1]))
+                spy1=str(list(lup[2]))
+                spz1=str(list(lup[0]))
                 spl=str(lup[0].size)
+                spx=self.DogStr(spx1)
+                spy=self.DogStr(spy1)
+                spz=self.DogStr(spz1)
                 apu=DogStr.replace('$num',str(chDog)).replace('$rx1',spx).replace('$ry1',spy).replace('$rz1',spz).replace('$l1',spl)
                 cumdog=cumdog+apu
                 print(chDog)
@@ -175,6 +180,12 @@ def _DogOutput_fired(self):
         with open(self.DogOutC, 'w+',encoding='utf-8') as f:
             f.write(self.DogC)
         pass
+
+    def DogStr(self,strindog):
+        luf=len(strindog)
+        cudog='{'+strindog[1:luf-1]+'}'
+        return cudog
+        pass
 
 
 

diff --git a/dogui3.py b/dogui3.py
@@ -0,0 +1,91 @@
+#!/usr/bin/env python3
+
+import numpy as npy
+import numpy as np
+from mayavi import mlab
+from mayavi.core.ui.mayavi_scene import MayaviScene
+from mayavi.tools.mlab_scene_model import MlabSceneModel
+from traits.api import *
+from traits.api import Enum, HasTraits, Int, Range, Str
+from traitsui.api import *
+from traitsui.api import Item, RangeEditor, View
+from traitsui.menu import *
+from tvtk.pyface.scene_editor import SceneEditor
+
+
+
+
+
+class DogPlotres(HasTraits):
+
+
+    scene = Instance(MlabSceneModel, ())
+    ImportDog =  Button("导入")
+    DogImpor=Directory('./res/')
+
+    DogName=Str('p0_0.npy')
+
+    MeenDog=Enum(['点图','等值面'])
+    Num=0
+
+    DogNext =  Button("下一个")
+
+    DogNext2 =  Button("上一个")
+
+    DogNum =Int(15)
+
+    DogNummax =Float(2)
+    DogNummin =Float(-2)
+
+    Plot =  Button("绘图")
+
+
+
+    dog=Item(name='scene',editor=SceneEditor(scene_class=MayaviScene),resizable=True,height=600,width=800)
+
+    view=View(HGroup(dog,VGroup(Item(name='DogImpor',label='导入目录'),Item(name='ImportDog',label='导入'),Item(name='DogName',label='文件名'),Item(name='MeenDog',label='类别'),Item(name='DogNum',label='等值面数量'),Item(name='DogNummax',label='等值面最大值'),Item(name='DogNummin',label='等值面最小值'),Item(name='DogNext',label='下一个'),Item(name='DogNext2',label='上一个')),'Plot',show_labels=False),resizable = True,title=u"TraitsUI",buttons=OKCancelButtons)
+
+    def _ImportDog_fired(self):
+        self.LsbDog=np.load(str(self.DogImpor)+str(self.DogName))
+
+        pass
+
+    def _Plot_fired(self, parameter_list):
+        self.ploot()
+
+        pass
+    def _DogNext_fired(self):
+        self.Num=self.Num+1
+        self.DogName='p0_'+str(self.Num)+'.npy'
+        self.LsbDog=np.load(str(self.DogImpor)+str(self.DogName))
+        self.ploot()
+
+        pass
+    def _DogNext2_fired(self):
+        self.Num=self.Num-1
+        self.DogName='p0_'+str(self.Num)+'.npy'
+        self.LsbDog=np.load(str(self.DogImpor)+str(self.DogName))
+        self.ploot()
+        pass
+    def ploot(self):
+        if self.MeenDog=='点图':
+            a,b,c=self.LsbDog.shape
+            mlab.clf()
+            mlab.contour3d(self.LsbDog[2:int(a-4),2:int(b-4),2:int(c-4)],contours=int(self.DogNum),vmax=float(self.DogNummax),extent=[0,a,0,b,0,c],vmin=float(self.DogNummin),transparent=True)
+            mlab.axes(xlabel='x', ylabel='y', zlabel='z')
+            mlab.colorbar()
+        else:
+            mlab.clf()
+            fileddog=mlab.pipeline.scalar_field(self.LsbDog)
+            mlab.pipeline.volume(fileddog)
+            mlab.axes(xlabel='x', ylabel='y', zlabel='z')
+            mlab.colorbar()
+        pass
+
+
+
+    pass
+
+
+app = DogPlotres(resizable=True)
+app.configure_traits()   
diff --git a/plotdog.py b/plotdog.py
@@ -0,0 +1,14 @@
+
+import numpy as np
+from mayavi import mlab
+
+ddd= np.load('p0_0.npy')
+# ddd[1,1:81,1:81]=6
+# ddd[:,1,:]=6
+fileddog=mlab.pipeline.scalar_field(ddd)
+# mlab.pipeline.volume(fileddog)
+# 
+print(ddd.shape)
+mlab.contour3d(ddd,contours=20,transparent=True)
+mlab.colorbar()
+mlab.show()
diff --git a/test/cat.stl b/test/cat.stl
diff --git a/cat2.stl → test/cat2.stl b/cat2.stl → test/cat2.stl
diff --git a/cat3.stl → test/cat3.stl b/cat3.stl → test/cat3.stl
diff --git a/cat4.stl → test/cat4.stl b/cat4.stl → test/cat4.stl
diff --git a/cat5.stl → test/cat5.stl b/cat5.stl → test/cat5.stl
diff --git a/test/cat6.stl b/test/cat6.stl
diff --git a/test/cat7.stl b/test/cat7.stl
diff --git a/test/cat8.stl b/test/cat8.stl
diff --git a/test/cat9.stl b/test/cat9.stl
diff --git a/test/dddog.ini b/test/dddog.ini
@@ -0,0 +1,19 @@
+[计算规模]
+长度 = 128
+宽度 = 128
+高度 = 128
+计算时间 = 5000
+[初始值]
+自动初值 = 是
+p0路径 = 
+vx路径 =
+vy路径 = 
+vz路径 =
+z0路径 = out4.npy
+m = 4
+[中间过程]
+添加方式 =./outc4.c
+
+[结果输出]
+输出频率 = 100
+输出路径 = ./res/
diff --git a/moban.c → test/moban.c b/moban.c → test/moban.c
diff --git a/outc.c → test/outc.c b/outc.c → test/outc.c
diff --git a/test/outc2.c b/test/outc2.c
@@ -0,0 +1,21 @@
+
+void usercode(float * p1,float * vx,float * vy,float * vz,float *z0,int xar, int yar,int zar,float m,int i)
+{
+
+
+    static int x2[]={0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2};
+    static int y2[]={10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 17, 17, 17, 17, 17, 17, 17, 17, 17, 17, 17, 17, 17, 17, 18, 18, 18, 18, 18, 18, 18, 18, 18, 18, 18, 18, 18, 18, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 10, 10, 10, 10, 10, 10, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 15, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 16, 17, 17, 17, 17, 17, 17, 17, 17, 17, 17, 17, 17, 17, 18, 18, 18, 18, 18, 18, 18, 18, 18, 18, 18, 18, 18, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 19, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 20, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 21, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 22, 10, 10, 10, 10, 10, 10, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 12, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 13, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 14, 15, 15, 15, 15, 15, 15, 15, 15, 15, 16, 16, 16, 16, 16, 16, 16, 16, 16, 17, 17, 17, 17, 17, 17, 17, 18, 18};
+    static int z2[]={17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 17, 18, 19, 20, 21, 22, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 17, 18, 19, 20, 21, 22, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 21, 22, 23, 24, 25, 26, 27, 28, 29, 21, 22, 23, 24, 25, 26, 27, 28, 29, 23, 24, 25, 26, 27, 28, 29, 25, 26};
+    int longdog=425;
+    for (int icd = 0; icd  < longdog; icd ++)
+    {
+       p1[x2[icd]*xar+y2[icd]*yar+z2[icd]*zar]=sin(0.1*i);
+        //vx[x2[icd]*xar+y2[icd]*yar+z2[icd]*zar]=
+        //vy[x2[icd]*xar+y2[icd]*yar+z2[icd]*zar]=
+        //vz[x2[icd]*xar+y2[icd]*yar+z2[icd]*zar]=
+        //z0[x2[icd]*xar+y2[icd]*yar+z2[icd]*zar]=
+    }
+
+
+}
+
diff --git a/test/outc3.c b/test/outc3.c
diff --git a/test/outc4.c b/test/outc4.c
@@ -0,0 +1,21 @@
+
+void usercode(float * p1,float * vx,float * vy,float * vz,float *z0,int xar, int yar,int zar,float m,int i)
+{
+
+
+    static int x3[]={31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 31, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 32, 33, 33, 33, 33, 33, 33, 33, 33, 33, 33, 33, 33, 33, 33, 33, 33, 33, 33, 33, 33, 33, 34, 34, 34, 34, 34, 34, 34, 34, 34, 34, 34, 34, 34, 34, 34, 34, 34, 34, 34, 34, 34, 35, 35, 35, 35, 35, 35, 35, 35, 35, 35, 35, 35, 35, 35, 35, 35, 35, 35, 35, 35, 35, 36, 36, 36, 36, 36, 36, 36, 36, 36, 36, 36, 36, 36, 36, 36, 36, 36, 36, 36, 36, 36, 37, 37, 37, 37, 37, 37, 37, 37, 37, 37, 37, 37, 37, 37, 37, 37, 37, 37, 37, 37, 37, 38, 38, 38, 38, 38, 38, 38, 38, 38, 38, 38, 38, 38, 38, 38, 38, 38, 38, 38, 38, 38, 39, 39, 39, 39, 39, 39, 39, 39, 39, 39, 39, 39, 39, 39, 39, 39, 39, 39, 39, 39, 39, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 40, 41, 41, 41, 41, 41, 41, 41, 41, 41, 41, 41, 41, 41, 41, 41, 41, 41, 41, 41, 41, 41, 42, 42, 42, 42, 42, 42, 42, 42, 42, 42, 42, 42, 42, 42, 42, 42, 42, 42, 42, 42, 42, 43, 43, 43, 43, 43, 43, 43, 43, 43, 43, 43, 43, 43, 43, 43, 43, 43, 43, 43, 43, 43, 44, 44, 44, 44, 44, 44, 44, 44, 44, 44, 44, 44, 44, 44, 44, 44, 44, 44, 44, 44, 44, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 46, 46, 46, 46, 46, 46, 46, 46, 46, 46, 46, 46, 46, 46, 46, 46, 46, 46, 46, 46, 46, 47, 47, 47, 47, 47, 47, 47, 47, 47, 47, 47, 47, 47, 47, 47, 47, 47, 47, 47, 47, 47, 48, 48, 48, 48, 48, 48, 48, 48, 48, 48, 48, 48, 48, 48, 48, 48, 48, 48, 48, 48, 48, 49, 49, 49, 49, 49, 49, 49, 49, 49, 49, 49, 49, 49, 49, 49, 49, 49, 49, 49, 49, 49, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 50, 51, 51, 51, 51, 51, 51, 51, 51, 51, 51, 51, 51, 51, 51, 51, 51, 51, 51, 51, 51, 51};
+    static int y3[]={27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47};
+    static int z3[]={7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, 7};
+    int longdog=441;
+    for (int icd = 0; icd  < longdog; icd ++)
+    {
+        p1[x3[icd]*xar+y3[icd]*yar+z3[icd]*zar]=sin(0.5*i);
+        //vx[x3[icd]*xar+y3[icd]*yar+z3[icd]*zar]=
+        //vy[x3[icd]*xar+y3[icd]*yar+z3[icd]*zar]=
+        //vz[x3[icd]*xar+y3[icd]*yar+z3[icd]*zar]=
+        //z0[x3[icd]*xar+y3[icd]*yar+z3[icd]*zar]=
+    }
+
+
+}
+