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