Actual source code: ex197.c
1: static char help[] = "Test MatMultHermitianTranspose() and MatMultHermitianTransposeAdd().\n\n";
3: #include <petscmat.h>
5: int main(int argc,char **args)
6: {
7: Mat A,B,C;
8: Vec x,y,ys;
9: PetscInt i,j;
10: PetscScalar v;
11: PetscBool flg;
13: PetscInitialize(&argc,&args,(char*)0,help);
14: MatCreate(PETSC_COMM_WORLD,&A);
15: MatSetSizes(A,PETSC_DECIDE,PETSC_DECIDE,2,2);
16: MatSetType(A,MATAIJ);
17: MatSetFromOptions(A);
18: MatSetUp(A);
20: i = 0; j = 0; v = 2.0;
21: MatSetValues(A,1,&i,1,&j,&v,INSERT_VALUES);
22: i = 0; j = 1; v = 3.0 + 4.0*PETSC_i;
23: MatSetValues(A,1,&i,1,&j,&v,INSERT_VALUES);
24: i = 1; j = 0; v = 5.0 + 6.0*PETSC_i;
25: MatSetValues(A,1,&i,1,&j,&v,INSERT_VALUES);
26: i = 1; j = 1; v = 7.0 + 8.0*PETSC_i;
27: MatSetValues(A,1,&i,1,&j,&v,INSERT_VALUES);
28: MatAssemblyBegin(A,MAT_FINAL_ASSEMBLY);
29: MatAssemblyEnd(A,MAT_FINAL_ASSEMBLY);
31: /* Create vectors */
32: VecCreate(PETSC_COMM_WORLD,&y);
33: VecSetSizes(y,PETSC_DECIDE,2);
34: VecSetFromOptions(y);
35: VecDuplicate(y,&ys);
36: VecDuplicate(y,&x);
38: i = 0; v = 10.0 + 11.0*PETSC_i;
39: VecSetValues(x,1,&i,&v,INSERT_VALUES);
40: i = 1; v = 100.0 + 120.0*PETSC_i;
41: VecSetValues(x,1,&i,&v,INSERT_VALUES);
42: VecAssemblyBegin(x);
43: VecAssemblyEnd(x);
45: MatMultHermitianTranspose(A,x,y);
46: VecView(y,PETSC_VIEWER_STDOUT_WORLD);
47: MatMultHermitianTransposeAdd(A,x,y,ys);
48: VecView(ys,PETSC_VIEWER_STDOUT_WORLD);
50: MatHermitianTranspose(A,MAT_INITIAL_MATRIX,&B);
51: MatCreateHermitianTranspose(A,&C);
52: MatMultHermitianTransposeEqual(B,C,4,&flg);
54: MatMultHermitianTransposeAddEqual(B,C,4,&flg);
56: MatDestroy(&C);
57: MatDestroy(&B);
59: MatDestroy(&A);
61: VecDestroy(&x);
62: VecDestroy(&y);
63: VecDestroy(&ys);
64: PetscFinalize();
65: return 0;
66: }
68: /*TEST
70: build:
71: requires: complex
72: test:
74: test:
75: suffix: 2
76: nsize: 2
78: TEST*/