I NEED TO CODE???

                    /* Put all coefficients to zero */
                    for (A = 0; A < 4; A++)
                    {
                        for (B = 0; B < 4; B++)
                        {
                            for (qq = 0; qq < 7; qq++)
                            {
                                for (ii = 0; ii < 3; ii++)
                                {
                                    for (ji = 0; ji < 3; ji++)
                                    {
                                        for (ki = 0; ki < 3; ki++)
                                        {
                                            coefg[A][B][ii][ji][ki][qq] = 0.0;
                                        }
                                    }
                                }
                            }
                        }
                    }

Since arrays with more than one dimension of do not get stored in one contigous region of memory, memset wouldn't work because it only works on contigous regions of memory. No doubt it would segfault as memset ran off the space allocated to the pointer coefg.

tl;dr:
>>4 fails.