Method LaplaceTiming

solve/nptests.cpp:301–364

  void NumProcTestJS::LaplaceTiming()
  {
    cout << "laplace timing" << endl;
    const MeshAccess & ma = pde.GetMeshAccess();

    clock_t starttime, endtime;
    starttime = clock();

    LaplaceIntegrator<2,FE_Trig1> lpi (new ConstantCoefficientFunction(1));
    ElasticityIntegrator<2> eli (new ConstantCoefficientFunction(1),
				 new ConstantCoefficientFunction(0.2));

    SourceIntegrator<2> source (new ConstantCoefficientFunction(1));


    //    BlockBilinearFormIntegrator<LaplaceIntegrator<2> > blpi(lpi, 3);
    //    BlockBilinearFormIntegrator<LaplaceIntegrator<2,FE_Trig1> > blpi(lpi, 3);

    BilinearFormIntegrator * bfi = &lpi;
    LinearFormIntegrator * lfi = &source;

    /*
    ElasticityIntegrator eli (2, new ConstantCoefficientFunction(1),
			      new ConstantCoefficientFunction(0.2));
    */
    FE_Trig1 trig;
    ElementTransformation eltrans;
    //    DenseMatrix elmat(trig.GetNDof());
    LocalHeap locheap (100000);

    int n = 1000000;

    ngbla::FlatMatrix<double> elmat;
    ngbla::FlatVector<double> elvec(3, locheap);


    void * heapp = locheap.GetPointer();

    for (int i = 1; i <= n; i++)
      {
	locheap.CleanUp(heapp);
	ma.GetElementTransformation (1, eltrans, locheap);
	bfi->CalcElementMatrix(trig, eltrans, elmat, locheap);
	//	elvec.Assign (lfi->CalcElementVector (trig, eltrans, locheap));
	// AssembleLaplace2 (eltrans, elmat, &locheap);
      }

    cout << "elmat = " << elmat << endl;
    cout << "elvec = " << elvec << endl;

    endtime = clock();
    cout << "sec per 1E6 els = " << 1e6 / n * double(endtime - starttime)/CLOCKS_PER_SEC << endl;

    /*
      Laplace, trig1: 169000
      Laplace, tet1:  141000

      Elast, trig1:   103000