[Getdp] Surface probability density .pro

Thu Jun 23 19:30:41 CEST 2011

Hi Folks,

I'm managing to illicit some non-standard behavior from getdp, particularly from PETSC:

the call is getdp test.pro -msh test.msh -pre testr and the error is 

[0]PETSC ERROR: Caught signal number 11 SEGV: Segmentation Violation, probably memory access 
out of range

I'm using this on MacOsX (10.6.7) and the test problem in the demo folder of both the 32 bit and 64 
bit stable OR snapshot build works. (i.e. not getdp's fault).  The problem is clearly between the 
keyboard and the chair.

I'm trying to calculate a probability at each point, and then integrate it within a region.  The 
probability density p(x, y) = exp(-(x+y)).

My problem file is quite happily digested by getdp, just cant compute the presolver info.

Any ideas?

Kind regards,

Mike Grant.

--------------test.pro-------------------
Group {
	Space = Region[101]; Building = Region[102];	
	Problem = Region[ {Space, Building} ];

}

Function {
	prob[] = Exp(-1*($X + $Y);
}

FunctionSpace {
{ Name gspace; Type Scalar;
	BasisFunction {
{Name sn; NameOfCoef vn; Function BF_Node;
Support Problem; Entity NodesOf[All]; }
	}
	GlobalQuantity {
	{Name gquantity; Type AssociatedWith; NameOfCoef vn;}
	}
}
}

Formulation {
{ Name testf; Type FemEquation;
	Quantity {
	{ Name gquantity; Type Global; NameOfSpace gspace; }
	}
	Equation {
		GlobalTerm {  [ prob[], {gquantity} ]  ; In Problem; }
	}
}
}

Resolution {
{ Name testr;
	System {
	{Name gsystem; NameOfFormulation testf;}
	}
	Operation {
	Generate [gsystem]; Solve[gsystem]; SaveSolution[gsystem];
	}
}
}

PostProcessing {
{ Name test; NameOfFormulation testf;
	Quantity {
	{Name gquantity; Value {Local{ [ {gquantity} ] ; In Problem; }}}
	}
}
}

PostOperation {
{Name post_out; NameOfPostProcessing test;
	Operation {
		Print [gquantity, OnElementsOf Problem, File "map.pos"];
	}
}
} 

<html><p><font face = "verdana" size = "0.8" color = "navy">This communication is intended for the addressee only. It is confidential. If you have received this communication in error, please notify us immediately and destroy the original message. You may not copy or disseminate this communication without the permission of the University. Only authorized signatories are competent to enter into agreements on behalf of the University and recipients are thus advised that the content of this message may not be legally binding on the University and may contain the personal views and opinions of the author, which are not necessarily the views and opinions of The University of the Witwatersrand, Johannesburg. All agreements between the University and outsiders are subject to South African Law unless the University agrees in writing to the contrary.</font></p></html>