# James Tursa

Interested in external applications, mex & engine applications with Fortran, C, C++, etc, particularly in the areas of speed improvement and memory management efficiencies.

**Statistics**

RANK**17**

of 250.280

REPUTATION**12.602**

CONTRIBUTIONS

**17** Questions

**4.038** Answers

ANSWER ACCEPTANCE **
41.18%
**

VOTES RECEIVED**2.271**

RANK**164** of 17.386

REPUTATION**7.217**

AVERAGE RATING**4.60**

CONTRIBUTIONS**32** Files

DOWNLOADS **125**

ALL TIME DOWNLOADS**67595**

CONTRIBUTIONS**0 Posts**

CONTRIBUTIONS**0** Public Channels

AVERAGE RATING

CONTRIBUTIONS**0** Highlights

AVERAGE NO. OF LIKES

**Content Feed**

how to read and write half precision arrays via mexFunction

Very unfortunately, MATLAB has implemented the half precision data type as an opaque classdef type instead of a simple numeric t...

1 día ago | 0

how to write mex function when input is a matlab function

This can be done, but it will not be efficient. The basic problem is that C/C++ doesn't understand anything about MATLAB functio...

1 día ago | 0

Is it possible to process the sparse matrix faster with vectorization instead of for loop?

You need to re-evaluate how you are doing things. In the first place, you have the Euclidean calculation wrong. It should be thi...

8 días ago | 0

How to exchange data between two C++ MEX files

Your current scheme is not foolproof. How does the 1st mex function know when it is safe to unlock and clear memory? How does th...

9 días ago | 0

| accepted

How to solve Lc=y without backslash operator?

First, I am assuming there is a typo and the system you are solving is Ly = c, not Lc=y. Second, you made a good start by writi...

11 días ago | 0

I am using euler's method to solve a differential equation, but when I run the code it doesn't plot.

Take a look at these lines: t0=0; %start time t1= 500; %finish time dt = 100000; Your stepsize is much larger than the total...

13 días ago | 1

| accepted

Accessing struct using C library always NULL

Typically, it is best just to post your code so we can see exactly what you are doing, instead of posting a description of your ...

13 días ago | 1

| accepted

Errors while trying to setup equation for root finding.

Did you mean multiply by the "a"? x_pdo = z_pdo/(1 + a*(k_pdo - 1)); x_water = z_water/(1 + a*(k_water - 1)); x_glycerol = z_...

14 días ago | 0

How do I access a field from a function's input?

Assuming you are passing in a character string for freq, use this syntax: output = myStruct.(freq);

16 días ago | 0

| accepted

How can I obtain the T and Y for R Runge Kutta method?

The first version has an error. This line: k3=h*feval(f,T(j),Y(j)); should be this instead: k3=h*feval(f,T(j)+h/2,Y(j)+k2/2);...

21 días ago | 0

How to solve array indices error?

x1 = 0*ft; : distance1=sqrt((x1-x2)^2+(y1-y2)^2+(z1-z2)^2); : sp = distance1(x1,y1,z1,dx,dyv(i),dz); % incident di...

alrededor de 1 mes ago | 0

Solving ODE using Euler Method and 4th Order Runge Kutta Method

You pretty much have the Euler scheme worked out, so I will help you with a vector formulation. Take this code: for i=1:n ...

alrededor de 1 mes ago | 0

| accepted

Problem finding sum of array using vectorization

The first element of x is 1, so that element produces 1/(1-1) = 1/0 = inf in the second line. You need to rewrite that second li...

alrededor de 1 mes ago | 0

Multiplication of large matrix with its transpose

The fastest way is to simply write it as A * A', because MATLAB will see that the operands are the same and call a special multi...

alrededor de 1 mes ago | 0

Solving ODE in MATLAB using Runge-Kutta method of order 4

Why aren't you using this to update z: z(i+1) = z(i) + (1/6)*(l1+(2*l2)+(2*l3)+l4);

alrededor de 1 mes ago | 0

Why is the string type not implemented as standard type?

All of the standard full numeric types as well as char and logical are implemented as simple rectangular data arrays. The string...

alrededor de 1 mes ago | 1

how to print randomly selected column?

Shouldn't that be size(data,2)? Also, generally you should be using string comparison functions for the tests, not the == opera...

alrededor de 2 meses ago | 0

NaN shows up in matrix when solving my ODEs, how can I fix this?

"All of my matrices I am solving for have NaN after the third time line. How do I solved this issue?" Set up a breakpoint somew...

alrededor de 2 meses ago | 1

| accepted

Mean of multidimensional array

Doesn't mean(mdar,3) do what you want?

alrededor de 2 meses ago | 0

| accepted

Function Error: Array indices must be positive integers or logical values.

When you have parentheses appear on the left hand side of an assignment like this: v(ro,ri) = (sqrt(m/(ri)))*((1/sqrt(R))-((sqr...

alrededor de 2 meses ago | 1

| accepted

MXCREATEDOUBLEMATRIX memory overflow error

Rather later for an Answer, but here goes anyway: You should never use literal integers for API calls, because you can't be sur...

alrededor de 2 meses ago | 0

Convert unsigned integer to decimal

Not sure what the real question is here. Does one of these do what you want? int1 = uint32(hex2dec('9d900000')); int1 = uint32...

alrededor de 2 meses ago | 0

convert integer to bytes

Depending on what you are doing downstream with this, you can use the typecast( ) function. E.g., bytes = typecast(my_integer,'...

alrededor de 2 meses ago | 0

| accepted

Selecting multiple different rows in a 3D matrix

You can use a vector for the indexing. E.g., x = matrix(:,:,[2:16,24:48]); Does that, or something like that, do what you wan...

alrededor de 2 meses ago | 0

Using a variable vector in a loop

Suppose you store the column vectors in a matrix called V. Then vk would simply be V(:,k). And if the uj are stored in a matrix...

alrededor de 2 meses ago | 1

| accepted

Fortran code in Matlab

I'm trying to figure out the purpose of the backslash \. If this was the last character in the format statement then I would ass...

alrededor de 2 meses ago | 0

Assigning mxUint16 variable to plhs[0]

You cannot attach native C/C++ memory to an mxArray. The MATLAB API "mxSetEtc" functions check to see that the memory addresses ...

alrededor de 2 meses ago | 0

| accepted

How to separate a vector into two different vectors?

Others have already pointed out better ways of doing this. But to answer your question as to why your current code is not workin...

alrededor de 2 meses ago | 1

Nested "For Loop" Array Size

To keep the indexing the same, just use the top value from the loops. E.g., DataArray = nan(300,50,100,numel(inputValsY3)); No...

alrededor de 2 meses ago | 0

| accepted