ext2_print.c

LastWave

/*..........................................................................*/
/*                                                                          */
/*      L a s t W a v e    P a c k a g e 'dwtrans2d' 2.1                    */
/*                                                                          */
/*      Copyright (C) 1998-2002  E.Bacry, J.Fraleu, J.Kalifa, E. Le Pennec, */
/*                         W.L. Hwang , S.Mallat, S.Zhong                   */
/*      emails : [e-mail]                             */
/*               [e-mail]                               */
/*               [e-mail]                               */
/*               [e-mail]                             */
/*               [e-mail]                               */
/*               [e-mail]                                   */
/*               [e-mail]                               */
/*                                                                          */
/*..........................................................................*/
/*                                                                          */
/*      This program is a free software, you can redistribute it and/or     */
/*      modify it under the terms of the GNU General Public License as      */
/*      published by the Free Software Foundation; either version 2 of the  */
/*      License, or (at your option) any later version                      */
/*                                                                          */
/*      This program is distributed in the hope that it will be useful,     */
/*      but WITHOUT ANY WARRANTY; without even the implied warranty of      */
/*      MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the       */
/*      GNU General Public License for more details.                        */
/*                                                                          */
/*      You should have received a copy of the GNU General Public License   */
/*      along with this program (in a file named COPYRIGHT);                */
/*      if not, write to the Free Software Foundation, Inc.,                */
/*      59 Temple Place, Suite 330, Boston, MA  02111-1307  USA             */
/*                                                                          */
/*..........................................................................*/



#include "lastwave.h"
#include "extrema2d.h"  

void W2_point_print(EXT2 ext)
{
  Printf("x = %3d, y = %3d, mag = %6.2f, arg = %5.2f\n",
	 ext->x, ext->y, ext->mag, ext->arg/M_PI);
}
          

/*********************************************************/
/** The functions deal with points print out. ************/
/*********************************************************/
 
void C_Point2Print(char **argv)
{
  WTRANS2 wtrans=NULL;
  int level;
  int row_min, row_max, col_min, col_max;
  int i, j;
  EXT2 ext;
  EXTLIS2 extlis;
  int ncol, nrow;
  EXT2 *values;
  int zap, nextline_flag = NO;
  char car;

  argv = ParseArgv(argv,tWTRANS2_,NULL,&wtrans,tINT,&level,-1);
  if (wtrans ==NULL)   wtrans= GetWtrans2Cur();
 
 
  if (!INRANGE(1, level, wtrans->extrep->noct))
   Errorf("Invalid level \n");

  extlis = wtrans->extrep->array[level];

  ncol = extlis->ncol;
  nrow = extlis->nrow;
  values = (EXT2 *)extlis->first;

  row_min = 0;
  row_max = nrow - 1;
  col_min = 0;
  col_max = ncol - 1;

  while (car=ParseOption(&argv)) {
    switch (car) {
    case 'r' :
      argv = ParseArgv(argv,tINT,&row_min,tINT,&row_max,-1);
      if (row_min > row_max) {zap= row_min; row_min= row_max; row_max= zap;}
      break;
    case 'c' :
      argv = ParseArgv(argv,tINT,&col_min,tINT,&col_max,-1); 
      if (col_min > col_max) {zap= col_min; col_min= col_max; col_max= zap;}
      break;
    default :
      ErrorOption(car);
    }
  }

  NoMoreArgs(argv);
  /* check whether valid row and col, if not put into boundary */
  if (INRANGE(0,row_min,nrow-1) == NO)
    row_min = 0;
  if (INRANGE(0,row_max,nrow-1) == NO)
    row_max = extlis->nrow - 1;
  if (INRANGE(0,col_min,ncol-1) == NO)
    col_min = 0;
  if (INRANGE(0,col_max,ncol-1) == NO)
    col_max = ncol - 1;

  for (i = row_min; i <= row_max; i++)
    for (j = col_min; j <= col_max; j++)
      if (ext = values[i * ncol + j]) {
     	W2_point_print(ext);
      }
  if (nextline_flag == NO) Printf("\n");  
}

/*********************************************/
/***    propagation chain print out        ***/
/*********************************************/

static void W2_prop_print(EXT2 ext)
{
  for (ext; ext; ext = ext->coarser)
     W2_point_print(ext);
  Printf("\n");
}

static void W2_propagation_print(EXTLIS2 extlis,int row_min,int row_max,int col_min,int col_max)
{
  int i, j, k, ncol;
  EXT2 *values;
  EXT2 ext;

  ncol = extlis->ncol;
  values = extlis->first;

  for (i = row_min; i <= row_max; i++)
    for (j = col_min; j <= col_max; j++) {
      k = i * ncol + j;
      if (ext = values[k]) W2_prop_print(ext);
    } 
}     




/*************************************************************/
/* print out the points information at the propagation curve */
/* defined within block -r .. and -c ...                     */
/*************************************************************/

void C_Point2ProPrint(char **argv)
{
  EXTLIS2 extlis;
  EXTREP2 extrep;
  WTRANS2 wtrans=NULL;
  int row_min, row_max, col_min, col_max;
  int zap;
  int nrow, ncol;
  int level;
  char car;

  argv = ParseArgv(argv,tWTRANS2_,NULL,&wtrans,tINT,&level,-1);
  if (wtrans ==NULL)   wtrans= GetWtrans2Cur();
 
 
  if (!INRANGE(1, level, wtrans->extrep->noct))
   Errorf("Invalid level \n");


  extrep = wtrans->extrep;
  extlis = extrep->array[level];
  nrow = extlis->nrow;
  ncol = extlis->ncol;

  row_min = 0;
  row_max = nrow - 1;
  col_min = 0;
  col_max = ncol - 1;

  while (car=ParseOption(&argv)) {
    switch (car) {
    case 'r' :
      argv = ParseArgv(argv,tINT,&row_min,tINT,&row_max,-1);
      if (row_min > row_max) {zap= row_min; row_min= row_max; row_max= zap;}
      break;
    case 'c' :
      argv = ParseArgv(argv,tINT,&col_min,tINT,&col_max,-1); 
      if (col_min > col_max) {zap= col_min; col_min= col_max; col_max= zap;}
      break;
    default :
      ErrorOption(car);
    }
 }   
  /* check whether valid row and col, if not put into boundary */
  if (INRANGE(0,row_min,nrow-1) == NO)
    row_min = 0;
  if (INRANGE(0,row_max,nrow-1) == NO)
    row_max = extlis->nrow - 1;
  if (INRANGE(0,col_min,ncol-1) == NO)
    col_min = 0;
  if (INRANGE(0,col_max,ncol-1) == NO)
    col_max = ncol - 1;

  W2_propagation_print(extlis, row_min, row_max, col_min, col_max);

}