matlab - Octave error when trying to integrate: "integrand F must return a single, real-valued vector of the same size as the input" -

-

i've confirmed issue in attempting compute glonball_2, don't see why it's producing error. function funct_2 should returning real numbers. "same size input" part, don't think i'm interpreting error message correctly. funct_2 takes 2 arguments , returns 1 value, funct_1, caused no error, that's not it. explain this? code:

function glonball = glonball(alpha) epsilon = 1; rmin = 0; rmax = 1; thmin = 0; thmax = 2*pi; funct_1 = @(r,th) (2).*(abs(alpha).^2-1).^(2)./abs(1-conj(alpha).*r.*exp(i.*th)).^4; glonball_1 = dblquad(funct_1, rmin, rmax, thmin, thmax); funct_2 = @(r,th) (abs(alpha.-r.*exp(i.*th))/abs(1-conj(alpha).*r.*exp(i.*th))).^2; glonball_2 = dblquad(funct_2, rmin, rmax, thmin, thmax); glonball = .5*glonball_1 + (1/((4).*epsilon^2))*glonball_2; end

here's error receive when attempting run input argument between 0 , 1 (the area of interest):

error: quadcc: integrand f must return single, real-valued vector of same size input error: called from: error:   c:\software\octave-3.6.4\share\octave\3.6.4\m\general\dblquad.m @ line  72, column 10 error:    @ line -1, column -1 error: quadcc: integrand f must return single, real-valued vector error:   c:\software\octave-3.6.4\share\octave\3.6.4\m\general\dblquad.m @ line  65, column 5 error:   c:\software\octave-3.6.4\glonball.m @ line 10, column 13

the .- unnecessary mention in comment above, doesn't cause errors in octave. give syntax error in matlab. similarly, / instead of ./ what's giving primary error.

after changing both of function ran fine in both octave 4.0.0 , matlab 2015a.

glonball(1) ans =     1.5708

Comments

Post a Comment

Popular posts from this blog

python - Healpy: From Data to Healpix map -

-
Image
i have data grid rows represent theta (0, pi) , columns represent phi (0, 2*pi) , f(theta,phi) density of dark matter @ location. wanted calculate power spectrum , have decided use healpy. what can not understand how format data healpy use. if provide code (in python obvious reasons) or point me tutorial, great! have tried hand @ doing following code: #grid dimensions nrows*ncols (subject change) theta = np.linspace(0, np.pi, num=grid.shape[0])[:, none] phi = np.linspace(0, 2*np.pi, num=grid.shape[1]) nside = 512 print "pixel area: %.2f square degrees" % hp.nside2pixarea(nside, degrees=true) pix = hp.ang2pix(nside, theta, phi) healpix_map = np.zeros(hp.nside2npix(nside), dtype=np.double) healpix_map[pix] = grid but, when try execute code power spectrum. specifically, : cl = hp.anafast(healpix_map[pix], lmax=1024) i error: typeerror: bad number of pixels if point me tutorial or edit code great. more specifications: data in 2d np array , can change numrows/...
Read more

c - Bitwise operation with (signed) enum value -

-
i using enumerator values flags: typedef enum { = 0x00, b = 0x01u, // u has no influence, expected c = 0x02u, // u has no influence, expected ... } enum_name; volatile unsigned char* reg = someaddress; *reg |= b; according misra-c:2004 bitwise operations shall not done signed type. unfortunately, compiler iar use signed int (or short or char) underlying type of enums, , option can find relates size, not signedness ("--enum-is-int"). according iar c/c++ development guide arm , pages 169 , 211, can define type of enum s if enable iar language extensions ( -e command-line option, or project > options > c/c++ compiler > language > allow iar extensions in ide). in particular, should define "sentinel" value, make sure compiler chooses correct type. prefers signed types, , uses smallest possible integer type, sentinel should largest positive integer corresponding unsigned integer type can describe. example, typedef enum ...
Read more

xslt - Unnest parent nodes by child node -

-
in xml there items 0-n attributes , item should copied each attribute new item 1 attribute. given xml this: <?xml version="1.0" encoding="utf-8" ?> <items> <item> <name>a</name> <attributes> <attribute> <key>attribute1</key> <value>1</value> </attribute> </attributes> </item> <item> <name>b</name> </item> <item> <name>c</name> <attributes> <attribute> <key>attribute1</key> <value>5</value> </attribute> <attribute> <key>attribute2</key> <value>2</value> </attribute> <attribute> <key>attri...
Read more