1 2 static char help[] = "Makes a simple histogram.\n"; 3 4 #include <petscsys.h> 5 #include <petscdraw.h> 6 7 int main(int argc,char **argv) 8 { 9 PetscDraw draw; 10 PetscDrawHG hist; 11 PetscDrawAxis axis; 12 PetscErrorCode ierr; 13 int n = 20,i,x = 0,y = 0,width = 400,height = 300,bins = 8; 14 PetscInt w = 400,h = 300,nn = 20,b = 8,c = PETSC_DRAW_GREEN; 15 int color = PETSC_DRAW_GREEN; 16 const char *xlabel,*ylabel,*toplabel; 17 PetscReal xd; 18 PetscBool flg; 19 20 xlabel = "X-axis Label"; toplabel = "Top Label"; ylabel = "Y-axis Label"; 21 22 ierr = PetscInitialize(&argc,&argv,NULL,help);if (ierr) return ierr; 23 ierr = PetscOptionsGetInt(NULL,NULL,"-width",&w,NULL);CHKERRQ(ierr); 24 ierr = PetscOptionsGetInt(NULL,NULL,"-height",&h,NULL);CHKERRQ(ierr); 25 ierr = PetscOptionsGetInt(NULL,NULL,"-n",&nn,NULL);CHKERRQ(ierr); 26 ierr = PetscOptionsGetInt(NULL,NULL,"-bins",&b,NULL);CHKERRQ(ierr); 27 ierr = PetscOptionsGetInt(NULL,NULL,"-color",&c,NULL);CHKERRQ(ierr); 28 ierr = PetscOptionsHasName(NULL,NULL,"-nolabels",&flg);CHKERRQ(ierr); 29 width = (int) w; height = (int)h; n = (int)nn; bins = (int) b; color = (int) c; 30 if (flg) { xlabel = NULL; ylabel = NULL; toplabel = NULL; } 31 32 ierr = PetscDrawCreate(PETSC_COMM_WORLD,0,"Title",x,y,width,height,&draw);CHKERRQ(ierr); 33 ierr = PetscDrawSetFromOptions(draw);CHKERRQ(ierr); 34 ierr = PetscDrawHGCreate(draw,bins,&hist);CHKERRQ(ierr); 35 ierr = PetscDrawHGSetColor(hist,color);CHKERRQ(ierr); 36 ierr = PetscDrawHGGetAxis(hist,&axis);CHKERRQ(ierr); 37 ierr = PetscDrawAxisSetColors(axis,PETSC_DRAW_BLACK,PETSC_DRAW_RED,PETSC_DRAW_BLUE);CHKERRQ(ierr); 38 ierr = PetscDrawAxisSetLabels(axis,toplabel,xlabel,ylabel);CHKERRQ(ierr); 39 /*ierr = PetscDrawHGSetFromOptions(hist);CHKERRQ(ierr);*/ 40 41 for (i=0; i<n; i++) { 42 xd = (PetscReal)(i - 5); 43 ierr = PetscDrawHGAddValue(hist,xd*xd);CHKERRQ(ierr); 44 } 45 ierr = PetscDrawHGDraw(hist);CHKERRQ(ierr); 46 ierr = PetscDrawHGSave(hist);CHKERRQ(ierr); 47 48 ierr = PetscDrawHGDestroy(&hist);CHKERRQ(ierr); 49 ierr = PetscDrawDestroy(&draw);CHKERRQ(ierr); 50 ierr = PetscFinalize(); 51 return ierr; 52 } 53 54 55 56 /*TEST 57 58 build: 59 requires: x 60 61 test: 62 output_file: output/ex1_1.out 63 64 TEST*/ 65