C++ passing array pointer address in function
<p>The code below uses an octahedron as a starting 3D shape, ( I found it on the web) and after some tweaking I can't still make it work.</p>

<pre><code>#include &lt;stdlib.h&gt;
#include &lt;stdio.h&gt;
#include &lt;string.h&gt;
#include &lt;stddef.h&gt; /* must include for the offsetof macro */
#include &lt;GL/glew.h&gt;
#include &lt;GL/glfw.h&gt;
#include &lt;glm/glm.hpp&gt;
#include &lt;glm/gtc/matrix_transform.hpp&gt;
#include &lt;glm/gtc/type_ptr.hpp&gt;
#include &lt;stdlib.h&gt;
#include &lt;vector&gt;
#include &lt;iostream&gt;

using namespace std; 

struct XYZ {
    GLdouble x;
    GLdouble y;
    GLdouble z;
};

struct FACET3 {
    XYZ p1;
    XYZ p2;
    XYZ p3;
};


void Normalise(XYZ *p_input)
{
    double magnitude = 0;
    magnitude = sqrt((p_input-&gt;x * p_input-&gt;x )+ (p_input-&gt;y * p_input-&gt;y) + (p_input-&gt;z * p_input-&gt;z));
    p_input-&gt;x = p_input-&gt;x / magnitude;
    p_input-&gt;y = p_input-&gt;y / magnitude;
    p_input-&gt;z = p_input-&gt;z / magnitude;

}

int CreateNSphere(FACET3 *f,int iterations)
{

   int i,it;
   double a;
   XYZ p[6] = {0,0,1,  0,0,-1,  -1,-1,0,  1,-1,0,  1,1,0, -1,1,0};
   XYZ pa,pb,pc;
   int nt = 0,ntold;

   /* Create the level 0 object */
   a = 1 / sqrt(2.0);
   for (i=0;i&lt;6;i++) {
      p[i].x *= a;
      p[i].y *= a;
   }
   f[0].p1 = p[0]; f[0].p2 = p[3]; f[0].p3 = p[4];
   f[1].p1 = p[0]; f[1].p2 = p[4]; f[1].p3 = p[5];
   f[2].p1 = p[0]; f[2].p2 = p[5]; f[2].p3 = p[2];
   f[3].p1 = p[0]; f[3].p2 = p[2]; f[3].p3 = p[3];
   f[4].p1 = p[1]; f[4].p2 = p[4]; f[4].p3 = p[3];
   f[5].p1 = p[1]; f[5].p2 = p[5]; f[5].p3 = p[4];
   f[6].p1 = p[1]; f[6].p2 = p[2]; f[6].p3 = p[5];
   f[7].p1 = p[1]; f[7].p2 = p[3]; f[7].p3 = p[2];
   nt = 8;

   if (iterations &lt; 1)
      return(nt);

   /* Bisect each edge and move to the surface of a unit sphere */
   for (it=0;it&lt;iterations;it++) {
      ntold = nt;
      for (i=0;i&lt;ntold;i++) {
         pa.x = (f[i].p1.x + f[i].p2.x) / 2;
         pa.y = (f[i].p1.y + f[i].p2.y) / 2;
         pa.z = (f[i].p1.z + f[i].p2.z) / 2;
         pb.x = (f[i].p2.x + f[i].p3.x) / 2;
         pb.y = (f[i].p2.y + f[i].p3.y) / 2;
         pb.z = (f[i].p2.z + f[i].p3.z) / 2;
         pc.x = (f[i].p3.x + f[i].p1.x) / 2;
         pc.y = (f[i].p3.y + f[i].p1.y) / 2;
         pc.z = (f[i].p3.z + f[i].p1.z) / 2;
         Normalise(&amp;pa);
         Normalise(&amp;pb);
         Normalise(&amp;pc);
         f[nt].p1 = f[i].p1; f[nt].p2 = pa; f[nt].p3 = pc; nt++;
         f[nt].p1 = pa; f[nt].p2 = f[i].p2; f[nt].p3 = pb; nt++;
         f[nt].p1 = pb; f[nt].p2 = f[i].p3; f[nt].p3 = pc; nt++;
         f[i].p1 = pa;
         f[i].p2 = pb;
         f[i].p3 = pc;
      }
   }

   return(nt);
}


int main()
{
    FACET3 facet[8];
    int facets = CreateNSphere(facet, 2);
    printf(" Result: %d", facets);
} 
</code></pre>

<p>When I try to compile it, it returns:</p>

<p>"Segmentation Error"</p>

<p>gdb debugger says:</p>

<pre><code>Program received signal SIGSEGV, Segmentation fault.
0x00000000004011df in CreateNSphere(FACET3*, int) ()
</code></pre>

<p>any help?</p>

<p>after compiling with "-g" option "gbd" returned this:</p>

<pre><code>(gdb) run
Starting program: /home/alexander/Graphics/Code-Handouts/tut_3/a.out 

Program received signal SIGSEGV, Segmentation fault.
0x00000000004011df in CreateNSphere (f=0x7fffffffdeb0, iterations=2) at sphere_model.cpp:86
86           f[nt].p1 = pb; f[nt].p2 = f[i].p3; f[nt].p3 = pc; nt++;
(gdb) where
#0  0x00000000004011df in CreateNSphere (f=0x7fffffffdeb0, iterations=2) at sphere_model.cpp:86
#1  0x0000000000401303 in main () at sphere_model.cpp:100
(gdb) 
</code></pre>
