function result = parseSoapResponse(response)
%parseSoapResponse Convert the response from a SOAP server into MATLAB types.
%   parseSoapResponse(RESPONSE) converts RESPONSE, a string returned by a SOAP
%   server, into a cell array of approprate MATLAB datatypes.
%
%   Example:
%
%   m = createSoapMessage( ...
%       'urn:xmethodsBabelFish', ...
%       'BabelFish', ...
%       {'en_it','Matthew likes MATLAB.'}, ...
%       {'translationmode','sourcedata'}, ...
%       repmat({'{http://www.w3.org/2001/XMLSchema}string'},1,2));
%
%   response = callSoapService( ...
%       'http://services.xmethods.net:80/perl/soaplite.cgi', ...
%       'urn:xmethodsBabelFish#BabelFish', ...
%        m);
%
%   results = parseSoapResponse(response)
%
%   See also createClassFromWsdl, createSoapMessage, callSoapService.

% Matthew J. Simoneau, June 2003
% $Revision: 1.1.6.5.2.1 $  $Date: 2005/01/22 17:23:28 $
% Copyright 1984-2004 The MathWorks, Inc.

% Parse the text into a DOM.
d = org.apache.xerces.parsers.DOMParser;
try
    d.parse(org.xml.sax.InputSource(java.io.StringReader(response)));
catch
    error('The server''s response was not valid XML:\n\n%s',char(response));
end
response = d.getDocument;

% Extract the SOAP envelope.
envelope = getFirstChildElementNode(response);
if ~isequal(getLocalPart(char(envelope.getNodeName)),'Envelope')
    % Started with a header.  The envelope must be next.
    envelope = getNextSiblingElementNode(envelope);
end

% Extract the SOAP body.
soapBody = getFirstChildElementNode(envelope);
if ~isequal(getLocalPart(char(soapBody.getNodeName)),'Body')
    % Started with a header.  The body must be next.
    soapBody = getNextSiblingElementNode(soapBody);
end

% Pull the result node out of the SOAP Body.
resultNode = getFirstChildElementNode(soapBody);

% If the server returned a fault, turn it into an error.
if strcmp(char(resultNode.getLocalName),'Fault')
    faultstringNode = response.getElementsByTagName('faultstring').item(0);
    error('SOAP Fault: %s',char(faultstringNode.getTextContent))
end

% Build hash table of references.
references = java.util.Hashtable;
refs = response.getElementsByTagName('multiRef');
for i = 1:refs.getLength
    ref = refs.item(i-1);
    references.put(ref.getAttribute('id'),ref);
end

% Convert XML to a structure of appropriate MATLAB values.
resultChild = getFirstChildElementNode(resultNode);
if isempty(resultChild)
    result = [];
else
    result = convert(resultChild,references);
end


%===============================================================================
function s = convert(d,references)

nChildren = d.getChildNodes.getLength;
nextChild = d.getFirstChild;
if (nChildren == 0)
    href = d.getAttribute('href');
    if href.length == 0
        % Real data, but empty value.
        s = convertType(...
            '', ...
            char(d.getAttribute('xsi:type')));
    else
        href = href.substring(1,href.length);
        s = convert(references.get(href),references);
    end
elseif (nextChild.getNodeType == nextChild.TEXT_NODE) && (nChildren == 1)
    % We've reached real data.  Convert it.
    s = convertType(...
        char(nextChild.getNodeValue), ...
        char(d.getAttribute('xsi:type')));
else
    s = [];
    while ~isempty(nextChild)
        if (nextChild.getNodeType ~= nextChild.TEXT_NODE)
            name = char(nextChild.getNodeName);
            name = name(max([find(name == ':') 0])+1:end);
            value = convert(nextChild,references);
            if isfield(s,name)
                if isstruct(value)
                    % We'd like to do something simple, like
                    %   s.(name)(end+1,1) = value;
                    % but this will error if the structures don't match.
                    % Pad out both structures with the missing fields and [].
                    existingFields = fieldnames(s.(name));
                    newFields = fieldnames(value);
                    fieldsToAdd = setdiff(newFields,existingFields);
                    for iFieldsToAdd = 1:length(fieldsToAdd)
                        s.(name)(1).(fieldsToAdd{iFieldsToAdd}) = [];
                    end
                    fieldsToAdd = setdiff(existingFields,newFields);
                    for iFieldsToAdd = 1:length(fieldsToAdd)
                        value(1).(fieldsToAdd{iFieldsToAdd}) = [];
                    end
                    s.(name) = [s.(name);value];
                elseif ischar(value)
                    if ischar(s.(name))
                        s.(name) = {s.(name);value};
                    else
                        s.(name){end+1,1} = value;
                    end
                else
                    s.(name)(end+1,1) = value;
                end
            else
                s.(name) = value;
            end
        end
        nextChild = nextChild.getNextSibling;
    end
    % The CHAR below is to handle the empty case.
    if ~isempty(regexp(char(d.getAttribute('xsi:type')),':Array$'))
        if length(fieldnames(s)) == 1
            s = s.(name);
        else
            error('Non-unique elements.')
        end
    end
end

%===============================================================================
function value = convertType(value,type)
% Convert each primitive type to a MATLAB-friendly type.

% TODO: Handle namespaces properly.
type = type(max([find(type == ':') 0])+1:end);

% If "type" is empty, make sure it is a CHAR.
type = char(type);
switch type
    case 'string'
        % Leave it alone.
    case 'boolean'
        switch value
            case 'false'
                value = false;
            case 'true'
                value = true;
            otherwise
                error('not true or false');
        end
    case {'int','integer','decimal','double'}
        value = str2double(value);
    case 'float'
        value = single(str2double(value));
    case {'base64','base64Binary'}
        value = org.apache.axis.encoding.Base64.decode(value);
    otherwise
        %disp(sprintf('Unknown type: "%s".',type))
end


%===============================================================================
function child = getFirstChildElementNode(node)
child = node.getFirstChild;
while ~isempty(child) && (child.getNodeType ~= child.ELEMENT_NODE)
    child = child.getNextSibling;
end

%===============================================================================
function sibling = getNextSiblingElementNode(node)
sibling = node.getNextSibling;
while ~isempty(sibling) && (sibling.getNodeType ~= sibling.ELEMENT_NODE)
    sibling = sibling.getNextSibling;
end

%===============================================================================
function s = getLocalPart(s)
s = regexprep(s,'.*[}:]','');