[Getdp] Surface probability density

Mon Jul 11 13:31:21 CEST 2011

Hi Mike,

in the case you haven't solved your problem yet, here is an example for integrating your probability function in certain regions (-> see attached files).
It is based on your test.pro file with some modifications. The result fits perfectly the analytical solution (calculated with Mathcad -> See Probability.JPG) I hope it is helpful.

Please let me know if you have further questions.

Kind regards,
Michael

-----Original Message-----
From: getdp-bounces at ace20.montefiore.ulg.ac.be [mailto:getdp-bounces at ace20.montefiore.ulg.ac.be] On Behalf Of getdp-request at ace20.montefiore.ulg.ac.be
Sent: Friday, June 24, 2011 12:00 PM
To: getdp at geuz.org
Subject: getdp Digest, Vol 99, Issue 6

Send getdp mailing list submissions to
	getdp at geuz.org

To subscribe or unsubscribe via the World Wide Web, visit
	http://www.geuz.org/mailman/listinfo/getdp
or, via email, send a message with subject or body 'help' to
	getdp-request at geuz.org

You can reach the person managing the list at
	getdp-owner at geuz.org

When replying, please edit your Subject line so it is more specific than "Re: Contents of getdp digest..."

Today's Topics:

   1. Surface probability density .pro (Michael Grant)

----------------------------------------------------------------------

Message: 1
Date: Thu, 23 Jun 2011 19:30:41 +0200
From: "Michael Grant" <m.grant at ee.wits.ac.za>
To: getdp at geuz.org
Subject: [Getdp] Surface probability density .pro
Message-ID: <20110623172504.M7641 at ee.wits.ac.za>
Content-Type: text/plain;	charset="iso-8859-1"

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>

------------------------------

_______________________________________________
getdp mailing list
getdp at geuz.org
http://www.geuz.org/mailman/listinfo/getdp

End of getdp Digest, Vol 99, Issue 6
************************************
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Probability.JPG
Type: image/jpeg
Size: 15557 bytes
Desc: Probability.JPG
URL: <http://www.geuz.org/pipermail/getdp/attachments/20110711/e67e1997/attachment.JPG>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Probability2.geo
Type: application/octet-stream
Size: 615 bytes
Desc: Probability2.geo
URL: <http://www.geuz.org/pipermail/getdp/attachments/20110711/e67e1997/attachment.geo>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Probability2.pro
Type: application/octet-stream
Size: 2494 bytes
Desc: Probability2.pro
URL: <http://www.geuz.org/pipermail/getdp/attachments/20110711/e67e1997/attachment.pro>