Fixed Laplace damping option for TE-problem

Author: daniel-koehn

par/visu/create_STF.ipynb

@@ -44,9 +44,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {
-    "collapsed": true
-   },
+   "metadata": {},
    "outputs": [],
    "source": [
     "# TD parameters\n",
@@ -126,9 +124,7 @@
   {
    "cell_type": "code",
    "execution_count": null,
-   "metadata": {
-    "collapsed": true
-   },
+   "metadata": {},
    "outputs": [],
    "source": [
     "plt.rc('font',**{'family':'sans-serif','sans-serif':['Helvetica']})\n",

src/TE/RHS_source_TE.c

@@ -19,7 +19,7 @@ void RHS_source_TE(struct waveAC *waveAC, float ** srcpos, int ishot){
 	complex float omega;
 
 	amp = vector(1,2);
-	omega = 2.0 * M_PI * ((*waveAC).freq - S * I);
+	omega = 2.0 * M_PI * ((*waveAC).freq + S * I);
 	mu0 = 4.0 * M_PI * 1e-7;
 
 	/* calculate discrete source position */

src/TE/RHS_source_TE_hess.c

@@ -17,7 +17,7 @@ void RHS_source_TE_hess(struct waveAC *waveAC, int ** recpos, int trace){
         float ampr, ampi, mu0;
 	complex float omega;
 
-	omega = 2.0 * M_PI * ((*waveAC).freq - S * I);
+	omega = 2.0 * M_PI * ((*waveAC).freq + S * I);
 
         /* initialize source vector */
         for (i=0;i<NXNY;i++){

src/TE/ass_grad_TE.c

@@ -19,8 +19,8 @@ void ass_grad_TE(struct fwiTE *fwiTE, struct waveAC *waveAC, struct matTE *matTE
 	complex float wien, Omega, Omega2;
 
 	wien = 	(*waveAC).stfr + (*waveAC).stfi * I;	
-	Omega = 2.0 * M_PI * (*waveAC).freq - (I * S);
-	Omega2 = cpowf(((2.0*M_PI*(*waveAC).freq) - (I * S)),2.0);
+	Omega = 2.0 * M_PI * ((*waveAC).freq + (I * S));
+	Omega2 = cpowf(2.0*M_PI*((*waveAC).freq + (I * S)),2.0);
 
 	/* printf("abs(wien) = %e + i %e \n",creal(cabsf(wien)),cimag(cabsf(wien))); */
 

src/TE/hessian_TE.c

@@ -54,8 +54,8 @@ void hessian_TE(struct fwiTE *fwiTE, struct waveAC *waveAC, struct PML_AC *PML_A
 		(*waveAC).freq = (*waveAC).stage_freq[nfreq];
 
 		/* set squared complex angular frequency*/
-		Omega = 2.0 * M_PI * (*waveAC).freq - (I * S);		
-                Omega2 = cpowf(((2.0*M_PI*(*waveAC).freq) - (I * S)),2.0);
+		Omega = 2.0 * M_PI * ((*waveAC).freq + (I * S));		
+                Omega2 = cpowf(2.0*M_PI*((*waveAC).freq + (I * S)),2.0);
 
 		/* define PML damping profiles */
 		pml_pro(PML_AC,waveAC);

src/TE/hessian_shin_TE.c

@@ -54,8 +54,8 @@ void hessian_shin_TE(struct fwiTE *fwiTE, struct waveAC *waveAC, struct PML_AC *
 		(*waveAC).freq = (*waveAC).stage_freq[nfreq];
 
 		/* set squared complex angular frequency*/
-		Omega = 2.0 * M_PI * (*waveAC).freq - (I * S);
-		Omega2 = cpowf(((2.0*M_PI*(*waveAC).freq) - (I * S)),2.0);
+		Omega = 2.0 * M_PI * ((*waveAC).freq + (I * S));
+		Omega2 = cpowf((2.0*M_PI*((*waveAC).freq + (I * S))),2.0);
 
 		/* define PML damping profiles */
 		pml_pro(PML_AC,waveAC);

src/TE/init_A_TE_9p_pml.c

@@ -51,8 +51,8 @@ void init_A_TE_9p_pml(struct PML_AC *PML_AC, struct matTE *matTE, struct waveAC
         count=0;    /* count non-zero elements */
 
 	/* set squared complex angular frequency*/
-	Omega = 2.0 * M_PI * ((*waveAC).freq - S * I);
-	Omega2 = cpowf(2.0 * M_PI * ((*waveAC).freq - S * I),2.0);
+	Omega = 2.0 * M_PI * ((*waveAC).freq + S * I);
+	Omega2 = cpowf(2.0 * M_PI * ((*waveAC).freq + S * I),2.0);
 
 	/* define magnetic permeability mu0 */
 	mu0 = 4.0 * M_PI * 1e-7;

src/cp_grad_frame.c

@@ -207,8 +207,8 @@ void cp_grad_frame(float ** A){
 		    }
 	     }
 
-	    sprintf(modfile,"grad_damp_PML.bin");
-	    writemod(modfile,absorb_coeff,3); 
+	    /*sprintf(modfile,"grad_damp_PML.bin");
+	    writemod(modfile,absorb_coeff,3); */
 
 	    free_vector(coeff,1,NPML);	
 	    free_matrix(absorb_coeff,1,NY,1,NX);