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: mmax_f.h,v 2.0 2003/02/22 15:23:25 judd Exp $ */
  #include"VU_mprintm_f.include"
  static void mmax_f(void){
     printf("\n*******\nTEST mmax_f\n\n");
     { vsip_scalar_f data1[] = {-1, 2, 0, -5, -6,  3.4, -3.4, 5.6, -.3};
       vsip_scalar_f data2[] = {-2, 1, 0, -3, -7, -3.4, -3.5, 5.6, -.2};
       vsip_scalar_f ans[]   = {-1, 2, 0, -3, -6,  3.4, -3.4, 5.6, -.2};
       vsip_block_f  *block1 = vsip_blockbind_f(data1,9,VSIP_MEM_NONE);
       vsip_block_f  *block2 = vsip_blockbind_f(data2,9,VSIP_MEM_NONE);
       vsip_block_f  *ans_block = vsip_blockbind_f(ans,9,VSIP_MEM_NONE);
       vsip_mview_f  *a = vsip_mbind_f(block1,0,3,3,1,3);
       vsip_mview_f  *b = vsip_mbind_f(block2,0,3,3,1,3);
       vsip_mview_f  *ansm = vsip_mbind_f(ans_block,0,3,3,1,3);
       vsip_mview_f  *c = vsip_mcreate_f(3,3,VSIP_ROW,VSIP_MEM_NONE);
       vsip_mview_f  *chk = vsip_mcreate_f(3,3,VSIP_COL,VSIP_MEM_NONE);
       vsip_blockadmit_f(block1,VSIP_TRUE);
       vsip_blockadmit_f(block2,VSIP_TRUE);
       vsip_blockadmit_f(ans_block,VSIP_TRUE);
       printf("test out of place, compact, user \"a\", vsipl \"b\"\n");
       vsip_mmax_f(a,b,c);
       printf("vsip_mmax_f(a,b,c)\n matrix a\n");VU_mprintm_f("8.6",a);
       printf("matrix b\n");VU_mprintm_f("8.6",b);
       printf("matrix c\n");VU_mprintm_f("8.6",c);
       printf("expected answer to 4 decimal digits\n");VU_mprintm_f("8.4",ansm);
  
       vsip_msub_f(c,ansm,chk); vsip_mmag_f(chk,chk); vsip_mclip_f(chk,.0001,.0001,0,1,chk);
       if(vsip_msumval_f(chk) > .5)
           printf("error\n");
       else
           printf("correct\n");
  
       printf("check mmax_f b,c inplace\n");
       vsip_mmax_f(a,b,b);
       vsip_msub_f(b,ansm,chk); vsip_mmag_f(chk,chk); vsip_mclip_f(chk,.0001,.0001,0,1,chk);
       if(vsip_msumval_f(chk) > .5)
           printf("error\n");
       else
           printf("correct\n");
  
       vsip_malldestroy_f(a);
       vsip_malldestroy_f(b);
       vsip_malldestroy_f(c);
       vsip_malldestroy_f(chk);
       vsip_malldestroy_f(ansm);
     }
     return;
  }
  