Date: 2024-04-26T10:23:23+00:00

Actual source code: test1.c
slepc-3.21.1 2024-04-26
Report Typos and Errors

      1: /*
      2:    - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
      3:    SLEPc - Scalable Library for Eigenvalue Problem Computations
      4:    Copyright (c) 2002-, Universitat Politecnica de Valencia, Spain

      6:    This file is part of SLEPc.
      7:    SLEPc is distributed under a 2-clause BSD license (see LICENSE).
      8:    - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
      9: */

     11: static char help[] = "Test RG interface functions.\n\n";

     13: #include <slepcrg.h>

     15: #define NPOINTS 10
     16: #define NVERTEX 7

     18: int main(int argc,char **argv)
     19: {
     20:   RG             rg;
     21:   PetscInt       i,inside,nv;
     22:   PetscBool      triv;
     23:   PetscReal      re,im,a,b,c,d;
     24:   PetscScalar    ar,ai,cr[NPOINTS],ci[NPOINTS],vr[NVERTEX],vi[NVERTEX],*pr,*pi;

     26:   PetscFunctionBeginUser;
     27:   PetscCall(SlepcInitialize(&argc,&argv,(char*)0,help));
     28:   PetscCall(RGCreate(PETSC_COMM_WORLD,&rg));

     30:   /* ellipse */
     31:   PetscCall(RGSetType(rg,RGELLIPSE));
     32:   PetscCall(RGIsTrivial(rg,&triv));
     33:   PetscCheck(triv,PETSC_COMM_WORLD,PETSC_ERR_PLIB,"Region should be trivial before setting parameters");
     34:   PetscCall(RGEllipseSetParameters(rg,1.1,2,0.1));
     35:   PetscCall(RGSetFromOptions(rg));
     36:   PetscCall(RGIsTrivial(rg,&triv));
     37:   PetscCheck(!triv,PETSC_COMM_WORLD,PETSC_ERR_PLIB,"Region should be non-trivial after setting parameters");
     38:   PetscCall(RGView(rg,NULL));
     39:   PetscCall(RGViewFromOptions(rg,NULL,"-rg_ellipse_view"));
     40:   re = 0.1; im = 0.3;
     41: #if defined(PETSC_USE_COMPLEX)
     42:   ar = PetscCMPLX(re,im);
     43: #else
     44:   ar = re; ai = im;
     45: #endif
     46:   PetscCall(RGCheckInside(rg,1,&ar,&ai,&inside));
     47:   PetscCall(PetscPrintf(PETSC_COMM_WORLD,"Point (%g,%g) is %s the region\n",(double)re,(double)im,(inside>=0)?"inside":"outside"));

     49:   PetscCall(RGComputeBoundingBox(rg,&a,&b,&c,&d));
     50:   PetscCall(PetscPrintf(PETSC_COMM_WORLD,"The bounding box is [%g,%g]x[%g,%g]\n",(double)a,(double)b,(double)c,(double)d));

     52:   PetscCall(PetscPrintf(PETSC_COMM_WORLD,"Contour points: "));
     53:   PetscCall(RGComputeContour(rg,NPOINTS,cr,ci));
     54:   for (i=0;i<NPOINTS;i++) {
     55: #if defined(PETSC_USE_COMPLEX)
     56:     re = PetscRealPart(cr[i]);
     57:     im = PetscImaginaryPart(cr[i]);
     58: #else
     59:     re = cr[i];
     60:     im = ci[i];
     61: #endif
     62:     PetscCall(PetscPrintf(PETSC_COMM_WORLD,"(%.3g,%.3g) ",(double)re,(double)im));
     63:   }
     64:   PetscCall(PetscPrintf(PETSC_COMM_WORLD,"\n"));

     66:   /* interval */
     67:   PetscCall(RGSetType(rg,RGINTERVAL));
     68:   PetscCall(RGIsTrivial(rg,&triv));
     69:   PetscCheck(triv,PETSC_COMM_WORLD,PETSC_ERR_PLIB,"Region should be trivial before setting parameters");
     70:   PetscCall(RGIntervalSetEndpoints(rg,-1,1,-0.1,0.1));
     71:   PetscCall(RGSetFromOptions(rg));
     72:   PetscCall(RGIsTrivial(rg,&triv));
     73:   PetscCheck(!triv,PETSC_COMM_WORLD,PETSC_ERR_PLIB,"Region should be non-trivial after setting parameters");
     74:   PetscCall(RGView(rg,NULL));
     75:   PetscCall(RGViewFromOptions(rg,NULL,"-rg_interval_view"));
     76:   re = 0.2; im = 0;
     77: #if defined(PETSC_USE_COMPLEX)
     78:   ar = PetscCMPLX(re,im);
     79: #else
     80:   ar = re; ai = im;
     81: #endif
     82:   PetscCall(RGCheckInside(rg,1,&ar,&ai,&inside));
     83:   PetscCall(PetscPrintf(PETSC_COMM_WORLD,"Point (%g,%g) is %s the region\n",(double)re,(double)im,(inside>=0)?"inside":"outside"));

     85:   PetscCall(RGComputeBoundingBox(rg,&a,&b,&c,&d));
     86:   PetscCall(PetscPrintf(PETSC_COMM_WORLD,"The bounding box is [%g,%g]x[%g,%g]\n",(double)a,(double)b,(double)c,(double)d));

     88:   PetscCall(PetscPrintf(PETSC_COMM_WORLD,"Contour points: "));
     89:   PetscCall(RGComputeContour(rg,NPOINTS,cr,ci));
     90:   for (i=0;i<NPOINTS;i++) {
     91: #if defined(PETSC_USE_COMPLEX)
     92:     re = PetscRealPart(cr[i]);
     93:     im = PetscImaginaryPart(cr[i]);
     94: #else
     95:     re = cr[i];
     96:     im = ci[i];
     97: #endif
     98:     PetscCall(PetscPrintf(PETSC_COMM_WORLD,"(%.3g,%.3g) ",(double)re,(double)im));
     99:   }
    100:   PetscCall(PetscPrintf(PETSC_COMM_WORLD,"\n"));

    102:   /* polygon */
    103: #if defined(PETSC_USE_COMPLEX)
    104:   vr[0] = PetscCMPLX(0.0,2.0);
    105:   vr[1] = PetscCMPLX(1.0,4.0);
    106:   vr[2] = PetscCMPLX(2.0,5.0);
    107:   vr[3] = PetscCMPLX(4.0,3.0);
    108:   vr[4] = PetscCMPLX(5.0,4.0);
    109:   vr[5] = PetscCMPLX(6.0,1.0);
    110:   vr[6] = PetscCMPLX(2.0,0.0);
    111: #else
    112:   vr[0] = 0.0; vi[0] = 1.0;
    113:   vr[1] = 0.0; vi[1] = -1.0;
    114:   vr[2] = 0.6; vi[2] = -0.8;
    115:   vr[3] = 1.0; vi[3] = -1.0;
    116:   vr[4] = 2.0; vi[4] = 0.0;
    117:   vr[5] = 1.0; vi[5] = 1.0;
    118:   vr[6] = 0.6; vi[6] = 0.8;
    119: #endif
    120:   PetscCall(RGSetType(rg,RGPOLYGON));
    121:   PetscCall(RGIsTrivial(rg,&triv));
    122:   PetscCheck(triv,PETSC_COMM_WORLD,PETSC_ERR_PLIB,"Region should be trivial before setting parameters");
    123:   PetscCall(RGPolygonSetVertices(rg,NVERTEX,vr,vi));
    124:   PetscCall(RGSetFromOptions(rg));
    125:   PetscCall(RGIsTrivial(rg,&triv));
    126:   PetscCheck(!triv,PETSC_COMM_WORLD,PETSC_ERR_PLIB,"Region should be non-trivial after setting parameters");
    127:   PetscCall(RGView(rg,NULL));
    128:   PetscCall(RGViewFromOptions(rg,NULL,"-rg_polygon_view"));
    129:   re = 5; im = 0.9;
    130: #if defined(PETSC_USE_COMPLEX)
    131:   ar = PetscCMPLX(re,im);
    132: #else
    133:   ar = re; ai = im;
    134: #endif
    135:   PetscCall(RGCheckInside(rg,1,&ar,&ai,&inside));
    136:   PetscCall(PetscPrintf(PETSC_COMM_WORLD,"Point (%g,%g) is %s the region\n",(double)re,(double)im,(inside>=0)?"inside":"outside"));

    138:   PetscCall(RGComputeBoundingBox(rg,&a,&b,&c,&d));
    139:   PetscCall(PetscPrintf(PETSC_COMM_WORLD,"The bounding box is [%g,%g]x[%g,%g]\n",(double)a,(double)b,(double)c,(double)d));

    141:   PetscCall(PetscPrintf(PETSC_COMM_WORLD,"Contour points: "));
    142:   PetscCall(RGComputeContour(rg,NPOINTS,cr,ci));
    143:   for (i=0;i<NPOINTS;i++) {
    144: #if defined(PETSC_USE_COMPLEX)
    145:     re = PetscRealPart(cr[i]);
    146:     im = PetscImaginaryPart(cr[i]);
    147: #else
    148:     re = cr[i];
    149:     im = ci[i];
    150: #endif
    151:     PetscCall(PetscPrintf(PETSC_COMM_WORLD,"(%.3g,%.3g) ",(double)re,(double)im));
    152:   }
    153:   PetscCall(PetscPrintf(PETSC_COMM_WORLD,"\n"));

    155:   /* check vertices */
    156:   PetscCall(RGPolygonGetVertices(rg,&nv,&pr,&pi));
    157:   PetscCheck(nv==NVERTEX,PETSC_COMM_WORLD,PETSC_ERR_PLIB,"Wrong number of vertices: %" PetscInt_FMT,nv);
    158:   for (i=0;i<nv;i++) {
    159: #if !defined(PETSC_USE_COMPLEX)
    160:     if (pr[i]!=vr[i] || pi[i]!=vi[i])
    161: #else
    162:     if (pr[i]!=vr[i])
    163: #endif
    164:        SETERRQ(PETSC_COMM_WORLD,PETSC_ERR_PLIB,"Vertex number %" PetscInt_FMT " does not match",i);
    165:   }

    167:   PetscCall(PetscFree(pr));
    168: #if !defined(PETSC_USE_COMPLEX)
    169:   PetscCall(PetscFree(pi));
    170: #endif
    171:   PetscCall(RGDestroy(&rg));
    172:   PetscCall(SlepcFinalize());
    173:   return 0;
    174: }

    176: /*TEST

    178:    test:
    179:       suffix: 1
    180:       requires: !complex

    182:    test:
    183:       suffix: 1_complex
    184:       requires: complex

    186:    test:
    187:       suffix: 2
    188:       requires: !complex
    189:       args: -rg_ellipse_view draw:tikz:ellipse.tikz -rg_interval_view draw:tikz:interval.tikz -rg_polygon_view draw:tikz:polygon.tikz
    190:       filter: cat - ellipse.tikz interval.tikz polygon.tikz
    191:       requires: !single

    193:    test:
    194:       suffix: 2_complex
    195:       requires: complex !single
    196:       args: -rg_ellipse_view draw:tikz:ellipse.tikz -rg_interval_view draw:tikz:interval.tikz -rg_polygon_view draw:tikz:polygon.tikz
    197:       filter: cat - ellipse.tikz interval.tikz polygon.tikz

    199: TEST*/