text: /* Created RJudd */
  /* SPAWARSYSCEN D857 */
  /**********************************************************************
  // For TASP VSIPL Documentation and Code neither the United States    /
  // Government, the United States Navy, nor any of their employees,    /
  // makes any warranty, express or implied, including the warranties   /
  // of merchantability and fitness for a particular purpose, or        /
  // assumes any legal liability or responsibility for the accuracy,    /
  // completeness, or usefulness of any information, apparatus,         /
  // product, or process disclosed, or represents that its use would    /
  // not infringe privately owned rights                                /
  **********************************************************************/
  /* $Id: mlge_d.h,v 2.0 2003/02/22 15:23:25 judd Exp $ */
  static void mlge_d(){
      printf("\n******\nTEST mlge_d\n");
      {
          vsip_scalar_d data1[]= {1,  2,  3,   4,  5,   6,   7,   8,   9};
          vsip_scalar_d data2[]= {1., 2., 3.3, 4., 5.0, 5.6, 7.7, 7.8, 9.0};
          vsip_scalar_d ans[] =  {1,  1,  0,   1,  1,   1,   0,   1,   1};
          vsip_mview_d *m1 = vsip_mbind_d(
                    vsip_blockbind_d(data1,9,VSIP_MEM_NONE),0,3,3,1,3);
          vsip_mview_d *m2 = vsip_mbind_d(
                    vsip_blockbind_d(data2,9,VSIP_MEM_NONE),0,3,3,1,3);
          vsip_mview_d *ma = vsip_mbind_d(
                    vsip_blockbind_d(ans,9,VSIP_MEM_NONE),0,3,3,1,3);
          vsip_mview_d *m3 = vsip_mcreate_d(3,3,VSIP_ROW,VSIP_MEM_NONE);
          vsip_mview_bl  *bl = vsip_mcreate_bl(3,3,VSIP_ROW,VSIP_MEM_NONE);
          vsip_mview_d *chk = vsip_mcreate_d(3,3,VSIP_ROW,VSIP_MEM_NONE);
          vsip_blockadmit_d(vsip_mgetblock_d(m1),VSIP_TRUE);
          vsip_blockadmit_d(vsip_mgetblock_d(m2),VSIP_TRUE);
          vsip_blockadmit_d(vsip_mgetblock_d(ma),VSIP_TRUE);
          printf("call vsip_mlge_d(a,b,c)\n");
          printf("a =\n");VU_mprintm_d("8.6",m1);
          printf("b =\n");VU_mprintm_d("8.6",m2);
          vsip_mlge_d(m1,m2,bl);
          vsip_mcopy_bl_d(bl,m3);
          printf("c =\n");VU_mprintm_d("8.6",m3);
          printf("\nright answer =\n"); VU_mprintm_d("6.4",ma);
          vsip_msub_d(ma,m3,chk); vsip_mmag_d(chk,chk); vsip_mclip_d(chk,0,.0001,0,1,chk);
          if(fabs(vsip_msumval_d(chk)) > .5)
              printf("error\n");
          else
              printf("correct\n");
         
          vsip_malldestroy_d(m1); vsip_malldestroy_d(m2);
          vsip_malldestroy_d(m3); vsip_malldestroy_d(ma);
          vsip_malldestroy_d(chk);
          vsip_malldestroy_bl(bl);
      }                                                       
        
      return;
  }
  