Home My Page Projects Code Snippets Project Openings diderot
Summary Activity Tracker Tasks SCM

SCM Repository

[diderot] Diff of /branches/fem/examples/mip_d3s/init.c
ViewVC logotype

Diff of /branches/fem/examples/mip_d3s/init.c

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

revision 3233, Tue Aug 11 15:07:12 2015 UTC revision 3234, Fri Aug 14 13:35:44 2015 UTC
# Line 13  Line 13 
13  #include "teem/air.h"  #include "teem/air.h"
14    
15    
16    void Nrrd2(Nrrd *nData, float  stepSize)
17    {
18        fprintf(stderr,"\n ******** Saving Results *********\n\n");
19        FILE *fp;
20        fp=fopen("b.txt","w");
21    
22        float *p = (float *)(nData->data);
23        unsigned int iPoints = nData->axis[0].size;
24        unsigned int jPoints = nData->axis[1].size;
25        unsigned int kPoints = nData->axis[2].size;
26        unsigned int lPoints = nData->axis[3].size;
27    
28        // for (uint32_t i = 0;  i < iPoints;  i++) {
29        for (uint32_t j = 0;  j < jPoints;  j++) {
30            for (uint32_t k = 0;  k < kPoints;  k++) {
31    
32                    //float pp0=p[(k+kPoints*(j))];
33                    //float pp1=p[(k+kPoints*(j+1))];
34               // float pp0=p[(k+kPoints+(j))];
35                //float pp1=p[(k+kPoints+(j+1))];
36                int fi=iPoints*(j+jPoints*(k));
37                int f2=1+2*(j+jPoints*(k));
38                float pp0=p[fi];
39                float pp1=p[f2];
40    
41              //  float pp0=p[(j+jPoints*(k))];
42               // float pp1=p[(j+jPoints*(k+1))];
43                //float pp0=p[(j+jPoints+k)];
44                //float pp1=p[(j+jPoints+(k+1))]
45               // float pp0=p[(k+j+jPoints)];
46                //float pp1=p[(k+1+j+jPoints)];
47              //  float pp0=p[(k*j+jPoints)];
48                //float pp1=p[((k+1)*j+jPoints)];
49    
50                    fprintf(fp,"\n F[%f][%f]=%f %f @ %d %d",(j*stepSize),(k*stepSize),pp0,pp1,fi,f2);
51    
52            }
53        }
54        // }
55        fprintf(stderr,"\n  %d-x-%d-x-%d -x-%d-x-%d  x%d*****************\n",iPoints,jPoints,kPoints,nData->axis[3].size,nData->axis[4].size,nData->axis[5].size);
56    
57    }
58    
59    
60    
61    
62    void Nrrd3(Nrrd *nData, float  stepSize)
63    {
64        fprintf(stderr,"\n ******** Saving Results *********\n\n");
65        FILE *fp;
66        fp=fopen("a.txt","w");
67    
68        float *p = (float *)(nData->data);
69        unsigned int iPoints = nData->axis[0].size;
70        unsigned int jPoints = nData->axis[1].size;
71        unsigned int kPoints = nData->axis[2].size;
72        unsigned int lPoints = nData->axis[3].size;
73    
74       // for (uint32_t i = 0;  i < iPoints;  i++) {
75            for (uint32_t j = 0;  j < jPoints;  j++) {
76                for (uint32_t k = 0;  k < kPoints;  k++) {
77                        for (uint32_t l = 0;  l < lPoints;  l++) {
78    
79                            int f1=iPoints*(j+jPoints*(k+kPoints*(l)));
80    
81                            float pp0=p[f1];
82                            float pp1=p[f1+1];
83                            float pp2=p[f1+2];
84                            float pp3=p[f1+3];
85    
86    
87                    fprintf(fp,"\n F[%f][%f][%f]=%f %f %f %f",(j*stepSize),(k*stepSize),(l*stepSize),pp0,pp1,pp2,pp3);
88                        }
89                }
90            }
91       // }
92        fprintf(stderr,"\n  %d-x-%d-x-%d -x-%d-x-%d  x%d*****************\n",iPoints,jPoints,kPoints,nData->axis[3].size,nData->axis[4].size,nData->axis[5].size);
93    
94    }
95    
96    
97    
98    
99  void WriteOutput (Nrrd *nData)  void WriteOutput (Nrrd *nData)
100  {  {
# Line 89  Line 171 
171      mip_d3s_InVarSet_imgResV(wrld,imgResV);      mip_d3s_InVarSet_imgResV(wrld,imgResV);
172      mip_d3s_InVarSet_rayStep(wrld,rayStep);      mip_d3s_InVarSet_rayStep(wrld,rayStep);
173    
     mip_d3s_InVarSet_camEye(wrld,camEye);  
     mip_d3s_InVarSet_camAt(wrld,camAt);  
     mip_d3s_InVarSet_camUp(wrld,camUp);  
     mip_d3s_InVarSet_camNear(wrld,camNear);  
     mip_d3s_InVarSet_camFar(wrld,camFar);  
     mip_d3s_InVarSet_camFOV(wrld,camFOV);  
174    
175      // nrrd for getting computational state      // nrrd for getting computational state
176      Nrrd *nData = nrrdNew();      Nrrd *nData = nrrdNew();
# Line 114  Line 190 
190          fprintf(stderr, "Error getting nrrd data: %s", mip_d3s_GetErrors(wrld));          fprintf(stderr, "Error getting nrrd data: %s", mip_d3s_GetErrors(wrld));
191          exit(1);          exit(1);
192      }      }
193        Nrrd3(nData,  1.0);
194      WriteOutput(nData);      WriteOutput(nData);
195      mip_d3s_Shutdown (wrld);      mip_d3s_Shutdown (wrld);
196    

Legend:
Removed from v.3233  
changed lines
  Added in v.3234

root@smlnj-gforge.cs.uchicago.edu
ViewVC Help
Powered by ViewVC 1.0.0