printdebuginfo.c

二维SIMPLEC算法

#include "global_var.h" 
void printdebuginfo()
{   int i,j;
	if (!lprintdebuginfo) return;
	//***************************************for debug***********************
      printf("\n");
	  printf("nf=%d iter=%d\n, before solve",nf,iter);
	  for (j=1;j<=m1;j++)
	  {
		  for (i=1;i<=l1;i++) printf("T=%6.3f ",t(i,j));
		  printf("\n");
	  }
         printf("\n");

     for (j=1;j<=m1;j++)
	  {
		  for (i=1;i<=l1;i++) printf("v=%6.3f ",v(i,j));
		  printf("\n");
	  }
         printf("\n");

    for (j=1;j<=m1;j++)
	  {
		  for (i=1;i<=l1;i++) printf("u=%6.3f ",u(i,j));
		  printf("\n");
	  }
         printf("\n");

	  for (j=2;j<=m2;j++)
	   {
		  for (i=2;i<=l2;i++) 	  printf("ap=%6.3f ",ap[i][j]); 
		  printf("\n");
	   }
	     printf("\n");
     for (j=2;j<=m2;j++)
	  {
		  for (i=2;i<=l2;i++) 	  printf("aim=%6.3f ",aim[i][j]); 
		  printf("\n");
	   }
	     printf("\n");
	   for (j=2;j<=m2;j++)
	    {
		  for (i=2;i<=l2;i++) 	  printf("aip=%6.3f ",aip[i][j]); 
		  printf("\n");
	   }
	     printf("\n");
	  for (j=2;j<=m2;j++)
	   {
	 	  for (i=2;i<=l2;i++) 	  printf("ajm=%6.3f ",ajm[i][j]); 
		  printf("\n");
	   }
	     printf("\n");
	  for (j=2;j<=m2;j++)
	   {
	 	  for (i=2;i<=l2;i++) 	  printf("ajp=%6.3f ",ajp[i][j]); 
		  printf("\n");
	   }
	     printf("\n");
	  for (j=2;j<=m2;j++)
	   {
		  for (i=2;i<=l2;i++) 	  printf("con=%6.3f ",con[i][j]); 
		  printf("\n");
	   }
	     printf("\n");
   //***************************************************************************
}